離散傅立葉變換及譜分析

1、數(shù)字信號處理實驗實驗二、離散傅立葉變換及譜分析學院：信息工程學院班級：電子101班姓名：*學號：*一、實驗目的掌握離散傅里葉變換的計算機實現(xiàn)方法。檢驗實序列傅里葉變換的性質(zhì)。掌握計算序列的循環(huán)卷積的方法。學習用DFT對連續(xù)信號和時域離散信號進行譜分析的方法，了解可能出現(xiàn)的分 析誤差，以便在實際中正確應用DFT。二、實驗內(nèi)容1.實現(xiàn)序列的離散傅里葉變換并對結果進行分析。(自己選擇序列，要求包括復 序列，實序列，實偶序列，實奇序列，虛奇序列)本例檢驗實序列的性質(zhì) DFTxec(n)=ReX(k)DFTxoc(n)=ImX(k)(1)設x(n)=10*(0.8).”n(0=nNerror(N mus

2、t be = the length of x1) endif length(x2)Nerror(N must be = the length of x2) endx1=x1 zeros(1,N-length(x1);x2=x2 zeros(1,N-length(x2);m=0:1:N-1;x2=x2(mod(-m,N)+1);H=zeros(N,N);for n=1:1:NH(n,:)=cirshftt(x2,n-1,N);endy=x1*H實驗結果：y =92102860110148160142(3)補零序列的離散傅立葉變換序列X(n) R5(n)，寫出序列的傅立葉變換程序和將原序列補零到2

3、0長序列的DFT。n=0:4;x=ones(1,5);k=0:999;w=(pi/500)*k;X=x*(exp(-j*pi/500).”(n*k);Xe=abs(X);subplot(3,2,1);stem(n,x);ylabel(x(n);subplot(3,2,2);plot(w/pi,Xe);ylabel(|X(ejw)|);N=20;x=ones(1,5),zeros(1,N-5);n=0:1:N-1;X=dft(x,N);magX=abs(X);k=(0:length(magX)-1)*N/length(magX);subplot(3,2,3);stem(n,x);ylabel(x

4、(n);subplot(3,2,4);stem(k,magX);axis(0,20,0,5);ylabel(|X(k)|);(4)高密度譜和高分辨率譜之間的區(qū)別。補零到40長：M=10;n=0:M-1;x=2*cos(0.35*pi*n)+cos(0.5*pi*n);N=40;x=x,zeros(1,N-10);Y=dft(x,N);k1=0:1:N-1;w1=2*pi/N*k1;subplot(2,1,2);stem(w1/pi,abs(Y);title(信號的頻譜);40個有效采樣點：M=40;n=0:M-1;x=2*cos(0.35*pi*n)+cos(0.5*pi*n);subplot(2,1,1);stem(n,x);title(沒有足夠采樣點的信號);Y=dft(x,M);k1=0:M-1;w1=2*pi/M*k1;subplot(2,1,2);stem(w1/pi,abs(Y);title(信號的頻譜);實驗說明:高密度是在原有序列后