課程設計---iir數字低通濾波器的matlab實現

上傳人：奔*** IP屬地：河北更新時間：2024-03-01 格式：doc 頁數：17 大?。?05.00KB 人氣指數：12 舉報 版權申訴

已閱讀1頁，還剩16頁未讀，繼續(xù)免費閱讀

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、　　IIR數字低通濾波器的MATLAB實現　　目錄　　一、引言4　?。ㄒ唬┰O計目的4　?。ǘ┰O計思想4</b

2、>　　二、IIR數字濾波器的MATLAB實現　　三、實驗步驟13　　四、實驗結果及截圖14　　結果14　　參考文獻15&l

3、t;/b>　　附錄16　　一、引言　　設計目的　　（1）熟悉Matlab在信號處理方面的應用　?。?）深刻了

4、解頻域上的頻譜特性。（3）了解數字濾波器的設計。（4）學會比較濾波前后語音信號的差異。（5）利用GUI工具做一個簡潔大方的圖形用戶界面。　　二、IIR數字低通濾波器的MATLAB實現　?。ㄒ唬㊣IR數字濾波器的傳遞函數及特點　　設IIR濾波器的輸入序列為x(n)，則IIR濾波器的輸入序列x(n)與輸出序列y

5、(n)之間的關系可以用下面的方程式表示：　　其中，和是濾波器的系數，其中中至少有一個非零。與之相對應的差分方程為：　　由傳遞函數可以發(fā)現無限常單位沖激響應濾波器有如下特點：　　單位沖激響應h(n)是無限長的。 　　系統(tǒng)傳遞函數H(z)在有限z平面上有極點存在。</p

6、>　　結構上存在著輸出到輸入的反饋，也就是結構上是遞歸型的。　?。ǘ㊣IR數字濾波器的設計與實現　　IIR數字濾波器的設計有多種方法，如頻率變換法、數字域直接設計以及計算輔助設計等。下面只介紹頻率變換設計法。首先考慮由模擬低通濾波器到數字低通濾波器的轉換，其基本的設計過程如下:　　將

7、數字濾波器的技術指標轉換為模擬濾波器的技術指標；　　設計模擬濾波器G(S)；　　將G(S)轉換成數字濾波器H(Z)；　　在低通濾波器的設計基礎上，可以得到數字高通、帶通、帶阻濾波器的設計流程如下：　　給定數字濾波器的設計要求（高通、帶阻、帶通）；&

8、lt;p>　　轉換為模擬（高通、帶阻、帶通）濾波器的技術指標；　　轉換為模擬低通濾波器的指標；　　設計得到滿足第三步要求的低通濾波器傳遞函數；　　通過頻率轉換得到模擬（高通、帶阻、帶通）濾波器；　　變換為數字（高通、帶阻、帶通）濾波器。<

9、;p>　　在matlab中設計IIR濾波器的方法及其它們所用到的函數如表6－1所示。　　表6－1 matlab中設計IIR濾波器的方法列表　?。ㄈ?雙線性變換法　　與沖激響應不變法比較，雙線性變換的主要優(yōu)點是靠頻率的非線性關系得到S平面與Z平面的單值一一對應關系，整個值對應于單位圓一周。所以從模擬傳遞函數

10、可直接通過代數置換得到數字濾波器的傳遞函數。　　【實例6-4】設計一個截止頻率為200Hz的數字低通濾波器，采用頻率為1000Hz。　　解：源代碼如下：　　[z,p,k]=buttap(3);　　[b,a]=zp2tf(z,p,k);

11、　　[bt,at]=lp2lp(b,a,200*2*pi);　　[bz,az]=bilinear(bt,at,1000);　　freqz(bz,az,512,1000)　　結果如圖6－4所示：　　bz =<

12、/b>　　0.0753 0.2259 0.2259 0.0753　　az =　　1.0000 -0.8266 0.5154 -0.0865　　【實例6-5】基于Butterworth模擬濾波器原型

13、，使用雙線性轉換設計數字濾波器，其中參數指標為：　　通帶截止頻率：通帶波動值：　　阻帶截止頻率：阻帶波動值：　　解：首先確定濾波器的階數N，同時根據確定＝0.5。　　接著使用bilinear進行雙線性轉換，最后繪制在頻域上的各種圖像，其源代碼如下：

14、　　％數字濾波器指標　　wp=0.2*pi;　　ws=0.3*pi;　　Rp=1;　　As=15;<p&g

