Simulink系統(tǒng)仿真課程設(shè)計(jì)

1、信息系統(tǒng)仿真課程設(shè)計(jì)課程設(shè)計(jì)報(bào)告題 目： 信息系統(tǒng)課程設(shè)計(jì)仿真 院 （系）： 信息科學(xué)與技術(shù)工程學(xué)院 專業(yè)班級(jí)： 通信工程 1003 學(xué)生姓名： 學(xué) 號(hào)： 指導(dǎo)教師： 吳莉 朱忠敏 2012 年 1 月 14 日至2012 年 1 月 25 日華中科技大學(xué)武昌分校制信息系統(tǒng)仿真課程設(shè)計(jì)任務(wù)書一、設(shè)計(jì)（調(diào)查報(bào)告/論文）題目MATLAB仿真設(shè)計(jì)：（1）自編程序?qū)崿F(xiàn)動(dòng)態(tài)卷積；（2）用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器。Simulink仿真設(shè)計(jì)：（1）低通抽樣定理的實(shí)現(xiàn)；（2）抽樣量化編碼器的設(shè)計(jì)。二、設(shè)計(jì)（調(diào)查報(bào)告/論文）主要內(nèi)容MATLAB仿真設(shè)計(jì)：（1）自編程序?qū)崿F(xiàn)動(dòng)態(tài)卷積：動(dòng)態(tài)演示序列f1=u(

2、n)(0n11)和f2=(0n16)的卷積和的過(guò)程。（2）用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器：采用雙線性變換法設(shè)計(jì)一個(gè)巴特沃斯數(shù)字低通濾波器，要求：wp0.2，Rp1 dB；ws0.3，As15 dB，濾波器采樣頻率Fs100 Hz。自編程序?qū)崿F(xiàn)上述濾波器的設(shè)計(jì)，滿足各個(gè)參數(shù)要求。Simulink仿真設(shè)計(jì)：（1）低通抽樣定理的實(shí)現(xiàn)：建立Simulink動(dòng)態(tài)仿真模型，設(shè)置模塊參數(shù)與系統(tǒng)仿真參數(shù)，進(jìn)行系統(tǒng)仿真。對(duì)仿真結(jié)果進(jìn)行調(diào)試、分析，得出結(jié)論。（2）抽樣量化編碼器的設(shè)計(jì)：設(shè)計(jì)一個(gè)PCM編碼器，當(dāng)輸入一個(gè)電平，通過(guò)該編碼器得到相應(yīng)的編碼輸出。三、原始資料電子課件，仿真參考例子。四、要求的設(shè)計(jì)（調(diào)查

3、/論文）成果對(duì)于MATLAB仿真設(shè)計(jì)部分：（1）自編卷積程序要出現(xiàn)動(dòng)態(tài)卷積的過(guò)程，并且結(jié)果正確；（2）雙線性變換法設(shè)計(jì)的IIR濾波器應(yīng)滿足參數(shù)要求，設(shè)計(jì)結(jié)果合理、正確。對(duì)于Simulink仿真設(shè)計(jì)部分：（1） 建立Simulink仿真模型完成低通抽樣定理的實(shí)現(xiàn)；（2） 建立Simulink仿真模型仿真完成PCM編碼器的設(shè)計(jì)，調(diào)試。當(dāng)輸入一個(gè)電平，應(yīng)得出正確的編碼輸出。五、進(jìn)程安排講述MATLAB語(yǔ)言及Simulink仿真技術(shù)。（2.5天）布置課設(shè)任務(wù)。（0.5天）Simulink仿真設(shè)計(jì)。（2.5天）MATLAB仿真設(shè)計(jì)。（2.5天）驗(yàn)收課程設(shè)計(jì)結(jié)果。（1天）撰寫課程設(shè)計(jì)報(bào)告。（1天）六、主要

4、參考資料1 李賀冰等.Simulink通信仿真教程.北京：國(guó)防工業(yè)出版社，2006.2 劉舒帆等.數(shù)字信號(hào)處理實(shí)驗(yàn)（MATLAB版）.西安：西安電子科技大學(xué)出版社，2008.3 黨宏社.信號(hào)與系統(tǒng)實(shí)驗(yàn)（MATLAB版）.西安：西安電子科技大學(xué)出版社，2007.4 郭文彬等.通信原理基于Matlab的計(jì)算機(jī)仿真.北京：北京郵電大學(xué)出版社，2006. 指導(dǎo)教師（簽名）： 20 年 月 日 目錄 摘要 .5一、 Simulink仿真設(shè)計(jì) .6 1.1 低通抽樣定理 .6 1.2 抽樣量化編碼 .9二、MATLAB仿真設(shè)計(jì) .12 2.1、自編程序?qū)崿F(xiàn)動(dòng)態(tài)卷積 .12 2.1.1 編程分析 .12 2

