窗函數(shù)數(shù)字信號課設(shè)

.\河北科技大學(xué)課程設(shè)計(jì)報(bào)告學(xué)生姓名： 學(xué)號：專業(yè)班級： 電子信息工程課程名稱： 數(shù)字信號處理課程設(shè)計(jì)學(xué)年學(xué)期 2013——2014學(xué)年第二學(xué)期指導(dǎo)教師： 張秀清2014年6月.\課程設(shè)計(jì)成績評定表學(xué)生姓名專業(yè)班級設(shè)計(jì)題目指導(dǎo)教師評語

學(xué) 號 成績起止時(shí)間 2014-6-9至2014-6-121.窗函數(shù)設(shè)計(jì)低通濾波器2.用哈明窗設(shè)計(jì)FIR帶通數(shù)字濾波器指導(dǎo)教師：年 月 日.\目 錄窗函數(shù)設(shè)計(jì)低通濾波器1.1設(shè)計(jì)目的……………………1精品文檔放心下載1.2設(shè)計(jì)原理推導(dǎo)與計(jì)算………1感謝閱讀1.3設(shè)計(jì)內(nèi)容與要求……………2精品文檔放心下載1.4設(shè)計(jì)源程序與運(yùn)行結(jié)果……………………3謝謝閱讀1.5思考題……………………10感謝閱讀用哈明窗設(shè)計(jì)FIR帶通數(shù)字濾波器2.1設(shè)計(jì)要求……………………14精品文檔放心下載2.2設(shè)計(jì)原理和分析……………14感謝閱讀2.3詳細(xì)設(shè)計(jì)……………………15精品文檔放心下載2.4調(diào)試分析及運(yùn)行結(jié)果……………………15感謝閱讀2.5心得體會…………………17謝謝閱讀參考文獻(xiàn)………………………17謝謝閱讀.\1.窗函數(shù)設(shè)計(jì)低通濾波器1.1設(shè)計(jì)目的熟悉設(shè)計(jì)線性相位數(shù)字濾波器的一般步驟。掌握用窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器的原理和方法。精品文檔放心下載熟悉各種窗函數(shù)的作用以及各種窗函數(shù)對濾波器特性的影響。感謝閱讀學(xué)會根據(jù)指標(biāo)要求選擇合適的窗函數(shù)。1.2設(shè)計(jì)原理推導(dǎo)與計(jì)算如果所希望的濾波器的理想的頻率響應(yīng)函數(shù)為Hej，則其對應(yīng)的單位脈感謝閱讀d沖響應(yīng)為1Hjjd（4.1）d2d窗函數(shù)設(shè)計(jì)法的基本原理是設(shè)計(jì)設(shè)計(jì)低通FIR數(shù)字濾波器時(shí)，一般以理想低通濾波特性為逼近函數(shù)Hej，即感謝閱讀ej，N1Hc,其中d0，2c1jjd1cejejdsinna22ddnac用有限長單位脈沖響應(yīng)序列hn逼近hn。由于hn往往是無限長序列，精品文檔放心下載d d而且是非因果的，所以用窗函數(shù)n將hn截?cái)?，并進(jìn)行加權(quán)處理，得到：謝謝閱讀d（4.2）hnhdnnhn就作為實(shí)際設(shè)計(jì)的FIR數(shù)字濾波器的單位脈沖響應(yīng)序列，其頻率響應(yīng)函精品文檔放心下載Hej為jN1jn（4.3）Hehnen0式中，N為所選窗函數(shù)n的長度。用窗函數(shù)法設(shè)計(jì)的濾波器性能取決于窗函數(shù)n的類型及窗口長度N的取感謝閱讀.\值。設(shè)計(jì)過程中，要根據(jù)對阻帶最小衰減和過渡帶寬度的要求選擇合適的窗函數(shù)類型和窗口長度N。各種類型的窗函數(shù)可達(dá)到的阻帶最小衰減和過渡帶寬度見謝謝閱讀表（一）。旁瓣峰值過阻帶最小過渡帶帶窗函數(shù)寬幅度/dB渡帶寬衰減/dB2N4π矩形窗-13/N-120.98π三角形窗-25/N-252.18π漢寧窗-31/N-443.18π哈明窗-41/N-533.312不萊克曼窗-57π/N-745.5凱塞窗10(=7.865)-57π/N-805表（一） 各種窗函數(shù)的基本參數(shù)這樣選定窗函數(shù)類型和長度N之后，求出單位脈沖響應(yīng)hnhdn?n，并按照式（4.3）求出Hej。Hej是否滿足要求，如果Hej不滿足要求，精品文檔放心下載則要重新選擇窗函數(shù)類型和長度N，再次驗(yàn)算，直至滿足要求。感謝閱讀1.3設(shè)計(jì)內(nèi)容與要求（一）設(shè)計(jì)要求：學(xué)會計(jì)算濾波器各項(xiàng)性能指標(biāo)及如何來滿足給定的指標(biāo)要求。精品文檔放心下載用MATLAB語言編程實(shí)現(xiàn)給定指標(biāo)要求的濾波器的設(shè)計(jì)。感謝閱讀熟悉MATLAB語言，獨(dú)立編寫程序。設(shè)計(jì)低通FIR濾波器的指標(biāo)：.\通帶最大波動Rp0.25dB,，p0.2阻帶最小衰減A50dB,，0.3ss（二）、設(shè)計(jì)內(nèi)容：1.熟悉各種窗函數(shù),在MATLAB命令窗下瀏覽各種窗函數(shù)，繪出（或打印）所看到的窗函數(shù)圖。精品文檔放心下載2.編寫計(jì)算理想低通濾波器單位抽樣響應(yīng)hd(n)的m函數(shù)文件ideal.m。精品文檔放心下載編寫計(jì)算N階差分方程所描述系統(tǒng)頻響函數(shù)H(ej)的m函數(shù)文件fr.m。感謝閱讀4.根據(jù)指標(biāo)要求選擇窗函數(shù)的形狀與長度N。（至少選擇兩種符合要求的窗函數(shù)及其對應(yīng)的長度）。精品文檔放心下載5.編寫.m程序文件，通過調(diào)用ideal.m和fr．m文件，計(jì)算你設(shè)計(jì)的實(shí)際低通FIR濾波器的單位抽樣響應(yīng)h(n)和頻率響應(yīng)H(ej)，打印在頻率區(qū)間[O，π]上的幅頻響應(yīng)特性曲線H(ej)~，幅度用分貝表示。謝謝閱讀6.驗(yàn)證所設(shè)計(jì)的濾波器是否滿足指標(biāo)要求。1.4設(shè)計(jì)的源程序及運(yùn)行結(jié)果：1、利用MATLAB窗口觀察各種窗函數(shù)：w(n)