15、t;　　％將數字濾波器指標反轉變化為模擬濾波器的參數　　T=1;　　fs=1/T;　　omegap=(2/T)*tan(wp/T);　　omegas=(2/T)*tan(ws/T);</p&

16、gt;　　ep=sqrt(10^(Rp/10)-1);　　Ripple=sqrt(1/(1+ep*ep));　　Attn=1/(10^(As/20));　?。utterworth原型模擬濾波器的設計　　[cs,ds]=afd_butt(omegap

17、,omegas,Rp,As);　　％雙線性變換　　[b,a]= bilinear(cs,ds,T,fs);　?。ヮl域圖像的繪制　　freqz(b,a);

18、　程序運行后，產生4階的butterworth數字濾波器，頻率響應如圖6－5所示的波形。　　圖6－5 4階Butterworth數字濾波器　　實驗步驟　　錄入一段語音信號，設置采樣率和采樣點數。　　輸入信號的頻譜分析

19、　　要求首先畫出輸入信號的時域波形；然后對該信號進行快速傅立葉變換，得到信號的頻譜特性。　　設計數字低通濾波器并畫出數字低通濾波器的頻率響應。　　用濾波器對信號進行濾波，比較濾波前后語音信號的波形及頻譜。　　利用Matlab GUI工具做一個圖形用戶界

20、面，窗口簡潔大方。　　實驗結果及截圖　　結果　　總結　　通過這次課程設計，使自己對通信原理中的信號的變換，特別是傅立葉變換及其實現過程有了更深的了解。我對ma

21、tlab 軟件也有了熟練的了解，對工程項目設計積累的一定經驗。起初在面對一個陌生時會有忌憚和不知所措，而通過這次訓練基本可以通過help 指令和doc 指令來自學函數用法。在程序調試的時候經常會出現錯誤，耐心調試找到錯誤并搞清楚為什么會出錯則是作為一個程序員的一項基本素質。　　我們團隊分工明確，成員們勤勞能干，奠定我們高效地完成課程設計的基礎。而這期間，我們也確實遇到了很多困難，但我們都攻克了。

22、總之，我們在這次課程設計中學到了很多。　　參考文獻：　　【1】通信原理第六版樊昌信曹麗娜編著國防工業(yè)出版社　　【2】信號與系統(tǒng)引論鄭君里應啟楊為理編著高等教育出版社　　【3】 MATLAB及在電子信息類課程

23、中的應用第二版唐向宏岳恒立鄭雪峰編著電子工業(yè)出版社　　附錄：　　源程序　　function varargout = idea2(varargin)　　% IDEA2 M-file f

24、or idea2.fig　　% IDEA2, by itself, creates a new IDEA2 or raises the existing　　% singleton*.　　%　　% H =

