数字信号处理实验三离散时间傅里叶变换DTFT及IDTFT

处理实验三一、实验目的：(2)熟悉应用DTFT对典型信号进行频谱分析的方法。(3)掌握用MATLAB进行离散时间傅里叶变换及其逆变换的方法。二、实验内容： (1)自己生成正弦序列(如矩形序列，正弦序列，指数序列等)，对其进行频谱分析，观察其时域波形和频域的幅频特性。记矩形序列：程序：M=10;N=2*M+1;T=;n=-4*M:4*M;x=[zeros(1,3*M),ones(1,N),zeros(1,3*M)];w=[-15::15]+1e-10;X=sin*N*w*T)./sin*w*T);subplot(1,3,1);stem(n,x,'.');axis([-20,20,,]),gridonxlabel('n'),title('(a)序列幅度')subplot(1,3,2),plot(w,X),gridonxlabel('\Omega'),title('(b)幅频特性')subplot(1,3,3),plot(w,X),gridonv=axis;axis([-pi/T,pi/T,v(3),v(4)]);xlabel('\Omega'),title('(c)横轴放大后幅频特性')setgcfcolor'w')正弦序列：程序:x=sin(n*pi);w=(pi/100)*k;X=x*(exp(-j*pi/100)).^(n'*k);stem(n,x,'.k');title('x(n)=sin(πn)');plot(w/pi,magX,'.k');title('X(e^jw)幅度谱');plot(w/pi,angX,'.k');title('X(e^jw)相位谱');n=-10:10;x=sin(n*pi);w=(pi/100)*k;X=x*(exp(-j*pi/100)).^(n'*k);stem(n,x,'.k');title('x(n)=sin(πn)');plot(w/pi,magX,'.k');title('X(e^jw)幅度谱');plot(w/pi,angX,'.k');title('X(e^jw)相位谱');波形如下:指数序列：程序：n=-5:5;x=.^n;k=-200:200;w=(pi/100)*k;X=x*(exp(-j*pi/100)).^(n'*k);magX=abs(X);angX=angle(X);subplot(2,1,1);plot(w/pi,magX,'k');grid;axis([-2,2,0,15])xlabel('frequencyinunitsof\pi');ylabel('|x|')gtext('MagnitdePart')subplot(2,1,2);plot(w/pi,angX,'k')/pi,grid;axis([-2,2,-4,4])xlabel('frequencyinunitsof\pi');ylabel('radians\pi')gtext('AnglePart');2.对于理想的低通，高通滤波器，用IDTFT求出它的逆变换所对应列。记录实验中观察到的现象，绘出相应的时域序列曲线。要求滤波器的截至频率可由用户在MATLAB界面自行输程序：wc=*pi;n=[-10:10]+1e-10;hd=sin(n*wc)./(n*pi);bplotplot([-pi,-wc,-wc,wc,wc,pi],[0,0,1,1,0,0])xlabel('频率(1/秒)');ylabel('幅度');axis([-pi,pi,,]),gridonsubplot(1,2,2);stem(n,hd),gridonxlabel('n');ylabel('序列');axis([-10,10,*wc,*wc])setgcfco