數(shù)字信號(hào)處理實(shí)驗(yàn) 基于Matlab

《數(shù)字信號(hào)處理實(shí)驗(yàn) 基于Matlab》由會(huì)員分享，可在線閱讀，更多相關(guān)《數(shù)字信號(hào)處理實(shí)驗(yàn) 基于Matlab（64頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、數(shù)字信號(hào)處理實(shí)驗(yàn)（基于Matlab）王春興物理與電子科學(xué)學(xué)院2003.5數(shù)字信號(hào)處理實(shí)驗(yàn)（基于Matlab）第一章 概述信號(hào)的定義及表達(dá)信號(hào)：帶有信息的任何物理量本課程討論的信號(hào)：時(shí)間函數(shù) 連續(xù)時(shí)間信號(hào) CT 離散時(shí)間信號(hào)DT 模擬信號(hào)到數(shù)字信號(hào) 采樣（時(shí)間離散化） 量化（取值離散化）（P.7 圖1.6）數(shù)字信號(hào)的表達(dá)：二進(jìn)制數(shù)組 -二進(jìn)制編碼 數(shù)制轉(zhuǎn)換： 10進(jìn)制-2進(jìn)制 MSB和LSB 符號(hào)數(shù)的表達(dá)：符號(hào)位（補(bǔ)碼） 數(shù)字信號(hào)處理分析方式： 理論分析設(shè)計(jì) 采用DT信號(hào)進(jìn)行 （本課程內(nèi)容）實(shí)時(shí)電路處理 采用數(shù)字信號(hào)進(jìn)行 （數(shù)字電路內(nèi)容）數(shù)字信號(hào)的特點(diǎn)及應(yīng)用特點(diǎn)：數(shù)字信號(hào)與模擬信號(hào)的比較抗干擾性

2、強(qiáng)、精度高、容易存儲(chǔ)、可靈活處理與計(jì)算機(jī)系統(tǒng)兼容數(shù)字技術(shù)的應(yīng)用領(lǐng)域語音技術(shù)（傳輸、識(shí)別、合成）圖象處理（靜止圖象、移動(dòng)圖象、三維動(dòng)畫）地波分析（地震探測(cè)、地質(zhì)探礦）諧振分析（高層建筑、橋梁、機(jī)翼等）自動(dòng)控制 實(shí)時(shí)檢測(cè)本課程的主要教學(xué)安排：主要內(nèi)容：（50學(xué)時(shí)）數(shù)字信號(hào)的頻率分析：定義、變換與計(jì)算 (22學(xué)時(shí)) 頻率定義、CTFS與DTFS、CTFT與DTFT、DFT與FFT數(shù)字信號(hào)的頻率處理：濾波器設(shè)計(jì) （26學(xué)時(shí)） LTI系統(tǒng)分析理想濾波器與低階數(shù)字濾波；FIR濾波器設(shè)計(jì)、IIR濾波器設(shè)計(jì) 數(shù)字濾波器結(jié)構(gòu)與誤差分析教材與參考書： 教材： Digital Signal Processing s

3、pectral computation and filter design （Third edition） (美) Chi-Tsong Chen 電子工業(yè)出版社 2002版 數(shù)字信號(hào)處理基礎(chǔ)（加）Joyce Van de Vegte 著 侯正信 王國(guó)安 等譯 電子工業(yè)出版社2003版信號(hào)與系統(tǒng)計(jì)算機(jī)練習(xí)利用MATLAB John R.Buck 劉樹棠 譯 西安交通大學(xué)出版社 2000版 主要工具：MATLAB ：信號(hào)波形圖、頻譜計(jì)算與分析濾波器設(shè)計(jì)及系統(tǒng)頻率特性分析考核方式：平時(shí)作業(yè) 30%考試（筆試、操作） 70% MATLAB的基本應(yīng)用方法命令窗口（Command window）的使用：輸

4、入各類變量或函數(shù)名稱，按回車即得到當(dāng)前變量或函數(shù)值；輸入各類命令，按回車即得到該命令執(zhí)行結(jié)果；若需要輸入多行命令或程序，各行間用“；”間隔；M文件的編制與調(diào)試執(zhí)行打開空白文件或已經(jīng)有的文件，進(jìn)行程序文件的輸入編輯；各行間用“；”間隔；一行中“%”以后內(nèi)容為注釋部分，不影響程序執(zhí)行；程序編制完畢后，可以按“F5”鍵保存執(zhí)行，注意根據(jù)屏幕提示建立文件名稱；如果出現(xiàn)錯(cuò)誤，可在命令窗口看到錯(cuò)誤類型及位置，根據(jù)錯(cuò)誤檢測(cè)信息對(duì)程序進(jìn)行調(diào)試；MATLAB命令及函數(shù)信號(hào)的表達(dá)方式及作圖在MATLAB中，任何變量或函數(shù)均表現(xiàn)為向量，任何向量的元素編號(hào)均從1開始；序列（向量）表達(dá)方式 設(shè)定坐標(biāo)向量n和信號(hào)向量x；