25、IDEA2 returns the handle to a new IDEA2 or the handle to　　% the existing singleton*.　　%　　% IDEA2('CALLBACK',hObject,even

26、tData,handles,...) calls the local　　% function named CALLBACK in IDEA2.M with the given input arguments.　　%　　% IDEA2('Proper

27、ty','Value',...) creates a new IDEA2 or raises the　　% existing singleton*. Starting from the left, property value pairs are　　% applied to the GUI before id

28、ea2_OpeningFunction gets called. An　　% unrecognized property name or invalid value makes property application　　% stop. All inputs are passed to idea2_OpeningFcn via v

29、arargin.　　%　　% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one　　% instance to run (singleton)"

31、dify the response to help idea2　　% Last Modified by GUIDE v2.5 13-Jun-2012 16:27:25　　% Begin initialization code - DO NOT EDIT　　gui_Singleton = 1;

32、　　gui_State = struct('gui_Name', mfilename, ...　　'gui_Singleton', gui_Singleton, ...　　'gui_OpeningFcn', @idea2_OpeningFcn, ...

33、　　'gui_OutputFcn', @idea2_OutputFcn, ...　　'gui_LayoutFcn', [] , ...　　'gui_Callback', []);　　if nargin && ischar(

34、varargin{1})　　gui_State.gui_Callback = str2func(varargin{1});　　end　　if nargout　　[varargout{1:nargout}] = gui_mainfcn(

35、gui_State, varargin{:});　　else　　gui_mainfcn(gui_State, varargin{:});　　end　　% End initialization co

36、de - DO NOT EDIT　　% --- Executes just before idea2 is made visible.　　function idea2_OpeningFcn(hObject, eventdata, handles, varargin)　　% This function has no

37、 output args, see OutputFcn.　　% hObject handle to figure　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with

38、 handles and user data (see GUIDATA)　　% varargin command line arguments to idea2 (see VARARGIN)　　% Choose default command line output for idea2　　handles.ou

39、tput = hObject;　　% Update handles structure　　guidata(hObject, handles);　　global SampleRate　　SampleRate=44100;　　glo

40、bal SecLength;　　SecLength=60; %因為wav文件是十分鐘的，不能全部讀入，所以將其分為10段，一段60秒　　global SecNum;　　SecNum=600/SecLength;　　global CurrentSecti

41、on　　% UIWAIT makes idea2 wait for user response (see UIRESUME)　　% uiwait(handles.figure1);　　% --- Outputs from this function are returned to the command line

42、.　　function varargout = idea2_OutputFcn(hObject, eventdata, handles) 　　% varargout cell array for returning output args (see VARARGOUT);　　% hObject handl

43、e to figure　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (see GUIDATA)　　% Get defa

44、ult command line output from handles structure　　varargout{1} = handles.output;　　% --- Executes on button press in pushbutton1.　　function Insert_Callback(hObj

45、ect, eventdata, handles)　　% hObject handle to pushbutton1 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles st

46、ructure with handles and user data (see GUIDATA)　　global SecLength;　　global SecNum;　　global CurrentSection;　　global SampleRate;</p&g

47、t;　　global Myrate;　　global SoundData;　　global Mycode;　　SampleRate = Myrate;　　%Sam =Myrate;　　%disp(Sam);&

48、lt;/p>　　[FileName,PathName] = uigetfile('*.wav');　　data=wavread(strcat(PathName,FileName),[SecLength*(CurrentSection-1)*(SampleRate)+1,SecLength*CurrentSection*SampleRate]);<

49、/p>　　SoundData=data(1:10:length(data),1);　　%N= Mycode;　　n= 0:length(SoundData)-1;　　t= n/SampleRate;　　figure,plot(t,SoundDa

50、ta);　　%CurrentSection = CurrentSection + 1　　% --- Executes on button press in pushbutton2.　　function Sure_Callback(hObject, eventdata, handles)&l

51、t;p>　　% hObject handle to pushbutton2 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (se

52、e GUIDATA)　　global Myrate ;　　global Mycode;　　Mycode = str2num(get(handles.edit1,'string'));　　Myrate = str2num(get(findobj('t

53、ag','edit2'),'string'));　　%disp(Myrate);　　% --- Executes on button press in pushbutton3.　　function DFT_Callback(hObject, eventdata, handl

54、es)　　% hObject handle to pushbutton3 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles

55、and user data (see GUIDATA)　　global SampleRate;　　global SoundData;　　global Mycode;　　global y;　　N

56、 = Mycode;　　%disp(N);　　y = fft(SoundData,N);%進行fft變換　　mag=abs(y);%求幅值　　n= 1:N;f = n*SampleRate/length(y);%進行對應的頻率轉換&l

57、t;/p>　　figure;　　plot(f,mag);%做頻譜圖　　xlabel('頻率(Hz)');　　ylabel('幅值');　　title('頻譜'

58、);　　grid;　　%Py =2*(y.*conj(y))/N; %計算功率譜密度Py　　%figure;　　%plot(f,Py);　　%xla

59、bel('頻率(Hz)');　　%ylabel('功率譜密度');　　%title('x=120.0*cos(2*3.14*sf*t)功率譜密度');　　%grid;　　%point =

60、 get(handles.edit1,'string');　　%fs = SampleRate; 　　%m = length(SoundData); % Window length　　%n = pow2(nextpow2(m)); % Transform length<

61、;p>　　%y = fft(SoundData,N); % DFT%做n點的FFT，n為2的多少次冪，次冪值為大于或等于length(SoundData)的絕對值的2的最小次冪　　%n = (0:N-1) % Frequency range%圖上顯示的最高頻率為采樣頻率　　%power = y.*conj(y)/n; % Power of the

62、DFT%計算復共軛值　　%f = n*SampleRate/2;　　%figure,plot(f,y);　　%xlabel('Frequency (Hz)');　　%ylabel('DFT');

63、%title('{\bf 傅里葉變換}');　　%pause(1);　　% --- Executes on button press in pushbutton4.　　function DTpass_Callback(hObject, eventdata, handles)

64、　　% hObject handle to pushbutton4 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (

65、see GUIDATA)　　global y;　　global SampleRate;　　%%根據原始信號與噪聲信號頻譜圖，選擇巴特沃斯低通濾波器　　%%濾波器性能指標：取wp=1500*2*pi rad/s,ws=2400*2*pi

66、 rad/s,Rp=2,Rs=15　　%選擇濾波器的最小階數　　wp=1000*2*pi;ws=6000*2*pi;Rp=2;Rs=15;　　[N1,Wn]=buttord(wp,ws,Rp,Rs,'s');　　%創(chuàng)建butterworth模擬濾波器<

67、;/p>　　[Z,P,K]=buttap(N1);　　%把濾波器的零極點模型轉化為傳遞函數模型　　[Bap,Aap]=zp2tf(Z,P,K);　　%把模擬濾波器原型轉換成截止頻率為Wn的低通濾波器　　[b,a]=lp2lp(Bap,Aap,W

68、n);　　%用雙線性變換法實現模擬濾波器到數字濾波器的轉換　　[bz,az]=bilinear(b,a,SampleRate);　　%繪制頻率響應曲線　　[H,W]=freqz(bz,az);

69、　　figure(3);　　subplot(211),　　plot(W*SampleRate/(2*pi),abs(H)); grid 　　title('butterworth低通濾波器頻響圖');xlabel('頻率/Hz'); ylabel('頻率響應幅度'

70、;);　　%%%%%%%用所設計的濾波器對語音+噪聲信號y濾波%%%%%%%　　global Mycode;　　y1=filter(bz,az,y);　　n=(0:Mycode-1)*SampleRate/length(y1);

71、　subplot(212),plot(n,y1);title('濾波后信號波形圖');　　xlabel('頻率/Hz'); ylabel('頻率響應幅度');　　%sound(x,fs);%播放語音信號　　%sound(y,fs);%播放語音+噪聲信號</p

72、>　　%sound(y1,fs);%播放濾波后的信號　　% --- Executes on button press in pushbutton5.　　function quit_Callback(hObject, eventdata, handles)　　% hObject h

73、andle to pushbutton5 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (see GUIDATA)&

74、lt;p>　　close;　　msgbox('謝謝使用，再見！');　　function edit1_Callback(hObject, eventdata, handles)　　% hObject handle to edit1 (see GCB

75、O)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (see GUIDATA)　　% Hints: get(hObjec

76、t,'String') returns contents of edit1 as text　　% str2double(get(hObject,'String')) returns contents of edit1 as a double　　%A=get(hObject,'string');

77、　　%guidata(hObject,handles); %必須要加這一句　　% --- Executes during object creation, after setting all properties.　　function edit1_CreateFcn(hObject, eventdata, han

78、dles)　　% hObject handle to edit1 (see GCBO)　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles empty - handles not created

79、 until after all CreateFcns called　　% Hint: edit controls usually have a white background on Windows.　　% See ISPC and COMPUTER.　　if ispc</b

80、>　　set(hObject,'BackgroundColor','white');　　else　　set(hObject,'BackgroundColor',get(0,'defaultUicontrolBack

81、groundColor'));　　end　　function edit2_Callback(hObject, eventdata, handles)　　% hObject handle to edit2 (see GCBO)<p&g

82、t;　　% eventdata reserved - to be defined in a future version of MATLAB　　% handles structure with handles and user data (see GUIDATA)　　% Hints: get(hObject,'String') r

83、eturns contents of edit2 as text　　% str2double(get(hObject,'String')) returns contents of edit2 as a double　　% --- Executes during object creation, after setting a

84、ll properties.　　function edit2_CreateFcn(hObject, eventdata, handles)　　% hObject handle to edit2 (see GCBO)　　% eventdata reserved - to be defined in a fu

85、ture version of MATLAB　　% handles empty - handles not created until after all CreateFcns called　　% Hint: edit controls usually have a white background on Windows.&

86、lt;p>　　% See ISPC and COMPUTER.　　if ispc　　set(hObject,'BackgroundColor','white');　　else</p&

眾賞文庫> 全部分類> 畢業(yè)設計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

課程設計---iir數字低通濾波器的matlab實現

文檔簡介

溫馨提示

最新文檔

評論

課程設計---iir數字低通濾波器的matlab實現

文檔簡介

溫馨提示

最新文檔

評論

免費下載