利用 ARMA、AR、MA 模型,以及周期图等进行系统参数估计
[Copy to clipboard] [ - ]
CODE:
N=456;
B1=[1 0.3544 0.3508 0.1736 0.2401];
A1=[1 -1.3817 1.5632 -0.8843 0.4096];
w=linspace(0,pi,512);
H1=freqz(B1,A1,w);%产生信号的频域响应
Ps1=abs(H1).^2;
SPy11=0;%20 次 AR(4)
SPy12=0;%20 次 AR(8)
SPy13=0;%20 次 AR 周期图
SPy14=0;%20 次 ARMA(4,4)
SPy15=0;%20 次 ARMA(8,8)
VSPy11=0;%20 次 AR(4)
VSPy12=0;%20 次 AR(8)
VSPy13=0;%20 次 AR 周期图
VSPy14=0;%20 次 ARMA(4,4)
VSPy15=0;%20 次 ARMA(8,8)
for k=1:20
%采用自协方差法对 AR 模型参数进行估计%
%gA1:AR 模型的参数;gE1:激励白噪声的方差%
y1=filter(B1,A1,randn(1,N)).*[zeros(1,200),ones(1,256)];
[Py11,F]=pcov(y1,4,512,1);%AR(4)的估计%
[Py12,F]=pcov(y1,8,512,1);%AR(8)的估计%
[Py13,F]=periodogram(y1,[],512,1);
SPy11=SPy11+Py11;
SPy12=SPy12+Py12;
SPy13=SPy13+Py13;
VSPy11=VSPy11+abs(Py11).^2;
VSPy12=VSPy12+abs(Py12).^2;
VSPy13=VSPy13+abs(Py13).^2;
figure(1)
plot(w./(2*pi),Ps1,F,Py11);
legend('真实功率谱','20 次 AR(4)估计图');
hold on;
figure(2)
plot(w./(2*pi),Ps1,F,Py12);
legend('真实功率谱','20 次 AR(8)估计图');
hold on;
figure(3)
评论0