5、x和n為長(zhǎng)度相同的向量，向量的編號(hào)從1開始； n=-2 :0.1:2 坐標(biāo)向量可以直接逐點(diǎn)寫出：n=2 3 4 5 6 7； 也可以采用起點(diǎn)，終點(diǎn)和步長(zhǎng)的形式寫出：n=-2 :0.1:2 ； 信號(hào)向量可以直接逐點(diǎn)寫出：x=1 2 3 4 3 2； 也可以采用與n有關(guān)的函數(shù)運(yùn)算形式寫出： 例如： x=3*n x=exp(j*(pi/8)*n) 作圖： 采用stem（n，x） 作出離散圖形 DT信號(hào) 采用plot（n，x） 作出連續(xù)圖形（折線連接） CT信號(hào)作圖時(shí)主要通過合理設(shè)置n的范圍及步長(zhǎng)來保證變量坐標(biāo)的正確性；可以利用title，axis等函數(shù)為圖形設(shè)置說明和坐標(biāo)范圍；特別注意：作圖時(shí)必須保

6、證坐標(biāo)向量與信號(hào)向量長(zhǎng)度完全一致；0101：離散序列的作圖直接表現(xiàn)離散序列n=2 3 4 5 6 7;x=1 2 3 4 3 2;stem(n,x);0102：將圖形表現(xiàn)為連續(xù)曲線n=2 3 4 5 6 7;x=1 2 3 4 3 2;plot(n,x);0203：信號(hào)表現(xiàn)為坐標(biāo)向量的函數(shù)n=2 3 4 5 6 7;x=exp(j*(pi/8)*n);plot(n,x);0204：圖形說明和坐標(biāo)范圍的設(shè)置n=-20:0.5:20;x=exp(j*(pi/8)*n);plot(n,x),title(n=-20:0.5:20;x=exp(j*(pi/8)*n);plot(n,x);axis(-20

7、,20,-2,2);第二章 DTFS和CTFS-周期信號(hào)的頻率分量信號(hào)的時(shí)域表達(dá)形式：連續(xù)時(shí)間信號(hào)CT 離散時(shí)間信號(hào)DTDT信號(hào)由CT信號(hào)采樣得出CT信號(hào) 采樣信號(hào) DT信號(hào) 周期信號(hào)：每隔一段時(shí)間重復(fù)的信號(hào) 信號(hào)變化快慢的描述：周期T：信號(hào)重復(fù)的時(shí)間間隔 頻率：?jiǎn)挝粫r(shí)間內(nèi)信號(hào)重復(fù)的次數(shù) 頻域表達(dá)形式 人類接受的自然信號(hào)主要以頻率形式表達(dá)： 聲音 色彩信號(hào)分析的重要任務(wù)：從時(shí)域信號(hào)得出頻域信號(hào)頻率的定義-單頻率時(shí)間信號(hào)CT信號(hào)的頻率 與周期性密切相關(guān) 標(biāo)準(zhǔn)信號(hào) 上述信號(hào)均為周期信號(hào)，周期為T；頻率均為 ； 周期與頻率的關(guān)系 取值范圍 DT信號(hào)的頻率 (仿照CT信號(hào)定義)DT信號(hào)由CT信號(hào)采樣得

8、到；周期信號(hào)的采樣不一定為周期信號(hào)；對(duì)于DT信號(hào)，頻率與周期沒有直接對(duì)應(yīng)關(guān)系；仿照CT信號(hào)定義基本信號(hào) 上述信號(hào)頻率均為 要點(diǎn)：信號(hào)的時(shí)域性質(zhì)：非周期性頻率性質(zhì)：多重性 頻率范圍 取值范圍 Nyquist frequency range 與采樣頻率有關(guān)CT與DT信號(hào)頻率之間的關(guān)系：P. 29 fd與采樣頻率有關(guān)，位于范圍內(nèi)； 將高頻按周期折合到低頻（主值區(qū)）；例 p.27 2.1 2.2付氏級(jí)數(shù)及頻率分量一般周期時(shí)間信號(hào)的頻率：采用單頻率信號(hào)表達(dá)CTFS （連續(xù)時(shí)間付氏級(jí)數(shù)）定義式 （周期信號(hào)：周期為P） 存在條件：絕對(duì)可積，而且在一個(gè)周期內(nèi)間斷點(diǎn)和極值點(diǎn)有限； 頻率分量：CTFS coeff

9、icient / frequency component通常為復(fù)數(shù)： Magnitude phase phase的范圍 共軛對(duì)稱性：奇偶性 real : even : odd real and even real and even real and odd imaginary and odd 周期信號(hào)頻率分量的計(jì)算 （例 P.38） 要點(diǎn)：2.4 利用簡(jiǎn)單分解求2.5 利用公式求矩形脈沖的 2.6 利用公式求沖激串（sampling function）的周期信號(hào)的時(shí)間范圍： 雙邊信號(hào)一般周期信號(hào)的頻率范圍： 若對(duì)于所有的，則稱為帶限信號(hào)；頻率分量的意義：平均功率: 只與magnitude有關(guān)p

