數(shù)字信號處理課程設(shè)計--數(shù)字濾波器

1、<p><b>　　目錄</b></p><p>　　一、課程設(shè)計的目的…………………………………………………………………………………….1</p><p>　　二、數(shù)字濾波器的設(shè)計步驟………………………………………………………………………….2</p><p>　　2.2、IIR數(shù)字濾波器與FIR數(shù)字濾波器的區(qū)別………………….2&

2、lt;/p><p>　　三、IIR數(shù)字濾波器………………………………………………………………………………………..2</p><p>　　3.1、IIR數(shù)字濾波器的特點2</p><p>　　3.1.2、IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)3</p><p>　　3.1.3、借助成熟的模擬濾波器的成果3</p><p>　

3、　3.1.4、需加相位校準(zhǔn)網(wǎng)絡(luò)3</p><p>　　3.2、用雙線性法設(shè)計IIR數(shù)字濾波器3</p><p>　　3.3、巴特沃斯低通濾波器的設(shè)計4</p><p>　　3.4、巴特沃斯高通濾波器的設(shè)計4</p><p>　　3.4.1、巴特沃斯高通濾波器各參數(shù)圖形4</p><p>　　3.4.2、巴特沃

4、斯高通濾波器濾波效果圖5</p><p>　　四、FIR數(shù)字濾波器……………………………………………………………………………………….5</p><p>　　4.1、FIR濾波器的特點5</p><p>　　4.2、窗函數(shù)法設(shè)計FIR數(shù)字濾波器5</p><p>　　五、程序?qū)嵗创a………………………………………………………………………

5、………………..7</p><p>　　六、問題分析………………………………………………………………………………………………12</p><p>　　七、心得體會………………………………………………………………………………………………13</p><p>　　八、參考文獻(xiàn)………………………………………………………………………………………………13</p>

6、<p><b>　　一、課程設(shè)計的目的</b></p><p>　　數(shù)字濾波是數(shù)字信號分析中最重要的組成部分之一，與模擬濾波相比，它具有精度和穩(wěn)定性高、系統(tǒng)函數(shù)容易改變、靈活性強、便于大規(guī)模集成和可實現(xiàn)多維濾波等優(yōu)點。在信號的過濾、檢測和參數(shù)的估計等方面，經(jīng)典數(shù)字濾波器是使用最廣泛的一種線性系統(tǒng)。</p><p>　　本次課程設(shè)計是通過對常用數(shù)字濾波器的

7、設(shè)計和實現(xiàn)，掌握數(shù)字信號處理的工作原理及設(shè)計方法；熟悉用雙線性變換法設(shè)計 IIR 數(shù)字濾波器和用窗函數(shù)法設(shè)計FIR數(shù)字濾波器的原理與方法，掌握利用數(shù)字濾波器對信號進行濾波的方法，掌握數(shù)字濾波器的計算機仿真方法，并能夠?qū)υO(shè)計結(jié)果加以分析。</p><p>　　二、數(shù)字濾波器的設(shè)計步驟</p><p>　　2.1、不論是IIR濾波器還是FIR濾波器的設(shè)計都包括三個步驟：</p>

8、<p>　　按照實際任務(wù)的要求，確定濾波器的性能指標(biāo)。</p><p>　　用一個因果、穩(wěn)定的離散線性時不變系統(tǒng)的系統(tǒng)函數(shù)去逼近這一性能指標(biāo)。根據(jù)不同的要求可以用IIR系統(tǒng)函數(shù)，也可以用FIR系統(tǒng)函數(shù)去逼近。</p><p>　　利用有限精度算法實現(xiàn)系統(tǒng)函數(shù)，包括結(jié)構(gòu)選擇，字長選擇等。</p><p>　　2.2、IIR數(shù)字濾波器與FIR數(shù)字濾波器的區(qū)別&

9、lt;/p><p>　　2.2.1、單位響應(yīng)</p><p>　　IIR數(shù)字濾波器單位響應(yīng)為無限脈沖序列，而FIR數(shù)字濾波器單位響應(yīng)為有限的； </p><p>　　FIR濾波器，也就是“非遞歸濾波器”，沒有引入反饋。這種濾波器的脈沖響應(yīng)是有限的。 </p><p>　　2.2.2、幅頻特性</p><p>　　IIR數(shù)字

