exa8_8.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 47 行

% exa080501_sigdec.m ,for example 8.5.1 
% to dectect signal from noise

clear;
% 产生数据
N=128;
n=0:N-1;
f1=0.10;
f2=0.19;
A1=1;
A2=3;
s1=A1*cos(2*pi*f1*n);
s2=A2*cos(2*pi*f2*n);
s=s1+s2;
noise=rand(1,N);
x=s+noise;

% 绘制噪声波形
subplot(321)
plot(noise)
title('纯噪声')

% 绘制噪声信号波形
subplot(322)
plot(x)
title('噪声+信号')

% 绘制噪声双谱平面图
subplot(323)
bspecn=bispeci(noise,10);
title('')

% 绘制信号双谱平面图
subplot(324)
bspecx=bispeci(x,10);
title('')

% 绘制噪声双谱三维图
subplot(325)
abs_bspecn=abs(bspecn);
mesh(abs_bspecn)

% 绘制噪声信号双谱三维图
subplot(326)
abs_bspecx=abs(bspecx);
mesh(abs_bspecx)