10、hase的影響：對(duì)信號(hào)波形的影響 對(duì)視頻有影響，對(duì)音頻沒有影響DTFS（離散時(shí)間付氏級(jí)數(shù)） 定義式 頻率分量討論： 的周期性 只有N個(gè)獨(dú)立系數(shù) 共軛對(duì)稱性：奇偶性 m的取值范圍：應(yīng)使 的范圍為時(shí)移的影響： 線性相位變化例 p.4854 2.8 利用公式直接計(jì)算系數(shù) N=3 2.9 利用公式直接計(jì)算系數(shù) N=42.10 時(shí)移的作用時(shí)移不影響DTFS系數(shù)的幅度，只在系數(shù)中加入線性相位；利用MATLAB計(jì)算頻率分量DTFS系數(shù)的FFT計(jì)算重要函數(shù) c=(1/N)*fft(x) x=N*ifft(c)應(yīng)用要點(diǎn)：x和c的序列都為N個(gè)元素，下標(biāo)排列都為 1 .N，分別對(duì)應(yīng)于離散時(shí)間n 0.N-1和離散頻率

11、 m 0.N-1；DTFS可以利用任何一個(gè)周期進(jìn)行計(jì)算；在采用FFT的時(shí)候，輸入數(shù)據(jù)必須從n=0到N-1；利用shift（）函數(shù)（P.56）可以使頻率向量排布在對(duì)稱區(qū)間內(nèi)：N為odd時(shí)，下標(biāo)0在正中，N為even時(shí)，下標(biāo)0偏左；利用m=ceil(-(N-1)/2):ceil(N-1/2)可以得到對(duì)應(yīng)的橫坐標(biāo)；注意：c為周期序列，周期為N；例2.8題的求解 N=3;T=0.5; x=-2 1 0.6;例2.9題的求解 N=4;T=1; x= 2.5 -0.4 1 -2;%program 2.2N=3;T=0.5; x=-2 1 0.6;D=2*pi/(N*T);X=fft(x/N);m=ceil

