離散傅里葉變換判斷方法

判斷系統(tǒng)是否為最小相位系統(tǒng)的簡(jiǎn)單方法是：如果兩個(gè)系統(tǒng)的傳遞函數(shù)分子和分母的最高次數(shù)都分別是m，n，則頻率ω趨于無(wú)窮時(shí)，兩個(gè)系統(tǒng)的對(duì)數(shù)幅頻曲線斜率均為－20（n－m）dB/dec但對(duì)數(shù)相頻曲線卻不同：最小相位系統(tǒng)趨于－90°（n－m），而非最小相位系統(tǒng)卻不這樣。

在工程實(shí)際中經(jīng)常遇到的模擬信號(hào)xn(t),其頻譜函數(shù)Xn(jΩ)也是連續(xù)函數(shù)，為了利用DFT對(duì)xn(t)進(jìn)行譜分析，對(duì)xn(t)進(jìn)行時(shí)域采樣得到x(n)= xn(nT),再對(duì)x(n)進(jìn)行DFT，得到X(k)則是x(n)的傅里葉變換X(ejω)在頻率區(qū)間[0,2π]上的N點(diǎn)等間隔采樣，這里x(n)和X(k)都是有限長(zhǎng)序列

然而，傅里葉變換理論證明，時(shí)間有限長(zhǎng)的信號(hào)其頻譜是無(wú)限寬的，反之，弱信號(hào)的頻譜有限寬的則其持續(xù)時(shí)間將為無(wú)限長(zhǎng)，因此，按采樣定理采樣時(shí)，采樣序列應(yīng)為無(wú)限長(zhǎng)，這不滿足DFT的條件。實(shí)際中，對(duì)于頻譜很寬的信號(hào)，為防止時(shí)域采樣后產(chǎn)生‘頻譜混疊’，一般用前置濾波器濾除幅度較小的高頻成分，使信號(hào)的帶寬小于折疊頻率；同樣對(duì)于持續(xù)時(shí)間很長(zhǎng)的信號(hào)，采樣點(diǎn)數(shù)太多也會(huì)導(dǎo)致存儲(chǔ)和計(jì)算困難，一般也是截取有限點(diǎn)進(jìn)行計(jì)算。上述可以看出，用DFT對(duì)模擬信號(hào)進(jìn)行譜分析，只能是近似的，其近似程度取決于信號(hào)帶寬、采樣頻率和截取長(zhǎng)度

模擬信號(hào)xn(t)的傅里葉變換對(duì)為

X(jΩ)={-∞, ∞}x(t)*exp^-jΩt dt

x（t）=1/2π{-∞, ∞} X(JΩ)*e^jΩt dΩ

用DFT方法計(jì)算這對(duì)變換對(duì)的方法如下：

（a）對(duì)xn(t)以T為間隔進(jìn)行采樣，即xn(t)|t=nT= xa(nT)= x(n),由于

t→nT,dt→T, {-∞, ∞}→∑n={-∞, ∞}

因此得到

X(jΩ)≈∑n={-∞, ∞}x(nT)*exp^-jΩnT*T

x（nT）≈1/2π{0, Ωs} X(JΩ)*e^jΩnT Dω

（b）將序列x(n)= xn(t)截?cái)喑砂蠳個(gè)抽樣點(diǎn)的有限長(zhǎng)序列

X(jΩ)≈T∑n={0,N-1}x(nT)*exp^-jΩnT*T

由于時(shí)域抽樣，抽樣頻率為fs=1/T,則頻域產(chǎn)生以fs為周期的周期延拓，如果頻域是帶限信號(hào)，則有可能不產(chǎn)生頻譜混疊，成為連續(xù)周期頻譜序列，頻譜的周期為fs=1/T

（c）為了數(shù)值計(jì)算，頻域上也要抽樣，即在頻域的一個(gè)周期中取N個(gè)樣點(diǎn)，fs=NF0,每個(gè)樣點(diǎn)間隔為F0，頻域抽樣使頻域的積分式變成求和式，而在時(shí)域就得到原來(lái)已經(jīng)截?cái)嗟碾x散時(shí)間序列的周期延拓，時(shí)間周期為T(mén)0=1/F0。因此有

Ω→kΩ0，dΩ→Ω0，{-∞, ∞} dΩ→∑n={-∞, ∞}Ω0

T0=1/F0=N/fs=NT

Ω0=2ΠF0

Ω0T=Ω0/fs=2π/N

X(jkΩ0)≈T∑n={0,N-1}x(nT)*exp^-jkΩ0nT

