91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 2KB
    文件類型: .m
    金幣: 1
    下載: 0 次
    發布日期: 2021-06-04
  • 語言: Matlab
  • 標簽:

資源簡介

遞推極大似然參數辨識法MATLAB程序 clear all%清理工作間變量 close all%關閉所有圖形 clc%清屏 %%%% M序列、噪聲信號產生%%%% L=1200;%四位移位積存器產生的M序列的周期 y1=1;y2=1;y3=1;y4=0;%四個移位積存器的輸出初始值 for i=1:L; x1=xor(y3,y4);%第一個移位積存器的輸入信號 x2=y1;%第二個移位積存器的輸入信號 x3=y2;%第三個移位積存器的輸入信號 x4=y3;%第四個移位積存器的輸入信號 y(i)=y4;%第四個移位積存器的輸出信號,幅值"0"和"1" if y(i)>0.5,u(i)=-1;%M序列的值為"1"時,辨識的輸入信號取“-1” else u(i)=1;%M序列的值為"0"時,辨識的輸入信號取“1” end y1=x1;y2=x2;y3=x3;y4=x4;%為下一次的輸入信號作準備 end ------

資源截圖

代碼片段和文件信息

clear?all%清理工作間變量
close?all%關閉所有圖形
clc%清屏
%%%%?M序列、噪聲信號產生%%%%
L=1200;%四位移位積存器產生的M序列的周期
y1=1;y2=1;y3=1;y4=0;%四個移位積存器的輸出初始值
for?i=1:L;
????x1=xor(y3y4);%第一個移位積存器的輸入信號
????x2=y1;%第二個移位積存器的輸入信號
????x3=y2;%第三個移位積存器的輸入信號
????x4=y3;%第四個移位積存器的輸入信號
????y(i)=y4;%第四個移位積存器的輸出信號,幅值“0“和“1“
????if?y(i)>0.5u(i)=-1;%M序列的值為“1“時辨識的輸入信號取“-1”
????else?u(i)=1;%M序列的值為“0“時辨識的輸入信號取“1”
????end
????y1=x1;y2=x2;y3=x3;y4=x4;%為下一次的輸入信號作準備
end
v=randn(11200);?%產生一組1200個正態分布的隨機噪聲
V=0;?%計算噪聲方差
for?i=1:1200
????V=V+v(i)*v(i);
end
????V1=V/1200?%計算噪聲方差結束
%%%%?M序列、噪聲信號產生結束?%%%
%%%遞推極大似然算法開始%%%
z(2)=0;z(1)=0;
for?k=3:1200?%根據v和u計算z(k)
????z(k)=-0.5*z(k-1)+0.2*z(k-2)+1.2*u(k-1)+0.3*u(k-2)+1.0*v(k)-1.0*v(k-1)+0.8*v(k-2);
end
????theta1=0.001*ones(61);p0=eye(66);
????zf(1)=0;zf(2)=0;vf(1)=0;vf(2)=0;uf(1)=0;uf(2)=0;
for?k=3:1200?%遞推迭代計算
????h=[-z(k-1);-z(k-2);u(k-1);u(k-2);v(k-1);v(k-2)];
????hf=h;
????K=p0*hf*inv(hf‘*p0*hf+1);
????p=[eye(66)-K*hf‘]*p0;
????v(k)=z(k)-h‘*th

評論

共有 條評論

相關資源