clc;
clear?all;
close?all;

%產生仿真所需信號
len?=?256;
M?=?20;
sigma?=?1;
v?=?normrnd（0sigma1len*M）;
B1?=?[1?0.3544?0.3508?0.1736?0.2401];
A1?=?[1?-1.3817?1.5632?-0.8843?0.4906];
%?B1?=?[1?1.5857?0.9604?0?0];
%?A1?=?[1?-1.6408?2.2044?-1.4808?0.8145];
s?=?filter（B1A1v）;

%求真實的頻譜
n?=?-len/2?:?len/2-1;
w?=?（2*pi/len）?*?n;
z?=?exp（1i*w）;
Px?=?（B1（1）?+?B1（2）./z?+?B1（3）./（z.^2）?+?B1（4）./（z.^3）?+?B1（5）./（z.^4））.^2;
Px?=?Px./（A1（1）?+?A1（2）./z?+?A1（3）./z.^2?+?A1（4）./z.^3?+?A1（5）./z.^4）.^2;
Px?=?10*log10（Px）;

%用ARMA模型法進行功率譜估計
p?=?4;
q?=?4;
fs?=?2;
f?=?linspace（-fs/2fs/22048）;
Sum_Px?=?zeros（1length（f））;
figure;
for?j?=?1?:?M
????z?=?iddata（（s（（1+（j-1）*len）:（len+（j-1）*len）