利用matlab編寫dtft,觀察群延時

1、一、實驗目的一、實驗目的1.自行編寫功能函數(shù)實現(xiàn)DTFT的變化，并運用該函數(shù)實現(xiàn)之后的功能。2.利用編寫的DTFT函數(shù)，觀察函數(shù)的離散傅里葉變換，并在三維空間中觀察函數(shù)的頻響特性。3.利用MATLAB設計奇對稱系統(tǒng)觀察群延時對輸入信號的影響。二、實驗條件二、實驗條件PC機，MATLAB7.03、實驗內(nèi)容實驗內(nèi)容1）編寫功能函數(shù)mydtft。根據(jù)公式，將h[n]用樣值向量h和位置向量n描述，???????????nsnjffenhH?2]

2、[)(則DTFT可以用一個向量和一個矩陣的相乘來實現(xiàn)。函數(shù)如下：functionH=mydtft(homga)L=length(h)n=[0:L1]M=exp(1inomga)H=hM2）做出DTFT，由于MATLAB無法計算連續(xù)變量，只能在某個區(qū)間范圍內(nèi)將omga賦值為很密的、長度很長的向量來近似為連續(xù)變量。由可知，當fs夠高時，h[n]可近似于h(t)。???????fethHtj????2][)(fs=50T=2A=5P=2pif

3、s40omga=pi8:P:pi8s=[zeros(1fs)Aones(1flo(Tfs)1)]%構(gòu)造門函數(shù)S=mydtft(somga)plot(omgaabs(S))gridonfigureplot3(omgareal(S)imag(S))gridonfigureP=2pifs40omga=0.2pi:P:pi2.2S=mydtft(somga)plot(omgaabs(S))figuresd=[sss]S=mydtft(sdomg