用matlab求idtft,数字信号处理实验讲义（修订版）

数字信号处理实验讲义

五、参考程序

1、给出一矩形序列，并且给出其频谱曲线。 程序如下：

M=4;N=2*M+1;T=0.5;n=-4*M:4*M;

x=[zeros(1,3*M),ones(1,N),zeros(1,3*M)]; % 给出输入序列 w=[-15:0.1:15]+1e-10;

X=sin(0.5*N*w*T)./sin(0.5*w*T); % 给出频谱序列 subplot(1,3,1),stem(n,x,'.') % 画出输入序列 axis([-20,20,-0.1,1.1]),grid on xlabel('n'),title('(a) 序列幅度')

subplot(1,3,2),plot(w,X),grid on % 画出频谱序列 xlabel('\\Omega'),title('(b) 幅频特性')

subplot(1,3,3),plot(w,X),grid on % 改变横轴比例，画出频谱序列 v=axis;axis([-pi/T,pi/T,v(3),v(4)])

xlabel('\\Omega'),title('(c) 横轴放大后幅特性')

set(gcf,'color','w') % 置图形背景色为白

2、 图示理想低通滤波器的频率特性曲线，并用IDTFT求出它的反变换。 程序如下： clear,close all

wc=input('wc= (范围从0到pi)'); % 低通滤波器的截止频率 n=[-10:10]+1e-10;

10

数字信号处理实验讲义

hd=sin(n*wc