10、濾波器幅頻特性精度很高，不是線性相位的，可以應(yīng)用于對相位信息不敏感的音頻信號上；FIR數(shù)字濾波器的幅頻特性精度較之于IIR數(shù)字濾波器低，但是線性相位，就是不同頻率分量的信號經(jīng)過fir濾波器后他們的時間差不變，這是很好的性質(zhì)。 </p><p>　　2.2.3、實時信號處理</p><p>　　FIR數(shù)字濾波器是有限的單位響應(yīng)也有利于對數(shù)字信號的處理，便于編程，用于計算的時延也小，這對實時的

11、信號處理很重要。</p><p>　　三、IIR數(shù)字濾波器</p><p>　　3.1、IIR數(shù)字濾波器的特點</p><p>　　IIR濾波器有以下幾個特點： </p><p>　　3.1.1、封閉函數(shù)</p><p>　　IIR數(shù)字濾波器的系統(tǒng)函數(shù)可以寫成封閉函數(shù)的形式。 </p><p>

12、　　3.1.2、IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)</p><p>　　IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準(zhǔn)型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。由于運算中的舍入處理，使誤差不斷累積，有時會產(chǎn)生微弱的寄生振蕩。 </p><p>　　3.1.3、借助成熟的模擬濾波器的成果</p

13、><p>　　IIR數(shù)字濾波器在設(shè)計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設(shè)計數(shù)據(jù)或圖表可查，其設(shè)計工作量比較小，對計算工具的要求不高。在設(shè)計一個IIR數(shù)字濾波器時，我們根據(jù)指標(biāo)先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。 </p><p>　　3.1.4、需加相位校準(zhǔn)網(wǎng)絡(luò)</p><p>

14、;　　IIR數(shù)字濾波器的相位特性不好控制，對相位要求較高時，需加相位校準(zhǔn)網(wǎng)絡(luò)。</p><p>　　3.2、用雙線性法設(shè)計IIR數(shù)字濾波器</p><p>　　脈沖響應(yīng)不變法的主要缺點是產(chǎn)生頻率響應(yīng)的混疊失真。這是因為從S平面到Ｚ平面是多值的映射關(guān)系所造成的。為了克服這一缺點，可以采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到-π/T～π/T之間，再用z=esT轉(zhuǎn)換到Z平面上。也就

15、是說，第一步先將整個S平面壓縮映射到S1平面的-π/T～π/T一條橫帶里；第二步再通過標(biāo)準(zhǔn)變換關(guān)系z=es1T將此橫帶變換到整個Z平面上去。這樣就使S平面與Z平面建立了一一對應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖1所示。</p><p>　　雙線性變換的映射關(guān)系</p><p>　　3.3、巴特沃斯低通濾波器的設(shè)計</p><p>　　

16、為了便于觀察濾波的效果，本次課程設(shè)計中我選用的信號只包括一高一低兩個頻率值，這樣在DFT頻譜分析圖中便可以清楚地觀察到頻率的分布及經(jīng)過各類濾波器后的效果。</p><p>　　所用信號：t=0:0.001:2;</p><p>　　y=100*sin(2*pi*2*t)+100*sin(2*pi*20*t);</p><p><b>　　f=80;<

17、/b></p><p>　　n=0:1/f:2;</p><p>　　yn=100*sin(2*pi*2*n)+ 100*sin(2*pi*20*n);</p><p>　　由信號經(jīng)過DFT頻譜分析之后的圖中我們可以讀出歸一化之后的頻率值，這樣便可以取出適當(dāng)?shù)膮?shù)設(shè)計濾波器了，由128點的DFT頻譜分析圖中可以看出，設(shè)計低通濾波器時應(yīng)取wp=0.2*pi,ws

18、=0.4*pi,和設(shè)計好的低通濾波器一起便可以對原信號進行濾波了。</p><p>　　3.4、巴特沃斯高通濾波器的設(shè)計</p><p>　　3.4.1、巴特沃斯高通濾波器各參數(shù)圖形</p><p>　　3.4.2、巴特沃斯高通濾波器濾波效果圖</p><p>　　四、FIR數(shù)字濾波器</p><p>　　4.1、FI

19、R濾波器的特點</p><p>　　有限長單位沖激響應(yīng)（FIR）濾波器有以下特點： </p><p>　　(1) 系統(tǒng)的單位沖激響應(yīng)h (n)在有限個n值處不為零； </p><p>　　(2) 系統(tǒng)函數(shù)H(z)在|z|>0處收斂，極點全部在z = 0處（因果系統(tǒng)）； </p><p>　　(3) 結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入