12、(-(N-1)/2):ceil(N-1)/2);w=m*D;subplot(2,1,1),stem(w,abs(shift(X),title(a);subplot(2,1,2),stem(w,angle(shift(X)*180/pi),title(b);CTFS系數(shù)的FFT計(jì)算問題：x為連續(xù)信號(hào)；m取值范圍為無限大區(qū)間； 方案：在一個(gè)周期內(nèi)取N點(diǎn)對(duì)x進(jìn)行采樣（離散化）； 求出DTFS系數(shù)-周期序列； 取主值范圍內(nèi)的序列即為對(duì)應(yīng)CTFS系數(shù)； 當(dāng)x為帶限信號(hào)時(shí)，在滿足采樣定理?xiàng)l件下可以得出準(zhǔn)確的CTFS系數(shù)；例 2.4題的求解 2.12%program 2.2P=2*pi/0.3;N=11;T

13、=P/N;D=2*pi/P;n=0:N-1;x=-1.2+0.8*sin(0.6*n*T)-1.6*cos(1.5*n*T);X=fft(x/N);m=ceil(-(N-1)/2):ceil(N-1)/2);w=m*D;subplot(2,1,1),stem(w,abs(shift(X),title(a);subplot(2,1,2),stem(w,angle(shift(X)*180/pi),title(b); 對(duì)于帶限信號(hào)，在滿足采樣定理的條件下，不同大小的N值（采樣數(shù)量）得到的幅頻分量相同； 若采樣周期不夠小，則將產(chǎn)生頻率混疊失真；例2.3 頻率混疊的影響例2.5題的求解通過改變采樣點(diǎn)數(shù)

14、量N，可以比較混疊的影響大小%program 2.4 hold on;N=42;P=4;T=P/N;D=2*pi/P;q=floor(1/T);x=ones(1,q+1) zeros(1,N-2*q-1) ones(1,q);X=fft(x/N);m=ceil(-(N-1)/2):ceil(N-1)/2);stem(m*D,shift(X),b,fill);對(duì)于在足夠逼近條件下，magnetude可以得到足夠良好近似值；（能量逼近）；但計(jì)算出的phase不會(huì)得出良好近似值（通常不采用）； 關(guān)鍵：N的選?。ㄗ銐虼笠垣@得良好近似，足夠小以減少運(yùn)算量） 根據(jù)設(shè)定的精確度，求出最小的N： 令N=2n，

15、選定n的特定值，再逐1增加；重復(fù)計(jì)算在Nyquist范圍內(nèi)的系數(shù)差，直到系數(shù)差小于設(shè)定值為止；例2.6 根據(jù)設(shè)定最大誤差，自動(dòng)選取最小采樣點(diǎn)數(shù)量，并求出滿足要求的頻譜例2.5題的求解% program 2.6a=1;b=100;P=4;D=2*pi/P;beta=1;while bbeta N1=2a;T1=P/N1;q1=floor(1/T1); x1=ones(1,q1+1) zeros(1,N1-2*q1-1) ones(1,q1); X1=fft(x1/N1); N2=2*N1;T2=P/N2;q2=floor(1/T2); x2=ones(1,q2+1) zeros(1,N2-2*q

16、2-1) ones(1,q2); X2=fft(x2/N2); m1p=0:N1/2; d=max(abs(abs(X1(m1p+1)-abs(X2(m1p+1); mm=max(abs(X1(m1p+1); b=d/mm*100; a=a+1;endN2,bm=-N2/2+1:N2/2;stem(m*D,abs(shift(X2);對(duì)此程序作少數(shù)改動(dòng)可以得到對(duì)其他信號(hào)的計(jì)算：P.73 例2.14 平均功率計(jì)算 CT信號(hào) DT信號(hào) 例2.15 分別采用時(shí)域序列和頻譜序列求信號(hào)平均功率例2.4題的信號(hào)功率%program 2.8P=2*pi/0.3;N=11;T=P/N;n=0:N-1;x=-1

17、.2+0.8*sin(0.6*n*T)-1.6*cos(1.5*n*T);P1=sum(x.2)*T/PX=fft(x/N);P2=sum(abs(X).2)例2.16 采用頻譜序列求信號(hào)平均功率例2.5題的信號(hào)功率%program 2.9N=1024;D=2*pi/4;x=ones(1,257) zeros(1,511) ones(1,256);X=fft(x/N);mu=floor(5/D);m=2:mu+1;p=abs(X(1)2+2*sum(abs(X(m).2)第三章 CTFT和DTFT-一般信號(hào)的頻譜實(shí)際信號(hào)都是非周期信號(hào)（非雙邊信號(hào)）周期信號(hào)對(duì)應(yīng)于離散頻率分量（離散頻譜）非周期信

18、號(hào)對(duì)應(yīng)于連續(xù)頻譜CTFT 連續(xù)時(shí)間信號(hào)的付氏變換 定義式 頻譜存在條件：絕對(duì)可積，而且在一個(gè)周期內(nèi)間斷點(diǎn)和極值點(diǎn)有限； 例 3.1 eatu(t) 的頻譜P.85 圖3.1例3.2 wa（t）的頻譜：時(shí)域窗口函數(shù)P.87 圖3.2 主瓣：高2a，寬2/a 旁瓣：寬/a 零點(diǎn)：n /a 窗口寬度與主瓣/旁瓣寬度成反比； 窗口的時(shí)移不改變幅頻特性，只引入線性相位；例3.3 模擬理想低通濾波器 ：頻域窗口函數(shù) P.89 圖3.3CT周期信號(hào)的頻譜 步驟：將CT周期信號(hào)先展開為CTFS，再進(jìn)行逐項(xiàng)變換； 若干常用信號(hào)的頻譜（P.91 圖3.4） 例 3.4 沖激串的CTFT 頻譜的性質(zhì) 連續(xù)有界：若絕

19、對(duì)可積，則有界并連續(xù)；奇偶性： P.94 表3.1時(shí)間實(shí)函數(shù)-幅頻偶、相頻奇時(shí)間實(shí)偶-頻譜實(shí)偶時(shí)間實(shí)奇-頻譜虛奇時(shí)移與頻移： 時(shí)移引入線性相位，頻移對(duì)應(yīng)復(fù)指數(shù)調(diào)制時(shí)間尺度變換： 時(shí)間壓縮對(duì)應(yīng)頻譜擴(kuò)展Parsevals relation能量的頻率分布 能量只與幅頻特性有關(guān)，時(shí)移不影響信號(hào)能量分布周期信號(hào)與非周期信號(hào)的能量對(duì)比周期信號(hào)能量為無限大，平均功率為，非零功率只存在于離散頻率點(diǎn)；絕對(duì)可積的非周期信號(hào)能量為，在任何特定頻率點(diǎn)能量為零，能量分布于頻率區(qū)間上；連續(xù)時(shí)間信號(hào)截?cái)鄬?duì)于頻譜的影響只有極少數(shù)信號(hào)可以求出頻譜的解析表達(dá)式；絕大多數(shù)實(shí)際信號(hào)只能采用數(shù)值方式求解頻譜；對(duì)無限長(zhǎng)時(shí)間連續(xù)信號(hào)，在實(shí)

20、際計(jì)算時(shí)必須考慮截?cái)嗖㈦x散化；時(shí)域截?cái)嗄Ｐ停阂源翱诤瘮?shù)乘以時(shí)間函數(shù) 時(shí)域乘積對(duì)應(yīng)于頻域卷積： 其中 單頻率信號(hào)的截?cái)嘈Ч?（P.104 圖3.8）使單頻率展寬，出現(xiàn)主瓣（高L=2a、寬4/L）和旁瓣（高0.2L、寬2/L）；對(duì)于有限帶寬信號(hào)，截?cái)鄬?dǎo)致帶外泄露（能量）和紋波現(xiàn)象；L越小，上述效應(yīng)越顯著;對(duì)于連續(xù)信號(hào)，增大L可以將上述效應(yīng)削弱到可以忽略的程度；例3.5 的頻譜：采用不同寬度的窗口截?cái)啵唬≒.104 圖3.8）Gibbs現(xiàn)象用付氏變換表達(dá)時(shí)間函數(shù)時(shí)，當(dāng)頻譜信號(hào)含有不連續(xù)點(diǎn)時(shí)，頻譜的紋波將會(huì)變窄并靠近該點(diǎn)，但紋波不會(huì)隨L的無限增大而消失，而是趨于一個(gè)常量（寬度無限小，高度約為不連續(xù)變化

21、量的9%）；（P.105 圖3.9）將頻譜變換為時(shí)間函數(shù)時(shí)存在相同的現(xiàn)象；采用矩形窗口截?cái)嘈盘?hào)必然出現(xiàn)Gibbs現(xiàn)象。DTFT 離散時(shí)間信號(hào)的付氏變換定義式 頻譜 具有周期性 可以表達(dá)為任何一個(gè)周期上的積分 主值區(qū)域：例 3.6 anun 的頻譜 P.109 圖3.11例 3.7 離散時(shí)間矩形窗口函數(shù)wdn的頻譜N=2M+1 P.110 圖3.12與連續(xù)時(shí)間窗口類似，離散時(shí)間窗口的頻譜在一個(gè)周期內(nèi)也產(chǎn)生主瓣與旁瓣，頻譜零點(diǎn)位于處；主瓣高度為N；旁瓣只有N-2個(gè)； 例3.8 數(shù)字理想低通濾波器 （P.112 圖3.13）DT周期信號(hào)的頻譜 標(biāo)準(zhǔn)信號(hào) 時(shí)域沖激串對(duì)應(yīng)于頻域沖激串一般周期信號(hào)將DT周

22、期信號(hào)先展開為DTFS，再進(jìn)行逐項(xiàng)變換； 截?cái)嗟挠绊憣?duì)無限長(zhǎng)時(shí)間DT信號(hào)進(jìn)行計(jì)算時(shí)，必須進(jìn)行截?cái)?；時(shí)域截?cái)嗄Ｐ停阂源翱诤瘮?shù)乘以時(shí)間函數(shù) 頻域?qū)?yīng)卷積： 其中 N=2M+1截?cái)嘈Ч?P.116 圖3.14使單頻率展寬，出現(xiàn)主瓣（高L=2a、寬4/L）和旁瓣（高0.2L、寬2/L）；對(duì)于有限帶寬信號(hào)，截?cái)鄬?dǎo)致帶外泄露（能量）和紋波現(xiàn)象；L越小，上述效應(yīng)越顯著。矩形窗截?cái)啾厝粚?dǎo)致Gibbs現(xiàn)象 P.117 圖3.15連續(xù)信號(hào)的離散化連續(xù)信號(hào)不能進(jìn)行數(shù)值計(jì)算，必須離散化為離散信號(hào)；連續(xù)信號(hào)離散化過程稱為采樣過程；Nyquist 采樣定理 理想采樣：利用沖激串相乘使連續(xù)時(shí)間信號(hào)離散化 時(shí)域離散化 頻域周

