19-3.m

高阶谱分析包

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模型估计的双谱三维图')