Description: hf1=zeros(1,N)
hf1(f<=B1&f>=-B1)=1
y1=hf1.*fftshift(y)
figure
subplot(211)
plot(f,abs(fftshift(y))/N,f,hf1, r-- )
axis([-fs/2,fs/2,0,1.2])
title( 原信号频谱 )
xlabel( f /Hz )
subplot(212)
plot(f,abs(y1)/N)
title( 低通滤波，截止频率90Hz )
xlabel( f /Hz )
To Search:
File list (Check if you may need any files):
example4_6.m
Download