18-2.m

高阶谱分析包

load qpc；
%间接双谱估计；
bspec=bispeci(zmat,21,64,0,'unbiased',128,1);
%直接双谱估计;
bsped=bispecd(zmat,128,1,64,0);
%直接互双谱估计
dbspec=bispecdx(zmat,zmat,zmat,128,3,64,0);