數(shù)字信號處理實驗三離散時間傅里葉變換DTFT及IDTFT.

1、數(shù)字信號處理實驗三離散時間傅里葉變換DTFT及IDTFT一、實驗目的：(1) 通過本實驗，加深對DTFT和IDFT的理解。(2) 熟悉應用DTFT對典型信號進行頻譜分析的方法。(3) 掌握用MATLAB進行離散時間傅里葉變換及其逆變換的方法。二、實驗內容： （1）自己生成正弦序列（如矩形序列，正弦序列，指數(shù)序列等），對其進行頻譜分析，觀察其時域波形和頻域的幅頻特性。記錄實驗中觀察到的現(xiàn)象，繪出相應的時域序列和幅頻特性曲線。矩形序列：程序： M=10;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:

2、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 onxlabel('n'),title('(a)序列幅度')subplot(1,3,2),plot(w,X),grid onxlabel('Omega'),title('(b)幅頻特性')subplot(1,3,3),plot(w,X),grid onv=axis;axis(-pi/T,pi/T,v(3),v(4);xlabe

3、l('Omega'),title('(c)橫軸放大后幅頻特性')set(gcf,'color','w')正弦序列：程序:n=-10:10;x=sin(n*pi);k=-200:200;w=(pi/100)*k;X=x*(exp(-j*pi/100).(n'*k);magX=abs(X);angX=angle(X);subplot(3,1,1);stem(n,x,'.k');title('x(n)=sin(n)');subplot(3,1,2);plot(w/pi,magX,'.k&#

4、39;);title('X(ejw)幅度譜');subplot(3,1,3);plot(w/pi,angX,'.k');title('X(ejw)相位譜');n=-10:10;x=sin(n*pi);k=-200:200;w=(pi/100)*k;X=x*(exp(-j*pi/100).(n'*k);magX=abs(X);angX=angle(X);subplot(3,1,1);stem(n,x,'.k');title('x(n)=sin(n)');subplot(3,1,2);plot(w/pi,mag

5、X,'.k');title('X(ejw)幅度譜');subplot(3,1,3);plot(w/pi,angX,'.k');title('X(ejw)相位譜');波形如下:指數(shù)序列：程序：n=-5:5;x=(-0.1).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('

6、;frequency in units ofpi');ylabel('x')gtext('Magnitde Part')subplot(2,1,2);plot(w/pi,angX,'k')/pi,grid;axis(-2,2,-4,4)xlabel('frequency in units ofpi');ylabel('radianspi')gtext('Angle Part');2.對于理想的低通，高通濾波器，用 IDTFT 求出它的逆變換所對應得離散時間序列。記錄實驗中觀察到的現(xiàn)象，繪出相

7、應的時域序列曲線。要求濾波器的截至頻率可由用戶在 MATLAB 界面自行輸入。 程序： wc=0.5*pi;n=-10:10+1e-10;hd=sin(n*wc)./(n*pi);subplot(1,2,1);plot(-pi,-wc,-wc,wc,wc,pi,0,0,1,1,0,0)xlabel('頻率(1/秒)');ylabel('幅度');axis(-pi,pi,-0.1,1.1),grid onsubplot(1,2,2);stem(n,hd),grid onxlabel('n');ylabel('序列');axis(-10,10,-0.1*wc,0.4*wc)set(gcf,'color','w')三、思考題 離散時間