5、.1.2自編matlab程序： .13 2.1.3 仿真圖形 .132.1.4仿真結(jié)果分析 .15 2.2用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器 .15 2.2.1雙線性變換法的基本知識(shí) .15 2.2.2采用雙線性變換法設(shè)計(jì)一個(gè)巴特沃斯數(shù)字低通濾波器 .16 2.2.3自編matlab程序 .16 2.2.4 仿真波形 .17 2.2.5仿真結(jié)果分析 .17三、總結(jié) .19四、參考文獻(xiàn) .19五、課程設(shè)計(jì)成績(jī) .20摘要 Matlab 是一種廣泛應(yīng)用于工程設(shè)計(jì)及數(shù)值分析領(lǐng)域的高級(jí)仿真平臺(tái)。它功能強(qiáng)大、簡(jiǎn)單易學(xué)、編程效率高，目前已發(fā)展成為由MATLAB語(yǔ)言、MATLAB工作環(huán)境、MATLAB圖形處

6、理系統(tǒng)、MATLAB數(shù)學(xué)函數(shù)庫(kù)和MATLAB應(yīng)用程序接口五大部分組成的集數(shù)值計(jì)算、圖形處理、程序開發(fā)為一體的功能強(qiáng)大的系統(tǒng)。本次課程設(shè)計(jì)主要包括MATLAB和SIMULINKL 兩個(gè)部分。首先利用SIMULINKL 實(shí)現(xiàn)了連續(xù)信號(hào)的采樣及重構(gòu)，通過(guò)改變抽樣頻率來(lái)實(shí)現(xiàn)過(guò)采樣、等采樣、欠采樣三種情況來(lái)驗(yàn)證低通抽樣定理，繪出原始信號(hào)、采樣信號(hào)、重構(gòu)信號(hào)的時(shí)域波形圖。然后利用SIMULINKL 實(shí)現(xiàn)抽樣量化編碼，首先用一連續(xù)信號(hào)通過(guò)一個(gè)抽樣量化編碼器按照A律13折線進(jìn)量化行，觀察其產(chǎn)生的量化誤差，其次利用折線近似的PCM編碼器對(duì)一連續(xù)信號(hào)進(jìn)行編碼。最后利用MATLAB進(jìn)行仿真設(shè)計(jì)，通過(guò)編程，在編程環(huán)境

7、中對(duì)程序進(jìn)行調(diào)試，實(shí)現(xiàn)動(dòng)態(tài)卷積以及雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器。 本次課程設(shè)計(jì)加深理解和鞏固通信原理、數(shù)字信號(hào)處理課上所學(xué)的有關(guān)基本概念、基本理論和基本方法，并鍛煉分析問(wèn)題和解決問(wèn)題的能力。一、Simulink仿真設(shè)計(jì)1.1 低通抽樣定理輸入信號(hào)為頻率為10Hz的正弦波，觀察對(duì)于同一輸入信號(hào)在不同的抽樣頻率時(shí)，恢復(fù)信號(hào)的不同波形形態(tài)。(1) 當(dāng)抽樣頻率大于信號(hào)頻率的兩倍，（eg:30Hz） 建立模型； 圖1-1 抽樣仿真框圖 圖1-3 參數(shù)設(shè)置圖1-2 F=30HZ抽樣及恢復(fù)波形 (2) 改變抽樣頻率為信號(hào)頻率的兩（eg:20Hz）圖1-4 F=20HZ的抽樣及恢復(fù)波形 (3) 改變抽樣頻

