最美情侣中文字幕电影,在线麻豆精品传媒,在线网站高清黄,久久黄色视频

歡迎光臨散文網(wǎng) 會(huì)員登陸 & 注冊(cè)

數(shù)字信號(hào)處理復(fù)習(xí)總結(jié)-最終版

2023-10-23 18:43 作者:答案資料  | 我要投稿

緒論:本章介紹數(shù)字信號(hào)處理課程的基本概念。

0.1信號(hào)、系統(tǒng)與信號(hào)處理

1.信號(hào)及其分類

信號(hào)是信息的載體,以某種函數(shù)的形式傳遞信息。這個(gè)函數(shù)可以是時(shí)間域、頻率域或其它域,但最基礎(chǔ)的域是時(shí)域。

分類:

周期信號(hào)/非周期信號(hào)

確定信號(hào)/隨機(jī)信號(hào)

能量信號(hào)/功率信號(hào)

連續(xù)時(shí)間信號(hào)/離散時(shí)間信號(hào)/數(shù)字信號(hào)

按自變量與函數(shù)值的取值形式不同分類:

2.系統(tǒng)

系統(tǒng)定義為處理(或變換)信號(hào)的物理設(shè)備,或者說(shuō),凡是能將信號(hào)加以變換以達(dá)到人們要求的各種設(shè)備都稱為系統(tǒng)。

3.信號(hào)處理

信號(hào)處理即是用系統(tǒng)對(duì)信號(hào)進(jìn)行某種加工。包括:濾波、分析、變換、綜合、壓縮、估計(jì)、識(shí)別等等。所謂“數(shù)字信號(hào)處理”,就是用數(shù)值計(jì)算的方法,完成對(duì)信號(hào)的處理。

0.2 數(shù)字信號(hào)處理系統(tǒng)的基本組成

數(shù)字信號(hào)處理就是用數(shù)值計(jì)算的方法對(duì)信號(hào)進(jìn)行變換和處理。不僅應(yīng)用于數(shù)字化信號(hào)的處理,而且


也可應(yīng)用于模擬信號(hào)的處理。以下討論模擬信號(hào)數(shù)字化處理系統(tǒng)框圖。

(1)前置濾波器

將輸入信號(hào)xa(t)中高于某一頻率(稱折疊頻率,等于抽樣頻率的一半)的分量加以濾除。

(2)A/D變換器

在A/D變換器中每隔T秒(抽樣周期)取出一次xa(t)的幅度,抽樣后的信號(hào)稱為離散信號(hào)。在A/D變換器中的保持電路中進(jìn)一步變換為若干位碼。

(3)數(shù)字信號(hào)處理器(DSP)

(4)D/A變換器

按照預(yù)定要求,在處理器中將信號(hào)序列x(n)進(jìn)行加工處理得到輸出信號(hào)y(n)。由一個(gè)二進(jìn)制碼流產(chǎn)生一個(gè)階梯波形,是形成模擬信號(hào)的第一步。

(5)模擬濾波器

把階梯波形平滑成預(yù)期的模擬信號(hào);以濾除掉不需要的高頻分量,生成所需的模擬信號(hào)ya(t)。

0.3 數(shù)字信號(hào)處理的特點(diǎn)

(1)靈活性。(2)高精度和高穩(wěn)定性。(3)便于大規(guī)模集成。(4)對(duì)數(shù)字信號(hào)可以存儲(chǔ)、運(yùn)算、系統(tǒng)可以獲得高性能指標(biāo)。

0.4 數(shù)字信號(hào)處理基本學(xué)科分支

數(shù)字信號(hào)處理(DSP)一般有兩層含義,一層是廣義的理解,為數(shù)字信號(hào)處理技術(shù)——DigitalSignalProcessing,另一層是狹義的理解,為數(shù)字信號(hào)處理器——DigitalSignalProcessor。

0.5 課程內(nèi)容

該課程在本科階段主要介紹以傅里葉變換為基礎(chǔ)的“經(jīng)典”處理方法,包括:(1)離散傅里葉變換及其快速算法。(2)濾波理論(線性時(shí)不變離散時(shí)間系統(tǒng),用于分離相加性組合的信號(hào),要求信號(hào)頻譜占據(jù)不同的頻段)。

在研究生階段相應(yīng)課程為“現(xiàn)代信號(hào)處理”(AdvancedSignalProcessing)。信號(hào)對(duì)象主要是隨機(jī)信號(hào),主要內(nèi)容是自適應(yīng)濾波(用于分離相加性組合的信號(hào),但頻譜占據(jù)同一頻段)和現(xiàn)代譜估計(jì)。

簡(jiǎn)答題:

1.按自變量與函數(shù)值的取值形式是否連續(xù)信號(hào)可以分成哪四種類型?

2.相對(duì)模擬信號(hào)處理,數(shù)字信號(hào)處理主要有哪些優(yōu)點(diǎn)?

3.?dāng)?shù)字信號(hào)處理系統(tǒng)的基本組成有哪些?


第一章:本章概念較多,需要理解和識(shí)記的內(nèi)容較多,學(xué)習(xí)時(shí)要注意。

1.1 離散時(shí)間信號(hào)

1.離散時(shí)間信號(hào)的定義

離散時(shí)間信號(hào)是指一個(gè)實(shí)數(shù)或復(fù)數(shù)的數(shù)字序列,它是整數(shù)自變量n的函數(shù),表示為x(n)。一般由模擬信號(hào)等間隔采樣得到:。時(shí)域離散信號(hào)有三種表示方法:1)用集合符號(hào)表示 2)用公式表示 3)用圖形表示

2.幾種基本離散時(shí)間信號(hào)(記住定義)

(1)單位采樣序列

(2)單位階躍序列

(3)矩形序列

(4)實(shí)指數(shù)序列

(5)正弦序列

ω是正弦序列數(shù)字域的頻率,單位是弧度。

對(duì)連續(xù)信號(hào)中的正弦信號(hào)進(jìn)行采樣,可得正弦序列。設(shè)連續(xù)信號(hào)為,它的采樣值為,因此(重點(diǎn))

這個(gè)式子具有一般性,它反映了由連續(xù)信號(hào)采樣得到的離散序列,其數(shù)字頻率與模擬頻率的一般關(guān)系。另外需要說(shuō)明的是,ω的單位為弧度,Ω的單位為弧度/秒。本書中,我們一律以ω表示數(shù)字域頻率,而以Ω及f表示模擬域頻率。

例:已知采樣頻率FT = 1000Hz, 則序列x(n) = cos(0.4πn) 對(duì)應(yīng)的模擬頻率為 ( 400π ) 弧度/s。

說(shuō)明:本題旨在理解數(shù)字頻率與模擬頻率之間的關(guān)系:。

(6)復(fù)指數(shù)序列

復(fù)指數(shù)序列是以余弦序列為實(shí)部、正弦序列為虛部所構(gòu)成的一個(gè)復(fù)數(shù)序列。

(7)周期序列(重點(diǎn))

所有

存在一個(gè)最小的正整數(shù)

,滿足:

,則稱序列

是周期序列,周期為

。(注意:按此定義,模擬信號(hào)是周期信號(hào),采用后的離散信號(hào)未必是周期的)

例:正弦序列

的周期性:
當(dāng)

,

為整數(shù)時(shí),

,即為周期性序列。周期

,式中,

、

限取整數(shù),且

的取值要保證

是最小的正整數(shù)。

可分幾種情況討論如下:(1)當(dāng)

為整數(shù)時(shí),只要

,

就為最小正整數(shù),即周期為

。(2)當(dāng)

不是整數(shù),而是一個(gè)有理數(shù)時(shí),設(shè)

,式中,

、

是互為素?cái)?shù)的整數(shù)(互為素?cái)?shù)就是兩個(gè)數(shù)沒(méi)有公約數(shù)),取

,則

,即周期為

。(3)當(dāng)

是無(wú)理數(shù)時(shí),則任何

皆不能使

為正整數(shù),這時(shí),正弦序列不是周期性的。

例:X(n) = cos(0.4πn)的基本周期為( 5 )。

[說(shuō)明]基本周期的定義即計(jì)算公式:,其中N和k均為整數(shù),N為基本周期(使得N為最小整數(shù)時(shí)k取值)。本題ω = 0.4π,代入上式得到:。

3.信號(hào)運(yùn)算

(1)加法:兩個(gè)信號(hào)之和 由同序號(hào)的序列值逐點(diǎn)對(duì)應(yīng)相加得到。

(2)乘法:兩個(gè)信號(hào)之積 由同序號(hào)的序列值逐點(diǎn)對(duì)應(yīng)相乘得到。

(3)移位:當(dāng),序列右移(稱為延時(shí));當(dāng),序列左移(稱為超前)。

(4)翻轉(zhuǎn):

(5)尺度變換:或,其中M和N都是正整數(shù)。

當(dāng)時(shí),序列是通過(guò)取x(n)的每第M個(gè)采樣形成,這種運(yùn)算稱為下采樣。對(duì)于序列,定義如下這種運(yùn)算稱為上采樣。

4.信號(hào)分解(重點(diǎn))

任一信號(hào)x(n)可表示成單位脈沖序列的移位加權(quán)和

簡(jiǎn)記為

1.2 時(shí)域離散系統(tǒng)

時(shí)域離散系統(tǒng)定義

1 線性系統(tǒng)(重點(diǎn))

判定公式
若=,=則

2 時(shí)不變系統(tǒng)(重點(diǎn))

判定公式:y(n)=T[x(n)] y(n-)=T[x(n-)]

例:判斷下列系統(tǒng)是否為線性、時(shí)不變系統(tǒng)。(重點(diǎn))

(1);

(2);

(1)令:輸入為

,輸出為

故該系統(tǒng)是時(shí)不變系統(tǒng)。

故該系統(tǒng)是線性系統(tǒng)。

(2)

令:輸入為

,輸出為

,因?yàn)?/p>

故系統(tǒng)是時(shí)不變系統(tǒng)。又因?yàn)?/p>

因此系統(tǒng)是非線性系統(tǒng)。

3 線性時(shí)不變系統(tǒng)(LTI或者LSI系統(tǒng))輸入與輸出之間關(guān)系(重點(diǎn)):

y(n)==x(n)*h(n)

重點(diǎn):線性離不變系統(tǒng)的輸出等于輸入序列和該系統(tǒng)的單位脈沖響應(yīng)的卷積

【說(shuō)明】離散時(shí)間LTI系統(tǒng)的單位沖激響應(yīng)h(n)為系統(tǒng)對(duì)單位沖激序列δ(n)的零狀態(tài)響應(yīng)。

單位沖激響應(yīng)的概念非常重要。在時(shí)域,LTI系統(tǒng)可以由其單位沖激響應(yīng)h(n)唯一確定,因此,我們常常用單位沖激響應(yīng)描述 LTI 系統(tǒng)。在這種情況下, LTI 系統(tǒng)的輸入輸出關(guān)系可以由卷積運(yùn)算描述:y(n)==x(n)*h(n)

物理意義: 卷積和運(yùn)算具有顯式意義,即可以用來(lái)確定系統(tǒng)的輸出。如果系統(tǒng)確定,則其單位沖激響應(yīng)是唯一的。由此,可求系統(tǒng)對(duì)任意輸入的響應(yīng)。

注意:計(jì)算卷積和的關(guān)鍵是求和區(qū)間的確定。因此,常常需要繪制序列x(m) 和h(n-m)的圖形。利用序列x(m) 和h(n-m)的圖形可助我們方便地確定求和區(qū)間。

卷積的求解方法(重點(diǎn)):

線性卷積是一種非常重要的一種運(yùn)算,對(duì)它的求解,一般我們采用作圖法。線性卷積滿足交換律,設(shè)兩序列長(zhǎng)度分別是N和M,線性卷積后序列的長(zhǎng)度為N+M-1。

卷積的計(jì)算過(guò)程包括翻轉(zhuǎn)、移位、相乘、相加四個(gè)過(guò)程。

