資源描述:
《實(shí)驗(yàn)五序列的傅里葉變換》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、數(shù)字信號處理實(shí)驗(yàn)五序列的傅里葉變換一、實(shí)驗(yàn)?zāi)康?1.充分熟悉exp函數(shù)的使用;2.熟悉序列傅里葉變換的計(jì)算;3.能夠畫出結(jié)果的圖形。二、實(shí)驗(yàn)步驟:1.用help查找exp函數(shù)的使用情況;2.編輯并生成函數(shù)FourierTran.m(帶長度的卷積和)function[X,magX,angX]=FourierTran(x,n,dot)%計(jì)算離散序列的付立葉變換%[X,magX,angX]=FourierTran(x,n)%或[X,magX,angX]=FourierTran(x,n,dot)ifnargin<3%nargin是輸入變量的個(gè)數(shù)dot=600;endk=-dot:dot;
2、w=(pi/dot)*k;X=x*(exp(-j).^(n'*w));%做傅里葉變換magX=abs(X);angX=angle(X);subplot(211);plot(w/pi,magX);xlabel('頻率(單位{pi})');ylabel('
3、X(e^{jomega})
4、');title('幅頻特性');subplot(212);plot(w/pi,angX/pi);xlabel('頻率(單位{pi})');ylabel('弧度/{pi}');title('相頻特性');3.編輯并生成函數(shù)sigfourier.m(序列卷積和)x=[1,1,1,1,1];n=0:
5、4;FourierTran(x,n);4.分析X=x*(exp(-j).^(n'*w))語句中n',如果用n替代,結(jié)果會(huì)有什么不同?為什么?5.當(dāng)dot=8時(shí),計(jì)算X幅頻特性的值,并X的幅頻特性的值一一列出來。1.改變x,n,dot的值兩次,畫出兩次傅里葉變換的結(jié)果。三、實(shí)驗(yàn)結(jié)果:1、分析X=x*(exp(-j).^(n'*w))語句中n',如果用n替代,結(jié)果會(huì)有什么不同?為什么?替換前:替換后:會(huì)出錯(cuò),因?yàn)閚'是行列式轉(zhuǎn)置,用n替代將不能正常計(jì)算行列式相乘,因?yàn)閮尚辛惺降牧袛?shù)與行數(shù)不相等。2、當(dāng)dot=8時(shí),計(jì)算X幅頻特性的值,并X的幅頻特性的值一一列出來。1、改變x,n,do
6、t的值兩次,畫出兩次傅里葉變換的結(jié)果。當(dāng)x=[1,2,3,4,3,2,1];n=0:6;FourierTran(x,n,10);當(dāng)x=[6,2,3,4,3,4,1,7];n=0:7;FourierTran(x,n,20);4、思考題:X=x*(exp(-j).^(n'*w));分析該語句如何實(shí)現(xiàn)傅里葉變換,舉例說明。