exa8_1.m

新版Matlab 7辅助信号处理技术与应用

% exa080301_bispecd.m ,for example 8.3.1 
% to test bispecd.m 

clear;
% 信号的产生
N=256;
n=0:N-1;
f1=0.15;
fai1=pi/2;
f2=0.25;
fai2=pi/3;
f3=0.3;
fai3=pi/4;
en=randn(size(n));
s1=cos(2*pi*f1*n+fai1);
s2=cos(2*pi*f2*n+fai2);
s3=cos(2*pi*f3*n+fai3);
x=s1+s2+s3+en;

% 双谱估计平面图
[bspec,waxis]=bispecd(x,16);
title('双谱平估计面图')

% 双谱估计三维图
bspec_abs=abs(bspec);
figure
mesh(bspec_abs)
title('双谱估计三维图')