exm046_2.m

matlab的初级教程

%exm046_2.m

%（1）构造受噪声污染的信号
clear,randn('state',0)
t=linspace(0,10,512);
y=3*sin(5*t)-6*cos(9*t)+5*randn(size(t));							%<3>
disp('按任意键继续！绘制带噪声的原始信号。')
pause
clf,shg
plot(t,y)  

disp('按任意键继续！绘制幅频曲线。')
pause
Y=fft(y);
Ts=t(2)-t(1)				%时间信号的采样周期
Ws=2*pi/Ts;				%时间信号的采样频率
Wn=Ws/2					%Nyquest频率
w=linspace(0,Wn,length(t)/2);		%半采样频率中相应的刻度
Ya=abs(Y(1:length(t)/2));

shg,plot(w,Ya)  

disp('按任意键继续！绘制局部放大的幅频曲线。')
pause
ii=find(w<=20);
shg
plot(w(ii),Ya(ii))
grid,xlabel('Frequency       Rad/s')