%巴特利特窗w=bartlett(20);subplot(3,2,1);plot(w);stem(w,'y');%'y'表示黃色%stem表示以離散圖輸出title('巴特利特床窗');xlabel('n');%橫坐標(biāo)為nylabel('w(n)');%縱坐標(biāo)為%布萊克曼窗w=blackman(20);

subplot(3,2,2);plot(w);stem(w,'b');%'b'表示藍(lán)色title('布萊克曼窗');xlabel('n');ylabel('w(n)');%矩形窗w=boxcar(20);subplot(3,2,3);plot(w);stem(w,'r');title('矩形窗');.\xlabel('n');title('漢寧窗');ylabel('w(n)');xlabel('n');%海明窗ylabel('w(n)');w=hamming(20);%凱澤窗plot(w);beta=5.6533;stem(w,'m');%'m'表示紫色w=kaiser(20,beta);title('海明窗');subplot(3,2,6);xlabel('n');plot(w);ylabel('w(n)');stem(w,'k');%'k'表示黑色%漢寧窗title('凱澤w=hanning(20);窗,beta=5.6533');subplot(3,2,5);xlabel('n');plot(w);ylabel('w(n)');stem(w,'g');%'g'表示綠色巴特利特床窗布萊克曼窗11))n0.5n0.5ww000510152005101520nn矩形窗海明窗11))n0.5n0.5ww000510152005101520nn漢寧窗凱澤窗,beta=5.653311))n0.5n0.5ww000510152005101520nn常用窗函數(shù)的圖形2、理想低通濾波器單位抽樣響應(yīng)hd(n)的m函數(shù)文件ideal.m。謝謝閱讀.\functionhd=ideal(wc,M)%理想低通濾波器計(jì)算%hd為0到M-1之間的理想脈沖響應(yīng)%wc為截止頻率%M為理想濾波器的長度alpha=(M-1)/2;n=0:M-1;m=n-alpha+eps;hd=sin(wc*m)./(pi*m);3、N階差分方程所描述的系統(tǒng)頻響函數(shù)的m函數(shù)文件fr.m。感謝閱讀function[db,mag,pha,gfd,w]=fr(b,a)精品文檔放心下載%求解系統(tǒng)響應(yīng)%db為相位振幅（db)%mag為絕對振幅%pha為相位響應(yīng)%grd為群延時(shí)%w為頻率樣本點(diǎn)矢量%b為Ha(z)分析多項(xiàng)式系數(shù)(對FIR而言，b=h)精品文檔放心下載%a為Hz(z)分母多項(xiàng)式系數(shù)(對FIR而言，a=1)精品文檔放心下載[H,w]=freqz(b,a,1000,'whole');感謝閱讀H=(H(1:501))';w=(w(1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));精品文檔放心下載pha=angle(H);gfd=grpdelay(b,a,w);4、實(shí)際低通濾波器FIR：%用海明窗設(shè)計(jì)低通濾波器wp=0.2*pi;.\ws=0.3*pi;tr_width=ws-wp;disp(['海明窗設(shè)計(jì)低通濾波器參數(shù):']);M=ceil(6.6*pi/tr_width)+1;感謝閱讀disp(['濾波器的長度為',num2str(M)]);感謝閱讀n=0:M-1;wc=(ws+wp)/2;%理想LPF的截止頻率感謝閱讀hd=ideal(wc,M);w_ham=(hamming(M))';h=hd.*w_ham;[db,mag,pha,gfd,w]=fr(h,[1]);謝謝閱讀delta_w=2*pi/1000;Rp=-(min(db(1:1:wp/delta_w+1)));%求出實(shí)際通帶波動感謝閱讀disp(['實(shí)際帶通波動為',num2str(Rp)]);感謝閱讀As=-round(max(db(ws/delta_w+1:1:501)));%求出最小阻帶衰減謝謝閱讀disp(['最小阻帶衰減為-',num2str(As)],’db’);感謝閱讀%繪圖subplot(1,1,1)subplot(2,6,1)stem(n,hd);title('理想沖擊響應(yīng)');axis([0M-1-0.10.3]);ylabel('hd(n)');subplot(2,6,2)stem(n,w_ham);title('海明窗');axis([0M-101.1]);ylabel('w(n)');subplot(2,6,7).\stem(n,h);title('實(shí)際沖激響應(yīng)');axis([0M-1-0.10.3]);xlabel('n');ylabel('h(n)');subplot(2,6,8)plot(w/pi,db);title('幅度響應(yīng)（db）');axis([01-10010]);grid;xlabel('以pi為單位的頻率');ylabel('分貝數(shù)');圖（1）海明窗設(shè)計(jì)的FIR海明窗設(shè)計(jì)低通濾波器參數(shù):濾波器的長度為67實(shí)際帶通波動為0.03936.\最小阻帶衰減為-52db%用布萊克曼窗設(shè)計(jì)低通濾波器wp=0.2*pi;ws=0.3*pi;tr_width=ws-wp;disp(['布萊克曼窗設(shè)計(jì)低通濾波器的參數(shù):']);謝謝閱讀M=ceil(11.0*pi/tr_width)+1;謝謝閱讀disp(['濾波器的長度為',num2str(M)]);謝謝閱讀n=0:M-1;%理想LPF的截止頻率wc=(ws+wp)/2;hd=ideal(wc,M);w_bla=(blackman(M))';h=hd.*w_bla;[db,mag,pha,gfd,w]=fr(h,[1]);謝謝閱讀delta_w=2*pi/1000;Rp=-(min(db(1:1:wp/delta_w+1)));%求出實(shí)際通帶波動謝謝閱讀disp(['實(shí)際帶通波動為',num2str(Rp)]);精品文檔放心下載As=-round(max(db(ws/delta_w+1:1:501)));%求出最小阻帶衰減謝謝閱讀disp(['最小阻帶衰減-',num2str(As)],’db’);精品文檔放心下載%繪圖subplot(2,6,3)stem(n,hd);title('理想沖擊響應(yīng)');axis([0M-1-0.10.3]);ylabel('hd(n)');subplot(2,6,4)stem(n,w_bla);title('布萊克曼窗');axis([0M-101.1]);.\ylabel('w(n)');subplot(2,6,9)stem(n,h);title('實(shí)際沖激響應(yīng)');axis([0M-1-0.10.3]);xlabel('n');ylabel('h(n)');subplot(2,6,10)plot(w/pi,db);title('幅度響應(yīng)（db）');axis([01-10010]);grid;xlabel('以pi為單位的頻率');ylabel('分貝數(shù)');圖（2）布萊克曼窗設(shè)計(jì)的FIR布萊克曼窗設(shè)計(jì)低通濾波器的參數(shù):感謝閱讀.\濾波器的長度為111實(shí)際帶通波動為0.0033304最小阻帶衰減為-73db5、技術(shù)指標(biāo)比較：（1）海明窗設(shè)計(jì)低通濾波器參數(shù):濾波器的長度為67實(shí)際帶通波動為0.03936最小阻帶衰減為-52db（2）布萊克曼窗設(shè)計(jì)低通濾波器的參數(shù):濾波器的長度為111實(shí)際帶通波動為0.0033304最小阻帶衰減為-73db在相同的技術(shù)指標(biāo)下用布萊克曼窗設(shè)計(jì)的低通濾波器實(shí)際帶通波動實(shí)際帶通波動最小，最小阻帶衰減，濾波器的長度最大；海明窗和凱澤窗最小阻帶衰減差不多，濾波器的長度頁差不多，但是海明窗實(shí)際波動小于凱澤窗；所以用布萊克曼窗用設(shè)計(jì)的FIR最逼近理想單位沖擊響應(yīng)。這三個(gè)窗設(shè)計(jì)的低通濾波器都符合要求。謝謝閱讀1.5思考題：1.設(shè)計(jì)線性相位數(shù)字濾波器的一般步驟：（1）、給定所要求的頻率響應(yīng)函數(shù)Hd(ejw)以及技術(shù)指標(biāo)阻帶衰減和;謝謝閱讀（2）、求hd=IDTFT[Hd(ejw)];由過渡帶帶寬及帶阻最小衰減的要求，利用六種常見的窗函數(shù)基本參數(shù)的比精品文檔放心下載較表或凱澤窗設(shè)計(jì)FIR濾波器的經(jīng)驗(yàn)公式，選擇窗函數(shù)的形狀及N的大?。ㄒ话憔肺臋n放心下載N要通過幾次試探而最后確定）；（3）求得所設(shè)計(jì)的FIR濾波器的單位抽樣響應(yīng)：h(n)=hd(n)w(n),N=0,1,2,3…,N-1精品文檔放心下載（4）、求Hd(ejw)=DIFT[h(n)],校驗(yàn)是否滿足設(shè)計(jì)要求，如果不滿足，則重謝謝閱讀新設(shè)計(jì)。.\窗函數(shù)有哪些指標(biāo)要求?對給定指標(biāo)要求的低通濾波器，理論計(jì)算所需窗函數(shù)的長度N。精品文檔放心下載答:窗函數(shù)的指標(biāo)要求：主瓣寬度，旁瓣峰值。海明窗設(shè)計(jì)低通濾波器的長精品文檔放心下載度為67，凱澤窗設(shè)計(jì)低通濾波器的長度為60，布萊克曼窗設(shè)計(jì)低通濾波器的長度為111。感謝閱讀用窗函數(shù)法設(shè)計(jì)FIR濾波器，濾波器的過渡帶寬度和阻帶衰減與哪些因素有關(guān)？精品文檔放心下載答：過渡帶寬度與窗函數(shù)的形狀和窗的寬度有關(guān)；阻帶衰減只有窗函數(shù)的精品文檔放心下載形狀決定，不受N的影響。4、計(jì)算理想帶通濾波器單位抽樣響應(yīng)hd(n)時(shí)取N為奇數(shù)和N為偶數(shù)有沒精品文檔放心下載有區(qū)別？你計(jì)算時(shí)所用的方法是僅適合于N為奇數(shù)或偶數(shù)的一種還是兩種都可精品文檔放心下載以用？答：以海明窗為例說明：理想沖擊響應(yīng)N=670.30.2)0.1(h0-0.10102030405060理想沖擊響應(yīng)N=680.30.2)0.1(h0-0.10102030405060理想沖擊響應(yīng)N=640.30.2)0.1(h0-0.10102030405060濾波器的長度為 67實(shí)際帶通波動為 0.03936最小阻帶衰減為 52.\濾波器的長度為68實(shí)際帶通波動為0.036424最小阻帶衰減為53濾波器的長度為64實(shí)際帶通波動為0.068677最小阻帶衰減為45取奇偶數(shù)時(shí)，低通濾波器的幅度函數(shù)是不同的，如上圖所示（海明窗設(shè)計(jì)低通濾波器N取奇偶數(shù)時(shí)的圖），通過比較，當(dāng)濾波器的長度大于技術(shù)指標(biāo)要求的長度時(shí)，選擇偶數(shù)也是滿足要求的，所以海明窗奇偶都滿足。同理，N取奇偶數(shù)，選擇其他的窗函數(shù)也滿足設(shè)計(jì)低通濾波器的指標(biāo)要求。感謝閱讀比較所選窗長N相同但窗形狀不同對濾波器設(shè)計(jì)結(jié)果的影響以及選同一種窗函數(shù)但窗長N不同時(shí)對濾波器設(shè)計(jì)結(jié)果的影響，將結(jié)論寫在報(bào)告中。感謝閱讀海明窗布萊克曼窗凱澤窗以海明窗設(shè)計(jì)的低通濾波器的長度為67濾波器的長度為67濾波器的長度為67濾波器的長度N=67，實(shí)際帶通波動0.03936實(shí)際帶通波動0.29523實(shí)帶通波動0.047919用布萊克曼窗和凱澤最小阻帶衰減為52最小阻帶衰減為30最小阻帶衰減為51窗設(shè)計(jì)低通濾波器線性斜率beta=4.5513以凱澤窗設(shè)計(jì)的低通濾波器的長度為60濾波器的長度為60濾波器的長度為60濾波器的長度N=60，實(shí)際帶通波動為0.12166實(shí)際帶通波動為0.47079實(shí)際帶通波動0.053711用布萊克曼窗和海明最小阻帶衰減為39最小阻帶衰減為26最小阻帶衰減為51窗窗設(shè)計(jì)低通濾波器線性斜率beta=4.5513.\以布萊克曼窗設(shè)計(jì)的濾波器的長度為111濾波器的長度為111濾波器的長度為111低通濾波器的長度實(shí)際帶通波動為0.02894實(shí)際帶通波動0.0033304實(shí)際帶通波動0.034426N=111，用凱澤窗和海最小阻帶衰減為59最小阻帶衰減為73最小阻帶衰減為58明窗窗設(shè)計(jì)低通濾波線性斜率beta=4.5513器結(jié)論 ：1、當(dāng)以海明窗設(shè)計(jì)的低通濾波器的長度N，用布萊克曼窗和凱澤窗設(shè)計(jì)低通濾波器，即N=67時(shí)，三個(gè)窗函數(shù)設(shè)計(jì)的低通濾波器中布萊克曼窗設(shè)計(jì)的不滿足要求，凱澤窗設(shè)計(jì)的帶通波動大于海明窗設(shè)計(jì)的低通濾波器帶通波動，最小阻帶衰減小于海明窗設(shè)計(jì)的低通濾波器。凱澤窗設(shè)計(jì)的結(jié)果沒有用海明窗設(shè)計(jì)結(jié)果好。感謝閱讀當(dāng)N取67、60、111時(shí)，用海明窗設(shè)計(jì)的低通濾波器的性能不同，通過上表比較，N也越大性能越好，滿足要求，當(dāng)N<67時(shí)不滿足設(shè)計(jì)技術(shù)指標(biāo)。謝謝閱讀2、以凱澤窗設(shè)計(jì)的低通濾波器的長度N，用布萊克曼窗和海明窗窗設(shè)計(jì)低通濾波器，即N=60時(shí)，海明窗設(shè)計(jì)的低通濾波器的最小阻帶衰減不滿足指標(biāo)要求，布萊克曼窗設(shè)計(jì)的低通濾波器的最小阻帶衰減和帶通波動都不滿足要。謝謝閱讀當(dāng)N取67、60、111時(shí)，用布萊克曼窗設(shè)計(jì)的低通濾波器的性能不同，通過上表比較，N也越大性能越好，滿足要求，當(dāng)N<111時(shí)不滿足設(shè)計(jì)技術(shù)指標(biāo)。感謝閱讀3、以布萊克曼窗設(shè)計(jì)的低通濾波器的長度N，用凱澤窗和海明窗窗設(shè)計(jì)低通濾波器，即N=111時(shí)，三個(gè)都滿足要求。但是用布萊克曼窗函數(shù)設(shè)計(jì)的低通濾波器的性能更好。謝謝閱讀當(dāng)N取67、60、111時(shí)，用凱澤窗設(shè)計(jì)的低通濾波器的性能不同，通過上表比較，N也越大性能越好，滿足要求，當(dāng)N>60時(shí)都滿足設(shè)計(jì)技術(shù)指標(biāo)。謝謝閱讀.\、用哈明窗設(shè)計(jì)FIR帶通數(shù)字濾波器2.1設(shè)計(jì)要求感謝閱讀針對一個(gè)含有5Hz、15Hz和30Hz的混合正弦波信號x(t)，設(shè)計(jì)濾波系統(tǒng)，濾除精品文檔放心下載5Hz和30Hz的正弦分量，阻帶的最小衰減不小于50dB。精品文檔放心下載① 確定采樣頻率fs、采樣長度N，對x(t)進(jìn)行采樣；感謝閱讀②選擇合適的窗函數(shù)，設(shè)計(jì)FIR帶通數(shù)字濾波器；③通過計(jì)算機(jī)仿真對濾波器的性能進(jìn)行分析。2.2設(shè)計(jì)原理分析和設(shè)計(jì)此題要求也是濾波，但是要求濾掉5Hz和30Hz的信號，而保留15Hz的信號，因此必須采用帶通濾波器。要求阻帶的最小衰減不低于50dB，查表可知可選哈明窗。精品文檔放心下載分析題目可知：f1=5Hz，T1=1/5s；f2=15Hz，T2=1/15s；f3=30Hz，T3=1/30s；則采樣時(shí)長tT至少應(yīng)為0.2s，取tT=1s；感謝閱讀采樣頻率fs≥2fc；可取fs=150；則采樣間隔T=1/fs；所以采樣長度N=tT/T=150；數(shù)字技術(shù)指標(biāo)可取謝謝閱讀wp1=2*pi*14/fs;wp2=2*pi*23/fs;ws1=2*pi*7/fs;ws2=2*pi*27/fs;濾波器設(shè)計(jì)好后，利用卷積和可得輸出信號y(n)=x(n)*h(n)，卷積和長度N=N1+N2-1;謝謝閱讀.\2.3詳細(xì)設(shè)計(jì)源程序代碼：tT=1; %采樣時(shí)長fs=150; %采樣頻率wp1=2*pi*14/fs;wp2=2*pi*20/fs;ws1=2*pi*7/fs;ws2=2*pi*27/fs;trwidth=wp1-ws1; %過渡帶帶寬N1=ceil(8*pi/trwidth) %計(jì)算濾波器階次感謝閱讀n=0:N1-1;wc1=(ws1+wp1)/2;alpha=(N1-1)/2;謝謝閱讀