20、的反饋，但有些結(jié)構(gòu)中（例如頻率抽樣結(jié)構(gòu)）也包含有反饋的遞歸部分。 </p><p>　　設(shè)FIR濾波器的單位沖激響應(yīng)h (n)為一個N點序列，0 ≤ n ≤ N —1，則濾波器的系統(tǒng)函數(shù)為 </p><p>　　H(z)=∑h(n)*z^-n</p><p>　　4.2、窗函數(shù)法設(shè)計FIR數(shù)字濾波器</p><p>　　用于產(chǎn)生窗函數(shù)的MAT

21、LAB文件有如下8個：</p><p>　?。?）bartlett.m (三角窗)</p><p>　　（2）blackman.m (布萊克曼窗)</p><p>　?。?）boxcar.m (矩形窗)</p><p>　?。?）hamming.m (海明窗)</p><p>　?。?）hanning.m (

22、漢寧窗)</p><p>　?。?）triang.m (三角窗)</p><p>　　（7）chebwin.m (切比雪夫窗)</p><p>　?。?）kaiser.m (凱澤窗)</p><p>　　2．FIR數(shù)字濾波器的文件</p><p>　　以下MATLAB文件是用于FIR數(shù)字濾波器的文件：</p

23、><p><b>　　（1）fir1.m</b></p><p>　　本文件采用窗函數(shù)法設(shè)計FIR數(shù)字濾波器，其調(diào)用格式是</p><p>　　1）b=fir1（N,Wn）</p><p>　　2）b=fir1（N,Wn,’high’）</p><p>　　3）b=fir1（N,Wn,’stop’）&

24、lt;/p><p>　　式中N為濾波器的階次，因此濾波器的長度為N+1；Wn是通帶截止頻率，其值在0~1之間，1對應(yīng)抽樣頻率的一半；b是設(shè)計好的濾波器系數(shù)。</p><p>　　對于格式（1），若Wn是一標(biāo)量，則可用來設(shè)計低通濾波器；若Wn是的向量，則用來設(shè)計帶通濾波器；若Wn是的向量，則可用來設(shè)計帶濾波器，此時，格式將變?yōu)椋?lt;/p><p>　　b=fir1（N,Wn

25、,’DC-1’）或b=fir1（N,Wn,’DC-0’）</p><p>　　其中，前者保證第一個帶為通帶，后者保證第一個帶為阻帶。</p><p>　　格式（2）用來設(shè)計高通濾波器；格式（3）用來設(shè)計帶阻濾波器。</p><p>　　值得注意是：在上述所有格式中，若不指定窗函數(shù)的類型，則fir1自動選擇漢明窗。</p><p><b&

26、gt;　?。?）fir2.m</b></p><p>　　本文件采用窗函數(shù)法設(shè)計具有任意幅頻特性的FIR數(shù)字濾波器。其調(diào)用格式是</p><p>　　b=fir1（N,F,M）</p><p>　　其中F是頻率向量，其值在0~1之間，M是與F相對應(yīng)的所希望的幅頻響應(yīng)。不指定窗函數(shù)的類型，則自動選擇漢明窗。</p><p>　?。?

27、）remez.m</p><p>　　本文件用來設(shè)計采用切比雪夫最佳一致逼近FIR數(shù)字濾波器。同時，還可以用來設(shè)計希爾伯特變換器和差分器。其調(diào)用格式是</p><p>　　1）b=remez（N,F,A）</p><p>　　2）b=remez（N,F,A,W）</p><p>　　3）b=remez（N,F,A,W,’hilbert’）&l

28、t;/p><p>　　4）b=remez（N,F,A,W,’differentiator’）</p><p>　　其中，N是給定的濾波器的階次；b是設(shè)計的濾波器的系數(shù)，其長度為N+1；F是頻率向量，其值在0~1之間；A是對應(yīng)F的各頻段上的理想幅頻響應(yīng)；W是各頻段上的加權(quán)向量。</p><p>　　值得注意是：若b的長度為偶數(shù)，設(shè)計高通和帶阻濾波器時有可能出現(xiàn)錯誤，因此最