1)將和用和表示,畫出和這兩個(gè)序列;

2)選擇一個(gè)序列,并將其按時(shí)間翻轉(zhuǎn)形成序列;

3)將移位n,得到;

4)將和相同m的序列值對(duì)應(yīng)相乘后,再相加。

例:設(shè)

,

,

如圖1所示。求

的卷積

。(重點(diǎn))


圖1

解 方法一:用圖解法求卷積和。

(1) 將和用和表示(圖2中(a)、(b)圖)。


圖2 圖解法求卷積過(guò)程

(2) 將進(jìn)行反折,形成(圖2中(c)圖);將移位,得到(圖2中(d)、(e)、(f)圖)。

(3) 將和相同的序列值相乘,再相加,得到(圖2中(g)圖)。

再討論解析法求線性卷積。

用式

求解上式首先要根據(jù)和的非零值區(qū)間確定求和的上下限,的非零值區(qū)間為,的非零值區(qū)間為,或,由兩個(gè)非零值區(qū)間可得的取值區(qū)間為,它們的乘積的非零值區(qū)間應(yīng)滿足:

因此

當(dāng) 、時(shí),;

當(dāng) 時(shí),;

當(dāng) 時(shí),。

與圖解法結(jié)果一致。

y(n)用公式表示為

方法二:當(dāng)序列和的長(zhǎng)度分別為有限長(zhǎng)和時(shí),可采用“不進(jìn)位乘法”求兩序列線卷積。

如圖1所示:,

例:兩線性時(shí)不變系統(tǒng)級(jí)聯(lián),其單位取樣響應(yīng)分別為

,輸入為

,求系統(tǒng)的輸出

。


已知:

,

,

。

解:設(shè)第一個(gè)系統(tǒng)的輸出為

,則


因而輸出為


4. 系統(tǒng)因果性和穩(wěn)定性的判定(重點(diǎn))

1)穩(wěn)定系統(tǒng):有界的輸入產(chǎn)生的輸出也有界的系統(tǒng),即:若,則(記住!!)

線性移不變系統(tǒng)是穩(wěn)定系統(tǒng)的充要條件:(系統(tǒng)穩(wěn)定的充分必要條件是系統(tǒng)的單位脈沖響應(yīng)絕對(duì)可和)(記住!!)

或:其系統(tǒng)函數(shù)H(z)的收斂域包含單位圓 |z|=1(記住!!)


2)因果系統(tǒng):時(shí)刻的輸出只由時(shí)刻之前的輸入決定(記住!!)

線性移不變系統(tǒng)是因果系統(tǒng)的充要條件:(記住!!)因果系統(tǒng)的單位脈沖響應(yīng)必然是因果序列。(記住!!)

或:其系統(tǒng)函數(shù)H(z)的收斂域在某圓外部:即:|z|>Rx(記住!!)

3)穩(wěn)定因果系統(tǒng):同時(shí)滿足上述兩個(gè)條件的系統(tǒng)。

線性移不變系統(tǒng)是因果穩(wěn)定系統(tǒng)的充要條件:,(記住!!)

或:H(z)的極點(diǎn)在單位圓內(nèi)H(z)的收斂域滿足:(記住!!)

例:判斷線性時(shí)不變系統(tǒng)的因果性、穩(wěn)定性,并給出依據(jù)。(重點(diǎn))

(1);

(2);

:(1)只要,該系統(tǒng)就是因果系統(tǒng),因?yàn)檩敵鲋慌cn時(shí)刻的和n時(shí)刻以前的輸入有關(guān)。如果,則,因此系統(tǒng)是穩(wěn)定系統(tǒng)。

(2)如果,,因此系統(tǒng)是穩(wěn)定的。系統(tǒng)是非因果的,因?yàn)檩敵鲞€和x(n)的將來(lái)值有關(guān)。

注意:如果給出的是h(n),用上面要求記住的充要條件判斷!

例:設(shè)某線性時(shí)不變系統(tǒng)的單位取樣響應(yīng)為

a為實(shí)數(shù)),分析系統(tǒng)的因果性和穩(wěn)定性。(重點(diǎn))

解:討論因果性:

因?yàn)?/p>

時(shí),

,所以該系統(tǒng)是因果系統(tǒng)。

討論穩(wěn)定性:

∴ 當(dāng)

時(shí),系統(tǒng)是穩(wěn)定的;否則,系統(tǒng)不穩(wěn)定。

例:設(shè)某線性時(shí)不變系統(tǒng)的單位取樣響應(yīng)為

a為實(shí)數(shù)),分析系統(tǒng)的因果性和穩(wěn)定性。(重點(diǎn))

解:討論因果性:

因?yàn)?/p>

時(shí),

,所以該系統(tǒng)是非因果系統(tǒng)。

討論穩(wěn)定性:

∴ 當(dāng)時(shí),系統(tǒng)是穩(wěn)定的;否則,系統(tǒng)不穩(wěn)定。

1.3 線性常系數(shù)差分方程

1 差分方程定義

卷積和是一種LTI 系統(tǒng)的數(shù)學(xué)模型,一般情況下,我們可以用差分方程描述LTI系統(tǒng)的輸入輸出關(guān)系。

差分方程給出了系統(tǒng)響應(yīng)y[n]的內(nèi)部關(guān)系。為得到y(tǒng)[n]的顯式解,必須求解方程。

2差分方程求解

1經(jīng)典法 2遞推法 3變換域法(參見下章z域變換)(重點(diǎn))

例:設(shè)系統(tǒng)的差分方程為

,輸入序列為

,求輸出序列

解:一階差分方程需一個(gè)初始條件。

設(shè)初始條件為:






設(shè)初始條件改為:






該例表明,對(duì)于同一個(gè)差分方程和同一個(gè)輸入信號(hào),因?yàn)槌跏紬l件不同,得到的輸出信號(hào)是不相同的。

幾點(diǎn)結(jié)論(重點(diǎn))

(1)對(duì)于實(shí)際系統(tǒng),用遞推解法求解,總是由初始條件向n>0的方向遞推,是一個(gè)因果解。但對(duì)于差分方程,其本身也可以向n<0的方向遞推,得到的是非因果解。因此差分方程本身不能確定該系統(tǒng)是因果系統(tǒng)還是非因果系統(tǒng),還需要用初始條件進(jìn)行限制。

(2)一個(gè)線性常系數(shù)差分方程描述的系統(tǒng)不一定是線性非時(shí)變系統(tǒng),這和系統(tǒng)的初始狀態(tài)有關(guān)。如果系統(tǒng)是因果的,一般在輸入x(n)=0(n<n0)時(shí),則輸出y(n)=0(n<n0),系統(tǒng)是線性非時(shí)變系統(tǒng)。

1.4 模擬信號(hào)數(shù)字處理方法

1 模擬信號(hào)數(shù)字處理框圖

:模擬信號(hào)輸入

預(yù)濾波:目的是限制帶寬(一般使用低通濾波器)

1采樣:將信號(hào)在時(shí)間上離散化

A/DC:模/數(shù)轉(zhuǎn)換

2量化:將信號(hào)在幅度上離散化(量化中幅度值=采樣幅度值)

3編碼:將幅度值表示成二進(jìn)制位(條件)

數(shù)字信號(hào)處理:對(duì)信號(hào)進(jìn)行運(yùn)算處理

D/AC:數(shù)/模轉(zhuǎn)換(一般用采樣保持電路實(shí)現(xiàn):臺(tái)階狀連續(xù)時(shí)間信號(hào)在采樣時(shí)刻幅度發(fā)生跳變 )

平滑濾波:濾除信號(hào)中高頻成分(低通濾波器),使信號(hào)變得平滑

:輸入信號(hào)經(jīng)過(guò)處理后的輸出信號(hào)

2.連續(xù)信號(hào)的采樣

對(duì)連續(xù)信號(hào)進(jìn)行理想采樣,設(shè)采樣脈沖,則采樣輸出

在討論理想采樣后,信號(hào)頻譜發(fā)生的變化時(shí),可遵循下面的思路:

1)由;2)由;

3)根據(jù)頻域卷積定理,由計(jì)算出。

計(jì)算過(guò)程:

1)

2)周期信號(hào)可以用傅里葉級(jí)數(shù)展開,因此

其中系數(shù)

所以

其傅里葉變換

3)

因此,采樣后信號(hào)頻譜產(chǎn)生周期延拓,周期為Ωs,同時(shí)幅度為原來(lái)的1/T倍。這是一個(gè)非常重要的性質(zhì),應(yīng)熟練掌握。

3 時(shí)域抽樣定理(重點(diǎn))

一個(gè)限帶模擬信號(hào),若其頻譜的最高頻率為,對(duì)它進(jìn)行等間隔抽樣而得,抽樣周期為T,或抽樣頻率為;只有在抽樣頻率時(shí),才可由準(zhǔn)確恢復(fù)。

例:有一連續(xù)信號(hào)式中,(1)求出的周期。

(2)用采樣間隔對(duì)進(jìn)行采樣,試寫出采樣信號(hào)的表達(dá)式。

(3)求出對(duì)應(yīng)的時(shí)域離散信號(hào)(序列) ,并求出的周期。

解:(1)周期為

(2)

(3)x(n)的數(shù)字頻率ω=0.8π,故,因而周期N=5,所以x(n)=cos(0.8πn+π/2)

簡(jiǎn)答題:(重點(diǎn))

1. 是不是任意連續(xù)信號(hào)離散后,都可從離散化后的信號(hào)恢復(fù)出原來(lái)的信號(hào)?為什么?

2. 一個(gè)連續(xù)時(shí)間信號(hào)經(jīng)過(guò)理想采樣以后,其頻譜會(huì)產(chǎn)生怎樣的變化?在什么條件下,頻譜不會(huì)產(chǎn)生失真?

3. 說(shuō)明時(shí)域采樣定理的要點(diǎn)?

4. 離散信號(hào)頻譜函數(shù)的一般特點(diǎn)是什么?

5. 畫出模擬信號(hào)數(shù)字處理框圖。并說(shuō)明各部分的作用。

名詞解釋:(重點(diǎn))

1. 時(shí)域采樣定理

2. 線性系統(tǒng)、時(shí)不變系統(tǒng)、穩(wěn)定系統(tǒng)、因果系統(tǒng)


第二章:本章涉及信號(hào)及系統(tǒng)的頻域分析方法,概念較多,但很基礎(chǔ),學(xué)習(xí)時(shí)要注意。

2.1 序列的傅里葉變換的定義及性質(zhì)

1.定義

DTFT 是一個(gè)用來(lái)確定離散時(shí)間序列頻譜的重要數(shù)學(xué)工具。

物理意義:傅里葉變換是將對(duì)信號(hào)的時(shí)域分析轉(zhuǎn)換為對(duì)其在頻域的分析,便于研究問(wèn)題。

若序列滿足絕對(duì)可和條件

則其離散時(shí)間傅里葉變換(DiscreteTimeFourierTransform-DTFT:非周期序列的傅里葉變換)定義為

------(記住!!)

反變換定義為

------

傅里葉變換對(duì)

例:設(shè)

,求其序列傅里葉變換。(重點(diǎn))


當(dāng)

時(shí)


(2-5)

的幅度和相位隨

變化曲線如圖2.1所示。




圖2.1 R4(n)的幅度與相位曲線

例:試求如下序列的傅里葉變換:(重點(diǎn))

(1)

(2)

(3)

(4)

解:

(1)

(2)

(3) ,

(4)

=

2.性質(zhì)

1)周期性(重點(diǎn)): DTFT是關(guān)于ω的周期為2π的周期函數(shù)

2)線性(重點(diǎn)):設(shè)

,

,那么

3)時(shí)移特性(重點(diǎn))

4)頻移特性

5)時(shí)域卷積定理(重點(diǎn))

6)頻域卷積定理

7)帕斯瓦爾定理

時(shí)域總能量等于頻域一周期內(nèi)總能量。

7) 幅度頻譜為ω的偶函數(shù),相位頻譜為ω的奇函數(shù)。

