離散傅立葉變換及譜分析

1、數(shù)字信號(hào)處理實(shí)驗(yàn)實(shí)驗(yàn)二、離散傅立葉變換及譜分析學(xué)院：信息工程學(xué)院班級(jí)：電子101班姓名：*學(xué)號(hào)：*一、實(shí)驗(yàn)?zāi)康恼莆针x散傅里葉變換的計(jì)算機(jī)實(shí)現(xiàn)方法。檢驗(yàn)實(shí)序列傅里葉變換的性質(zhì)。掌握計(jì)算序列的循環(huán)卷積的方法。學(xué)習(xí)用DFT對(duì)連續(xù)信號(hào)和時(shí)域離散信號(hào)進(jìn)行譜分析的方法，了解可能出現(xiàn)的分 析誤差，以便在實(shí)際中正確應(yīng)用DFT。二、實(shí)驗(yàn)內(nèi)容1.實(shí)現(xiàn)序列的離散傅里葉變換并對(duì)結(jié)果進(jìn)行分析。(自己選擇序列，要求包括復(fù) 序列，實(shí)序列，實(shí)偶序列，實(shí)奇序列，虛奇序列)本例檢驗(yàn)實(shí)序列的性質(zhì) DFTxec(n)=ReX(k)DFTxoc(n)=ImX(k)(1)設(shè)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實(shí)驗(yàn)結(jié)果：y =92102860110148160142(3)補(bǔ)零序列的離散傅立葉變換序列X(n) R5(n)，寫(xiě)出序列的傅立葉變換程序和將原序列補(bǔ)零到2

3、0長(zhǎng)序列的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ū)別。補(bǔ)零到40長(zhǎng)：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(信號(hào)的頻譜);40個(gè)有效采樣點(diǎn)：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(沒(méi)有足夠采樣點(diǎn)的信號(hào));Y=dft(x,M);k1=0:M-1;w1=2*pi/M*k1;subplot(2,1,2);stem(w1/pi,abs(Y);title(信號(hào)的頻譜);實(shí)驗(yàn)說(shuō)明:高密度是在原有序列后