29、好保證b的長度為奇數(shù)，即N應(yīng)為偶數(shù)。</p><p>　?。?）remexord.m</p><p>　　本文件采用切比雪夫一致逼近設(shè)計FIR數(shù)字濾波器時所需要的濾波器階次。其調(diào)用格式是</p><p>　　[N,Fo,Ao,W]=remexord（F,A,DEV,Fs）</p><p>　　式中，F(xiàn)、A的含義同文件（3），是通帶和阻帶上的偏

30、差；該文件輸出的是符合要求的濾波器階次N、頻率向量Fo、幅度向量Ao和加權(quán)向量W。若設(shè)計者事先不能確定自己要設(shè)計的濾波器的階次，那么，調(diào)用remexord后，就可利用這一族參數(shù)再調(diào)用remez，即b=remez（N,Fo,Ao,W），從而設(shè)計出所需要的濾波器。因此，通常remez和remexord結(jié)合使用。</p><p>　　值得說明是：remexord給出的階次N有可能偏低，這時適當(dāng)增加N即可；另外，若N為奇

31、數(shù)，就可令其加1，使其變?yōu)榕紨?shù)，這樣b的長度為奇數(shù)。</p><p>　　（5）sgolay.m</p><p>　　本文件用來設(shè)計Savitzky-Golay平滑濾波器。其調(diào)用格式是</p><p>　　b=sgolay（k,f）</p><p>　　式中k是多項式的階次，f是擬合的雙邊點數(shù)。要求，且f為奇數(shù)。</p><

32、;p>　?。?）firls.m</p><p>　　本文件用最小平方法設(shè)計線性相位FIR數(shù)字濾波器?？稍O(shè)計任意給定的理想幅頻特性。</p><p>　?。?）fircls.m</p><p>　　用帶約束的最小平方法設(shè)計線性相位FIR數(shù)字濾波器?？稍O(shè)計任意給定的理想幅頻特性。</p><p>　?。?）fircls1.m</p&g

33、t;<p>　　用帶約束最小平方法設(shè)計線性相位FIR低通和高通濾波器。可設(shè)計任意給定的理想幅頻特性。</p><p>　?。?）firrcos.m</p><p>　　用來設(shè)計低通線性相位FIR數(shù)字濾波器，其過渡帶為余弦函數(shù)形狀。</p><p><b>　　五、程序?qū)嵗创a</b></p><p>　　1

34、、題目：A．試用MATLAB設(shè)計一巴特沃斯低通數(shù)字濾波器，要求通帶截止頻率Wp=30Hz，阻帶截止頻率Ws=35Hz，通帶衰減不大于0.5dB，阻帶衰減不小于40dB，抽樣頻率Fs=100Hz</p><p><b>　　fp = 30;</b></p><p><b>　　fs = 35;</b></p><p><

35、;b>　　Fs = 100;</b></p><p>　　wp = 2*pi*fp/Fs;</p><p>　　ws = 2*pi*fs/Fs;</p><p>　　wp = tan(wp/2);</p><p>　　ws = tan(ws/2);</p><p>　　[N, wn] = buttor

36、d(wp, ws, 0.5, 40, 's');</p><p>　　[z, p, k] = buttap(N);</p><p>　　[b, a] = zp2tf(z, p, k);</p><p>　　[B, A] = lp2lp(b, a, wp);</p><p>　　[bz, az] = bilinear(B, A,

37、 .5);</p><p>　　[h, w] = freqz(bz, az, 256, Fs);</p><p><b>　　figure</b></p><p>　　plot(w, abs(h))</p><p><b>　　grid on</b></p><p><

38、b>　　運行結(jié)果：</b></p><p>　　2、題目：A. 用窗函數(shù)法設(shè)計一個線性相位FIR低通濾波器，并滿足性能指標(biāo)；通帶邊界頻率，阻帶邊界頻率，阻帶衰減不小于40dB，通帶波紋不大于3dB。要求選擇漢寧窗。</p><p>　　wp =0.5*pi;</p><p>　　ws=0.66*pi;</p><p>　　w

39、delta =ws-wp;</p><p>　　N= ceil(8*pi/wdelta)</p><p>　　if rem(N,2)==0</p><p><b>　　N=N+1;</b></p><p><b>　　end</b></p><p><b>　　N

40、w =N;</b></p><p>　　wc =(wp+ws)/2;</p><p>　　n =0: N-1;</p><p>　　alpha =(N-1)/2;</p><p>　　m =n-alpha+0.00001;</p><p>　　hd =sin(wc*m)./(pi*m);</p>

