例9-26 .m

这是一个MATLAB7.0基础与提高例题的所有源码

t = 0:0.001:0.6;  
％产生原始信号
x = sin(2*pi*50*t)+sin(2*pi*120*t);
％对原始信号加入0平均的随机噪音干扰
y = x + 2*randn(size(t))
％绘制干扰后的信号频谱图，见图9-28所示
plot(1000*t(1:50),y(1:50))
title('Signal Corrupted with Zero.Mean Random Noise')
xlabel('time (milliseconds)')
％运行512点得快速傅立叶变换
Y = fft(y,512)
Pyy = Y.* conj(Y) / 512
%获得开始的257个数据点，其余的225个点不用
f = 1000*(0:256)/512
%绘制频谱图，见图9-29示
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')