8、率小于信號(hào)頻率的兩倍（eg:5Hz） 圖1-5 F=5HZ的抽樣及恢復(fù)波形 分析：1、由仿真圖可知當(dāng)抽樣頻率大于或等于被抽樣信號(hào)頻率的2倍時(shí)，抽樣輸出通過(guò)模擬低通濾波器能夠恢復(fù)出被抽樣信號(hào)；當(dāng)抽樣頻率小于被抽樣信號(hào)頻率的2 倍，模擬低通濾波器的輸出波形的形狀已失真，即不能恢復(fù)出原始信號(hào)。 2、驗(yàn)證抽樣定理，即一個(gè)頻帶限制在0f內(nèi)的連續(xù)信號(hào)m(t)，如果取樣速率Fs大于或等于2f，則可以由抽樣值序列m(nTs)無(wú)失真地重建原始信號(hào)m(t)。即一個(gè)頻帶限制在（0，）內(nèi)的時(shí)間連續(xù)信號(hào)，如果以T秒的間隔對(duì)它進(jìn)行等間隔抽樣，則將被所得到的抽樣值完全確定。3、采樣信號(hào)通過(guò)模擬低通濾波器后，將其截止頻率設(shè)置

9、為80HZ，為了接近理想低通濾波器，將其濾波階數(shù)設(shè)置的高一點(diǎn)，信號(hào)通過(guò)低通濾波器后，除了被濾除阻帶上的高頻分量外，幅度還會(huì)有衰減，因此還需通過(guò)一個(gè)增益放大器。4、低通巴特沃茲濾波器的階次要選擇適宜，階次過(guò)高或過(guò)，輸出的信號(hào)都有失真另外，低截止頻率不能過(guò)大，否則會(huì)使多余的諧波通過(guò)，也會(huì)導(dǎo)致復(fù)原的信號(hào)失真。5、當(dāng)輸入的連續(xù)信號(hào)仿真時(shí)長(zhǎng)也會(huì)影響信號(hào)的波形，仿真是長(zhǎng)越短越好，輸出的波形越好。由于計(jì)算機(jī)處理的信號(hào)本質(zhì)上不是模擬信號(hào)，為了更接近模擬信號(hào)，就將仿真時(shí)長(zhǎng)盡可能的減小，但不可能減少到零，減少到零后的波形也有一定程度上的失真。1.2 抽樣量化編碼 1） 用一個(gè)正弦信號(hào)通過(guò)一個(gè)抽樣量化編碼器后按照A

10、律13折線產(chǎn)生量化輸出信號(hào)，從示波器（Scope）上可以觀察到產(chǎn)生的量化誤差。正弦信號(hào)為幅度為1，頻率為1Hz的連續(xù)時(shí)間信號(hào)。 建立模型； 圖1-6 抽樣量化編碼仿真模塊框 圖1-7 參數(shù)設(shè)置 圖1-8 參數(shù)設(shè)置 圖1-9 量化波形 圖1-10 量化誤差分析：1、正弦信號(hào)通過(guò)一個(gè)抽樣量化編碼器后按照A律13折線產(chǎn)生量化輸出信號(hào)，從示波器（scope）上可以觀察到產(chǎn)生的量化誤差，為了比較量化之前和之后的正弦信號(hào)，正弦信號(hào)產(chǎn)生器和抽樣量化編碼器的第2個(gè)輸出端口的輸出信號(hào)通過(guò)一個(gè)復(fù)用器連接到示波器1（scope1）2、抽樣量化編碼器用于產(chǎn)生1A律3折線，它把正弦信號(hào)產(chǎn)生器產(chǎn)生的正弦信號(hào)轉(zhuǎn)換成量化信號(hào)

11、，并且計(jì)算這個(gè)過(guò)程中產(chǎn)生的量化噪聲。抽樣量化編碼的參數(shù)設(shè)置如下：quantization partition設(shè)置為 -1/2 -1/4 -1/8 -1/16 -1/32 -1/64 -1/128 0 1/128 1/64 1/32 1/16 1/8 1/4 1/2 1 3、示波器運(yùn)行結(jié)果，圖中平滑曲線是表示抽樣之前的正弦信號(hào)，折線表示通過(guò)抽樣量化編碼之后的信號(hào)?？梢钥吹?，抽樣量化之后的信號(hào)與原來(lái)的連續(xù)信號(hào)之間存在一定的量化誤差，示波器1運(yùn)行結(jié)果，從中可以看出抽樣量化產(chǎn)生的量化誤差。1.3 設(shè)計(jì)一個(gè)13折線近似的PCM編碼器，設(shè)計(jì)一個(gè)13折線近似的PCM編碼器使它能夠?qū)θ≈翟?1，1內(nèi)的歸一化信