41、<p>　　win =hanning(Nw)';</p><p>　　h=hd.*win;</p><p><b>　　b=h;</b></p><p>　　freqz(b,1,512)</p><p><b>　　運行結(jié)果：</b></p><p>　

42、　3、題目：用雙線性變換法設(shè)計一個數(shù)字巴特沃夫低通濾波器，在頻帶截止頻率0.2pi處的衰減不大于1dB，在阻帶截止頻率0.3pi處的衰減不小于15dB。</p><p><b>　　syms z;</b></p><p>　　H=1/10*(1+2/z+4/z^2+2/z^3+1/z^4)</p><p>　　h=iztrans(H)</

43、p><p>　　w=-2*pi:0.1:2*pi;</p><p>　　hf=1/10*(1+2./(exp(j*w))+4./(exp(j*w).^2)+2./(exp(j*w).^3)+1./(exp(j*w).^4));</p><p>　　subplot(131)</p><p>　　plot(w,abs(hf))</p>

44、<p>　　subplot(132)</p><p>　　plot(w,angle(hf))</p><p>　　subplot(133)</p><p>　　h=[0.1 0.2 0.4 0.2 0.1]</p><p><b>　　stem(h)</b></p><p><b&

45、gt;　　運行結(jié)果：</b></p><p>　　4、設(shè)某FIR數(shù)字濾波器的系統(tǒng)函數(shù)為，試求h(n)的表達(dá)式、的幅頻響應(yīng)和相頻響應(yīng)的表達(dá)式，并畫出對應(yīng)表達(dá)式的波形。</p><p>　　wp=0.2*pi;ws=0.3*pi;</p><p>　　Fs=4000;T=1/Fs; </p><p>　　OmegaP=(2/T)*ta

46、n(wp/2);</p><p>　　OmegaS=(2/T)*tan(ws/2);</p><p>　　rp=1;rs=15;as=15;</p><p>　　ripple=10^(-rp/20);attn=10^(-rs/20);</p><p>　　[n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s

47、9;);</p><p>　　[z,p,k]=Buttap(n);</p><p>　　[b,a]=zp2tf(z,p,k);</p><p>　　[bt,at]=lp2lp(b,a,wn);</p><p>　　[b,a]=bilinear(bt,at,Fs);</p><p>　　[db,mag,pha,grd,w

48、]=freqz_m(b,a);</p><p><b>　　%下面繪出各條曲線</b></p><p>　　subplot(2,2,1);plot(w/pi,mag);title('Magnitude Frequency幅頻特性');</p><p>　　xlabel('w(/pi)');ylabel('

49、|H(jw)|');</p><p>　　axis([0,1,0,1.1]);</p><p>　　set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);</p><p>　　set(gca,'YTickMode','manual'

50、;,'YTick',[0 attn ripple 1]);grid</p><p>　　subplot(2,2,2);plot(w/pi,db);title('Magnitude Frequency幅頻特性(db)');</p><p>　　xlabel('w(/pi)');ylabel('dB');</p>&

51、lt;p>　　axis([0,1,-30,5]);</p><p>　　set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);</p><p>　　set(gca,'YTickMode','manual','YTick',[-60 -as

52、-rp 0]);grid</p><p>　　subplot(2,2,3);plot(w/pi,pha/pi);title('Phase Frequency相頻特性');</p><p>　　xlabel('w(/pi)');ylabel('pha(/pi)');</p><p>　　axis([0,1,-1,1]);

53、</p><p>　　subplot(2,2,4);plot(w/pi,grd);title('Group Delay群延時');</p><p>　　xlabel('w(/pi)');ylabel('Sample');</p><p>　　axis([0,1,0,15]);</p><p>　

54、　set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1]);grid</p><p>　　freqz_m.m文件</p><p>　　function [db,mag,pha,grd,w]=freqz_m(b,a);</p><p>　　[H,w]=freqz(b,a,10

55、00,'whole');</p><p>　　H=(H(1:501))';w=(w(1:501))';</p><p>　　mag=abs(H);</p><p>　　db=20*log10((mag+eps)/max(mag));</p><p>　　pha=angle(H);</p><

