利用matlab編寫(xiě)dtft,觀察群延時(shí)

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

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

3、s40omga=pi8:P:pi8s=[zeros(1fs)Aones(1flo(Tfs)1)]%構(gòu)造門(mén)函數(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