離散LTI系統(tǒng)的頻域分析

數(shù)字信號處理設計說明書離散LTI系統(tǒng)的頻域分析起止日期：2012年12月丄乙日至2012年12月二8日學生姓名尹斌班級電信科2班成績指導教師（簽字）計算機與信息工程學院2012年12月28日天津城市建設學院課程設計任務書2012—2013學年第1學期計算機與信息工程學院 電子信息工程系電子信息科學與技術專業(yè)課程設計名稱： 數(shù)字信號處理 設計題目： 離散LTI系統(tǒng)的頻域分析 完成期限：自2012年12月_17_日至2012年12月28日共_2_周設計依據(jù)、要求及主要內(nèi)容：一．課程設計依據(jù)在掌握系統(tǒng)函數(shù)、頻率響應和系統(tǒng)分析基本知識的基礎上。通過MATLAB仿真，加深對系統(tǒng)分析方法的理解；鍛煉運用所學知識，獨立分析問題、解決問題的綜合能力。二．課程設計內(nèi)容系統(tǒng)方程分別是y(n)—L2y(n—1)+0.72y(n—2)=x(n—1)y(n)—1.6y(n—1)+y(n—2)=x(n—1)y(n)—2y(n—1)+1.36y(n—2)=x(n—1)分別求出系統(tǒng)的系統(tǒng)函數(shù)，并利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。利用Matlab中的函數(shù)freqz計算出以上三個系統(tǒng)的頻率響應H(e加)，畫出幅頻特性與相頻特性。三．課程設計要求要求獨立完成設計任務。課程設計說明書封面格式要求見《天津城市建設學院課程設計教學工作規(guī)范》附表1課程設計的說明書要求簡潔、通順，計算正確，圖紙表達內(nèi)容完整、清楚、規(guī)范。運行程序，觀察并保存程序運行結(jié)果，能夠?qū)\行結(jié)果進行結(jié)果分析。課設說明書要求：說明題目的設計原理和思路、采用方法及程序。詳細說明調(diào)試方法和調(diào)試過程，并給程序加注釋。給出程序運行結(jié)果，并對其進行說明和分析。指導教師(簽字)： 系/教研室主任(簽字)： 批準日期：2012年12月13日目錄TOC\o"1-5"\h\z\o"CurrentDocument"第1章設計任務及要求 1\o"CurrentDocument"設計任務 1\o"CurrentDocument"課設要求 1\o"CurrentDocument"第2章設計原理 2\o"CurrentDocument"2.1設計原理 2\o"CurrentDocument"設計思路 2\o"CurrentDocument"第3章設計實現(xiàn) 3\o"CurrentDocument"3.1系統(tǒng)函數(shù) 3\o"CurrentDocument"設計程序 3\o"CurrentDocument"第4章設計結(jié)果及分析 5\o"CurrentDocument"第5章設計總結(jié) 8\o"CurrentDocument"參考文獻 9附錄 10第1章設計任務及要求1.1設計任務求出系統(tǒng)y(n)-1.2y(n-1)+0.72y(n-2)=x(n-1)的系統(tǒng)函數(shù)，并利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(e)，畫出幅頻特性與相頻特性。求出系統(tǒng)y(n)-1.6y(n-1)+y(n-2)=x(n-1)的系統(tǒng)函數(shù)，并利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(ej?)，畫出幅頻特性與相頻特性。求出系統(tǒng)y(n)-2y(n-1)+1.36y(n-2)=x(n-1)的系統(tǒng)函數(shù)，并利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(ejo)，畫出幅頻特性與相頻特性。課設要求說明題目的設計原理和思路、采用方法及程序。詳細說明調(diào)試方法和調(diào)試過程，并給程序加注釋。給出程序運行結(jié)果，并對其進行說明和分析。第2章設計原理2.1設計原理首先要了解系統(tǒng)函數(shù)的定義，用單位脈沖響應h(n)可以表示線性時不變離散系統(tǒng)，這時y(n)=x(n)*h(n)兩邊取z變換：Y(z)=X(z)H(z)。H(z)=Y(Z)/X(Z)則定義為系統(tǒng)函數(shù)。信號頻域分析是采用傅里葉變換將時域信號x(n)變換為X(f),從而幫助人們從另一個角度了解信號的額特征。在認識到系統(tǒng)函數(shù)的定義后再了解頻域分析的幾個結(jié)論：頻域分析法是一種圖解分析法，頻率特性是系統(tǒng)的一種數(shù)學模型。系統(tǒng)頻率特性的三種圖形為極坐標圖、對數(shù)頻率特性圖和對數(shù)幅相圖。若系統(tǒng)傳遞函數(shù)的極點和零點均位于z平面的左半平面，該系統(tǒng)為最小相位系統(tǒng)。反之，若系統(tǒng)的傳遞函數(shù)具有位于有半平面的零點或極點或有純滯后環(huán)節(jié)，則稱為非最小相位系統(tǒng)。利用Nyquist穩(wěn)定判據(jù)，可用開環(huán)頻率特性判別閉環(huán)系統(tǒng)穩(wěn)定性。同時可用相角裕量和幅角裕量來反映系統(tǒng)的相對穩(wěn)定性。設計思路我先就單位脈沖相應、零極點分布、系統(tǒng)頻率響應、幅頻和相頻在Matlab中運行的函數(shù)找到，然后查找相關資料編寫程序，然后進行修改，最后得出結(jié)果。下面就分別介紹了其中的幾個函數(shù)。離散LTI系統(tǒng)的頻率響應：Matlab提供的freqz函數(shù)可計算系統(tǒng)的頻率響應，其一般調(diào)用形式為H=freqz(b，a，w)式子中b和a分別為H(jw)分子分子多項式和分母多項式的系數(shù)向量；w為需要計算的H*(jw)的頻率采樣點向量。Matlab提供的zplane函數(shù)可以直接求解H(z)的零極點分布，其調(diào)用形式為：zplane(b，a)式子中：b和a分別為系統(tǒng)函數(shù)H(z)分子多項式和分母多項式的系數(shù)向量，該函數(shù)的作用是在z平面上畫出單位圓及系統(tǒng)的零點和極點。第3章設計實現(xiàn)3.1系統(tǒng)函數(shù)y(n)-1.2y(n-1)+0.72y(n-2)=x(n-1)求出其系統(tǒng)函數(shù)為H(z)=z-i/(1-1.2z-i+0.72z-2)。y(n)-1.6y(n-1)+y(n-2)=x(n-1)求出其系統(tǒng)函數(shù)為H(z)=z-i/(1-1.6z-i+z-2)。y(n)-2y(n-1)+1.36y(n-2)=x(n-1)求出其系統(tǒng)函數(shù)為H(z)=z-1/(1-2z-1+1.36z-2)。設計程序利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。源程序如下：A=[1,-1.2,0.72];B=[1];%設置系統(tǒng)函數(shù)向量A和Bsubplot(2,2,1);zplane(B,A);%繪制零極點圖hn=impz(B,A,58); %求系統(tǒng)單位脈沖響應h(n)n=0:length(hn)-1;subplot(2,2,2);y=hn;stem(n,hn.'.');%stab:%系統(tǒng)穩(wěn)定性判定函數(shù)disp('系統(tǒng)的極點為:')P=roots(A)%求H(z)的極點，并顯示disp('系統(tǒng)極點模的最大值為:')M=max(abs(P))%求所有極點模的最大值，并顯示ifM<1disp('系統(tǒng)穩(wěn)定'),else,disp('系統(tǒng)不穩(wěn)定'),end利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(ejo)，畫出幅頻特性與相頻特性。A=[1,-1.2,0.72];B=[1];%設置系統(tǒng)函數(shù)向量A和B[H,w]=freqz(B,A);%計算頻率響應subplot(2,2,1);plot(w/pi,abs(H));%繪制幅頻響應曲線xlabel('\omega/\pi');ylabel('IH(eAjA\omega)l');subplot(2,2,2);plot(w/pi,angle(H));%繪制相頻響應曲線xlabel('\omega/\pi');ylabel('phi(\omega)');利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。源程序如下：A=[1,-1.6,1];B=[1];%設置系統(tǒng)函數(shù)向量A和Bsubplot(2,2,1);zplane(B,A);%繪制零極點圖hn=impz(B,A,58); %求系統(tǒng)單位脈沖響應h(n)n=0:length(hn)-1;subplot(2,2,2);y=hn;stem(n,hn.'.');%stab:%系統(tǒng)穩(wěn)定性判定函數(shù)disp('系統(tǒng)的極點為:')P=roots(A)%求H(z)的極點，并顯示disp('系統(tǒng)極點模的最大值為:')M=max(abs(P))%求所有極點模的最大值，并顯示ifMvldisp('系統(tǒng)穩(wěn)定'),else,disp('系統(tǒng)不穩(wěn)定'),end,利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(ej?)，畫出幅頻特性與相頻特性。源程序如下：A=[1,-1.6,1];B=[1];%設置系統(tǒng)函數(shù)向量A和B[H,w]=freqz(B,A);%計算頻率響應subplot(2,2,1);plot(w/pi,abs(H));%繪制幅頻響應曲線xlabel('\omega/\pi');ylabel('IH(eAjA\omega)l');subplot(2,2,2);plot(w/pi,angle(H));%繪制相頻響應曲線xlabel('\omega/\pi');ylabel('phi(\omega)');利用Matlab畫出系統(tǒng)的零極點分布圖和相應的單位脈沖響應h(n)的波形，判斷系統(tǒng)的穩(wěn)定性。源程序如下：A=[1,-2,1.36];B=[1];%設置系統(tǒng)函數(shù)向量A和Bsubplot(2,2,1);zplane(B,A);%繪制零極點圖hn=impz(B,A,58); %求系統(tǒng)單位脈沖響應h(n)n=0:length(hn)-1;subplot(2,2,2);y=hn;stem(n,hn.'.');%stab:%系統(tǒng)穩(wěn)定性判定函數(shù)disp('系統(tǒng)的極點為:')P=roots(A)%求H(z)的極點，并顯示disp('系統(tǒng)極點模的最大值為:')M=max(abs(P))%求所有極點模的最大值，并顯示ifM<1disp('系統(tǒng)穩(wěn)定'),else,disp('系統(tǒng)不穩(wěn)定'),end利用Matlab中的函數(shù)freqz計算出系統(tǒng)的頻率響應H(ej?)，畫出幅頻特性與相頻特性。源程序如下：A=[1,-2,1.36];B=[1];%設置系統(tǒng)函數(shù)向量A和B[H,w]=freqz(B,A);%計算頻率響應subplot(2,2,1);plot(w/pi,abs(H));%繪制幅頻響應曲線xlabel('\omega/\pi');ylabel('|H(eAjA\omega)|');subplot(2,2,2);plot(w/pi,angle(H));%繪制相頻響應曲線xlabel('\omega/\pi');ylabel('phi(\omega)');