8) X(ejω)的實(shí)部為ω的偶函數(shù), X(ejω) 的虛部為ω的奇函數(shù)。

對(duì)稱關(guān)系的總結(jié)(重點(diǎn))

如果x[n]為復(fù)數(shù)序列,其DTFT為 X(ejω),

(a) x[n]實(shí)部的DTFT為X(ejω)的共軛對(duì)稱部分-----------

(b) x[n]虛部的DTFT 為X(ejω)的反共軛對(duì)稱部分-----------

(c) x[n]的共軛對(duì)稱部分的DTFT為 X(ejω)的實(shí)部-----------

(d) x[n] 的反共軛對(duì)稱部分的DTFT 為X(ejω)的虛部-----------

如果實(shí)序列x[n] 的 DTFT 為X(ejω),

(e) x[n]的偶對(duì)稱部分的DTFT為X(ejω) 的實(shí)部,

-----------

(f) x[n]的奇對(duì)稱部分的DTFT為 X(ejω) 的虛部

-----------

例:設(shè)系統(tǒng)的單位取樣響應(yīng),輸入序列為,完成下面各題:(1)求出系統(tǒng)輸出序列;(2)分別求出、和的傅里葉變換。(重點(diǎn))

:(1)

(2)

2.2 時(shí)域離散信號(hào)的傅里葉變換與模擬信號(hào)傅里葉變換之間的關(guān)系:

式中

2.3 序列的Z變換

1 Z變換定義(重點(diǎn))

Z變換為離散時(shí)間信號(hào)與LTI系統(tǒng)分析的重要數(shù)學(xué)工具。給定一離散時(shí)間序列x(n),其z變換定義為: ------(記住!!)

其中,,。z變換存在情況下的Z變量取值范圍稱為收斂域(ROC)。

注意:Z變換+不同收斂域?qū)?yīng)不同收斂域的不同序列

序列(Z變換+收斂域)(重點(diǎn))

例:求以下序列的Z變換及收斂域:(重點(diǎn))

(1);

(2);

(3)

:(1)

(2)

(3)

[說(shuō)明]上題也可以改為求序列的傅立葉變換??梢岳谩?/strong>

2 Z變換和DTFT之間的關(guān)系(重點(diǎn))

DTFT 為單位圓上的z變換。數(shù)學(xué)表達(dá)為: ------ 記住并理解!

3. 序列特性與X(z)的收斂域ROC的關(guān)系。(重點(diǎn))

收斂區(qū)域要依據(jù)序列的性質(zhì)而定。同時(shí),也只有Z變換的收斂區(qū)域確定之后,才能由Z變換唯一地確定序列。

一般來(lái)來(lái)說(shuō),序列的Z變換的收斂域在Z平面上的一環(huán)狀區(qū)域:

總結(jié):a. ROC不包含任何極點(diǎn)。

b.有理 z變換的收斂域ROC由其極點(diǎn)界定。

c. 對(duì)于有限長(zhǎng)序列x[n],其z變換的收斂域ROC 為整個(gè)z-平面,可能在 z = 0 或z = ∞除外。只有序列為時(shí),收斂域是整個(gè)Z平面。

d. 對(duì)于右邊序列x[n],其 z變換的收斂域ROC由其離原點(diǎn)最遠(yuǎn)的極點(diǎn)確定,其形式為。

e. 對(duì)于左邊序列x[n], 其 z變換的收斂域ROC由其離原點(diǎn)最近的極點(diǎn)確定,其形式為。

f. 對(duì)于雙邊序列x[n], 其 z變換的收斂域ROC環(huán)狀收斂域,,其形式為公共收斂域。

4. Z反變換(重點(diǎn))

常用序列的Z變換(重點(diǎn)--記住!!)

逆變換

x,C:收斂域內(nèi)繞原點(diǎn)逆時(shí)針的一條閉合曲線

留數(shù)定理:

留數(shù)輔助定理:

利用部分分式展開:,然后利用定義域及常用序列的Z變換求解。(重點(diǎn))

基本要求:用部分分式展開法求z反變換。(重點(diǎn))

例:假設(shè) ,收斂域ROC 為,則 的z反變換為( )。(重點(diǎn))

說(shuō)明:本題要求掌握序列的時(shí)域特性域z變換收斂域之間的對(duì)應(yīng)關(guān)系。具體說(shuō),有限長(zhǎng)序列的z變換的ROC是怎樣的,右邊序列的z變換的ROC是怎樣的,因果序列的z變換的ROC是怎樣的,左邊序列的z變換的ROC是怎樣的,反因果序列的z變換的ROC是怎樣的。

典型序列的z變換表達(dá)式是否記住了?這兩個(gè)典型z變換對(duì),對(duì)求z變換或逆z變換非常重要。

例:已知

,試求與

對(duì)應(yīng)的所有可能的序列

。(重點(diǎn))

解:同一個(gè)Z變換函數(shù),收斂域不同,對(duì)應(yīng)的序列也不同。本題沒(méi)有給定收斂域,所以必須先確定收斂域。

有兩個(gè)極點(diǎn):,,因?yàn)槭諗坑蚩偸且詷O點(diǎn)為邊界,所以收斂域有以下三種情況:,,,三種收斂域?qū)?yīng)三種不同的原序列,分別討論如下:

(1) 對(duì)應(yīng)左邊序列 ∴

(2) 對(duì)應(yīng)雙邊序列 ∴

(3) 對(duì)應(yīng)右邊序列 ∴

例:設(shè)


,用部分分式展開法求逆Z變換。(重點(diǎn))


解:先去掉z的負(fù)冪次,以便于求解,將的分子分母同乘以,得:

將等式兩端同時(shí)除以z,得:

因而得:

由收斂域知,為右邊序列,得:

主要應(yīng)用于單階極點(diǎn)的序列。

5 Z變換的性質(zhì)

1線性性質(zhì)(重點(diǎn))

2序列的移位性質(zhì)(重點(diǎn))

3序列乘以指數(shù)序列的性質(zhì)(重點(diǎn))

4序列乘以n的ZT

5復(fù)共軛序列的ZT

6初值定理

7終值定理

8時(shí)域卷積定理(重點(diǎn))

設(shè)

9復(fù)卷積定理

10帕斯維爾定理

,

那么

2.4 離散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)及頻率響應(yīng)

1 系統(tǒng)函數(shù)定義(重點(diǎn))

一個(gè)線性時(shí)不變離散時(shí)間系統(tǒng)在時(shí)域中可以用它的單位取樣響應(yīng)來(lái)表征,即: 對(duì)等式兩邊取Z變換并根據(jù)時(shí)域卷積定理,有:

則: 一般稱為系統(tǒng)的系統(tǒng)函數(shù)(系統(tǒng)零狀態(tài)響應(yīng)的Z變換與輸入的Z變換之比),它表征了系統(tǒng)的復(fù)頻域特性。

2 系統(tǒng)函數(shù)與差分方程的關(guān)系

(給定差分方程,能計(jì)算其系統(tǒng)函數(shù),或給定系統(tǒng)函數(shù),能計(jì)算得到差分方程。) (重點(diǎn))

3 頻率響應(yīng)(重點(diǎn))

頻率響應(yīng)是一個(gè)重要的概念,根據(jù)頻率響應(yīng),可理解濾波。

頻率響應(yīng)定義為系統(tǒng)單位沖激響應(yīng)的DTFT:

(重點(diǎn))

其中, |H(ejω)| 稱為幅頻響應(yīng), 稱為相頻響應(yīng)。系統(tǒng)的頻率響應(yīng)是以2π為周期的ω的連續(xù)函數(shù),這一點(diǎn)和連續(xù)系統(tǒng)的頻率響應(yīng)是不同的,學(xué)習(xí)時(shí)應(yīng)加以注意。若h(n)為實(shí)數(shù),則系統(tǒng)的幅度響應(yīng)在區(qū)間內(nèi)是偶對(duì)稱的,而相位響應(yīng)是奇對(duì)稱的。

注意僅當(dāng)穩(wěn)定系統(tǒng)才有頻率響應(yīng)。頻率響應(yīng)H(ejω)可根據(jù)DTFT 與z變換之間的關(guān)系簡(jiǎn)單得到:



穩(wěn)態(tài)響應(yīng)的求解

結(jié)論:

對(duì)于LTI 系統(tǒng),如果輸入為正弦序列x(n)=cos(ω0t+φ0), 則輸出響應(yīng)y(n)必為相同形式的正弦序列,但需在 ω=ω0的幅頻響應(yīng)|H(ejω)|進(jìn)行加權(quán),并通過(guò)相頻響應(yīng)

在 ω=ω0的值進(jìn)行移位,即:y[n]= |H(ejω0)|cos(ω0t+φ0+

)

例:假設(shè)實(shí)序列x[n]的DTFT 記為, 則其幅值 是關(guān)于ω的(偶函數(shù))。

說(shuō)明:還記得反復(fù)強(qiáng)調(diào)的一句話,實(shí)序列的DTFT的幅度、實(shí)部是關(guān)于頻率ω偶函數(shù),而相位和虛部則是關(guān)于頻率ω奇函數(shù)。

例:對(duì)于一LTI 離散時(shí)間系統(tǒng)其頻率響應(yīng),如果系統(tǒng)輸x(n) =, 響應(yīng)的穩(wěn)態(tài)輸出響應(yīng)y(n) = ( )。

說(shuō)明:將系統(tǒng)的頻率響應(yīng)寫成幅度相位表達(dá)式:,則輸出信號(hào)為:。這里由于給出了的具體表達(dá)式,所以需要分別計(jì)算出和之值。

4 用系統(tǒng)函數(shù)極點(diǎn)分布分析系統(tǒng)的因果性和穩(wěn)定性(重點(diǎn))

系統(tǒng)函數(shù):傳輸函數(shù)H(z) 為系統(tǒng)的單位沖激響應(yīng)h(n)的Z變換。

1)穩(wěn)定系統(tǒng):有界的輸入產(chǎn)生的輸出也有界的系統(tǒng),即:若,則

線性移不變系統(tǒng)是穩(wěn)定系統(tǒng)的充要條件:

或:其系統(tǒng)函數(shù)H(z)的收斂域包含單位圓 |z|=1(牢記此結(jié)論!)

2)因果系統(tǒng):時(shí)刻的輸出只由時(shí)刻之前的輸入決定

線性移不變系統(tǒng)是因果系統(tǒng)的充要條件:

或:其系統(tǒng)函數(shù)H(z)的收斂域在某圓外部:即:|z|>Rx(牢記此結(jié)論!)

3)穩(wěn)定因果系統(tǒng):同時(shí)滿足上述兩個(gè)條件的系統(tǒng)。

線性移不變系統(tǒng)是因果穩(wěn)定系統(tǒng)的充要條件:,

或:H(z)的極點(diǎn)在單位圓內(nèi)H(z)的收斂域滿足:(牢記此結(jié)論!)

例:.一因果LTI 離散時(shí)間系統(tǒng)的傳輸函數(shù), 則系統(tǒng)的單位沖激響應(yīng)為( 0.5nu(n) )。

說(shuō)明:根據(jù)傳遞函數(shù)求系統(tǒng)的單位沖激響應(yīng),其實(shí)就是將傳遞函數(shù)進(jìn)行逆z變換,但要注意系統(tǒng)的因果性如何。

例:因果IIR 離散時(shí)間LTI 系統(tǒng),其傳輸函數(shù),則系統(tǒng)( 穩(wěn)定)。

例:一FIR離散時(shí)間 LTI 系統(tǒng)總是( 穩(wěn)定)。

