【正文】 defaultUicontrolBackgroundColor39。amp。String39。 t=get(,39。,39。), get(0,39。)) returns contents of edit9 as a double 課程設(shè)計(jì) 第 36 頁(yè) 共 44 頁(yè) % Executes during object creation, after setting all properties. function edit9_CreateF(hObject, eventdata, handles) % hObject handle to edit9 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles empty handles not created until after all CreateFs called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc amp。 end function edit9_Callback(hObject, eventdata, handles) % hObject handle to edit9 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,39。BackgroundColor39。BackgroundColor39。String39。)。)) set(hObject,39。 isequal(get(hObject,39。) returns contents of edit7 as text % str2double(get(hObject,39。white39。defaultUicontrolBackgroundColor39。amp。String39。,39。), get(0,39。)) returns contents of edit5 as a double % Executes during object creation, after setting all properties. function edit5_CreateF(hObject, eventdata, handles) % hObject handle to edit5 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles empty handles not created until after all CreateFs called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc amp。 end 課程設(shè)計(jì) 第 33 頁(yè) 共 44 頁(yè) function edit5_Callback(hObject, eventdata, handles) % hObject handle to edit5 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,39。BackgroundColor39。BackgroundColor39。String39。)。)) set(hObject,39。 isequal(get(hObject,39。) returns contents of edit3 as text % str2double(get(hObject,39。white39。defaultUicontrolBackgroundColor39。amp。String39。,39。), get(0,39。)) returns contents of edit1 as a double % Executes during object creation, after setting all properties. function edit1_CreateF(hObject, eventdata, handles) % hObject handle to edit1 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles empty handles not created until after all CreateFs called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc amp。 function edit1_Callback(hObject, eventdata, handles) % hObject handle to edit1 (see GCBO) % eventdata reserved to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,39。 % Update handles structure guidata(hObject, handles)。 ischar(varargin{1}) = str2func(varargin{1})。gui_Callback39。gui_OutputF39。gui_Singleton39。s Tools menu. Choose GUI allows only one % instance to run (singleton). % % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help lyl % Last Modified by GUIDE 25Jun2020 20:22:38 % Begin initialization code DO NOT EDIT gui_Singleton = 1。Property39。 課程設(shè)計(jì) 第 27 頁(yè) 共 44 頁(yè) 參考資料： 1謝平，林洪彬 信號(hào)處理原理及應(yīng)用 北京：機(jī)械工程出版社 2劉波 .MATLAB信號(hào)處理 [M].北京：電子工業(yè)出版社， 2020。在做實(shí)驗(yàn)是一定要小心謹(jǐn)慎，也許一個(gè)小小的語(yǔ)法錯(cuò)誤和常量變量的定義的錯(cuò)誤就造成整個(gè)程序出現(xiàn)問(wèn)題，得不到所需的波形，導(dǎo)致實(shí)驗(yàn)結(jié)果不正確。在這個(gè) 課程設(shè)計(jì)中我們須利用 MATLAB 設(shè)計(jì)各種濾波器。 注： Wph 和 Wsh 兩輸入只有在帶通和帯阻時(shí)，輸入?yún)?shù)有效。 2，在信號(hào)輸入通道中輸入濾波器的參數(shù)。 課程設(shè)計(jì) 第 25 頁(yè) 共 44 頁(yè) 軟件功能及使用方法。 grid。 y=y1+y2。String39。 m=str2double(get(,39。 帶通濾波器容許一定頻率范圍信號(hào)通過(guò) , 但減弱 (或減少 )頻率低于於下限截止頻率和高于上限截止頻率的信號(hào)的通過(guò)。 得出幅頻特性如下圖： 課程設(shè)計(jì) 第 22 頁(yè) 共 44 頁(yè) 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 10 . 20 . 40 . 60 . 811 . 21 . 41 . 6巴特沃斯帯阻數(shù)字濾波器F r e q u e n c y / H z 小結(jié)： 經(jīng)觀察濾波器幅頻特性圖得， 巴特沃斯濾波器的特點(diǎn)是 通帶 內(nèi)的頻率響應(yīng)曲線(xiàn)最大限度平坦，沒(méi)有起伏，而在阻頻帶則逐漸下降為零。)。返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %調(diào)用脈沖響應(yīng)不變法 )。 %調(diào)用 butter函數(shù)確定巴特沃斯濾波器階數(shù) [B,A]=butter(N,Wc,39。 As= 11。 fs=1/T。Magnitude39。 xlabel(39。 %繪出巴特沃斯數(shù)字濾波器的幅頻特性曲線(xiàn) grid on。 Hz=freqz(D,C,W)。)。 %調(diào)用 butter函數(shù)確定巴特沃斯濾波器階數(shù) [B,A]=butter(N,Wc,39。 As= 10。 %采樣周期 fs=1/T。)。Frequency/Hz39。 title(39。 Hz=freqz(D,C,W)。 %調(diào)用 butter函數(shù)設(shè)計(jì)巴特沃斯濾波器 W=linspace(0,pi,400*pi)。high39。 %設(shè)置通帶最大最小衰減 [N,Wc]=buttord(Wp,Ws,Ap,As,39。 %采樣頻率 課程設(shè)計(jì) 第 17 頁(yè) 共 44 頁(yè) Wp= .*pi。)。Frequency/Hz39。 title(39。 % %調(diào)用 butter函數(shù)設(shè)計(jì)巴特沃斯濾波器 W=linspace(0,pi,400*pi)。)。 %設(shè)計(jì)歸一化通帶阻帶截止頻率 Ap=1。 數(shù)字低通的技術(shù)指標(biāo)為： Wp= rad Ap=1dB T=2s Ws= rad As=15dB 程序?yàn)椋? T=2。通帶內(nèi)，分母 Ω /Ω c1,相應(yīng) (Ω /Ω c)2N 隨 N 的增加而趨于 0， A(Ω 2)→1 ， 在過(guò)渡帶和阻帶， Ω /Ω c1， 隨 N 的增加， Ω e/Ω c1，所以 A(Ω2)快速下降。在振幅的對(duì)數(shù)對(duì)角頻率的波得圖上，從某一邊界角頻率開(kāi)始，振幅隨著角頻率的增加而逐步減少，趨向負(fù)無(wú)窮大。 Z=est 的映射關(guān)系反映的是 Ha(s)的周期嚴(yán)拓與 H(z)的關(guān)系，而不是 Ha(s)b 本身與 H(z)的關(guān)系，因此，使用脈沖響應(yīng)不變法時(shí)，從 Ha(s)到 H(z)并沒(méi)有一個(gè)由 S 平面到 Z 平面的簡(jiǎn)單代數(shù)映射關(guān)系，即沒(méi)有一個(gè) s=f(z)的代數(shù)關(guān)系式。 脈沖響應(yīng)不變法是從濾波器的脈沖響應(yīng)出發(fā)，使數(shù)字濾波器的單位脈沖響應(yīng)序 h(n)模仿模擬濾波去的沖擊響應(yīng) ha(t)，使 h(n)正好等于 ha(t)的采樣值，即 h(n)=ha(nT) （ 23） T 為采樣周期。模擬濾波器設(shè)計(jì)已經(jīng)有了一套相當(dāng)成熟的方法，它不但有完整的設(shè)計(jì)公式，而且還有較為完整的圖表供查詢(xún)，因此，充分利用這些已有的資源將會(huì)給數(shù)字濾波器的設(shè)計(jì)帶來(lái)很大方便， IIR 數(shù)字濾波器的設(shè)計(jì)步驟是： 課程設(shè)計(jì) 第 11 頁(yè) 共 44 頁(yè) (1)按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)； (2)根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器 H(s)； (3)在按一定規(guī)則將 H(s)轉(zhuǎn)換為 H(z)。以低通濾波器特性為例，頻率響應(yīng)有通帶、過(guò)渡帶及阻帶三個(gè)范圍。 優(yōu)點(diǎn)： ? 簡(jiǎn)化實(shí)現(xiàn)，用一個(gè)二階節(jié)，通過(guò)變換