wc2=(ws2+wp2)/2;%單位響應(yīng)的對稱中心

%計(jì)算3dB截止頻率m=n-alpha;hd=sin(wc2*m)./(pi*m)-sin(wc1*m)./(pi*m);%理想帶通濾波器的單位響應(yīng)wham=(hamming(N1))';精品文檔放心下載h=hd.*wham; %實(shí)際帶通濾波器單位響應(yīng)w=0:0.01:pi;H=freqz(h,1,w); %單位濾波器幅頻響應(yīng)謝謝閱讀T=1/fs; %對輸入信號進(jìn)行采樣N2=fs*tT;t=(0:N2-1)*T;x=sin(2*pi*5*t)+sin(2*pi*15*t)+sin(2*pi*30*t);感謝閱讀subplot(411) %繪制數(shù)字濾波器幅頻響應(yīng)plot(w/pi,abs(H));gridon;title('數(shù)字濾波器幅頻響應(yīng)');謝謝閱讀subplot(412) %繪制數(shù)字濾波器單位響應(yīng)stem(n,h,'.');gridon;title('數(shù)字濾波器單位響應(yīng)')精品文檔放心下載subplot(413) %繪制輸入信號plot(t,x);axis([00.4-44])gridon;title('輸入信號')subplot(414) %繪制輸出信號y=conv(x,h); %用卷積求輸出信號N=N1+N2-1; %計(jì)算卷積和長度tt=(0:N-1)*T;plot(tt,y);axis([00.8-11]).\grid