說(shuō)明:系統(tǒng)的穩(wěn)定性如何判斷?按照教材中的說(shuō)法,就是系統(tǒng)傳遞函數(shù)的收斂域如果包括“單位圓”,則系統(tǒng)是穩(wěn)定的。如果你熟悉了序列的z變換的ROC的性質(zhì),則此題不難回答。對(duì)于因果系統(tǒng)來(lái)說(shuō),其單位沖激響應(yīng)為因果序列,故其z變換的ROC一定是某圓外部的整個(gè)區(qū)域。而這個(gè)圓就位于離原點(diǎn)最遠(yuǎn)的極點(diǎn)上,所以,對(duì)于因果系統(tǒng),如果系統(tǒng)傳遞函數(shù)的全部極點(diǎn)都位于單位圓以內(nèi)的話,則系統(tǒng)是穩(wěn)定的。

對(duì)于FIR系統(tǒng),其單位沖激響應(yīng)是一個(gè)有限長(zhǎng)序列,其z變換的ROC為除了無(wú)窮遠(yuǎn)和原點(diǎn)之外的整個(gè)z平面,自然包括單位圓,所以FIR系統(tǒng)始終是穩(wěn)定的。

5 系統(tǒng)的頻率特性可由系統(tǒng)函數(shù)零點(diǎn)及極點(diǎn)確定

(式中,zk是極點(diǎn),zi是零點(diǎn);在極點(diǎn)處,序列x(n)的Z變換是不收斂的,因此收斂區(qū)域內(nèi)不應(yīng)包括極點(diǎn)。)

系統(tǒng)函數(shù)H(z)的極點(diǎn)位置主要影響頻響的峰值位置及尖銳程度,零點(diǎn)位置主要影響頻響的谷點(diǎn)位置及形狀。

例:設(shè)一階系統(tǒng)的差分方程為,,用幾何法分析其幅頻特性。(重點(diǎn))

解:對(duì)差分方程兩邊取Z變換,得:

系統(tǒng)函數(shù)為: ,極點(diǎn)為,零點(diǎn)為,如下圖左所示:






當(dāng)

時(shí),由于極點(diǎn)矢量長(zhǎng)度最短,幅頻特性出現(xiàn)峰值,隨著

的增加,幅度逐漸減小,當(dāng)

時(shí),由于極點(diǎn)矢量長(zhǎng)度最長(zhǎng),幅頻特性出現(xiàn)谷值,隨著

的增加,幅度逐漸增大,直到

時(shí),幅頻特性出現(xiàn)峰值,如上圖右所示。

簡(jiǎn)答題:(重點(diǎn))

1. 說(shuō)明有限長(zhǎng)序列、左邊序列、右邊序列、雙邊序列的概念和收斂域各是什么?

2. 說(shuō)明系統(tǒng)頻率響應(yīng)的概念?系統(tǒng)的頻率響應(yīng)和系統(tǒng)函數(shù)是什么關(guān)系?(單位圓上(

)的系統(tǒng)函數(shù)就是系統(tǒng)的頻率響應(yīng))

3. 說(shuō)明FIR系統(tǒng)為什么始終是穩(wěn)定的?

4. 怎樣在z域表示離散時(shí)間LTI 系統(tǒng)? 答案:傳輸函數(shù)H(z)表示離散時(shí)間LTI 系統(tǒng)。


第三章:DFT是為適應(yīng)計(jì)算機(jī)分析傅里葉變換規(guī)定的一種專門運(yùn)算,本章是數(shù)字信號(hào)處理課程的重點(diǎn)章節(jié)。

前言

信號(hào)處理中會(huì)遇到幾種信號(hào)形式:(1)連續(xù)周期信號(hào)(2)連續(xù)非周期信號(hào)(3)離散非周期信號(hào)(4)離散周期信號(hào)(重點(diǎn))

各種信號(hào)在時(shí)域和頻域之間總的來(lái)說(shuō)都是傅里葉變換,但具體形式及應(yīng)用是不同的。

1.連續(xù)周期信號(hào) —— 傅里葉級(jí)數(shù)(FS)

連續(xù)周期信號(hào)

可展開成傅里葉級(jí)數(shù):


(*)

式中,

,

的周期。

傅里葉級(jí)數(shù)的系數(shù)為:

幅度頻譜是指各次諧波的振幅隨頻率的變化關(guān)系,即:




2.連續(xù)非周期信號(hào) —— 傅里葉變換(FT)

連續(xù)非周期信號(hào)

的傅里葉變換為:


因?yàn)榉侵芷诳梢暈?/p>

,則離散頻譜間距

,則

變成

的連續(xù)函數(shù)。


3.離散非周期信號(hào) —— 序列的傅里葉變換(DTFT)


如果把序列看成連續(xù)時(shí)間信號(hào)的采樣,采樣間隔為

,則數(shù)字頻率

和模擬角頻率

的關(guān)系為

,且

,代入上式,得:



4.離散周期信號(hào) —— 離散傅里葉級(jí)數(shù)(DFS)

設(shè)

是周期為

的周期序列,即:


為任意整數(shù)

表3.1 四種傅里葉變換形式的歸納


一般規(guī)律:一個(gè)域的離散對(duì)應(yīng)另一個(gè)域的周期延拓, 一個(gè)域的連續(xù)必定對(duì)應(yīng)另一個(gè)域的非周期。(重點(diǎn))

3.1 離散傅里葉級(jí)數(shù)

1.周期序列的離散傅里葉級(jí)數(shù)(DFS)

說(shuō)明:離散傅里葉級(jí)數(shù)系數(shù),用DFS(Discrete Fourier Series)表示。

連續(xù)時(shí)間周期信號(hào)可以用傅里葉級(jí)數(shù)表示,離散周期序列也可以表示成傅里葉級(jí)數(shù)形式。

周期為N的復(fù)指數(shù)序列的基頻序列為

k次諧波序列為

由于,即,因而,離散傅里葉級(jí)數(shù)的所有諧波成分中只有N個(gè)是獨(dú)立的。因此在展開成離散傅里葉級(jí)數(shù)時(shí),我們只能取N個(gè)獨(dú)立的諧波分量,通常取k=0到(N-1),即

(*)

式中,1/N是習(xí)慣上采用的常數(shù),是k次諧波的系數(shù)。利用

將(*)式兩端同乘以,并對(duì)一個(gè)周期求和

由于

所以也是一個(gè)以N為周期的周期序列。因此,時(shí)域離散周期序列的離散傅里葉級(jí)數(shù)在頻域上仍然是一個(gè)周期序列。稱為離散傅里葉級(jí)數(shù)系數(shù),用DFS(Discrete Fourier Series)表示。

令,則

其中,符號(hào)DFS[.]表示離散傅里葉級(jí)數(shù)正變換,IDFS[.]表示離散傅里葉級(jí)數(shù)反變換。

例:設(shè),將以為周期進(jìn)行周期延拓,得到周期序列,求的DFS。

解:



其幅度特性為:


2.周期序列的傅里葉變換

思路:由

利用和DTFT的頻移特性,可得

傅里葉變換時(shí)域、頻域?qū)?yīng)關(guān)系:

根據(jù)序列的傅里葉變換和離散傅里葉級(jí)數(shù)頻域特性,再結(jié)合連續(xù)時(shí)間信號(hào)的傅里葉變換頻域特性,我們可以得出傅里葉變換時(shí)、頻域的一般對(duì)應(yīng)關(guān)系:連續(xù)→非周期,離散→周期這種對(duì)應(yīng)關(guān)系很重要,要求熟記(重點(diǎn))。

3.2 有限長(zhǎng)序列的離散傅立葉變換(DFT)

說(shuō)明:(Discrete Fourier Transform,DFT離散傅里葉變換)

1 定義(重點(diǎn))

,0≤

------(記住!!)


,0≤n≤

------記住!

其中,

應(yīng)當(dāng)注意,雖然和都是長(zhǎng)度為得有限長(zhǎng)序列,但他們分別是由周期序列和截取其主周期得到的,周期為的周期序列可以看成長(zhǎng)度為的有限長(zhǎng)序列周期延拓的結(jié)果。本質(zhì)上是做DFS或IDFS,所以不能忘記它們的隱含周期性。尤其是涉及其位移特性時(shí)更要注意。(重點(diǎn))

DFT的隱含周期性:(重點(diǎn))

例:設(shè)

,求

的4點(diǎn)DFT。(重點(diǎn))

解:


的4點(diǎn)離散傅里葉變換為:



為周期將

延拓成周期序列,得


其離散傅里葉級(jí)數(shù)為:



例:設(shè)

,求

的8點(diǎn)DFT。(重點(diǎn))

解:


的8點(diǎn)離散傅里葉變換為:



為周期將

延拓成周期序列,得


其離散傅里葉級(jí)數(shù)為:




由例可見,離散傅里葉變換的結(jié)果與變換區(qū)間長(zhǎng)度

的取值有關(guān)。

2 離散傅立葉變換與DTFT、Z變換的關(guān)系(重點(diǎn))

DFT的物理意義:X(k)為x(n)的傅里葉變換在區(qū)間上的等間隔采樣。為在Z平面單位圓上的點(diǎn)等間隔采樣。

3 時(shí)域分析

記住結(jié)論:時(shí)域抽樣對(duì)應(yīng)頻域的周期拓展,頻率抽樣對(duì)應(yīng)時(shí)域的以周期N的周期拓展。

這可以表述為如下公式:


3.3 離散傅里葉變換的基本性質(zhì)

1 線性性質(zhì)

2 循環(huán)移位性質(zhì)

設(shè)是長(zhǎng)度為的有限長(zhǎng)序列,則的點(diǎn)循環(huán)移位定義為():

循環(huán)移位的實(shí)現(xiàn)步驟:

3 循環(huán)卷積定理(重點(diǎn))

1)設(shè)序列h(n)和x(n)的長(zhǎng)度分別為NM。h(n)與x(n)的L點(diǎn)循環(huán)卷積定義為

式中,L稱為循環(huán)卷積區(qū)間長(zhǎng)度,L≥max[N,M]。

2) 循環(huán)卷積矩陣






特點(diǎn):

(1)第1行是序列{x(0),x(1),…,x(L-1)}的循環(huán)倒相序列。注意,如果x(n)的長(zhǎng)度M<L,則需要在x(n)末尾補(bǔ)L-M個(gè)零后,再形成第一行的循環(huán)倒相序列。

(2)第1行以后的各行均是前一行向右循環(huán)移1位形成的。

(3)矩陣的各主對(duì)角線上的序列值均相等。

循環(huán)卷積和線性卷積的區(qū)別

線性卷積:翻折—>乘加—>移位 :y(n)=x(n)*h(n)=∑h(k)x(n-k)

循環(huán)卷積:補(bǔ)零—>周期延拓—>翻折—>循環(huán)移位—>對(duì)應(yīng)值相加

例:計(jì)算下面給出的兩個(gè)長(zhǎng)度為4的序列h(n)與x(n)的4點(diǎn)和8點(diǎn)循環(huán)卷積。(重點(diǎn))

解:按照循環(huán)卷積矩陣寫出h(n)與x(n)的4點(diǎn)循環(huán)卷積矩陣形式為

h(n)與x(n)的8點(diǎn)循環(huán)卷積矩陣形式為

【補(bǔ)充】①計(jì)算h(n)與x(n)的線性卷積?②哪一種情況下計(jì)算的循環(huán)卷積結(jié)果就等于線性卷積?

【說(shuō)明】當(dāng)循環(huán)卷積區(qū)間長(zhǎng)度L大于等于y(n) = h(n)*x(n)的長(zhǎng)度時(shí),循環(huán)卷積結(jié)果就等于線性卷積。假設(shè)h(n)和x(n)都是有限長(zhǎng)序列,長(zhǎng)度分別是NM。循環(huán)卷積等于線性卷積的條件是LNM-1。(重點(diǎn))

3) 時(shí)域循環(huán)卷積定理

設(shè)h(n)和x(n)的長(zhǎng)度分別為N和M,其L點(diǎn)循環(huán)卷積為

L

則由DFT的循環(huán)卷積定理有

4 復(fù)共軛序列的DFT

性質(zhì):設(shè)是x(n)的復(fù)共軛序列,長(zhǎng)度為N,,

例:給定一16-點(diǎn)實(shí)序列x(n), 其 16-點(diǎn) DFT 記為X(k), 已知X(13)= 2 + j3,則X*(3) = ( 2 + j3 )。