DFT的一個(gè)重要特點(diǎn)就是隱含的周期性,從表面上看,離散傅里葉變換在時(shí)域和頻域都是非周期的,有限長(zhǎng)的序列,但實(shí)質(zhì)上DFT是從DFS引申出來(lái)的,它們的本質(zhì)是一致的,因此DTS的周期性決定DFT具有隱含的周期性?？梢詮囊韵氯齻€(gè)不同的角度去理解這種隱含的周期性

(1)從序列DFT與序列FT之間的關(guān)系考慮X(k)是對(duì)頻譜X(ejω)在[0,2π]上的N點(diǎn)等間隔采樣，當(dāng)不限定k的取值范圍在[0,N-1]時(shí)，那么k的取值就在[0,2π]以外，從而形成了對(duì)頻譜X(ejω)的等間隔采樣。由于X(ejω)是周期的，這種采樣就必然形成一個(gè)周期序列

(2)從DFT與DFS之間的關(guān)系考慮。X(k)= ∑n={0,N-1}x(n) WNexp^nk，當(dāng)不限定N時(shí)，具有周期性

(3)從WN來(lái)考慮，當(dāng)不限定N時(shí)，具有周期性

1.線性性質(zhì)

如果X1（n）和X2(N)是兩個(gè)有限長(zhǎng)序列，長(zhǎng)度分別為N1和N2，且Y(N)=AX1(N) BX2(N)

式中A,B為常數(shù)，取N=max[N1,N2]，則Y(N）的N點(diǎn)DFT為

Y(K)=DFT[Y(N)]=AX1(K) BX2(K), 0≤K≤N-1;

2.循環(huán)移位特性

設(shè)X(N)為有限長(zhǎng)序列，長(zhǎng)度為N，則X(N)地循環(huán)移位定義為

Y(N)=X((N M))下標(biāo)nR（N）

式中表明將X(N）以N為周期進(jìn)行周期拓延得到新序列X'(N)=X((N))下標(biāo)n，再將X'(N）左移M位，最后取主值序列得到循環(huán)移位序列Y(N)

（1）物理意義

設(shè)x(n)是長(zhǎng)度為N的有限長(zhǎng)序列，則其傅里葉變換，Z變換與離散傅里葉變換分別用以下三個(gè)關(guān)系式表示

X(e^jω)= ∑n={0,N-1}x(n) e^j-ωn

X(z)= ∑n={0,N-1}x(n)z^-n

X(k)= ∑n={0,N-1}x(n) e^-j2πkn/N

單位圓上的Z變換就是序列的傅里葉變換

離散傅里葉變換是x(n)的頻譜X(ejω)在[0,2π]上的N點(diǎn)等間隔采樣,也就是對(duì)序列頻譜的離散化,這就是DFT的物理意義.

離散傅里葉變換（DFT），是傅里葉變換在時(shí)域和頻域上都呈現(xiàn)離散的形式，將時(shí)域信號(hào)的采樣變換為在離散時(shí)間傅里葉變換（DTFT）頻域的采樣。在形式上，變換兩端（時(shí)域和頻域上）的序列是有限長(zhǎng)的，而實(shí)際上這兩組序列都應(yīng)當(dāng)被認(rèn)為是離散周期信號(hào)的主值序列。即使對(duì)有限長(zhǎng)的離散信號(hào)作DFT，也應(yīng)當(dāng)將其看作經(jīng)過(guò)周期延拓成為周期信號(hào)再作變換。在實(shí)際應(yīng)用中通常采用快速傅里葉變換以高效計(jì)算DFT。

（1）時(shí)域和頻域混疊

根據(jù)采樣定理，只有當(dāng)采樣頻率大于信號(hào)最高頻率的兩倍時(shí)，才能避免頻域混疊。實(shí)際信號(hào)的持續(xù)時(shí)間是有限的，因而從理論上來(lái)說(shuō)，其頻譜寬度是無(wú)限的，無(wú)論多 大的采樣頻率也不能滿足采樣定理。但是超過(guò)一定范圍的高頻分量對(duì)信號(hào)已沒(méi)有多大的影響，因而在工程上總是對(duì)信號(hào)先進(jìn)行低通濾波

另一方面，DFT得到的頻率函數(shù)也是離散的，其頻域抽樣間隔為F0，即頻率分辨力。為了對(duì)全部信號(hào)進(jìn)行采樣，必須是抽樣點(diǎn)數(shù)N滿足條件

N=T0/T=fs/F0

從以上兩個(gè)公式來(lái)看，信號(hào)最高頻率分量fc和頻率分辨力F0有矛盾。若要fc增加，則抽樣間隔T就要減小，而FS就要增加，若在抽樣點(diǎn)數(shù)N不變的情況下，必然是F0增加，分辨力下降。唯一有效的方法是增加記錄長(zhǎng)度內(nèi)的點(diǎn)數(shù)N，在fc和F0給定的條件下，N必須滿足

N>2fc/F0

（2）截?cái)嘈?yīng)

在實(shí)際中遇到的序列x(n)，其長(zhǎng)度往往是有限長(zhǎng)，甚至是無(wú)限長(zhǎng)，用DFT對(duì)其進(jìn)行譜分析時(shí)，必須將其截?cái)酁殚L(zhǎng)度為N的有限長(zhǎng)序列

