exa8_6.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 25 行

% exa080404_maspect.m ,for example 8.4.4 
% to test maorder.m,maest.m,and bispect.m 

clear;
% 装载数据
load ma1;

% 估计MA模型阶次
p=maorder(y);

% 计算MA模型参数
bvec=maest(y,p)

% 计算基于MA模型的双谱，并画出其平面图
ma=bvec;
ar=1;
subplot(211)
bisp=bispect(ma,ar,64);
title('基于MA模型估计的双谱平面图')

% 画出基于MA模型的双谱三维图
abs_bisp=abs(bisp);
subplot(212)
mesh(abs_bisp);
title('基于MA模型估计的双谱三维图')