exa8_2.m

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

% exa080302_bispeci.m ,for example 8.3.2 
% to test bispeci.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;

subplot(2,2,1)
plot(n/N,x)
title('原始信号')

% 信号的FFT变换
X=fft(x);
X=fftshift(X);
subplot(2,2,2)
plot((n-N/2)/N,abs(X))
title('傅立叶变换')

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

% 双谱三维图
bspec_abs=abs(bspec);
subplot(2,2,4)
mesh(bspec_abs)
title('双谱估计三维图')