23、期性復(fù)制恢復(fù)采樣信號(hào)的條件：1 帶限信號(hào) 存在最高頻率 2 采樣頻率（Nyquist rate）滿足上述條件時(shí)，可通過頻域的低通濾波（截止頻率為）分離出原始頻譜，恢復(fù)信號(hào)；該操作等效于時(shí)域的理想內(nèi)插恢復(fù)（在每個(gè)采樣點(diǎn)插入連續(xù)取樣函數(shù)）。頻率混疊問題不滿足條件2時(shí)，可能產(chǎn)生部分混疊（高頻），低頻信號(hào)可恢復(fù)；不滿足條件1時(shí)，總是產(chǎn)生混疊；當(dāng)CT信號(hào)能量有限時(shí)，可以增大T使得混疊影響足夠?。粚?duì)信號(hào)進(jìn)行預(yù)先低通濾波是消除混疊的有效方式。時(shí)限帶限理論任何非零連續(xù)信號(hào)都不可能即為時(shí)限又為帶限；利用數(shù)字技術(shù)處理連續(xù)信號(hào)時(shí)必然需要截?cái)啵厝划a(chǎn)生誤差；在一定誤差范圍內(nèi)，有限能量的連續(xù)信號(hào)可以近視看作為時(shí)限帶限信

24、號(hào)，并利用數(shù)字技術(shù)處理。音頻信號(hào)的實(shí)際處理過程模擬信號(hào)-采樣及零階保持-ADC-數(shù)字編碼信號(hào)數(shù)字編碼信號(hào)-DAC-零階保持信號(hào)-低通濾波-模擬信號(hào)DT周期信號(hào)的頻譜 標(biāo)準(zhǔn)信號(hào) 時(shí)域沖激串對(duì)應(yīng)于頻域沖激串一般周期信號(hào)將DT周期信號(hào)先展開為DTFS，再進(jìn)行逐項(xiàng)變換； 截?cái)嗟挠绊憣?duì)無限長(zhǎng)時(shí)間DT信號(hào)進(jìn)行計(jì)算時(shí)，必須進(jìn)行截?cái)?；時(shí)域截?cái)嗄Ｐ停阂源翱诤瘮?shù)乘以時(shí)間函數(shù) 頻域?qū)?yīng)卷積： 離散取樣函數(shù) N=2M+1截?cái)嘈Ч?P.1116 圖3.14使單頻率展寬，出現(xiàn)主瓣（高L=2a、寬4/L）和旁瓣（高0.2L、寬2/L）；對(duì)于有限帶寬信號(hào)，截?cái)鄬?dǎo)致帶外泄露（能量）和紋波現(xiàn)象；L越小，上述效應(yīng)越顯著。矩形窗截?cái)?/p>

25、必然導(dǎo)致Gibbs現(xiàn)象 P.1116 圖3.14Nyquist 采樣定理 采樣：利用沖激串相乘使連續(xù)時(shí)間信號(hào)離散化 恢復(fù)采樣信號(hào)的條件：3 帶限信號(hào) 存在最高頻率 4 采樣頻率（Nyquist rate）滿足上述條件時(shí)，可通過頻域的低通濾波（截止頻率為）分離出原始頻譜，恢復(fù)信號(hào)；該操作等效于時(shí)域的理想內(nèi)插恢復(fù)（在每個(gè)采樣點(diǎn)插入連續(xù)取樣函數(shù)）。頻率混疊問題不滿足條件2時(shí)，可能產(chǎn)生部分混疊（高頻），低頻信號(hào)可恢復(fù)；不滿足條件1時(shí)，總是產(chǎn)生混疊；當(dāng)CT信號(hào)能量有限時(shí)，可以增大T使得混疊影響足夠小；對(duì)信號(hào)進(jìn)行預(yù)先低通濾波是消除混疊的有效方式。時(shí)限帶限理論任何非零連續(xù)信號(hào)都不可能即為時(shí)限又為帶限；利用數(shù)