12、號(hào)樣值進(jìn)行編碼。（eg：當(dāng)輸入為843，輸出編碼應(yīng)為）圖1-11 13折線PCM編碼器模型模型 圖1-12 常數(shù)發(fā)生器模塊參數(shù) 圖1-13限幅器模塊參數(shù) 圖1-14 查表模塊參數(shù) 圖1-15 增益參數(shù) 圖1-16 延時(shí)器參數(shù) 圖1-17 絕對(duì)值參數(shù) 圖1-18 量化器參數(shù) 圖1-19 十進(jìn)制轉(zhuǎn)二進(jìn)制參數(shù) 分析：1由常數(shù)發(fā)生器模塊產(chǎn)生一個(gè)對(duì)2048個(gè)量化臺(tái)階的歸一化的電壓，連續(xù)信號(hào)通過(guò)限幅器將輸入信號(hào)的幅度限制在PCM編碼定義的范圍-1 ，1內(nèi)，將延遲的門限設(shè)置為0，用來(lái)判斷輸入電平的極性，其輸出作為PCM編碼的極性碼。樣值取絕對(duì)值后，通過(guò)查表模塊進(jìn)行13折線壓縮，由于以取絕對(duì)值，表格中的值只取

13、13折線第一象限中的九個(gè)分段點(diǎn)的坐標(biāo)，查表方法為內(nèi)插法，以便能精確近似，然后用增益模塊將樣值放大128倍，用間距為1 的Quantier 進(jìn)行四舍五入量化，最后將整數(shù)編碼為7位二進(jìn)制序列，作為PCM編碼的低7位。2、在設(shè)置增益時(shí)，增益為127或128時(shí)輸出的量化電平正確，二、MATLAB仿真設(shè)計(jì)2.1、自編程序?qū)崿F(xiàn)動(dòng)態(tài)卷積2.1.1 編程分析1）掌握離散卷積過(guò)程：序列倒置 移位 相乘 求和；2）自編程序?qū)崿F(xiàn)兩個(gè)序列的動(dòng)態(tài)卷積。題目：動(dòng)態(tài)演示序列f1=u(n)(0n11)和f2=(0n16)的卷積和的過(guò)程。2.1.2自編matlab程序：clear; clc; %清除數(shù)據(jù)區(qū)，指令區(qū)nf1=1:1

14、0;lf1=length(nf1); %f1序列長(zhǎng)度f(wàn)1=ones(1,lf1); %f1序列nf2=1:15;lf2=length(nf2); %f2長(zhǎng)度f(wàn)2=exp(-0.1.*nf2); %f2序列f3=fliplr(f2); %f2翻轉(zhuǎn)h=zeros(1,lf2+1),f2,zeros(1,lf1-1); %f2補(bǔ)零u=zeros(1,lf2+1),f1,zeros(1,lf2-1); %f1補(bǔ)零nf3=-lf2:lf1+lf2-1; %求卷積for k=0:lf1+lf2 %卷積的動(dòng)態(tài)顯示過(guò)程p=zeros(1,k),f3,zeros(1,lf1+lf2-k); %序列補(bǔ)零并且進(jìn)行循

