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

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

第4章設(shè)計(jì)結(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醱險(xiǎn)心? □E簽口￡□-t-HJ-srE-4DU_I■jfTwiHdnirTMndWfcfcw?>[jP-L.ia.T2];B-[L];?ubplxii：聳芻j].ip)mlBj趟;hrpjrci團(tuán)退人規(guī)；irfaLefiftkUmS-L,5Bfeplertt馮斗因;in旳IdW<〒辿的ffi±Aia]F-XK-tS^I喩黔段爲(wèi)隹的■■留l:過gl>〕)l￡BCLdaxpL*￡■注霸丘']』-Le.fthzpi：'岸開■工尋宅')..nrilP■DLflMD+DLMDDliDL4MD-D.HDD1電圧整白■的總±1■対i眉￡Ui dL^pl號(hào)湖用占有⑴>.ihc1戸>〕ii?LdLipI：' >B■L?,dirp[-^^7 '].,?nd話尊極點(diǎn)分布圖，單位脈沖響應(yīng)及穩(wěn)定性判斷話尊極點(diǎn)分布圖，單位脈沖響應(yīng)及穩(wěn)定性判斷予十■il- _ 一 * 亠. 一 —AtBdi:蘭kMrr&artreei(D^chupWr-ric*-Hdp >D?U*b甌&凸?4?□□?；胤l響應(yīng)和相頻響應(yīng)曲線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?點(diǎn)為i）P=F3iXi：fll°l°ldupl：'吊駛補(bǔ)百腥取?立昭為「）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叮爭(zhēng)站■電）?■XqtLgr離繇?gh-nd零極點(diǎn)分布圖，單位脈沖響應(yīng)及系統(tǒng)穩(wěn)定性判斷幅頻響應(yīng)和相頻響應(yīng)曲線*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-嚴(yán)-9匸aH三一?A-[L-E；].?];>[□;jijtplirfUj2jL).:p]wfflj的；遠(yuǎn)1平iBjAf6SJ;Yi=[l|liMiffthChii>-1；futfiliyt閔2j2);j=]WL；i!ri:ei^M.■「)；MtsbsdLivCf-LCiiXtSI如dL沖〔'耳牛出占l￡陽(yáng)■衣町切')Hkkf(atd1刃)3<K]dirp(初橢7ET』ii】mHLip[KWFBn1)□nd-疏的枇點(diǎn)為：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零極點(diǎn)分布圖，單位脈沖響應(yīng)及系統(tǒng)穩(wěn)定性判斷旳隹umL~Ug?皿1rts-Edc?葉kb4rrTeels匕卑卜如WfriswH?lp1D*日』k負(fù)召凸?4d□ □幅頻響應(yīng)和相頻響應(yīng)曲線第5章設(shè)計(jì)總結(jié)這學(xué)期我們學(xué)習(xí)了數(shù)字信號(hào)處理這門課程并開設(shè)了實(shí)驗(yàn)，同時(shí)在短學(xué)期中進(jìn)行了為期兩周的課程設(shè)計(jì)，目的是讓我們更深刻的理解數(shù)字信號(hào)處理這門課，并且把所學(xué)的理論與實(shí)踐相聯(lián)系。我拿到題目后，上網(wǎng)查閱老師給的參考資料，并查找了相關(guān)的書籍，了解了離散LTI系統(tǒng)的頻域分析的原理。參考資料所給內(nèi)容比較簡(jiǎn)單且Matlab的程序有缺省和錯(cuò)誤，所以我又學(xué)習(xí)數(shù)字信號(hào)處理中相關(guān)知識(shí)并查閱一些資料，學(xué)習(xí)Matlab的使用方法，最終寫出程序，得到了相應(yīng)的結(jié)果，并完成報(bào)告。這次課程設(shè)計(jì)，我學(xué)會(huì)了很多東西。熟悉了仿真軟件Matlab的使用方法，學(xué)會(huì)了Matlab程序的編寫，完成仿真。我還意識(shí)到英文水平的重要性，我們學(xué)習(xí)電信科專業(yè)，會(huì)需要了解國(guó)外先進(jìn)的技術(shù)，也需要閱讀經(jīng)典書籍和各種文獻(xiàn)、資料，因此需要扎實(shí)的英語閱讀和翻譯功底，尤其是相關(guān)的專業(yè)詞匯。此外，這次課程設(shè)