Y(n)=x(n).R_N(n)

根據(jù)頻率卷積定理

Y(e)=1/2Πx(e)*H(e)

|ω|<2π/N叫做主瓣，其余部分叫做旁瓣

（3）頻譜泄露

原序列x(n)的頻譜是離散譜線，經(jīng)截?cái)嗪笫姑扛V線都帶上一個(gè)辛格譜，就好像使譜線向兩邊延申，通常將這種是遇上的截?cái)鄬?dǎo)致頻譜展寬成為泄露，泄露使得頻譜變得模糊，分辨率降低

（4）譜間干擾

因截?cái)嗍怪髯V線兩邊形成許多旁瓣，引起不同分量間的干擾，成為譜間干擾，這不僅影響頻譜分辨率，嚴(yán)重時(shí)強(qiáng)信號(hào)的旁瓣可能湮滅弱信號(hào)的主譜線。

截?cái)嘈?yīng)是無(wú)法完全消除的，只能根據(jù)要求折中選擇有關(guān)參量。

（5）柵欄效應(yīng)

N點(diǎn)DFT是在頻率區(qū)間[0,2π]上對(duì)信號(hào)的頻譜進(jìn)行N點(diǎn)等間隔采樣，得到的是若干個(gè)離散點(diǎn)X(k)，且它們之限制為基頻F0的整數(shù)倍，這部好像在柵欄的一邊通過(guò)縫隙看另一邊的景象，只能在離散點(diǎn)的地方看到真實(shí)的景象，其余部分頻譜成分被遮攔，所以稱(chēng)為柵欄效應(yīng)。

減小柵欄效應(yīng)，可以在時(shí)域數(shù)據(jù)末端增加一些零值點(diǎn)，是一個(gè)周期內(nèi)的點(diǎn)數(shù)增加

（6）信號(hào)長(zhǎng)度的選擇

在時(shí)域內(nèi)對(duì)信號(hào)長(zhǎng)度的選擇會(huì)影響DFT運(yùn)算的正確性。實(shí)際的信號(hào)往往是隨機(jī)的，沒(méi)有確定的周期，因此在實(shí)際中，應(yīng)經(jīng)可能估計(jì)出幾個(gè)典型的、帶有一定周期性的信號(hào)區(qū)域進(jìn)行頻譜分析，然后在取其平均值，從而得到合理的結(jié)果。

1.C語(yǔ)言實(shí)現(xiàn)代碼

