exa8_26.m

有关matlab r2007的实用教程

t = 0:0.001:0.6;  
%产生原始信号
x = sin(2*pi*50*t)+sin(2*pi*120*t);
%对原始信号加入0平均的随机噪音干扰
subplot(2,2,1)
plot(1000*t(1:50),x(1:50))
title('（a）原始信号')
ylabel('u /v')
xlabel('t /ms')

yt = 2*randn(size(t));
subplot(2,2,2)
plot(1000*t(1:50),yt(1:50))
title('（b）零均值随机噪声')
ylabel('u /v')
xlabel('t /ms')

y = x + yt;
%绘制干扰后的信号频谱图
subplot(2,2,3)
plot(1000*t(1:50),y(1:50))
title('（c）零均值随机噪声干扰的信号')
ylabel('u /v')
xlabel('t /ms')

%运行512点得快速傅立叶变换
Y = fft(y,512);
Pyy = Y.* conj(Y) / 512;
%获得开始的257个数据点，其余的225个点不用
f = 1000*(0:256)/512;
%绘制频谱图
subplot(2,2,4)
plot(f,Pyy(1:257))
title('（d）fft变换频谱')
ylabel('u /v')
xlabel('f /Hz')