說(shuō)明:DFT的性質(zhì)。實(shí)序列的DFT的共軛對(duì)稱性:X(k) =X*(N-k),或X(N-k) =X*(k)。

5 DFT的共軛對(duì)稱性(重點(diǎn))

可總結(jié)出DFT的共軛對(duì)稱性質(zhì):如果序列x(n)的DFT為X(k),則x(n)的實(shí)部和虛部(包括j)的DFT分別為X(k)的共軛對(duì)稱分量和共軛反對(duì)稱分量;而x(n)的共軛對(duì)稱分量和共軛反對(duì)稱分量的DFT分別為X(k)的實(shí)部和虛部乘以j。

3.4 頻域采樣定理

離散傅里葉變換相當(dāng)于信號(hào)傅里葉變換的等間隔采樣,也就是說(shuō)實(shí)現(xiàn)了頻域的采樣,便于計(jì)算機(jī)計(jì)算。那么是否任一序列都能用頻域采樣的方法去逼近呢?這是一個(gè)很吸引人的問(wèn)題。

我們考慮一個(gè)任意的絕對(duì)可和的序列x(n),它的z變換為

如果對(duì)X(z)單位圓上進(jìn)行等距離采樣

現(xiàn)在要問(wèn),這樣采樣以后,信息有沒(méi)有損失?或者說(shuō),采樣后所獲得的有限長(zhǎng)序列xN(n)能不能代表原序列x(n)。

為了弄清這個(gè)問(wèn)題,我們從周期序列開始

由于

所以

也即是原非周期序列x(n)的周期延拓序列,其時(shí)域周期為頻域采樣點(diǎn)數(shù)N。在第一章我們看到,時(shí)域的采樣造成頻域的周期延拓,這里又對(duì)稱的看到,頻域采樣同樣造成時(shí)域的周期延拓。

因此,如果序列x(n)不是有限長(zhǎng)的,則時(shí)域周期延拓時(shí),必然造成混疊現(xiàn)象,因而一定會(huì)產(chǎn)生誤差。

對(duì)于長(zhǎng)度為M的有限長(zhǎng)序列,只有當(dāng)頻域采樣點(diǎn)數(shù)N大于或等于序列長(zhǎng)度M時(shí),才有

即可由頻域采樣值X(k)恢復(fù)出原序列x(n),否則產(chǎn)生時(shí)域混疊現(xiàn)象,這就是所謂的頻域采樣定理。(重點(diǎn))

內(nèi)插公式:

3.5 DFT的應(yīng)用舉例

1.用DFT計(jì)算線性卷積(重點(diǎn))

用循環(huán)(周期)卷積計(jì)算有限長(zhǎng)序列的線性卷積(重點(diǎn))

對(duì)周期要求:(N1、N2分別為兩個(gè)序列的長(zhǎng)度)(記住!!)

2.用DFT進(jìn)行譜分析的誤差問(wèn)題(重點(diǎn))

(1)混疊現(xiàn)象

利用DFT逼近連續(xù)時(shí)間信號(hào)的傅里葉變換,為避免混疊失真,按照抽樣定理的要求,采樣頻率至少是信號(hào)最高頻率的兩倍。

解決混疊問(wèn)題的唯一方法是保證采樣頻率足夠高。

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

任何帶限信號(hào)都是非時(shí)限的,任何時(shí)限信號(hào)都是非帶限的。實(shí)際問(wèn)題中遇到的離散時(shí)間序列可能是非時(shí)限的、無(wú)限長(zhǎng)序列,在對(duì)該序列利用DFT進(jìn)行處理時(shí),由于作DFT的點(diǎn)數(shù)總是有限的,因此就有一個(gè)必須將該序列截?cái)嗟膯?wèn)題。序列截?cái)嗟倪^(guò)程相當(dāng)于給該序列乘上一個(gè)矩形窗口函數(shù)RN(n)。如果原來(lái)序列的頻譜為,矩形窗函數(shù)的頻譜為,則截?cái)嗪笥邢揲L(zhǎng)序列的頻譜為

截?cái)嗪笮蛄械念l譜與原序列頻譜必然有差別,這種差別對(duì)譜分析的影響主要表現(xiàn)在如下兩個(gè)方面

頻譜泄露:由于矩形窗函數(shù)頻譜的引入,使卷積后的頻譜被展寬了,即的頻譜“泄露”到其它頻率處,稱為頻譜泄露。

在進(jìn)行DFT時(shí),由于取無(wú)限個(gè)數(shù)據(jù)是不可能的,所以序列的時(shí)域截?cái)嗍潜厝坏模孤妒请y以避免的。為了盡量減少泄露的影響,截?cái)鄷r(shí)要根據(jù)具體的情況,選擇適當(dāng)形狀的窗函數(shù),如漢寧窗或漢明窗等。

譜間干擾。 在主譜線兩邊形成很多旁瓣,引起不同頻率分量間的干擾(簡(jiǎn)稱譜間干擾),特別是強(qiáng)信號(hào)譜的旁瓣可能湮沒(méi)弱信號(hào)的主譜線,或者把強(qiáng)信號(hào)譜的旁瓣誤認(rèn)為是另一頻率的信號(hào)的譜線,從而造成假信號(hào),這樣就會(huì)使譜分析產(chǎn)生較大偏差。 

(3)柵欄效應(yīng)

由于DFT是有限長(zhǎng)序列的頻譜等間隔采樣所得到的樣本值,這就相當(dāng)于透過(guò)一個(gè)柵欄去觀察原來(lái)信號(hào)的頻譜,因此必然有一些地方被柵欄所遮擋,這些被遮擋的部分就是未被采樣到的部分,這種現(xiàn)象稱為柵欄效應(yīng)。由于柵欄效應(yīng)總是存在的,因而可能會(huì)使信號(hào)頻率中某些較大的頻率分量由于被“遮擋”而無(wú)法得到反映。此時(shí),通常在有限長(zhǎng)序列的尾部增補(bǔ)若干個(gè)零值,借以改變?cè)蛄械拈L(zhǎng)度。這樣對(duì)加長(zhǎng)的序列作DFT時(shí),由于點(diǎn)數(shù)增加就相當(dāng)于調(diào)整了原來(lái)柵欄的間隙,可以使原來(lái)得不到反映的那些較大的頻率分量落在采樣點(diǎn)上而得到反映。

產(chǎn)生原因說(shuō)明:由傅里葉變換理論知道,若信號(hào)持續(xù)時(shí)間有限長(zhǎng),則其頻譜無(wú)限寬;若信號(hào)的頻譜有限寬,則其持續(xù)時(shí)間必然為無(wú)限長(zhǎng)。(重點(diǎn))

如果用DFT分析連續(xù)信號(hào)的頻譜,在對(duì)連續(xù)信號(hào)采樣時(shí),無(wú)法滿足采樣定理,那么就會(huì)出現(xiàn)頻譜混疊現(xiàn)象。解決混疊問(wèn)題的唯一方法是保證采樣頻率足夠高。當(dāng)連續(xù)信號(hào)無(wú)限長(zhǎng)或很長(zhǎng)時(shí),在對(duì)連續(xù)信號(hào)采樣時(shí),采樣點(diǎn)數(shù)太多以致無(wú)法存儲(chǔ)和計(jì)算,需要將信號(hào)截?cái)?/strong>,這樣將導(dǎo)致頻譜的泄漏現(xiàn)象。為了盡量減少泄露的影響,截?cái)鄷r(shí)要根據(jù)具體的情況,選擇適當(dāng)形狀的窗函數(shù),如漢寧窗或漢明窗等。用DFT計(jì)算連續(xù)信號(hào)的頻譜只能得到采樣點(diǎn)上的頻譜,而不能看到整個(gè)頻譜,這種現(xiàn)象稱為柵欄效應(yīng)??梢酝ㄟ^(guò)增加點(diǎn)數(shù),因?yàn)辄c(diǎn)數(shù)增加就相當(dāng)于調(diào)整了原來(lái)柵欄的間隙,可以使原來(lái)得不到反映的那些較大的頻率分量落在采樣點(diǎn)上而得到反映。

3.用DFT進(jìn)行譜分析的參數(shù)選擇問(wèn)題(重點(diǎn))

對(duì)模擬信號(hào)頻譜的采樣間隔,稱之為頻率分辨率。

(1)在已知信號(hào)的最高頻率fc(即譜分析范圍)時(shí),為了避免頻率混疊現(xiàn)象,要求采樣頻率Fs滿足:Fs>2fc。

(2)采樣頻率Fs,采樣點(diǎn)數(shù)N,譜分辨率F=Fs/N,如果保持采樣點(diǎn)數(shù)N不變,要提高頻譜分辨率(減小F),就必須降低采樣頻率,采樣頻率的降低會(huì)引起譜分析范圍變窄和頻譜混疊失真。如維持Fs不變,為提高頻率分辨率可以增加采樣點(diǎn)數(shù)N。因?yàn)镹T=Tp,T=Fs-1,只有增加對(duì)信號(hào)的觀察時(shí)間Tp,才能增加N。

(3)采樣點(diǎn)數(shù)N>2fc/F

(4)最小記錄時(shí)間Tp≥1/F

例:用DFT對(duì)實(shí)信號(hào)進(jìn)行譜分析,要求頻率分辨率

,信號(hào)最高頻率為

,試確定以下參數(shù):(1)最小記錄時(shí)間;(2)最大取樣間隔;(3)最少采樣點(diǎn)數(shù);(4)若要求頻率分辨率提高一倍,求最少采樣點(diǎn)數(shù)。(重點(diǎn))

解 (1)

(2)

(3)

(4)


簡(jiǎn)答題:(重點(diǎn))

1. 一個(gè)序列的DFT與序列的傅里葉變換之間的關(guān)系是什么?

2. 序列的DTFT和序列的z變換間的關(guān)系是什么?序列的DFT和序列的Z變換間的關(guān)系是什么?

3. 有限長(zhǎng)序列

的長(zhǎng)度為M,對(duì)其進(jìn)行頻域采樣,不失真的條件是什么??jī)蓚€(gè)有限長(zhǎng)序列

,對(duì)它們進(jìn)行線性卷積,結(jié)果用

表示,

的長(zhǎng)度是多少?如果進(jìn)行循環(huán)卷積,那么什么時(shí)候線性卷積和循環(huán)卷積的結(jié)果相等?

4. 用DFT進(jìn)行譜分析帶來(lái)哪些誤差問(wèn)題?采取什么措施可以減少這些誤差?

5. 時(shí)域采樣定理的要點(diǎn)是什么?頻域采樣定理的要點(diǎn)是什么?


第四章:快速傅里葉變換并不是一種新的變換,它是離散傅里葉變換的一種快速算法。

4.1 直接計(jì)算DFT的問(wèn)題及改進(jìn)的途徑

直接計(jì)算DFT,需要

次復(fù)數(shù)乘法,

次復(fù)數(shù)加法。

直接計(jì)算離散傅里葉變換,由于計(jì)算量近似正比于N2,顯然對(duì)于很大的N值,直接計(jì)算離散傅里葉變換要求的算術(shù)運(yùn)算量非常大。(重點(diǎn))我們可以利用系數(shù)WNnk的特性來(lái)改善離散傅里葉變換的計(jì)算效率。

(1)

的對(duì)稱性

(2)

的周期性

利用

的對(duì)稱性和周期性,將大點(diǎn)數(shù)的DFT分解成若干個(gè)小點(diǎn)數(shù)的DFT,F(xiàn)FT正是基于這個(gè)基本思路發(fā)展起來(lái)的。(重點(diǎn))

說(shuō)明:快速傅里葉變換FFT(Fast Fourier Transform)

分類:按時(shí)間抽?。―IT)算法和按頻率抽取(DIF)算法。

4.2 基2FFT的算法原理和FFT運(yùn)算特點(diǎn)

1)數(shù)據(jù)要求:

