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

處理實驗三一、實驗?zāi)康模?2)熟悉應(yīng)用DTFT對典型信號進行頻譜分析的方法。(3)掌握用MATLAB進行離散時間傅里葉變換及其逆變換的方法。二、實驗內(nèi)容： (1)自己生成正弦序列(如矩形序列，正弦序列，指數(shù)序列等)，對其進行頻譜分析，觀察其時域波形和頻域的幅頻特性。記矩形序列：程序：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)相位譜');波形如下:指數(shù)序列：程序：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求出它的逆變換所對應(yīng)列。記錄實驗中觀察到的現(xiàn)象，繪出相應(yīng)的時域序列曲線。要求濾波器的截至頻率可由用戶在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