intDFT(intdir,intm,double*x1,double*y1)
{
longi,k;
doublearg;
doublecosarg,sinarg;
double*x2=NULL,*y2=NULL;
x2=malloc(m*sizeof(double));
y2=malloc(m*sizeof(double));
if(x2==NULL||y2==NULL)return(FALSE);
for(i=0;i
 

離散余弦變換(DCT for Discrete Cosine Transform)是與傅里葉變換相關(guān)的一種變換，它類(lèi)似于離散傅里葉變換(DFT for Discrete Fourier Transform),但是只使用實(shí)數(shù)。離散余弦變換相當(dāng)于一個(gè)長(zhǎng)度大概是它兩倍的離散傅里葉變換，這個(gè)離散傅里葉變換是對(duì)一個(gè)實(shí)偶函數(shù)進(jìn)行的（因?yàn)橐粋€(gè)實(shí)偶函數(shù)的傅里葉變換仍然是一個(gè)實(shí)偶函數(shù)），在有些變形里面需要將輸入或者輸出的位置移動(dòng)半個(gè)單位(DCT有8種標(biāo)準(zhǔn)類(lèi)型，其中4種是常見(jiàn)的)。

最常用的一種離散余弦變換的類(lèi)型是下面給出的第二種類(lèi)型，通常我們所說(shuō)的離散余弦變換指的就是這種。它的逆，也就是下面給出的第三種類(lèi)型，通常相應(yīng)的被稱(chēng)為"反離散余弦變換"，"逆離散余弦變換"或者"IDCT"。

有兩個(gè)相關(guān)的變換，一個(gè)是離散正弦變換(DST for Discrete Sine Transform),它相當(dāng)于一個(gè)長(zhǎng)度大概是它兩倍的實(shí)奇函數(shù)的離散傅里葉變換；另一個(gè)是改進(jìn)的離散余弦變換(MDCT for Modified Discrete Cosine Transform),它相當(dāng)于對(duì)交疊的數(shù)據(jù)進(jìn)行離散余弦變換。

離散余弦變換，尤其是它的第二種類(lèi)型，經(jīng)常被信號(hào)處理和圖像處理使用，用于對(duì)信號(hào)和圖像(包括靜止圖像和運(yùn)動(dòng)圖像)進(jìn)行有損數(shù)據(jù)壓縮。這是由于離散余弦變換具有很強(qiáng)的"能量集中"特性:大多數(shù)的自然信號(hào)(包括聲音和圖像)的能量都集中在離散余弦變換后的低頻部分，而且當(dāng)信號(hào)具有接近馬爾科夫過(guò)程(Markov processes)的統(tǒng)計(jì)特性時(shí)，離散余弦變換的去相關(guān)性接近于K-L變換(Karhunen-Loève 變換--它具有最優(yōu)的去相關(guān)性)的性能。

例如，在靜止圖像編碼標(biāo)準(zhǔn)JPEG中，在運(yùn)動(dòng)圖像編碼標(biāo)準(zhǔn)MJPEG和MPEG的各個(gè)標(biāo)準(zhǔn)中都使用了離散余弦變換。在這些標(biāo)準(zhǔn)制中都使用了二維的第二種類(lèi)型離散余弦變換，并將結(jié)果進(jìn)行量化之后進(jìn)行熵編碼。這時(shí)對(duì)應(yīng)第二種類(lèi)型離散余弦變換中的n通常是8，并用該公式對(duì)每個(gè)8x8塊的每行進(jìn)行變換，然后每列進(jìn)行變換。得到的是一個(gè)8x8的變換系數(shù)矩陣。其中(0,0)位置的元素就是直流分量，矩陣中的其他元素根據(jù)其位置表示不同頻率的交流分量。

一個(gè)類(lèi)似的變換, 改進(jìn)的離散余弦變換被用在高級(jí)音頻編碼(AAC for Advanced Audio Coding)，Vorbis 和 MP3 音頻壓縮當(dāng)中。

離散余弦變換也經(jīng)常被用來(lái)使用譜方法來(lái)解偏微分方程，這時(shí)候離散余弦變換的不同的變量對(duì)應(yīng)著數(shù)組兩端不同的奇/偶邊界條件。

離散系數(shù)是衡量資料中各觀測(cè)值離散程度的一個(gè)統(tǒng)計(jì)量。當(dāng)進(jìn)行兩個(gè)或多個(gè)資料離散程度的比較時(shí)，如果度量單位與平均數(shù)相同，可以直接利用標(biāo)準(zhǔn)差來(lái)比較。如果單位和（或）平均數(shù)不同時(shí)，比較其離散程度就不能采用標(biāo)準(zhǔn)差，而需采用標(biāo)準(zhǔn)差與平均數(shù)的比值（相對(duì)值）來(lái)比較 ：

離散系數(shù)通?？梢赃M(jìn)行多個(gè)總體的對(duì)比，通過(guò)離散系數(shù)大小的比較可以說(shuō)明不同總體平均指標(biāo)（一般來(lái)說(shuō)是平均數(shù)）的代表性或穩(wěn)定性大小。一般來(lái)說(shuō)，離散系數(shù)越小，說(shuō)明平均指標(biāo)的代表性越好；離散系數(shù)越大，平均指標(biāo)的代表性越差。

離散系數(shù)只對(duì)由比率標(biāo)量計(jì)算出來(lái)的數(shù)值有意義。舉例來(lái)說(shuō)，對(duì)于一個(gè)氣溫的分布，使用開(kāi)爾文或攝氏度來(lái)計(jì)算的話并不會(huì)改變標(biāo)準(zhǔn)差的值，但是溫度的平均值會(huì)改變，因此使用不同的溫標(biāo)的話得出的變異系數(shù)是不同的。也就是說(shuō)，使用區(qū)間標(biāo)量得到的變異系數(shù)是沒(méi)有意義的。

離散系數(shù)在概率論的許多分支中都有應(yīng)用，比如說(shuō)在更新理論、排隊(duì)理論和可靠性理論中。在這些理論中，指數(shù)分布通常比正態(tài)分布更為常見(jiàn)。

由于指數(shù)分布的標(biāo)準(zhǔn)差等于其平均值，所以它的離散系數(shù)等于一。離散系數(shù)小于一的分布，比如愛(ài)爾朗分布稱(chēng)為低差別的 ，而離散系數(shù)大于一的分布，如超指數(shù)分布則被稱(chēng)為高差別的。