2)計(jì)算效率(乘法運(yùn)算次數(shù):,加法計(jì)算次數(shù):NM )(復(fù)數(shù)運(yùn)算)

(DFT運(yùn)算:乘法運(yùn)算次數(shù):,加法計(jì)算次數(shù):)(復(fù)數(shù)運(yùn)算)

對(duì)于算法原理,要求能夠看懂分解流圖。

1 時(shí)域抽取法如下:(Decimation In Time, DIT–FFT)

設(shè)序列x(n)長(zhǎng)度為N,且滿足N=2M,M為正整數(shù)。按n的奇偶把x(n)分解為兩個(gè)N/2點(diǎn)的子序列:

則x(n)的DFT為


所以

X(k)又可以寫為

上式將N點(diǎn)DFT分解為兩個(gè)N/2點(diǎn)的DFT運(yùn)算,運(yùn)算過(guò)程如下圖示

利用蝶形運(yùn)算求解。

運(yùn)算量:由按時(shí)間抽取的FFT流圖可見:

每級(jí)都由個(gè)蝶形單元構(gòu)成,因此每級(jí)都需要次復(fù)數(shù)乘法和次復(fù)數(shù)加法。這樣,級(jí)運(yùn)算共需要:

復(fù)數(shù)乘法

復(fù)數(shù)加法

直接計(jì)算DFT需要

復(fù)數(shù)乘法:

復(fù)數(shù)加法:

以乘法為例,對(duì)FFT算法與直接DFT算法的運(yùn)算量進(jìn)行比較:





8

64

12

5.4



128

16384

448

36.6



1024

1048576

5120

204.8



可以看出:當(dāng)越大時(shí),F(xiàn)FT算法的優(yōu)越性越突出。


DIT-FFT算法與DFT運(yùn)算量的比較

直接計(jì)算DFT與FFT算法的計(jì)算量之比為

N越大,F(xiàn)FT的優(yōu)點(diǎn)越為明顯

說(shuō)明:掌握給定點(diǎn)數(shù)的基2 DIT-FFT蝶形圖

8點(diǎn)DFT的完整FFT流圖:


2頻域抽樣法(Decimation In Frequency, DIF–FFT)

將長(zhǎng)度為N=2M的序列x(n)前后對(duì)半分開, 其N點(diǎn)DFT可表示為

按k的奇偶可將X(k)分為兩部分

k取偶數(shù)時(shí)

k取奇數(shù)時(shí)

令 得到

注:DIT—FFT與DIF—FFT特點(diǎn)比較(重點(diǎn))

相同之處

(1)DIF與DIT兩種算法均為原位運(yùn)算

(2)DIF與DIT運(yùn)算量相同

所以,DIF與DIT是兩種等價(jià)的FFT算法。

不同之處

(1)DIF與DIT兩種算法結(jié)構(gòu)倒過(guò)來(lái)

DIF的輸入序列為自然順序,輸出為倒序排列,與DIT的正好相反。

(2)蝶形結(jié)構(gòu)不同

DIF的復(fù)數(shù)乘法只出現(xiàn)在減法之后,DIT則是先作復(fù)數(shù)乘法后作加減法。

簡(jiǎn)答題:(重點(diǎn))

1. 比較DIT—FFT與DIF—FFT特點(diǎn)。

2. 為什么要進(jìn)行FFT變換?FFT變換的基本思想是什么?說(shuō)明基-2FFT分成哪兩種算法。說(shuō)出它們的英文名稱和中文含義。

(DFT的計(jì)算在數(shù)字信號(hào)處理中非常有用,但是由于DFT的計(jì)算量較大,即使采用計(jì)算機(jī)也很難對(duì)問(wèn)題進(jìn)行實(shí)時(shí)處理,通過(guò)引入其快速算法FFT,使DFT的計(jì)算大大簡(jiǎn)化,運(yùn)算時(shí)間一般可縮短一、二個(gè)數(shù)量級(jí)?;舅悸罚篋FT的運(yùn)算量與

成正比;如果一個(gè)大點(diǎn)數(shù)

的DFT能分解為若干小點(diǎn)數(shù)DFT的組合,則顯然可以達(dá)到減小運(yùn)算量的效果。根據(jù)把長(zhǎng)序列分解為短序列的分解形式不同,基-2FFT算法基本上可以分為兩大類:按時(shí)間抽取法—decimation in time, DIT;按頻率抽取法——decimation in frequency, DIF))



第五章:本章主要掌握IIR和FIR兩種濾波器的基本網(wǎng)絡(luò)結(jié)構(gòu)。

5.1 基本單元結(jié)構(gòu)

一個(gè)數(shù)字網(wǎng)絡(luò)可以用差分方程表示,也可以用單位脈沖響應(yīng)來(lái)表示,也可以用系統(tǒng)函數(shù)來(lái)表示。但是對(duì)于研究這個(gè)系統(tǒng)的實(shí)現(xiàn)方法,即它的運(yùn)算結(jié)構(gòu)來(lái)說(shuō),用方框圖或信號(hào)流圖最直接。對(duì)于延時(shí)、乘以系數(shù)以及相加三種基本運(yùn)算來(lái)說(shuō),方框圖和信號(hào)流圖表示法如下圖所示。

以二階數(shù)字濾波器y(n)=b1y(n-1)+b2y(n-2)+ax(n)為例,它的方框圖和信號(hào)流圖如下圖所示。

一般來(lái)說(shuō),用方框圖表示數(shù)字濾波器,結(jié)構(gòu)明顯、直觀;而用信號(hào)流圖來(lái)表示,則簡(jiǎn)單、方便。

利用圖論中的轉(zhuǎn)置定理,可以把一個(gè)信號(hào)流圖轉(zhuǎn)化為另一個(gè)等價(jià)的信號(hào)流圖。

轉(zhuǎn)置定理如果將流圖中所有支路方向都顛倒或反向,并交換輸入x(n)和輸出y(n),則其特性保持不變,新流圖是原流圖的轉(zhuǎn)置形式。

例如,上圖中流圖的轉(zhuǎn)置形式如下圖(a)所示,但通常的習(xí)慣是將輸入x(n)畫在流圖的左邊,而輸出畫在流圖的右邊,這樣得到圖(b)所示的轉(zhuǎn)置結(jié)構(gòu)。

5.2 無(wú)限長(zhǎng)脈沖響應(yīng)基本網(wǎng)絡(luò)結(jié)構(gòu)

IIR(Infinite Impulse Response無(wú)限長(zhǎng)脈沖響應(yīng))濾波器具有以下特點(diǎn)(重點(diǎn)):?jiǎn)挝幻}沖響應(yīng)h(n)無(wú)限長(zhǎng);系統(tǒng)函數(shù)H(z)在有限z平面(0<|z|<∞)上有極點(diǎn)存在;結(jié)構(gòu)上存在從輸出到輸入的反饋,即結(jié)構(gòu)是遞歸型的。

1. 直接型

對(duì)應(yīng)的系統(tǒng)函數(shù)為:

直接型包括直接Ⅰ型和直接Ⅱ型,書本講授的為直接Ⅱ型。直接Ⅱ型的推導(dǎo),利用到線性移不變系統(tǒng),交換級(jí)聯(lián)子系統(tǒng)的次序,系統(tǒng)函數(shù)不變。

對(duì)于直接Ⅱ型,要求能夠直接由差分方程或系統(tǒng)函數(shù)繪出相應(yīng)的信號(hào)流圖,反之亦然。

特點(diǎn):便于理解,累積誤差大,運(yùn)算速度相對(duì)慢。(重點(diǎn))

2.級(jí)聯(lián)型

對(duì)應(yīng)的系統(tǒng)函數(shù)為:

把濾波器用若干二階子網(wǎng)絡(luò)級(jí)聯(lián)起來(lái)構(gòu)成,每個(gè)二階子網(wǎng)絡(luò)采用直接Ⅱ型結(jié)構(gòu)來(lái)實(shí)現(xiàn)。

特點(diǎn):級(jí)聯(lián)型結(jié)構(gòu)中每一個(gè)一階網(wǎng)絡(luò)決定一個(gè)零點(diǎn)、一個(gè)極點(diǎn),每一個(gè)二階網(wǎng)絡(luò)決定一對(duì)零點(diǎn)、一對(duì)極點(diǎn)。相對(duì)直接型結(jié)構(gòu),其優(yōu)點(diǎn)是調(diào)整方便,此外,運(yùn)算累積誤差較直接型小。(重點(diǎn))

3.并聯(lián)型

對(duì)應(yīng)的系統(tǒng)函數(shù)為:

特點(diǎn):每一個(gè)一階網(wǎng)絡(luò)決定一個(gè)實(shí)數(shù)極點(diǎn),每一個(gè)二階網(wǎng)絡(luò)決定一對(duì)共軛極點(diǎn),調(diào)整極點(diǎn)位置方便,但調(diào)整零點(diǎn)位置不如級(jí)聯(lián)型方便。運(yùn)算誤差不積累。運(yùn)算速度最高。(重點(diǎn))

例:已知IIR DF的系統(tǒng)函數(shù)為畫出直接I型、直接II型的結(jié)構(gòu)流圖。

解:先將化為的有理式

直接I型:

直接II型:

例:已知IIR DF的系統(tǒng)函數(shù)為畫出級(jí)聯(lián)型和并聯(lián)型的結(jié)構(gòu)流圖。

解:級(jí)聯(lián)型



并聯(lián)型

5.3 有限長(zhǎng)脈沖響應(yīng)基本網(wǎng)絡(luò)結(jié)構(gòu)

FIR(Finite Impulse Response有限長(zhǎng)脈沖響應(yīng))濾波器具有以下特點(diǎn)(重點(diǎn)):?jiǎn)挝幻}沖響應(yīng)h(n)有限長(zhǎng);系統(tǒng)函數(shù)H(z)在|z|>0處收斂,對(duì)因果系統(tǒng)而言,極點(diǎn)全部位于z=0處;結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu),沒(méi)有輸出到輸入的反饋。

FIR濾波器有以下幾種基本結(jié)構(gòu):直接型;級(jí)聯(lián)型。

1.直接型

特點(diǎn):直觀明了,便于理解,但不便于調(diào)整參數(shù)。(重點(diǎn))

2.級(jí)聯(lián)型

將H(z)因式分解得到

特點(diǎn):每一個(gè)一階因子控制一個(gè)零點(diǎn),每一個(gè)二階因子控制一對(duì)共軛極點(diǎn),調(diào)整零點(diǎn)位置比直接型方便,但H(z)中的系數(shù)比直接型多(近似3/2N),因而需要的乘法器多。(重點(diǎn))

例:已知FIR DF的系統(tǒng)函數(shù)為畫出直接型和級(jí)聯(lián)型的結(jié)構(gòu)流圖。(重點(diǎn))

解:直接型:

級(jí)聯(lián)型:

簡(jiǎn)答題:(重點(diǎn))

1. IIR和FIR濾波器的基本結(jié)構(gòu)形式有哪些?各自有什么特點(diǎn)?

2. IIR和FIR濾波器的基本特點(diǎn)是什么?

名詞解釋:(重點(diǎn))

1. IIR濾波器

2. FIR濾波器


第六章:本章講授了設(shè)計(jì)IIR濾波器常用的兩種設(shè)計(jì)方法——脈沖響應(yīng)不變法和雙線性變換法。

6.1 引言

1.數(shù)字濾波器的分類

(1)IIR和FIR數(shù)字濾波器(重點(diǎn))

這是根據(jù)濾波器的單位脈沖響應(yīng)h(n)的長(zhǎng)度是否有限來(lái)劃分的。若h(n)是一個(gè)長(zhǎng)度為M+1的有限長(zhǎng)序列,通常將此時(shí)的系統(tǒng)稱為有限長(zhǎng)單位脈沖響應(yīng)(FIR,F(xiàn)initeImpulseResponse)系統(tǒng)。

