課程設(shè)計基于MATLAB的FIR濾波器設(shè)計

1、基于matlab的fir濾波器設(shè)計課程設(shè)計報告專 業(yè)： 通信工程 班 級： 2008114 學(xué) 號： 200811403 姓 名： 李 石 指導(dǎo)教師： 許麗 劉明堂 2010 年 6 月 30 日目 錄一、 概述31、設(shè)計目的32、設(shè)計要求33、系統(tǒng)功能3二、濾波器原理3三、軟件設(shè)計流程51、語音信號的采集52、語音信號的頻譜分析53、用濾波器對語音信號進行分析54、回放語音信號6四、源程序代碼及整體波形6五、總結(jié)8 (要求：給出一級目錄，宋體加粗，四號字,1.5倍行距。)一、概述1）設(shè)計目的通過對課程設(shè)計任務(wù)的完成，使我們進一步鞏固數(shù)字信號處理的基本概念、理論、分析方法和實現(xiàn)方法；使掌握的基

2、本理論和分析方法方面的知識得到進一步擴展；能有效地將理論和實際緊密結(jié)合；增強學(xué)生軟件編程實現(xiàn)能力和解決實際問題的能力。2）設(shè)計要求1、認真獨立的完成課程設(shè)計。2、通過課程設(shè)計，加深對所學(xué)知識的理解和認識。3、仿真調(diào)試通過,達到預(yù)期效果。4、寫出設(shè)計報告。3）系統(tǒng)功能因?qū)φZ音信號加噪后如不進行濾波會產(chǎn)生很大的失真，濾波器的作用會把濾波器工作頻率以外的頻譜成分濾除掉，可較好的回復(fù)原語音信號，濾波器的作用在實際工作中的作用也如此，如要傳輸?shù)男盘柾ㄟ^信道后會疊加上高斯白噪聲，在解調(diào)之前可先通過bpf濾除bpf以外的噪聲。二、濾波器原理1）fir和iir數(shù)字濾波器的選擇數(shù)字濾波器根據(jù)其實現(xiàn)的網(wǎng)絡(luò)結(jié)構(gòu)或者

3、其沖激響應(yīng)函數(shù)的時域特性，可分為兩種，即無限沖擊響應(yīng)iir濾波器和有限沖擊響應(yīng)fir濾波器。iir濾波器設(shè)計方法有間接法和直接法，間接法是借助于模擬濾波器的設(shè)計方法進行的。其設(shè)計步驟是：先設(shè)計過渡模擬濾波器得到系統(tǒng)函數(shù)h（s），然后將h（s）按某種方法轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)h(z)。fir濾波器必須采用間接法，常用的方法有窗函數(shù)法、頻率采樣發(fā)和切比雪夫等波紋逼近法。對于線性相位濾波器，經(jīng)常采用fir濾波器。fir濾波器結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入的反饋，系統(tǒng)函數(shù)h (z)在極點處收斂，極點全部在z = 0處（因果系統(tǒng)），因而只能用較高的階數(shù)達到高的選擇性。fir數(shù)字濾波器的幅頻特

4、性精度較之于iir數(shù)字濾波器低，但是線性相位，就是不同頻率分量的信號經(jīng)過fir濾波器后他們的時間差不變，這是很好的性質(zhì)。fir數(shù)字濾波器是有限的單位響應(yīng)也有利于對數(shù)字信號的處理，便于編程，用于計算的時延也小，這對實時的信號處理很重要。 fir濾波器因具有系統(tǒng)穩(wěn)定，易實現(xiàn)相位控制，允許設(shè)計多通帶（或多阻帶）濾波器等優(yōu)點收到人們的青睞。iir濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。iir濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準(zhǔn)型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。同時，iir數(shù)字濾波器在設(shè)計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和

