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

  • 大小: 3KB
    文件類型: .m
    金幣: 1
    下載: 0 次
    發(fā)布日期: 2021-05-29
  • 語言: Matlab
  • 標(biāo)簽: 系統(tǒng)辨識??

資源簡介

采用最小二乘法,利用含有噪聲的輸入輸出數(shù)據(jù)進行系統(tǒng)辨識,得到系統(tǒng)傳遞函數(shù)

資源截圖

代碼片段和文件信息

clear?all;
close?all;

r=35;
Np=15;
N=(r+1)*Np;?%數(shù)據(jù)長度
%%?------生成M序列和輸入數(shù)據(jù)
P=4;%寄存器長度
M=[0?1?0?0?1];
u=zeros(1N);
a=1;?%M序列幅度
s=1;?%方波
for?k=1:N
????M(1)=M(2)+M(5);
????if?M(1)==2
????????M(1)=0;
????end?
????IM=M(1)+s;?%逆M序列
????if?IM==2?||?IM==0
????????u(k)=a;
????else
?????????u(k)=-a;???????
????end????
????for?i=P+1:-1:2?%移位
????????M(i)=M(i-1);
????end?
????s=~s;
end
%%?------生成白噪聲
A=179;
m=32768;
sigma=1;?%白噪聲標(biāo)準(zhǔn)差
v=zeros(1N);?%白噪聲
xi=11;
for?k=1:N
????ksai=0;
????for?i=1:12???
??????????xi=A*xi;
??????????xi=mod(xim);
??????????ksai=ksai+xi/m;
????end
????v(k)=sigma*(ksai-6.0);
end
%%?------生成噪聲
LAMBDA=0.1;?%噪聲標(biāo)準(zhǔn)差
e=zeros(1N);
for?k=3:N
????e(k)=1.5*e(k-1)-0.7*e(k-2)+LAMBDA*v(k);
end
%%?------生成輸出數(shù)據(jù)
y=zeros(1N);?%不含噪聲
for?k=3:N
????y(k)=1.5*y(k-1)-0.7*y(k-2)+u(k-1)+0.5*u(k-2);
end
z=y+e;?%含噪聲
%%?------遺忘因子法
Mu=0.98;?%遺忘因子
L=300;?%數(shù)據(jù)長度
EPSILON=1e-3;
a2=1e6;
n=1;?%?初始階次
nend=10;?%最大階次
J=zeros(1nend);?%損失函數(shù)
if?L==100
????talpha=3.09;?%閾值
else
????if?L==300
????????talpha=3.03;
????else
????????if?L==500
????????????talpha=3.01;
????????end
????end
end
flag=1;
t=talpha+1;
while?(flag?&&?n????THETA=EPSILON*ones(n+n1);?%參數(shù)向量初始值
????P0=a2*eye(n+n);?%協(xié)方差矩陣初始值
????h=zeros(n+n1);?%數(shù)據(jù)向量???
????THETA0=THETA(:1);
????for?k=n+1:L
???????h=[-z(k-1:-1:k-n)u(k-1:-1:k-n)]‘;?%輸入輸出數(shù)據(jù)
???????K=P0*h/(h‘*P0*h+Mu);
????

評論

共有 條評論

相關(guān)資源