如果系統(tǒng)函數(shù)的分母中除a0外,還有其它的ak不為零,則相應(yīng)的h(n)將是無(wú)限長(zhǎng)序列,稱這種系統(tǒng)為無(wú)限長(zhǎng)單位脈沖響應(yīng)(IIR,InfiniteImpulseResponse)系統(tǒng)。

(2)低通、高通、帶通、帶阻濾波器

注意:數(shù)字濾波器(DF)與模擬濾波器(AF)的區(qū)別

數(shù)字濾波器的頻率響應(yīng)都是以2π為周期的,濾波器的低通頻帶處于2π的整數(shù)倍處,而高頻頻帶處于π的奇數(shù)倍附近。

2.設(shè)計(jì)指標(biāo)描述

濾波器的指標(biāo)通常在頻域給出。數(shù)字濾波器的頻率響應(yīng)一般為復(fù)函數(shù),通常表示為

其中,稱為幅頻響應(yīng),稱為相頻響應(yīng)。對(duì)IIR數(shù)字濾波器,通常用幅頻響應(yīng)來(lái)描述設(shè)計(jì)指標(biāo),而對(duì)于線性相位特性的濾波器,一般用FIR濾波器設(shè)計(jì)實(shí)現(xiàn)。

IIR低通濾波器指標(biāo)描述:

——通帶截止頻率,——阻帶截止頻率,——通帶最大衰減,——阻帶最小衰減,——3dB通帶截止頻率

3.設(shè)計(jì)方法(重點(diǎn))

三步:(1)按照實(shí)際需要確定濾波器的性能要求。(2)用一個(gè)因果穩(wěn)定的系統(tǒng)函數(shù)去逼近這個(gè)性能要求。(3)用一個(gè)有限精度的算法去實(shí)現(xiàn)這個(gè)系統(tǒng)函數(shù)。

IIR濾波器常借助模擬濾波器理論來(lái)設(shè)計(jì)數(shù)字濾波器,(重點(diǎn))設(shè)計(jì)步驟為:先根據(jù)所給的濾波器性能指標(biāo)設(shè)計(jì)出相應(yīng)的模擬濾波器傳遞函數(shù)Ha(s)( butterworth濾波器設(shè)計(jì)法等,有封閉公式利用),然后由Ha(s)經(jīng)變換(脈沖響應(yīng)不變法或者雙線性變換法等)得到所需的數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。在變換中,一般要求所得到的數(shù)字濾波器頻率響應(yīng)應(yīng)保留原模擬濾波器頻率響應(yīng)的主要特性。為此要求:(重點(diǎn))

(1)因果穩(wěn)定的模擬濾波器必須變成因果穩(wěn)定的數(shù)字濾波器;

(2)數(shù)字濾波器的頻響應(yīng)模仿模擬濾波器的頻響。

6.2 脈沖響應(yīng)不變法、雙線性不變法設(shè)計(jì)IIR數(shù)字低通濾波器

設(shè)計(jì)數(shù)字濾波器可以按照技術(shù)要求先設(shè)計(jì)一個(gè)模擬低通濾波器,得到模擬低通濾波器的傳輸函數(shù)

,再按一定的轉(zhuǎn)換關(guān)系將

轉(zhuǎn)換成數(shù)字低通濾波器的系統(tǒng)函數(shù)

。其設(shè)計(jì)流程如圖所示6.1所示。

利用模擬濾波器的設(shè)計(jì)結(jié)果來(lái)求相應(yīng)的數(shù)字濾波器,可以用映射的方法來(lái)實(shí)現(xiàn),把s平面映射到z平面。這種由復(fù)變量s到復(fù)變量z之間的映射關(guān)系,必須滿足三點(diǎn)要求:

圖6.1 從模擬濾波器設(shè)計(jì)數(shù)字濾波器流程圖

因果穩(wěn)定的

映射成因果穩(wěn)定的

,即s平面的左半平面必須映射到z平面單位圓的內(nèi)部。(重點(diǎn))

的頻率響應(yīng)能模仿

的頻率響應(yīng),即s平面的虛軸必須映射到z平面的單位圓上。變換前后的濾波器在時(shí)域或頻域的主要特征(頻率響應(yīng)或單位沖激響應(yīng)等)應(yīng)盡可能相同或接近。將傳輸函數(shù)

從s平面轉(zhuǎn)換到z平面的方法有多種,主要有沖激不變法和雙線性變換法。

1 脈沖響應(yīng)不變法

設(shè)模擬濾波器的系統(tǒng)函數(shù)為

,相應(yīng)的單位沖擊響應(yīng)是

,

。LT[.]代表拉氏變換,對(duì)

進(jìn)行等間隔采樣,采樣間隔為T,得到

,將h(n)=

作為數(shù)字濾波器的單位脈沖響應(yīng),那么數(shù)字濾波器的系統(tǒng)函數(shù)

便是

變換。因此脈沖響應(yīng)不變法是一種時(shí)域逼近方法,它使

在采樣點(diǎn)上等于

。但是,模擬濾波器的設(shè)計(jì)結(jié)果是

,所以下面基于脈沖響應(yīng)不變法的思想,導(dǎo)出直接從

的轉(zhuǎn)換公式。

設(shè)模擬濾波器只有單階極點(diǎn),且分母多項(xiàng)式的階次高于多項(xiàng)式的階次,將用部分分式表示: 式中為的單階極點(diǎn)。將進(jìn)行逆拉氏變換,得到: 式中,是單位階躍函數(shù)。對(duì)進(jìn)行等間隔采樣,采樣間隔為,得到:

對(duì)上式進(jìn)行變換,得到數(shù)字濾波器的系統(tǒng)函數(shù),即

由于沖激響應(yīng)不變法(又稱標(biāo)準(zhǔn)z變換法)從s平面→z平面是通過(guò)的多值映射,并不是一種簡(jiǎn)單的一一對(duì)應(yīng)的代數(shù)映射關(guān)系。脈沖響應(yīng)不變法所得到的數(shù)字濾波器的頻率響應(yīng),不是簡(jiǎn)單的重現(xiàn)模擬濾波器的頻率響應(yīng),而是模擬濾波器頻率頻率響應(yīng)的周期延拓。所以希望所設(shè)計(jì)的濾波器是帶限的,否則頻率混疊現(xiàn)象會(huì)使設(shè)計(jì)出的數(shù)字濾波器在附近的頻率特性嚴(yán)重的偏離模擬濾波器在處的頻率特性,嚴(yán)重時(shí)使數(shù)字濾波器的指標(biāo)得不到滿足。

特點(diǎn)(重點(diǎn))

優(yōu)點(diǎn):1.頻率變換關(guān)系是線性的,即,如果不存在頻譜混疊現(xiàn)象,用這種方法設(shè)計(jì)的數(shù)字濾波器會(huì)很好地重現(xiàn)原模擬濾波器的頻響特性。2.數(shù)字濾波器的單位脈沖響應(yīng)完全模仿模擬濾波器的單位沖擊響應(yīng)波形,時(shí)域特性逼近好。缺點(diǎn):會(huì)產(chǎn)生不同程度的頻譜混疊失真,其適合用于帶限濾波器(低通、帶通濾波器)的設(shè)計(jì),不適合用于高通、帶阻濾波器的設(shè)計(jì)。

2 雙線性變換法

(這種映射法能保證使s平面與z平面建立單值對(duì)應(yīng),從而消除混疊現(xiàn)象)

將雙線性變換帶入,得

s平面的

與z平面的

成非線性正切關(guān)系。正是因?yàn)槿绱?,雙線性變換法消除了頻率混疊現(xiàn)象。但與此同時(shí)也直接影響數(shù)字濾波器頻響逼真的模仿模擬濾波器的頻響。這種頻率之間的非線性變換關(guān)系帶來(lái)了以下兩個(gè)問(wèn)題:

(1) 一個(gè)線性相位的模擬濾波器經(jīng)雙線性變換后得到非線性相位的數(shù)字濾波器;

(2) 如果模擬濾波器的頻響具有片斷常數(shù)特性,則轉(zhuǎn)換到z平面上,數(shù)字濾波器仍具有片斷常數(shù)特性,但特性轉(zhuǎn)折點(diǎn)頻率值與模擬濾波器特性轉(zhuǎn)折點(diǎn)頻率值成非線性關(guān)系。

這種頻率的畸變,可通過(guò)對(duì)頻率進(jìn)行預(yù)修正來(lái)予以校正。

特點(diǎn)(重點(diǎn))

優(yōu)點(diǎn):1.不產(chǎn)生頻域混疊現(xiàn)象2.雙線性變換法可由簡(jiǎn)單的代數(shù)公式將直接轉(zhuǎn)換成。缺點(diǎn):與之間的非線性關(guān)系是雙線性變換法的缺點(diǎn),是數(shù)字濾波器頻響曲線不能保真地模仿模擬濾波器的頻響曲線形狀。只適合片段常數(shù)特性的濾波器設(shè)計(jì)。

例:已知模擬濾波器傳輸函數(shù)為,設(shè),用脈沖響應(yīng)不變法和雙線性變換法將轉(zhuǎn)換為數(shù)字濾波器系統(tǒng)函數(shù)。

解:用脈沖響應(yīng)不變法(令)將轉(zhuǎn)換為數(shù)字濾波器系統(tǒng)函數(shù)。

。

用雙線性變換法將轉(zhuǎn)換為數(shù)字濾波器系統(tǒng)函數(shù)。

6.3 設(shè)計(jì)IIR數(shù)字濾波器的頻率變換法

前面只介紹了IIR低通數(shù)字濾波器的設(shè)計(jì)方法,但是在工程上常常要設(shè)計(jì)各種截止頻率的低通、高通、帶通和帶阻數(shù)字濾波器,這些數(shù)字濾波器的設(shè)計(jì)方法通常是在設(shè)計(jì)一個(gè)低通濾波器的基礎(chǔ)上采用頻率變換法把低通濾波器轉(zhuǎn)換成所要求的濾波器。因?yàn)閿?shù)字濾波器設(shè)計(jì)都是從低通濾波器開始的,所以其他類型濾波器的設(shè)計(jì)都要把給定的技術(shù)指標(biāo)轉(zhuǎn)換成相應(yīng)的模擬低通濾波器指標(biāo)進(jìn)行設(shè)計(jì),然后再轉(zhuǎn)換成相應(yīng)的數(shù)字濾波器指標(biāo)。

這種轉(zhuǎn)換方法可以在模擬域進(jìn)行,也可以在數(shù)字域完成。用兩種變換方法進(jìn)行各種頻率范圍的數(shù)字濾波器的設(shè)計(jì)過(guò)程可表示如下:

(a) 模擬頻率變換法設(shè)計(jì)流程

(b) 數(shù)字頻率變換法設(shè)計(jì)流程

圖6.2 頻率變換法設(shè)計(jì)數(shù)字濾波器流程

上面兩種方法中模擬低通濾波器的技術(shù)指標(biāo)可根據(jù)所設(shè)計(jì)的技術(shù)指標(biāo)轉(zhuǎn)換獲得。也可給定一組相應(yīng)指標(biāo)進(jìn)行設(shè)計(jì),然后以此濾波器為基礎(chǔ)進(jìn)行各種變換。設(shè)計(jì)好的濾波器必須求其頻率特性,進(jìn)行指標(biāo)驗(yàn)證。如果不滿足給定數(shù)字濾波器指標(biāo),則需要修改模擬指標(biāo)重新設(shè)計(jì)。

簡(jiǎn)答題:(重點(diǎn))

1. 分析采用雙線性變換法將該模擬濾波器轉(zhuǎn)換成數(shù)字濾波器,數(shù)字濾波器的頻率特性相對(duì)原模擬濾波器頻率特性是否有失真,為什么?

2. 是否可以用脈沖響應(yīng)不變法或雙線性變換法把模擬帶阻濾波器變換成數(shù)字帶阻濾波器,為什么?

