load ar1
%計算AR模型參數
ar(:,1)=arrcest(y,2,0,2,12,128);
ar(:,2)=arrcest(y,2,0,3,12,128);
ar(:,3)=arrcest(y,2,0,4,12,128);
ar(:,4)=arrcest(y,2,0,-3,12,128);
ar(:,5)=arrcest(y,2,0,-4,12,128);
%計算AR模型的雙譜，并畫出平面圖
ma=1
ar1=ar(:,1)'
subplot(211)
bisp=bispect(ma,ar,64)
title('基于AR模型估計的雙譜平面圖');
%畫出基于AR模型的雙譜三維圖
abs_bisp=abs(bisp)
subplot(212)
mesh(abs_bisp)
title('基于AR模型估計的雙譜三維圖')