56、p>　　grd=grpdelay(b,a,w);</p><p><b>　　運行結(jié)果：</b></p><p><b>　　六、問題分析</b></p><p>　　6.1、雙線性變換法分析：</p><p>　　與脈沖響應(yīng)不變法相比，雙線性變換的主要優(yōu)點：靠頻率的嚴(yán)重非線性關(guān)系得到S平面與

57、Z平面的單值一一對應(yīng)關(guān)系，整個jΩ軸單值對應(yīng)于單位圓一周，這個關(guān)系就是式 所表示的，其中ω和Ω為非線性關(guān)系。如圖圖中看到，在零頻率附近，Ω～ω接近于線性關(guān)系，Ω進一步增加時，ω增長變得緩慢， (ω終止于折疊頻率處)，所以雙線性變換不會出現(xiàn)由于高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象。</p><p>　　雙線性變換法的缺點：Ω與ω的非線性關(guān)系，導(dǎo)致數(shù)字濾波器的幅頻響應(yīng)相對于模擬濾波器的幅頻響應(yīng)有畸變，(使數(shù)字

58、濾波器與模擬濾波器在響應(yīng)與頻率的對應(yīng)關(guān)系上發(fā)生畸變)。例如，一個模擬微分器，它的幅度與頻率是線性關(guān)系，但通過雙線性變換后，就不可能得到數(shù)字微分器。另外，一個線性相位的模擬濾波器經(jīng)雙線性變換后，濾波器就不再有線性相位特性。</p><p><b>　　6.2、實驗觀察</b></p><p>　　在雙線性變換法中，模擬頻率與數(shù)字頻率不再是線性關(guān)系，所以一個線性相位模擬器

59、經(jīng)過雙線性變換后得到的數(shù)字濾波器不再保持原有的線性相位了。如以上實驗過程中，采用雙線性變化法設(shè)計的butter和cheby1數(shù)字濾波器，從圖中可以看到這種非線性關(guān)系。</p><p><b>　　七、心得體會</b></p><p>　　這次課程設(shè)計，由于自己數(shù)字信號處理的理論課程跟不上進度，對數(shù)字濾波 器的設(shè)計缺少認(rèn)識。而且需要運用 matlab 軟件，所以比較吃力

60、。設(shè)計過程，看 了一遍數(shù)字信號處理課程關(guān)于數(shù)字濾波器的設(shè)計的內(nèi)容， 再通過利用參考文獻(xiàn)與 網(wǎng)絡(luò)，完成了用 Matlab 進行數(shù)字信號處理課程設(shè)計。通過課程設(shè)計，加深了對 課堂抽象概念的理解,鞏固了課堂上所學(xué)的理論知識,并能很好地理解與掌握數(shù) 字信號處理中的基本概念、基本原理、基本分析方法。同時掌握編程方法和解決 實際問題的技巧。 與其他高級語言的程序設(shè)計相比， MATLAB 環(huán)境下可以更方便、 快捷地設(shè)計出具有嚴(yán)格線性相位的 FIR 濾

61、波器，節(jié)省大量的編程時間，提高編程 效率，且參數(shù)的修改也十分方便，還可以進一步進行優(yōu)化設(shè)計。相信隨著版本的 不斷提高，MATLAB 在數(shù)字濾波器技術(shù)中必將發(fā)揮更大的作用。同時，用 MATLAB 計算有關(guān)數(shù)字濾波器的設(shè)計參數(shù)，如 H(z)、h(n)等，對于數(shù)字濾波器的硬件實 現(xiàn)也提供了一條簡單而準(zhǔn)確的途徑和依據(jù)。</p><p><b>　　八、參考文獻(xiàn)</b></p><

62、p>　?。?）數(shù)字信號處理教程 程佩青 清華大學(xué)出版社</p><p>　?。?）MATLAB信號處理詳解 陳亞勇 人民郵電出版社</p><p>　　（3）應(yīng)用MATLAB實現(xiàn)信號分析和處理  張明照  科學(xué)出版社 </p><p>　?。?）數(shù)字信號處理及MATLAB實現(xiàn)  余成波 

63、清華大學(xué)出版社 </p><p>　?。?）MATLAB7.0在數(shù)字信號處理中的應(yīng)用  羅軍輝  機械工業(yè)出版社 </p><p>　?。?）MATLAB信號處理  劉波  電子工業(yè)出版社 </p><p>　?。?）Matlab信號處理與應(yīng)用  董長虹 &