(答案:可以用雙線性變換法把模擬帶阻濾波器變換成數(shù)字帶阻濾波器,但不能用脈沖響應(yīng)不變法把模擬帶阻濾波器變換成數(shù)字帶阻濾波器。因?yàn)槊}沖響應(yīng)不變法中從s平面到z平面的映射是多值映射,將模擬帶阻濾波器變成數(shù)字帶阻濾波器會(huì)存在頻率響應(yīng)的混疊失真,而雙線性變換法從中從s平面到z平面的映射是一一對(duì)應(yīng)的,不存在頻率響應(yīng)的混疊失真。)

3. 說(shuō)明脈沖響應(yīng)不變法和雙線性變換法各自的優(yōu)點(diǎn)和缺點(diǎn)和適應(yīng)范圍?

4. 說(shuō)明IIR數(shù)字濾波器間接設(shè)計(jì)法(由模擬濾波器到數(shù)字濾波器)的基本思想和設(shè)計(jì)要求。

名詞解釋:(重點(diǎn))

1. IIR數(shù)字濾波器間接設(shè)計(jì)法

2. 雙線性變換法
(雙線性變換法是為克服脈沖響應(yīng)不變法頻率混疊現(xiàn)象提出的,其思想是將s平面上的點(diǎn)經(jīng)過(guò)和二次映射到z平面。雙線性變換法的優(yōu)點(diǎn)是不會(huì)產(chǎn)生頻率混疊,但卻引入頻率的非線性失真。)

3. 脈沖響應(yīng)不變法
(脈沖響應(yīng)不變法思想是把模擬濾波器的單位沖激響應(yīng)經(jīng)過(guò)采樣作為相應(yīng)數(shù)字濾波器的單位沖激響應(yīng)。其

平面和

平面的映射關(guān)系為:

。沖激不變法的優(yōu)點(diǎn)是數(shù)字頻率和模擬頻率為線性關(guān)系即:

。缺點(diǎn)是會(huì)產(chǎn)生頻率混疊,適用低通設(shè)計(jì)。)


第七章:本章主要講授線性相位FIR濾波器常用的設(shè)計(jì)方法——窗函數(shù)法


7.1 線性相位FIR濾波器的特點(diǎn)(重點(diǎn))

特點(diǎn):可實(shí)現(xiàn)嚴(yán)格的線性相位特性、系統(tǒng)是穩(wěn)定的、因果的、階數(shù)較高

穩(wěn)定和線性相位特性是FIR濾波器最突出的優(yōu)點(diǎn)。

設(shè)FIR單位脈沖響應(yīng)

長(zhǎng)度為

,其傳輸函數(shù)

(重點(diǎn))

其中,——幅度特性(重點(diǎn)),純實(shí)數(shù),可正可負(fù),即

——相位特性(重點(diǎn))

線性相位是指是的線性函數(shù),即群時(shí)延(常數(shù))

1 線性相位FIR數(shù)字濾波器時(shí)域約束條件(重點(diǎn))

1第一類線性相位(嚴(yán)格線性相位特性)對(duì)h(n)的約束條件,要求 和滿足:

2第二類線性相位對(duì)h(n)的約束條件,要求 和滿足:

問(wèn):何為線性相位濾波器?FIR濾波器成為線性相位濾波器的充分條件是什么?(重點(diǎn))

答:線性相位的濾波器是指其相位函數(shù)與數(shù)字頻率成線性關(guān)系,即。

FIR濾波器成為線性相位的充分條件是:

①是實(shí)數(shù)。

②滿足以為中心的偶對(duì)稱或者奇對(duì)稱,即

2.線性相位FIR濾波器幅度特性的特點(diǎn)(重點(diǎn))

(1)h(n)=h(N-1-n),N為奇數(shù)——1型

(2)h(n)=h(N-1-n),N為偶數(shù)——2型

(3)h(n)=-h(N-1-n),N為奇數(shù)——3型

(4)h(n)=-h(N-1-n),N為偶數(shù)——4型



類型

h(n)

Hg(ω)



1型

h(n)=h(N-1-n),N為奇數(shù)

Hg(ω)關(guān)于ω=0、π、2π偶對(duì)稱



2型

h(n)=h(N-1-n),N為偶數(shù)

Hg(ω)關(guān)于ω=0、2π偶對(duì)稱,關(guān)于ω=π奇對(duì)稱



3型

h(n)=-h(N-1-n),N為奇數(shù)

Hg(ω)關(guān)于ω=0、π、2π奇對(duì)稱



4型

h(n)=-h(N-1-n),N為偶數(shù)

Hg(ω)關(guān)于ω=0、2π奇對(duì)稱,關(guān)于ω=π偶對(duì)稱



實(shí)際使用時(shí),一般來(lái)說(shuō),1型適合構(gòu)成低通、高通、帶通、帶阻濾波器;2型適合構(gòu)成低通、帶通濾波器;3型適合構(gòu)成帶通濾波器;4型適合構(gòu)成高通、帶通濾波器。(重點(diǎn))

舉例說(shuō)明(重點(diǎn)):

系統(tǒng)函數(shù):

差分方程:

因?yàn)闉槠鎸?duì)稱,N=6為偶數(shù)。所以是4型第二類線性相位、N為偶數(shù)情況

3.線性相位FIR濾波器零點(diǎn)分布特點(diǎn)

由于,線性相位FIR濾波器的系統(tǒng)函數(shù)的零點(diǎn)分布具有如下特點(diǎn):如果是的零點(diǎn),其倒數(shù)也必然是其零點(diǎn);因?yàn)槭菍?shí)序列,故的零點(diǎn)也必定共軛成對(duì),所以和也是其零點(diǎn)。只要知道其中一個(gè),其他三個(gè)也就確定了。所以,零點(diǎn)必是互為倒數(shù)的共軛對(duì)(重點(diǎn))

7.2窗函數(shù)法設(shè)計(jì)線性相位FIR濾波器的一般步驟

通過(guò)第6章可知IIR數(shù)字濾波器的設(shè)計(jì)的主要方法是先設(shè)計(jì)一個(gè)模擬低通濾波器,然后把它轉(zhuǎn)換成形式的數(shù)字濾波器。但是對(duì)于FIR濾波器來(lái)說(shuō),設(shè)計(jì)方法的關(guān)鍵要求之一就是保證線性相位條件。而IIR濾波器的設(shè)計(jì)方法中只對(duì)幅度特性進(jìn)行了設(shè)計(jì)故無(wú)法保證這一點(diǎn)。所以,F(xiàn)IR濾波器的設(shè)計(jì)需要采用完全不同的方法。FIR濾波器的設(shè)計(jì)方法主要有窗函數(shù)法、頻率采樣法、切比雪夫逼近法等。


FIR濾波器的窗函數(shù)法設(shè)計(jì)思想是:在保證線性相位條件的前提下,選擇合適的長(zhǎng)度,使其傳輸函數(shù)滿足技術(shù)指標(biāo)要求。用一個(gè)長(zhǎng)度為的序列替代無(wú)限長(zhǎng)作為實(shí)際設(shè)計(jì)的濾波器的的單位脈沖響應(yīng),其系統(tǒng)函數(shù)為這種設(shè)計(jì)思想稱為窗函數(shù)設(shè)計(jì)法。顯然在保證對(duì)稱性的前提下,窗函數(shù)長(zhǎng)度越長(zhǎng),則越接近。但是誤差是肯定存在的,這種誤差稱為截?cái)嗾`差。

對(duì)FIR濾波器的影響:由于加窗后無(wú)限長(zhǎng)的變?yōu)橛邢揲L(zhǎng)的,肯定會(huì)引起誤差,表現(xiàn)在頻域就是通常所說(shuō)的吉布斯(Gibbs)效應(yīng)。吉布斯效應(yīng)直接影響濾波的性能,導(dǎo)致通帶內(nèi)的平穩(wěn)性變差和阻帶衰減不能滿足技術(shù)指標(biāo)。通常濾波器設(shè)計(jì)都要求過(guò)渡帶越窄越好,阻帶衰減越大越好。所以設(shè)計(jì)濾波器的方法要使吉布斯效應(yīng)的影響降低到最小。調(diào)整窗口長(zhǎng)度N只能有效的控制過(guò)渡帶的寬度,并不能減少帶內(nèi)波動(dòng)以及增大阻帶衰減。要想減小吉布斯效應(yīng)的影響,增加是無(wú)法實(shí)現(xiàn)的。如果改變窗函數(shù)的形狀,使其幅度函數(shù)具有較低的旁瓣幅度,就可減小通帶、阻帶的波動(dòng),并加大阻帶衰減。但是這時(shí)主瓣將會(huì)加寬以包含更多的能量,故而將會(huì)增加過(guò)渡帶寬度。所以當(dāng)一定時(shí),減小波動(dòng)和減小過(guò)渡帶是一對(duì)矛盾。必須根據(jù)實(shí)際要求,選擇合適的窗函數(shù)以滿足波動(dòng)要求,然后選擇滿足過(guò)渡帶指標(biāo)。

設(shè)計(jì)步驟:(1)根據(jù)對(duì)阻帶衰減以及過(guò)渡帶的指標(biāo)要求,選擇窗函數(shù)的類型,并估計(jì)窗口長(zhǎng)度N;(2)構(gòu)造希望逼近的頻率響應(yīng)函數(shù),即(3)計(jì)算:如果給出待求濾波器的頻響函數(shù)為,那么在單位脈沖響應(yīng)作用下:;(4)加窗得到設(shè)計(jì)結(jié)果:

簡(jiǎn)答題(重點(diǎn)):

1. 用窗口法設(shè)計(jì)FIR數(shù)字濾波器時(shí),為了改善阻帶的衰減特性,窗函數(shù)形狀需要滿足的兩個(gè)標(biāo)準(zhǔn)是什么?

對(duì)窗口函數(shù)形狀要求:(1)盡量減少窗口頻譜中的旁瓣,使能量盡量集中在主瓣中,這樣可減少上沖和余振,提高阻帶衰減。(2)主瓣寬度盡量窄,以獲得較陡的過(guò)渡帶。

2. 什么是吉布斯現(xiàn)象(效應(yīng))?如何有效減少該效應(yīng)的影響?(有限長(zhǎng)的序列h(n)去代替無(wú)限長(zhǎng)的hd(n),肯定會(huì)引起誤差,表現(xiàn)在頻域就是通常所說(shuō)的吉布斯(Gibbs)效應(yīng)。該效應(yīng)引起過(guò)渡帶加寬以及通帶和阻帶內(nèi)的波動(dòng),尤其使阻帶的衰減小,從而滿足不了技術(shù)上的要求。減少該效應(yīng)只能從窗函數(shù)的形狀上找解決問(wèn)題的方法。構(gòu)造新的窗函數(shù)形狀,使其譜函數(shù)的主瓣包含更多的能量,相應(yīng)旁瓣幅度更小。旁瓣的減小可使通帶、阻帶波動(dòng)減小,從而加大阻帶衰減。)

3. 窗口法設(shè)計(jì)FIR數(shù)字濾波器,改變窗的寬度對(duì)濾波器的頻率特性有什么影響?

4. 簡(jiǎn)要說(shuō)明窗函數(shù)法設(shè)計(jì)線性相位FIR濾波器的一般步驟。

5. 綜述數(shù)字濾波器的兩個(gè)主要分類IIR和FIR濾波器的特點(diǎn)和設(shè)計(jì)方法的不同。




數(shù)字信號(hào)處理復(fù)習(xí)總結(jié)-最終版的評(píng)論 (共 條)

分享到微博請(qǐng)遵守國(guó)家法律
尼木县| 四会市| 海林市| 台州市| 西盟| 西平县| 阆中市| 蒙山县| 西昌市| 浠水县| 柳河县| 临西县| 侯马市| 谢通门县| 宣汉县| 焦作市| 拉萨市| 临沂市| 若尔盖县| 仁怀市| 蒙阴县| 通道| 陕西省| 塘沽区| 岐山县| 淮北市| 贵南县| 岳西县| 湘阴县| 信宜市| 托克托县| 汝城县| 青阳县| 福清市| 洪湖市| 高淳县| 嘉兴市| 如皋市| 德清县| 页游| 临夏县|