26、字技術(shù)處理連續(xù)信號(hào)時(shí)必然需要截?cái)?，必然產(chǎn)生誤差；在一定誤差范圍內(nèi)，有限能量的連續(xù)信號(hào)可以近視看作為時(shí)限帶限信號(hào)，并利用數(shù)字技術(shù)處理。音頻信號(hào)的實(shí)際處理過程模擬信號(hào)-采樣及零階保持-ADC-數(shù)字編碼信號(hào)數(shù)字編碼信號(hào)-DAC-零階保持信號(hào)-低通濾波-模擬信號(hào)第四章 DFT和FFT-頻譜的計(jì)算絕大多數(shù)信號(hào)不能采用解析方式進(jìn)行頻譜分析，只能采用數(shù)值計(jì)算方法計(jì)算頻譜；DFT 離散付氏變換 只對(duì)有限長(zhǎng)度的序列定義，由DTFT頻譜在一個(gè)周期內(nèi)進(jìn)行N點(diǎn)采樣得出定義式 特點(diǎn)：將有限信號(hào)的無限寬頻譜壓縮到一個(gè)周期內(nèi)（有限寬度）將離散信號(hào)的連續(xù)頻譜采樣為有限序列設(shè) 則 重要性質(zhì)時(shí)間序列和頻譜序列均以N為周期；可擴(kuò)展