5、橢圓濾波器等，在設(shè)計一個iir數(shù)字濾波器時，我們根據(jù)指標(biāo)先寫出模擬濾波器的公式，然后通過一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。本次設(shè)計中采用的是iir濾波器，iir數(shù)字濾波器的幅頻特性精度較高，且階數(shù)較fir低，同時，iir數(shù)字濾波器在設(shè)計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等。本次設(shè)計中接住了巴特沃斯濾波器的設(shè)計。2）fir濾波器的程序設(shè)計由頻譜圖可以看出信號能量集中在低頻部分，為了減少無用的高頻分量我設(shè)計了以下的fir低通濾波器將其濾除。該設(shè)計采用窗函數(shù)法來完成以下指標(biāo)的低通濾波器。通帶截頻fp=3000hz，阻帶截頻fs=3500hz，通帶最

6、大衰減0.3db，阻帶最小衰減50db，采樣頻率fs為信號源采樣頻率22050hz。源程序如下：fs=22050;%信號的采樣頻率fc1=1500;%基帶信號頻率wc1=2*pi*fc1/fs;%相對角頻率wp1=2*pi*1000/fs;ws1=2*pi*1200/fs;n1=ceil(200*2*pi/ws1-wp1);window=boxcar(n1+1); %長度為n1的矩形窗windowb1=fir1(n1,wc1/pi,window);figure(1);freqz(b1,1,512); %繪制濾波器的幅頻響應(yīng)圖title(低通濾波器的頻率響應(yīng));三、軟件設(shè)計流程1、語音信號的采集

7、 利用windows下的錄音機, 錄制一段自己的話音, 時間在幾秒內(nèi)。然后在matlab軟件平臺下, 利用函數(shù)wavread對語音信號進行采樣, 記住采樣頻率和采樣點數(shù)。x,fs,bits=wavread();2、語音信號的頻譜分析 首先畫出語音信號的時域波形， 然后對語音信號進行頻譜分析。 在matlab 中, 可以利用函數(shù)fft對信號進行快速傅里葉變換, 得到信號的頻譜特性。3、用濾波器對語音信號進行濾波 設(shè)計濾波器對采集的信號進行濾波, 比較濾波前后語音信號的波形及頻譜。要求自己確定濾波器的性能指標(biāo)，用沖激響應(yīng)不變法、雙線性變換法或窗函數(shù)法設(shè)計低通濾波器。4、回放語音信號 在matlab

8、中, 函數(shù)sound 可以對聲音進行回放，其調(diào)用格式: sound (x, fs, bits)；四、源程序代碼及整體波形1.程序代碼:fs=8000;fc1=1500;wc1=2*pi*fc1/fs;wp1=2*pi*1000/fs;ws1=2*pi*1200/fs;n1=ceil(200*2*pi/ws1-wp1);window=boxcar(n1+1); %長度為n1的矩形窗windowb1=fir1(n1,wc1/pi,window);figure(1);freqz(b1,1,512);title(低通濾波器的頻率響應(yīng));x,fs,bits=wavread(d:/音樂/錄音.wav);%

9、信號讀取n= length (x) ; %求出語音信號的長度%sound(x,fs,nbits); %回放語音信號x=fft(x,n);figure(2)%傅里葉變換subplot(211);plot(x);title(原始信號波形);subplot(212);plot(abs(x);title(原始信號頻譜)sound(x,fs);si=awgn(x,40)figure(3)subplot(211);plot(si)title(加噪語音信號的時域波形);y1=fft(si);.subplot(212)plot(abs(y1);title(加噪信號的頻譜);y=filter(b1,1,si); %濾波輸出figure(4)%繪制輸入信號subplot(2,1,1); plot(si),title(輸入信號)%繪制輸出信號subplot(2,1,2);plot(y) title(輸出信號)sound(y,fs);%輸入輸出頻譜x=fft(si);y=fft(y);figure(5);subplot(2,1,1);plot(abs(x);xlabel(頻率);ylabel(幅度);title(輸入信號頻譜);grid on;subplot(2,1,2);plot(abs(y);xlabel(頻率);ylabel(幅度);title(輸出信號頻譜)