第4章設計結(jié)果及分析"WTlAfi 亠 _ 0AnEdkDuhigDflcktapWfadmuHa^D“上■(ft小J算Bf曹如““対叱■:創(chuàng)￡j亦臨必上|hw臥鋤出啊啊呦rCvntiiiL

也D-

r，"呷

.4r*ixE]Flgurt1F3-?Ec*t蘭申[njmTmJjDjesL-tapwiidnwH?lp□q?qCvntiiiL

也D-

r，"呷

.4r*ixE]Flgurt1F3-?Ec*t蘭申[njmTmJjDjesL-tapwiidnwH?lp□q?qan醱險心? □E簽口￡□-t-HJ-srE-4DU_I■jfTwiHdnirTMndWfcfcw?>[jP-L.ia.T2];B-[L];?ubplxii：聳芻j].ip)mlBj趟;hrpjrci團退人規(guī)；irfaLefiftkUmS-L,5Bfeplertt馮斗因;in旳IdW<〒辿的ffi±Aia]F-XK-tS^I喩黔段爲隹的■■留l:過gl>〕)l￡BCLdaxpL*￡■注霸丘']』-Le.fthzpi：'岸開■工尋宅')..nrilP■DLflMD+DLMDDliDL4MD-D.HDD1電圧整白■的總±1■対i眉￡Ui dL^pl號湖用占有⑴>.ihc1戸>〕ii?LdLipI：' >B■L?,dirp[-^^7 '].,?nd話尊極點分布圖，單位脈沖響應及穩(wěn)定性判斷話尊極點分布圖，單位脈沖響應及穩(wěn)定性判斷予十■il- _ 一 * 亠. 一 —AtBdi:蘭kMrr&artreei(D^chupWr-ric*-Hdp >D?U*b甌&凸?4?□□?；胤l響應和相頻響應曲線4kMATLAB軒EdEDtbug 輕如Hdp非rfjfrJ&?iD*iri=<vl怦2尉_J」住je?jrti?rr^fedijwl=-]Rjre]..un?E?ArCLj-L.flJ];B-[L1;sijbplptUPSfD;:A）.e?jrti?rr^fedijwl=-]Rjre]..un?E?ArCLj-L.flJ];B-[L1;sijbplptUPSfD;:A）.Lt=twf：iB怎殖;rpDilitfiiftliEhn'!-];fiut<>ltd（2p基黔；Khn:玳旺快bna"I'i.S?rsb6：di;pC壬址的LE?點為i）P=F3iXi：fll°l°ldupl：'吊駛補百腥取?立昭為「）rtfruIzbai；F〕）討radixpC戦科宅h.isb沖]岸4不孫竜秫的M張ij.NM+D.flMDi0.MM-D.QMDi辛詭祝樸HI的ON■対—P=ra&i：—P=ra&i：ial^亦町I?鼻弭慢占誓則玄1[為:?〕V-aiu<aKcIP）>itKLdip叮爭站■電）?■XqtLgr離繇?gh-nd零極點分布圖，單位脈沖響應及系統(tǒng)穩(wěn)定性判斷幅頻響應和相頻響應曲線*WT1AS-X_ 旦&EdfDabi^D^kl^pWrdOh-申BD呼總■.由 書Id?*6?2[!ID云TL唧wwfc創(chuàng)』囪awlcij-dkTffk!AddTW5H?r#fe?UYffl日出田田田田田田田田EdfView Taohip4^hEe0>豈ndew申?fe?UYffl日出田田田田田田田田EdfView Taohip4^hEe0>豈ndew申?DOh f)?祖□□it□RflalPstLWQ-嚴-9匸aH三一?A-[L-E；].?];>[□;jijtplirfUj2jL).:p]wfflj的；遠1平iBjAf6SJ;Yi=[l|liMiffthChii>-1；futfiliyt閔2j2);j=]WL；i!ri:ei^M.■「)；MtsbsdLivCf-LCiiXtSI如dL沖〔'耳牛出占l￡陽■衣町切')Hkkf(atd1刃)3<K]dirp(初橢7ET』ii】mHLip[KWFBn1)□nd-疏的枇點為：F■I.DQM+■0.60GD1L.OOIM-O.dOOOi斥籠105511的R天■為IJ.LHZ-*srflcr叭r聽般旳-、dt-pi： 占羯郎ilHdT力J〕^BUi；lfar>>：l)■ifKLdaxpC ■Lkipdup['f1rT'>Pmd. A .. A .A2///\\\-x、一一-■2零極點分布圖，單位脈沖響應及系統(tǒng)穩(wěn)定性判斷旳隹umL~Ug?皿1rts-Edc?葉kb4rrTeels匕卑卜如WfriswH?lp1D*日』k負召凸?4d□ □幅頻響應和相頻響應曲線第5章設計總結(jié)這學期我們學習了數(shù)字信號處理這門課程并開設了實驗，同時在短學期中進行了為期兩周的課程設計，目的是讓我們更深刻的理解數(shù)字信號處理這門課，并且把所學的理論與實踐相聯(lián)系。我拿到題目后，上網(wǎng)查閱老師給的參考資料，并查找了相關的書籍，了解了離散LTI系統(tǒng)的頻域分析的原理。參考資料所給內(nèi)容比較簡單且Matlab的程序有缺省和錯誤，所以我又學習數(shù)字信號處理中相關知識并查閱一些資料，學習Matlab的使用方法，最終寫出程序，得到了相應的結(jié)果，并完成報告。這次課程設計，我學會了很多東西。熟悉了仿真軟件Matlab的使用方法，學會了Matlab程序的編寫，完成仿真。我還意識到英文水平的重要性，我們學習電信科專業(yè)，會需要了解國外先進的技術，也需要閱讀經(jīng)典書籍和各種文獻、資料，因此需要扎實的英語閱讀和翻譯功底，尤其是相關的專業(yè)詞匯。此外，這次課程設