27、為周期信號(hào)；可以從任何一個(gè)周期中計(jì)算；例 DFT和DTFT的對(duì)比計(jì)算DFT（離散序列）是DTFT（連續(xù)頻譜）的采樣P.136140 例4.1-4.3DFT與DTFS的關(guān)系：將有限長(zhǎng)時(shí)間序列看作無限長(zhǎng)時(shí)間序列的一個(gè)周期，就可以進(jìn)行DTFS變換 只相差常數(shù)N反變換的時(shí)域混疊問題當(dāng)時(shí)域序列為有限長(zhǎng)m時(shí)，作N點(diǎn)的DFTNm 不會(huì)產(chǎn)生混疊 DFT和DTFT反變換結(jié)果相同；Nbeta N1=2a; n1=0:N1-1; x1=0.9.n1;X1=fft(x1); N2=2*N1; n2=0:N2-1; x2=0.9.n2;X2=fft(x2); m1p=0:N1/2; d=max(abs(abs(X1(m

28、1p+1)-abs(X2(2*m1p+1); mm=max(abs(X1(m1p+1); b=d/mm*100;a=a+1;endN2,b連續(xù)時(shí)間信號(hào)的頻譜計(jì)算當(dāng)連續(xù)時(shí)間信號(hào)不能表達(dá)為閉合形式時(shí)，只能采用數(shù)值計(jì)算；連續(xù)信號(hào)的計(jì)算必須先在長(zhǎng)度為L(zhǎng)的區(qū)間內(nèi)經(jīng)時(shí)間采樣（周期T）成為N點(diǎn)序列，才能進(jìn)行計(jì)算；本節(jié)只考慮絕對(duì)可積的連續(xù)時(shí)間信號(hào)；正區(qū)間信號(hào) （t1%最大值）的計(jì)算： 為系統(tǒng)的最大極點(diǎn)，為序列長(zhǎng)度；有限時(shí)間后，可以認(rèn)為系統(tǒng)響應(yīng)進(jìn)入穩(wěn)態(tài)；在進(jìn)行系統(tǒng)分析時(shí)，通常可以不考慮暫態(tài)響應(yīng)。連續(xù)時(shí)間LTIL系統(tǒng)沖激響應(yīng)： 卷積形式： 微分方程形式：Laplace變換： 基本變換對(duì) P.241 表5.4重要性

29、質(zhì)Laplace變換與z-變換的關(guān)系系統(tǒng)函數(shù)對(duì)于穩(wěn)定因果系統(tǒng)，所有極點(diǎn)位于左半平面；穩(wěn)定性： 頻率響應(yīng)： 對(duì)于穩(wěn)定系統(tǒng)成立穩(wěn)態(tài)響應(yīng)與暫態(tài)響應(yīng)時(shí)間常數(shù)= （極點(diǎn)實(shí)部絕對(duì)值的倒數(shù)）MAX系統(tǒng)響應(yīng)與極點(diǎn)位置的關(guān)系P.244 圖5.12連續(xù)時(shí)間系統(tǒng)頻率響應(yīng)的計(jì)算與測(cè)量MATLAB函數(shù)b=b(1) b(2) b(M+1) ; a=a(1) a(2)a(N+1);H,w=freqs(b,a) 已知系統(tǒng)函數(shù)時(shí)，在區(qū)間內(nèi)自動(dòng)選擇200個(gè)頻率點(diǎn)計(jì)算頻率響應(yīng)；不知道系統(tǒng)函數(shù)時(shí)，只能通過測(cè)量得到頻率響應(yīng)，再求得系統(tǒng)函數(shù)；測(cè)量方案：1 使用標(biāo)準(zhǔn)正弦信號(hào)實(shí)行掃頻測(cè)量，采用頻譜分析儀或網(wǎng)絡(luò)分析儀進(jìn)行測(cè)量并直接得出系統(tǒng)函數(shù)

30、；2 加入沖激信號(hào)，采用FFT由輸出響應(yīng)計(jì)算輸出頻譜；輸出頻譜直接對(duì)應(yīng)于系統(tǒng)頻率響應(yīng)；第6章 數(shù)字濾波器理想濾波器與實(shí)際濾波器濾波器基本概念 系統(tǒng)的作用：將輸入頻譜改變?yōu)檩敵鲱l譜濾波：對(duì)于幅頻特性的改變 頻譜變形 頻率選擇能物理實(shí)現(xiàn)的濾波器的系數(shù)必定為實(shí)數(shù)，其幅頻響應(yīng)為偶，相頻響應(yīng)為奇；所以濾波器頻率響應(yīng)只需要在正頻率區(qū)間內(nèi)表達(dá)；名詞：通帶 阻帶 截止頻率 帶寬 低通 高通 帶通 帶阻 理想低通數(shù)字濾波器頻率響應(yīng) P.254 圖6.1理想低通；通帶內(nèi)：幅度為1， 線性相位（與時(shí)移對(duì)應(yīng)）；阻帶內(nèi)：幅度為0；沒有過渡帶；只容許低頻信號(hào)通過；通帶內(nèi)的信號(hào)通過時(shí)只發(fā)生延遲，不發(fā)生信號(hào)畸變；濾波器的實(shí)現(xiàn)

31、物理實(shí)現(xiàn)的必要條件因果系統(tǒng)：對(duì)于理想低通濾波器 理想低通濾波器為非因果系統(tǒng)，不能實(shí)現(xiàn)；只考慮幅頻特性時(shí)的逼近處理：1 采用長(zhǎng)延時(shí)，則在負(fù)時(shí)間區(qū)域內(nèi)，可以足夠接近于0；代價(jià)為系統(tǒng)延時(shí)增加；（圖6.4）2 采用截止頻率為的高階Butterworth濾波器；（圖6.5）實(shí)際濾波器參數(shù) （圖6.6）通帶波動(dòng) 阻帶波動(dòng) 過渡帶 通帶截止頻率 阻帶截止頻率群延遲范圍 MATLAB函數(shù)： gd,w=grpdelay (b,a,256) 對(duì)采用b，a向量描述的系統(tǒng)，在區(qū)間等分256點(diǎn)給出群延遲和對(duì)應(yīng)頻率； 容許范圍越小，濾波器結(jié)構(gòu)越復(fù)雜；對(duì)于FIR濾波器，總是具有線性相位，設(shè)計(jì)目標(biāo)為滿足指定的幅頻特性；對(duì)于I

32、IR濾波器，在設(shè)計(jì)時(shí)只能考慮滿足指定的幅頻特性，然后另行設(shè)計(jì)一個(gè)全通濾波器來改善其相頻特性；設(shè)計(jì)目標(biāo)：最簡(jiǎn)單的數(shù)字因果濾波器，滿足指定的幅頻特性系統(tǒng)函數(shù)為有理真分式，階數(shù)最小，系統(tǒng)穩(wěn)定；模擬濾波器基本概念P.260 圖6.8MATLAB函數(shù) H,w=freqz(bn,an) 對(duì)采用bn，an向量描述的系統(tǒng)，在區(qū)間等分512點(diǎn)給出頻率響應(yīng)和對(duì)應(yīng)頻率點(diǎn)； 一階數(shù)字濾波器單極點(diǎn)，無零點(diǎn)系統(tǒng)函數(shù) 系統(tǒng)頻率響應(yīng)(圖6.9 P. 262)討論：極點(diǎn)位置對(duì)于頻率響應(yīng)的影響a=0 全通，線性相位 全通濾波器（單位延遲元件）a0 低通a絕對(duì)值越大，變化越大；單極點(diǎn)，單零點(diǎn)系統(tǒng)函數(shù) 系統(tǒng)頻率響應(yīng)(圖6.10 P.

33、 263)討論：零點(diǎn)位置對(duì)于頻率響應(yīng)的影響b絕對(duì)值越大，頻率幅度越?。唬ㄗ饔门c極點(diǎn)相反）設(shè)計(jì)濾波器時(shí)，零點(diǎn)應(yīng)盡可能遠(yuǎn)離極點(diǎn)；低通濾波設(shè)計(jì)方案：帶寬確定，通帶截止頻率為 ： 通帶寬度（3-db帶寬） db：即： 結(jié)論： 帶寬等于極點(diǎn)到單位圓的距離例6.1 低通濾波 對(duì)連續(xù)時(shí)間信號(hào)的處理離散化：根據(jù)采樣定理選擇采樣周期T，避免頻率混疊； 頻率范圍： 歸一化：對(duì)頻率乘以T； 頻率范圍： 根據(jù)所需頻率選擇帶寬，確定極點(diǎn)值a；根據(jù)歸一化條件得出系統(tǒng)函數(shù) 利用MATLAB程序畫出輸入/輸出信號(hào)：6.2 (P.267) 圖6.13高通濾波設(shè)計(jì)方案 例6.2 高通濾波 對(duì)連續(xù)時(shí)間信號(hào)的處理步驟與例6.1相同

34、結(jié)果： 圖6.15數(shù)字高通濾波器的問題 由于采樣的限制，會(huì)引入對(duì)信號(hào)的低頻調(diào)制；這種調(diào)制效應(yīng)難以消除；二階數(shù)字濾波器雙極點(diǎn)，雙零點(diǎn)系統(tǒng)函數(shù)頻率響應(yīng)(P.270 圖6.16)特點(diǎn)：頻譜幅度隨極點(diǎn)到單位圓距離增大而減小，隨零點(diǎn)到單位圓距離增大而增大；極點(diǎn)靠近單位圓：出現(xiàn)幅頻峰值；零點(diǎn)靠近單位圓：出現(xiàn)幅頻谷值；適當(dāng)排布零點(diǎn)和極點(diǎn)的位置，可以得出不同的幅頻特性；低通濾波器 高通濾波器 帶通濾波器 數(shù)字諧振器 圖6.17例6.3 數(shù)字諧振器 連續(xù)信號(hào)的帶通濾波根據(jù)采樣定理選取采樣周期T；進(jìn)行頻率范圍的歸一化；確定帶寬范圍：上下截止頻率 ；利用標(biāo)準(zhǔn)公式求出歸一化常數(shù)：結(jié)果：P.273 圖6.18結(jié)果顯然優(yōu)

35、于例6.2的數(shù)字高通濾波器；阻塞濾波器零點(diǎn)位于，極點(diǎn)靠近零點(diǎn) 圖6.19例6.4 阻塞濾波器 濾除60Hz信號(hào)選取采樣周期T，得出歸一化頻率；帶入公式進(jìn)行設(shè)計(jì)；結(jié)果：圖6.20零點(diǎn)對(duì)濾波器特性的影響極點(diǎn)只能位于單位圓內(nèi)，但零點(diǎn)可以位于單位圓內(nèi)外的任意地方；零點(diǎn)與互為倒數(shù)根；它們分別位于單位圓內(nèi)外的對(duì)稱位置；在采用系統(tǒng)函數(shù)分析頻率特性時(shí)，將任何零點(diǎn)用其倒數(shù)根替代，其幅頻特性不受影響；因此，設(shè)計(jì)數(shù)字濾波器時(shí)，總可以假定所有零點(diǎn)都位于單位圓內(nèi)；將零點(diǎn)用其倒數(shù)根替代會(huì)改變系統(tǒng)的相頻特性；單位圓內(nèi)的零點(diǎn)稱為最小相位零點(diǎn)；若系統(tǒng)函數(shù)的零點(diǎn)全在單位圓內(nèi)，稱為最小相位系統(tǒng)函數(shù)；全通濾波器當(dāng)所有極點(diǎn)與零點(diǎn)互為倒

36、數(shù)根時(shí)，幅頻特性恒定為1；所有零點(diǎn)均位于單位圓外，相頻特性變化大；用于延遲單元或相位均衡器：提供延遲，改變系統(tǒng)相頻特性；梳狀濾波器阻塞濾波器的擴(kuò)展：具有等距阻塞點(diǎn)的濾波器；圖6.24正弦發(fā)生器在單位圓上有一個(gè)單極點(diǎn)：臨界穩(wěn)定，保持暫態(tài)； 第七章 FIR濾波器設(shè)計(jì)FIR濾波器基本特點(diǎn)沖激響應(yīng)長(zhǎng)度有限： 設(shè) ： ， ；N階FIR濾波器，長(zhǎng)度為N+1系統(tǒng)函數(shù) N個(gè)極點(diǎn)都位于原點(diǎn)，N個(gè)零點(diǎn)；頻率響應(yīng)計(jì)算MATLAB函數(shù) H,w=freqz (h,1) 給出區(qū)間上512個(gè)頻率點(diǎn)的頻率響應(yīng)設(shè)計(jì)目標(biāo)：求解有限沖激響應(yīng)，使其盡可能逼近要求的頻率響應(yīng)線性相位濾波器；圖7.1 P.295 對(duì)低通濾波器的設(shè)計(jì)要求設(shè)計(jì)思想對(duì)于絕大多數(shù)系統(tǒng)，希望實(shí)現(xiàn)的頻率響應(yīng)在區(qū)間具有奇偶對(duì)稱性；這種對(duì)稱性反映到實(shí)沖激響應(yīng)上，也可以使沖激響應(yīng)表達(dá)為偶序列或奇序列；對(duì)于穩(wěn)定的系統(tǒng)，當(dāng)N足夠大時(shí)，總有，因此可以將無限沖激響應(yīng)截?cái)酁橹缓琋項(xiàng)的有限沖激響應(yīng)，使得后者頻率響應(yīng)逼近于前者；若截?cái)鄥^(qū)間為對(duì)稱，則得到的有限沖激響應(yīng)具有奇偶對(duì)稱性，其頻率響應(yīng)可以表現(xiàn)為實(shí)函數(shù)；例：設(shè)無限沖激響應(yīng)序列為實(shí)偶序列；其頻率響應(yīng)可表示為