15、環(huán)移位y1=p.*u; %求乘積yk=sum(y1); %各點(diǎn)之和，即一個(gè)點(diǎn)卷積結(jié)果y(k+1)=yk; %將卷積結(jié)果存入y序列z=zeros(1,lf2-1),y,zeros(1,lf2+lf1-k); %y序列進(jìn)行補(bǔ)零后的卷積序列 % 繪圖subplot(5,1,1);stem(nf3,u,filled);title(u(n);ylabel(u(n);subplot(5,1,2);stem(nf3,h,filled);title(h(n);ylabel(h(n);subplot(5,1,3);stem(nf3,p,filled);title(翻轉(zhuǎn)移位);ylabel(h(n-k);subp

16、lot(5,1,4);stem(nf3,y1,filled);title(乘積);ylabel(u(n)*h(n);subplot(5,1,5);stem(nf3,z,filled);title(卷積結(jié)果);ylabel(y(n);if k=0 disp(開始卷積),pause;EndIf k=1f1+1 disp/(暫停),pause;Else pause(0.5)Endend2.1.3 仿真圖形圖2-2 卷積動(dòng)態(tài)2.1.4仿真結(jié)果分析:1、 卷積后的結(jié)果0 0 0.90484 1.7236 2.4644 3.1347 3.7412 4.29 4.7866 5.236 5.6425 6.01

17、04 5.4384 4.9209 3.6455 3.0967 2.6001 2.1508 1.7442 1.3763 1.0435 0/74226 0.46973 0.22313 2、 u（n）的長(zhǎng)度為10，h（n）的長(zhǎng)度為15，卷積后的長(zhǎng)度Ny=Nu+Nh-1=24.有圖分析知，y（n）的長(zhǎng)度為24點(diǎn)，滿足卷積定理。2.2用雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器 2.2.1雙線性變換法的基本知識(shí)1）雙線性變換法是將整個(gè)s平面映射到整個(gè)z平面，其映射關(guān)系為雙線性變換法克服了脈沖響應(yīng)不變法從s平面到z平面的多值映射的缺點(diǎn)，消除了頻譜混疊現(xiàn)象。但其在變換過(guò)程中產(chǎn)生了非線性的畸變，在設(shè)計(jì)IIR數(shù)字濾波器的

18、過(guò)程中需要進(jìn)行一定的預(yù)修正。2）設(shè)計(jì)中有關(guān)雙線性變換法的子函數(shù)。Bilinear功能：雙線性變換將s域(模擬域)映射到z域(數(shù)字域)的標(biāo)準(zhǔn)方法，將模擬濾波器變換成離散等效濾波器。調(diào)用格式：numd，dendbilinear(num，den，F(xiàn)s)；將模擬域傳遞函數(shù)變換為數(shù)字域傳遞函數(shù)，F(xiàn)s為取樣頻率。numd，dendbilinear(num，den，F(xiàn)s，F(xiàn)p)；將模擬域傳遞函數(shù)變換為數(shù)字域傳遞函數(shù)，F(xiàn)s為取樣頻率，F(xiàn)p為通帶截止頻率。zd，pd，kdbilinear(z，p，k，F(xiàn)s)；將模擬域零極點(diǎn)增益系數(shù)變換到數(shù)字域，F(xiàn)s為取樣頻率。zd，pd，kdbilinear(z，p，k，F(xiàn)s，

19、Fp)；將模擬域零極點(diǎn)增益系數(shù)變換到數(shù)字域，F(xiàn)s為取樣頻率，F(xiàn)p為通帶截止頻率。Ad，Bd，Cd，Ddbilinear(A，B，C，D，F(xiàn)s)；將模擬域狀態(tài)變量系數(shù)變換到數(shù)字域，F(xiàn)s為取樣頻率。3）用MATLAB雙線性變換法進(jìn)行IIR數(shù)字濾波器設(shè)計(jì)的步驟：(1)輸入給定的數(shù)字濾波器設(shè)計(jì)指標(biāo)；(2)根據(jù)公式進(jìn)行預(yù)修正，將數(shù)字濾波器指標(biāo)轉(zhuǎn)換成模擬濾波器設(shè)計(jì)指標(biāo)；(3)確定模擬濾波器的最小階數(shù)和截止頻率；(4)計(jì)算模擬低通原型濾波器的系統(tǒng)傳遞函數(shù)；(5)利用模擬域頻率變換法，求解實(shí)際模擬濾波器的系統(tǒng)傳遞函數(shù)；(6)用雙線性變換法將模擬濾波器轉(zhuǎn)換為數(shù)字濾波器。2.2.2采用雙線性變換法設(shè)計(jì)一個(gè)巴特沃斯

20、數(shù)字低通濾波器，要求：wp0.2，Rp1 dB；ws0.3，As15 dB，濾波器采樣頻率Fs100 Hz。2.2.3自編matlab程序wp=0.2*pi; %濾波器的通帶截止頻率ws=0.3*pi; %濾波器的阻帶截止頻率Rp=1;As=15; %輸入濾波器的通阻帶衰減指標(biāo)ripple=10(-Rp/20);Attn=10(-As/20); Fs=100;T=1/Fs;Omgp=(2/T)*tan(wp/2); %原型通帶頻率預(yù)修正 Omgs=(2/T)*tan(ws/2); %原型阻帶頻率預(yù)修正 n,Omgc=buttord(Omgp,Omgs,Rp,As,s);%計(jì)算階數(shù)n和3dB截止

21、頻率z0,p0,k0=buttap(n); %歸一化原型設(shè)計(jì)ba=k0*real(poly(z0); %求原型濾波器系數(shù)baa=real(poly(p0); %求原型濾波器系數(shù)aba1,aa1=lp2lp(ba,aa,Omgc);%變換為模擬低通濾波器系數(shù)b,abd,ad=bilinear(ba1,aa1,Fs) %用雙線性變換法求數(shù)字濾波器系數(shù)b,asos,g=tf2sos(bd,ad); %由直接型變換為級(jí)聯(lián)型H,w=freqz(bd,ad);dbH=20*log10(abs(H)+eps)/max(abs(H);%化為分貝值subplot(2,2,1),plot(w/pi,abs(H);

22、ylabel(|H|);title(幅度響應(yīng));axis(0,1,0,1.1);set(gca,XTickMode,manual,XTick,0,0.25,0.4,1); set(gca,YTickMode,manual,YTick,0,Attn,ripple,1);gridsubplot(2,2,2),plot(w/pi,angle(H)/pi);ylabel(phi);title(相位響應(yīng));axis(0,1,-1,1);set(gca,XTickMode,manual,XTick,0,0.25,0.4,1);set(gca,YTickMode,manual,YTick,-1,0,1);g

23、ridsubplot(2,2,3),plot(w/pi,dbH);title(幅度響應(yīng)(dB);ylabel(dB);xlabel(頻率(pi);axis(0,1,-40,5);set(gca,XTickMode,manual,XTick,0,0.25,0.4,1);set(gca,YTickMode,manual,YTick,-50,-15,-1,0);gridsubplot(2,2,4),zplane(bd,ad);2.2.4 仿真波形2.2.5仿真結(jié)果分析1、濾波器階數(shù)n=6, 3dB截止頻率Omgc=76.623.分子系數(shù)bd = 0.0007 0.0044 0.0111 0.0148

24、 0.0111 0.0044 0.0007分母系數(shù)ad 1.0000 -3.1836 4.6222 -3.7795 1.8136 -0.4800 0.0544二階環(huán)節(jié)級(jí)聯(lián)結(jié)構(gòu)的分子分母系數(shù)矩陣sosSos=1 2.0004 1 1 -0.90437 0.21552 1 2.018 1.0184 1 -1.0106 0.35827 1 1.9816 0.98194 1 -1.2686 0.70513增益g=0.2、 由頻率特性曲線可知，設(shè)計(jì)結(jié)果在通帶截止頻率處滿足Rp1dB ，As15dB 的設(shè)計(jì)指標(biāo)要求，系統(tǒng)的極點(diǎn)全部在單位圓內(nèi)，是一個(gè)穩(wěn)定的系統(tǒng)，由N=6可知，設(shè)計(jì)的巴特沃斯數(shù)字低通濾波器是一

25、個(gè)6階的系統(tǒng)，原型H（s）在s=-處有六個(gè)零點(diǎn)，映射到z=-1處。3、 由幅度相應(yīng)圖可看出，通帶有富余量，這是由于用最小阻帶頻率指標(biāo)計(jì)算的結(jié)果。4、計(jì)算過(guò)程：Wp=0.2 Rp=1dB Ws=0.3 As=15dB Omgp=(2/T)*tan(wp/2)=200tan0.1=64.9839Omgs=(2/T)*tan(ws/2)=200tan0.15=101.9051sp= Omgs/Omgp=1.568Ksp= =10.875N=lgksp/lgsp=5.3056N取6Omgc=Omgs（100.1as-1)-1/2N=76.6229Gp=1/(p2+0.5176p+1)(p2+1.414

26、2p+1)(p2+1.9319p+1)P=S/Omgc s=2/T(1-Z-1)(1+Z1)H(z)= ( 0.0007+ 0.0044z-1+ 0.0111z-2+ 0.0148z -3+ 0.0111z-4+ 0.0044z-5 +0.0007z-6 )/(1.0000 -3.1836z-1 +4.6222z-2 -3.7795z -3 +1.8136z-4 -0.4800z-5 +0.0544z-6)本題所設(shè)計(jì)的濾波器是一個(gè)六階濾波器，由于MATLAB精度有限，表現(xiàn)的sos有誤差，把它們調(diào)整到精確值后，濾波器的級(jí)連型系統(tǒng)函數(shù)應(yīng)該為H（z）=0.（1+z-1)-6 /（1-0.90437 z-10.21552z-2）(1-1.0106 z-1+0.35827 z-2)(1 -1.2686 z-1+0