matlab【濾波器】程序

1、matlab濾波器程序wc=(400/1000)*pi;                             %求截止頻率w1=boxcar(81);              

2、;                  %窗函數(shù)w2=triang(81);w3=hamming(81);w4=hanning(81);w5=bartlett(81);w6=blackman(81);w7=chebwin(81,30);w8=kaiser(81,7.856);n=1:1:81;hd=sin(wc*(n-41)./(pi*(n-41);      

3、0;         %求h(d)hd(41)=wc/pi;h1=hd.*w1'                                    %加窗h2

4、=hd.*w2'h3=hd.*w3'h4=hd.*w4'h5=hd.*w5'h6=hd.*w6'h7=hd.*w7'h8=hd.*w8'mag1,rad=freqz(h1);                           %求幅頻特性曲線mag2,rad=freqz(h2);m

5、ag3,rad=freqz(h3);mag4,rad=freqz(h4);mag5,rad=freqz(h5);mag6,rad=freqz(h6);mag7,rad=freqz(h7);mag8,rad=freqz(h8);figure(1);                             

6、60;        %畫(huà)幅頻特性曲線plot(rad,20*log10(abs(mag1);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-80,0);title('利用矩形窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(2);plot(rad,20*log10(abs(mag2);xlabel('Normalized Frequency(rad)&

7、#39;);ylabel('Normaliaed Magnitude(dB)');axis(0,3,-80,0);title('利用三角窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(3);plot(rad,20*log10(abs(mag3);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-80,0);title('利用海明設(shè)計(jì)的數(shù)字濾波器');grid on;figure(4);plot(ra

8、d,20*log10(abs(mag4);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-80,0);title('利用漢寧窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(5);plot(rad,20*log10(abs(mag5);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-80

9、,0);title('利用巴特里特窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(6);plot(rad,20*log10(abs(mag6);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-100,0);title('利用布拉克曼窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(7);plot(rad,20*log10(abs(mag7);xlabel('Normalized Frequency(r

10、ad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-100,0);title('利用切比雪夫窗設(shè)計(jì)的數(shù)字濾波器');grid on;figure(8);plot(rad,20*log10(abs(mag8);xlabel('Normalized Frequency(rad)');ylabel('Normaliaed Magnitude(dB)');axis(0,3,-100,0);title('利用凱塞窗設(shè)計(jì)的數(shù)字濾波器');grid on;MATLAB參考

11、程序和仿真內(nèi)容%*%mode: 1-巴特沃斯低通；2-巴特沃斯高通；3-巴特沃斯帶通；4-巴特沃斯帶阻% 5-契比雪夫低通；6-契比雪夫高通；7-契比雪夫帶通；8-契比雪夫帶阻%fp1,fp2： 通帶截止頻率，當(dāng)高通或低通時(shí)只有fp1有效%fs1, fs2： 阻帶截止頻率，當(dāng)高通或低通時(shí)只有fs1有效%rp: 通帶波紋系數(shù)%as: 阻帶衰減系數(shù)%sample: 采樣率%h: 返回設(shè)計(jì)好的濾波器系數(shù)%*%functionb,a=iirfilt(mode,fp1,fp2,fs1,fs2,rp,as,sample)wp1=2*fp1/sample;wp2=2*fp2/sample;ws1=2*fs1

12、/sample;ws2=2*fs2/sample;%得到巴特沃斯濾波器的最小階數(shù)N和3bd頻率wnif mode<3N,wn=buttord(wp1,ws1,rp,as);elseif mode<5N,wn=buttord(wp1 wp2,ws1 ws2,rp,as);%得到契比雪夫?yàn)V波器的最小階數(shù)N和3bd頻率wnelseif mode<7N,wn=cheb1ord(wp1,ws1,rp,as);elseN,wn=cheblord(wp1 wp2,ws1 ws2,rp,as);end%得到濾波器系數(shù)的分子b和分母aif mode= =1b,a=butter(N,wn);en

13、dif mode= =2b,a=butter(N,wn,/high/);endif mode= =3b,a=butter(N,wn);endif mode= =4b,a=butter(N,wn,/stop/);endif mode= =5b,a=cheby1(N,rp,wn);endif mode= =6b,a=cheby1(N,rp,wn,/high/);endif mode= =7b,a=cheby1(N,rp,wn);endif mode= =8b,a=cheby1(N,rp,wn,/stop/);endset(gcf,/menubar/,menubar);freq_response=f

14、reqz(b,a);magnitude=20*log10(abs(freq_response);m=0:511;f=m*sample/(2*511);subplot(3,1,1);plot(f,magnitude);grid; %幅頻特性axis(0 sample/2 1.1*min(magnitude) 1.1*max(magnitude);ylabel('Magnitude');xlabel('Frequency->');phase=angle(freq_response);subplot(3,1,2);plot(f,phase);grid; %相頻特

15、性axis(0 sample/2 1.1*min(phase) 1.1*max(phase);ylabel('Phase');xlabel('Frequency->');h=impz(b,a,32); %32點(diǎn)的單位函數(shù)響應(yīng)t=1:32;subplot(3,1,3);stem(t,h);grid;axis(0 32 1.2*min(h) 1.1*max(h);ylabel('h(n)');xlabel('n->');基于MATLAB信號(hào)處理工具箱的數(shù)字濾波器設(shè)計(jì)與仿真摘要：傳統(tǒng)的數(shù)字濾波器的設(shè)計(jì)過(guò)程復(fù)雜，計(jì)算工作量大，

16、濾波特性調(diào)整困難，影響了它的應(yīng)用。本文介紹了一種利用matlab信號(hào)處理工具箱（signal processing toolbox）快速有效的設(shè)計(jì)由軟件組成的常規(guī)數(shù)字濾波器的設(shè)計(jì)方法。給出了使用matlab語(yǔ)言進(jìn)行程序設(shè)計(jì)和利用信號(hào)處理工具箱的fdatool工具進(jìn)行界面設(shè)計(jì)的詳細(xì)步驟。利用matlab設(shè)計(jì)濾波器，可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性調(diào)整參數(shù)，直觀簡(jiǎn)便，極大的減輕了工作量，有利于濾波器設(shè)計(jì)的最優(yōu)化。本文還介紹了如何利用matlab環(huán)境下的仿真軟件simulink對(duì)所設(shè)計(jì)的濾波器進(jìn)行模擬仿真。關(guān)鍵詞：數(shù)字濾波器 matlab fir iir引言：在電力系統(tǒng)微機(jī)保護(hù)和二次控制中，很多信號(hào)

17、的處理與分析都是基于對(duì)正弦基波和某些整次諧波的分析，而系統(tǒng)電壓電流信號(hào)（尤其是故障瞬變過(guò)程）中混有各種復(fù)雜成分，所以濾波器一直是電力系統(tǒng)二次裝置的關(guān)鍵部件【1】。目前微機(jī)保護(hù)和二次信號(hào)處理軟件主要采用數(shù)字濾波器。傳統(tǒng)的數(shù)字濾波器設(shè)計(jì)使用繁瑣的公式計(jì)算，改變參數(shù)后需要重新計(jì)算，在設(shè)計(jì)濾波器尤其是高階濾波器時(shí)工作量很大。利用matlab信號(hào)處理工具箱（signal processing toolbox）可以快速有效的實(shí)現(xiàn)數(shù)字濾波器的設(shè)計(jì)與仿真。1 數(shù)字濾波器及傳統(tǒng)設(shè)計(jì)方法數(shù)字濾波器可以理解為是一個(gè)計(jì)算程序或算法，將代表輸入信號(hào)的數(shù)字時(shí)間序列轉(zhuǎn)化為代表輸出信號(hào)的數(shù)字時(shí)間序列，并在轉(zhuǎn)化過(guò)程中，使信號(hào)按

18、預(yù)定的形式變化。數(shù)字濾波器有多種分類，根據(jù)數(shù)字濾波器沖激響應(yīng)的時(shí)域特征，可將數(shù)字濾波器分為兩種，即無(wú)限長(zhǎng)沖激響應(yīng)（iir）濾波器和有限長(zhǎng)沖激響應(yīng)（fir）濾波器。iir數(shù)字濾波器具有無(wú)限寬的沖激響應(yīng)，與模擬濾波器相匹配。所以iir濾波器的設(shè)計(jì)可以采取在模擬濾波器設(shè)計(jì)的基礎(chǔ)上進(jìn)一步變換的方法。fir數(shù)字濾波器的單位脈沖響應(yīng)是有限長(zhǎng)序列。它的設(shè)計(jì)問(wèn)題實(shí)質(zhì)上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖響應(yīng)的常數(shù)問(wèn)題，設(shè)計(jì)方法主要有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等。在對(duì)濾波器實(shí)際設(shè)計(jì)時(shí)，整個(gè)過(guò)程的運(yùn)算量是很大的。例如利用窗函數(shù)法【2】設(shè)計(jì)m階fir低通濾波器時(shí)，首先要根據(jù)（1）式計(jì)算出理想低通濾波器的單

19、位沖激響應(yīng)序列，然后根據(jù)（2）式計(jì)算出m個(gè)濾波器系數(shù)。當(dāng)濾波器階數(shù)比較高時(shí)，計(jì)算量比較大，設(shè)計(jì)過(guò)程中改變參數(shù)或?yàn)V波器類型時(shí)都要重新計(jì)算。                             （1）          

20、                              （2）設(shè)計(jì)完成后對(duì)已設(shè)計(jì)的濾波器的頻率響應(yīng)要進(jìn)行校核，要得到幅頻相頻響應(yīng)特性，運(yùn)算量也是很大的。我們平時(shí)所要設(shè)計(jì)的數(shù)字濾波器，階數(shù)和類型并不一定是完全給定的，很多時(shí)候都是要根據(jù)設(shè)計(jì)要求和濾波效果不斷的調(diào)整，以達(dá)到設(shè)計(jì)的最優(yōu)化。在這種情況下，濾波器的設(shè)

21、計(jì)就要進(jìn)行大量復(fù)雜的運(yùn)算，單純的靠公式計(jì)算和編制簡(jiǎn)單的程序很難在短時(shí)間內(nèi)完成設(shè)計(jì)。利用matlab強(qiáng)大的計(jì)算功能進(jìn)行計(jì)算機(jī)輔助設(shè)計(jì)，可以快速有效的設(shè)計(jì)數(shù)字濾波器，大大的簡(jiǎn)化了計(jì)算量，直觀簡(jiǎn)便。2數(shù)字濾波器的matlab設(shè)計(jì)2.1    fdatool界面設(shè)計(jì)2.1.1   fdatool的介紹fdatool（filter design & analysis tool）是matlab信號(hào)處理工具箱里專用的濾波器設(shè)計(jì)分析工具，matlab6.0以上的版本還專門(mén)增加了濾波器設(shè)計(jì)工具箱（filter design toolbox）。fdatoo

22、l可以設(shè)計(jì)幾乎所有的基本的常規(guī)濾波器，包括fir和iir的各種設(shè)計(jì)方法。它操作簡(jiǎn)單，方便靈活。fdatool界面總共分兩大部分，一部分是design filter，在界面的下半部，用來(lái)設(shè)置濾波器的設(shè)計(jì)參數(shù)，另一部分則是特性區(qū)，在界面的上半部分，用來(lái)顯示濾波器的各種特性。design filter部分主要分為：filter type（濾波器類型）選項(xiàng)，包括lowpass（低通）、highpass（高通）、bandpass（帶通）、bandstop（帶阻）和特殊的fir濾波器。design method（設(shè)計(jì)方法）選項(xiàng)，包括iir濾波器的butterworth（巴特沃思）法、chebyshev t

23、ype i（切比雪夫i型）法、 chebyshev type ii（切比雪夫ii型） 法、elliptic（橢圓濾波器）法和fir濾波器的equiripple法、least-squares（最小乘方）法、window（窗函數(shù)）法。filter order（濾波器階數(shù)）選項(xiàng)，定義濾波器的階數(shù)，包括specify order（指定階數(shù)）和minimum order（最小階數(shù)）。在specify order中填入所要設(shè)計(jì)的濾波器的階數(shù)（n階濾波器，specify ordern-1），如果選擇minimum order則matlab根據(jù)所選擇的濾波器類型自動(dòng)使用最小階數(shù)。frenquency spec

24、ifications選項(xiàng)，可以詳細(xì)定義頻帶的各參數(shù)，包括采樣頻率fs和頻帶的截止頻率。它的具體選項(xiàng)由filter type選項(xiàng)和design method選項(xiàng)決定，例如bandpass（帶通）濾波器需要定義fstop1（下阻帶截止頻率）、fpass1（通帶下限截止頻率）、fpass2（通帶上限截止頻率）、fstop2（上阻帶截止頻率），而lowpass（低通）濾波器只需要定義fstop1、fpass1。采用窗函數(shù)設(shè)計(jì)濾波器時(shí)，由于過(guò)渡帶是由窗函數(shù)的類型和階數(shù)所決定的，所以只需要定義通帶截止頻率，而不必定義阻帶參數(shù)。magnitude specifications選項(xiàng)，可以定義幅值衰減的情況。例

25、如設(shè)計(jì)帶通濾波器時(shí)，可以定義wstop1（頻率fstop1處的幅值衰減）、wpass（通帶范圍內(nèi)的幅值衰減）、wstop2（頻率fstop2處的幅值衰減）。當(dāng)采用窗函數(shù)設(shè)計(jì)時(shí)，通帶截止頻率處的幅值衰減固定為6db，所以不必定義。window specifications選項(xiàng)，當(dāng)選取采用窗函數(shù)設(shè)計(jì)時(shí)，該選項(xiàng)可定義，它包含了各種窗函數(shù)。2.1.2  帶通濾波器設(shè)計(jì)實(shí)例本文將以一個(gè)fir 濾波器的設(shè)計(jì)為例來(lái)說(shuō)明如何使用matlab設(shè)計(jì)數(shù)字濾波器：在小電流接地系統(tǒng)中注入83.3hz的正弦信號(hào)，對(duì)其進(jìn)行跟蹤分析，要求設(shè)計(jì)一帶通數(shù)字濾波器，濾除工頻及整次諧波，以便在非常復(fù)雜的信號(hào)中分離

26、出該注入信號(hào)。參數(shù)要求：96階fir數(shù)字濾波器，采樣頻率1000hz，采用hamming窗函數(shù)設(shè)計(jì)。本例中，首先在filter type中選擇bandpass（帶通濾波器）；在design method選項(xiàng)中選擇fir window（fir濾波器窗函數(shù)法），接著在window specifications選項(xiàng)中選取hamming；指定filter order項(xiàng)中的specify order95；由于采用窗函數(shù)法設(shè)計(jì)，只要給出通帶下限截止頻率fc1和通帶上限截止頻率fc2，選取fc170hz，fc284hz。設(shè)置完以后點(diǎn)擊design filter即可得到所設(shè)計(jì)的fir濾波器。通過(guò)菜單選項(xiàng)ana

27、lysis可以在特性區(qū)看到所設(shè)計(jì)濾波器的幅頻響應(yīng)、相頻響應(yīng)、零極點(diǎn)配置和濾波器系數(shù)等各種特性。設(shè)計(jì)完成后將結(jié)果保存為1.fda文件。在設(shè)計(jì)過(guò)程中，可以對(duì)比濾波器幅頻相頻特性和設(shè)計(jì)要求，隨時(shí)調(diào)整參數(shù)和濾波器類型，以便得到最佳效果。其它類型的fir濾波器和iir濾波器也都可以使用fdatool來(lái)設(shè)計(jì)。圖1 濾波器幅頻和相頻響應(yīng)（特性區(qū)）fig.1 magnitude response and phase response of the filter2.2 程序設(shè)計(jì)法在matlab中，對(duì)各種濾波器的設(shè)計(jì)都有相應(yīng)的計(jì)算振幅響應(yīng)的函數(shù)【3】，可以用來(lái)做濾波器的程序設(shè)計(jì)。上例的帶通濾波器可以用程序設(shè)計(jì)：c

28、=95;                         定義濾波器階數(shù)96階w1=2*pi*fc1/fs;w2=2*pi*fc2/fs;                %參數(shù)轉(zhuǎn)換，將模擬濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為

29、數(shù)字濾波器的技術(shù)指標(biāo)window=hamming(c+1);         %使用hamming窗函數(shù)h=fir1(c,w1/pi w2/pi,window);   使用標(biāo)準(zhǔn)響應(yīng)的加窗設(shè)計(jì)函數(shù)fir1freqz(h,1,512);                 數(shù)字濾波器頻率響應(yīng)在matlab環(huán)境下運(yùn)行該程序即可得到濾波器幅頻相頻響應(yīng)

30、曲線和濾波器系數(shù)h。篇幅所限，這里不再將源程序詳細(xì)列出。3    simulink仿真本文通過(guò)調(diào)用simulink中的功能模塊構(gòu)成數(shù)字濾波器的仿真框圖，在仿真過(guò)程中，可以雙擊各功能模塊，隨時(shí)改變參數(shù)，獲得不同狀態(tài)下的仿真結(jié)果。例如構(gòu)造以基波為主的原始信號(hào)，通過(guò)simulink環(huán)境下的digital filter design（數(shù)字濾波器設(shè)計(jì)）模塊導(dǎo)入2.1.2中fdatool所設(shè)計(jì)的濾波器文件1.fda。仿真圖和濾波效果圖如圖2所示。  圖2 simulink仿真圖及濾波效果圖fig.2 simulated connections and wa

31、veform可以看到經(jīng)過(guò)離散采樣、數(shù)字濾波后分離出了83.3hz的頻率分量（scope1）。之所以選取上面的疊加信號(hào)作為原始信號(hào)，是由于在實(shí)際工作中是要對(duì)已經(jīng)經(jīng)過(guò)差分濾波的信號(hào)進(jìn)一步做帶通濾波，信號(hào)的各分量基本同一致，可以反映實(shí)際的情況。本例設(shè)計(jì)的濾波器已在實(shí)際工作中應(yīng)用，取得了不錯(cuò)的效果。4   結(jié)論利用matlab的強(qiáng)大運(yùn)算功能，基于matlab信號(hào)處理工具箱（signal processing toolbox）的數(shù)字濾波器設(shè)計(jì)法可以快速有效的設(shè)計(jì)由軟件組成的常規(guī)數(shù)字濾波器，設(shè)計(jì)方便、快捷，極大的減輕了工作量。在設(shè)計(jì)過(guò)程中可以對(duì)比濾波器特性，隨時(shí)更改參數(shù)，以達(dá)

32、到濾波器設(shè)計(jì)的最優(yōu)化。利用matlab設(shè)計(jì)數(shù)字濾波器在電力系統(tǒng)二次信號(hào)處理軟件和微機(jī)保護(hù)中，有著廣泛的應(yīng)用前景。設(shè)計(jì)巴特沃斯濾波器的MATLAB函數(shù)2009-05-09 21:121. buttord（1）N,wc=buttord(wp，ws，p，s)       用于計(jì)算巴特沃斯數(shù)字濾波器的階數(shù)N和3dB截止頻率wc。       調(diào)用參數(shù)wp，ws分別為數(shù)字濾波器的通帶、阻帶截止頻率的歸一化值，要求：0wp1，0ws1。1表示數(shù)字頻率pi。  

33、0;    p，s分別為通帶最大衰減和組帶最小衰減（dB）。       當(dāng)wswp時(shí)，為高通濾波器；       當(dāng)wp和ws為二元矢量時(shí)，為帶通或帶阻濾波器，這時(shí)wc也是二元向量。       N，wc作為butter函數(shù)的調(diào)用參數(shù)。（2）N,c=buttord(p，s，p，s，s)        用于計(jì)算巴特沃斯

34、模擬濾波器的階數(shù)N和3dB截止頻率c。       p，s，c均為實(shí)際模擬角頻率。       說(shuō)明：buttord函數(shù)使用阻帶指標(biāo)計(jì)算3dB截止頻率，這樣阻帶會(huì)剛好滿足要求，而通帶會(huì)有富余。2.buttap（N）      z0，p0，k0=buttap(N)      用于計(jì)算N階巴特沃斯歸一化（3dB截止頻率c=1）模擬低通原型濾波器系統(tǒng)函數(shù)的零、極點(diǎn)和增益因子。

35、0;     說(shuō)明：如果要從零、極點(diǎn)模型得到系統(tǒng)函數(shù)的分子、分母多項(xiàng)式系數(shù)向量ba、aa，可調(diào)用                 B，A=zp2tf（z0，p0，k0）3.butter（1）b，a=butter（N，wc，ftype）       計(jì)算N階巴特沃斯數(shù)字濾波器系統(tǒng)函數(shù)分子、分母多項(xiàng)式的系數(shù)向量b、a。  

36、    調(diào)用參數(shù)N和wc分別為巴特沃斯數(shù)字濾波器的階數(shù)和3dB截止頻率的歸一化值（關(guān)于pi歸一化），一般是調(diào)用buttord（1）格式計(jì)算N和wc。       系數(shù)b、a是按照z-1的升冪排列。（2）B，A=butter（N，c，ftype，s）         計(jì)算巴特沃斯模擬濾波器系統(tǒng)函數(shù)的分子、分母多項(xiàng)式系數(shù)向量ba、aa。        

37、 調(diào)用參數(shù)N和c分別為巴特沃斯模擬濾波器的階數(shù)和3dB截止頻率(實(shí)際角頻率)，可調(diào)用buttord（2）格式計(jì)算N和c。        系數(shù)B、A按s的正降冪排列。        tfype為濾波器的類型：        ftype=high時(shí)，高通；c只有1個(gè)值。        ftype=stop時(shí)，帶阻阻；此時(shí)c

38、=cl,cu，分別為帶阻濾波器的通帶3dB下截止頻率和上截止頻率。       ftype缺省時(shí)：         若c只有1個(gè)值，則默認(rèn)為低通；        若c有2個(gè)值，則默認(rèn)為帶通；其通帶頻率區(qū)間cl < < cu。             

39、;    注意：所設(shè)計(jì)的帶通和帶阻濾波器系統(tǒng)函數(shù)是2N階。因?yàn)閹V波器相當(dāng)于N階低通濾波器與N階高通濾波器級(jí)聯(lián)。MATLAB 數(shù)字濾波器設(shè)計(jì) 作業(yè)2009-06-25 20:01其實(shí)并不是說(shuō)為了交作業(yè)才寫(xiě)的程序，而是自己在做作業(yè)的過(guò)程當(dāng)中學(xué)到了新的東西，這才是最有樂(lè)趣的！對(duì)信號(hào)處理課的總結(jié)： 信號(hào)處理是對(duì)各種信號(hào)處理的總稱，用處是從看似雜亂的信號(hào)中提取出有用的信息。 主要理論是（離散）傅立葉變換，使用時(shí)采用快速傅立葉變換，計(jì)算速度會(huì)快很多。 對(duì)信號(hào)快速傅立葉變換后會(huì)得到幅值頻譜圖（通常是分離的線，表示信號(hào)組成成分），從中我們可以得到需要的信息（比如是否包含特定頻

40、率的波）。 濾波器是對(duì)信號(hào)進(jìn)行處理，過(guò)濾掉我們不需要的波，比如說(shuō)噪聲，更具體的比如圖像的噪點(diǎn)的過(guò)濾。原理就是讓只讓特定頻率范圍的波通過(guò)，否則就衰減的比較嚴(yán)重，從而實(shí)現(xiàn)了指定頻率的分離。濾波器的幾個(gè)重要參數(shù)的含義。 MATLAB軟件在工程上確實(shí)是有這比較大的優(yōu)勢(shì)，基本上常用的函數(shù)都已經(jīng)內(nèi)置了，從理念上實(shí)現(xiàn)了傻瓜化。用GOOGLE翻譯成英文如下：Summary of the Signal Processing class: signal processing is a general term for a variety of signal processing, the use of clutt

41、er from the signal appears to extract useful information. The main theory is the (discrete) Fourier transform, the use of fast Fourier transform used to calculate the speed will be much faster. Fast Fourier Transform of the signal after the amplitude of frequency spectrum will be (usually the separa

42、tion line, said signal components), from which we can get the information needed (such as whether a particular wave frequency). filter for signal processing, we do not need to filter out the wave, for example, noise, more specific images such as noise filters. The principle is to allow only a specif

43、ic frequency range of the wave through, or on the attenuation of the more serious, in order to achieve the specified frequency separation. Filter the meaning of several important parameters. MATLAB software in the works are indeed there is such a big advantage, is basically a function of commonly us

44、ed have been built from the concept of a fool.(機(jī)器翻譯確實(shí)還有很長(zhǎng)的路要走)題目要求：有一個(gè)信號(hào)，經(jīng)AD轉(zhuǎn)換后得到數(shù)據(jù)文件dataX（采樣頻率 ，數(shù)據(jù)點(diǎn)數(shù) ）（見(jiàn)表），請(qǐng)用Matlab完成以下工作：（1）確定此信號(hào)有哪些主要頻率成分；（2）用雙性變換法設(shè)計(jì)一個(gè)濾波器，保留數(shù)據(jù)中頻率最低的頻率成份，去掉其余成份；繪出濾波器的頻率響應(yīng)特性；繪出濾波后數(shù)據(jù)的時(shí)域圖形、幅值譜；（3）設(shè)計(jì)一個(gè)FIR Butterworth數(shù)字濾波器，保留數(shù)據(jù)中第3高的頻率成分，繪出濾波器的頻率響應(yīng)特性；繪出濾波后數(shù)據(jù)的時(shí)域圖形、幅值譜；完整MATLAB 代碼如下：%信號(hào)檢

45、測(cè)MATLAB上機(jī)作業(yè)%winxos 2009-6-25%讀取數(shù)據(jù)文件，我選擇的是data11data=load('data11'); %read file%基本參數(shù)設(shè)置Nums=512; %數(shù)據(jù)data11文件，采樣頻率400Hz，數(shù)據(jù)個(gè)數(shù)512個(gè)Fs=400;T=Nums/Fs;%設(shè)置窗口大小scrsz = get(0,'ScreenSize'); %get screen sizefigure('Position',100 100 scrsz(3)-200 scrsz(4)-200,. % .is the Link-Line symbols&

46、#160;   'Name','信號(hào)檢測(cè)MATLAB示例 LiWei','NumberTitle','off'); %set window property,the size is less than screen size%繪制原始波形x=0:1/Fs:T-1/Fs; subplot(221); %第一張圖plot(x,data);xlabel('Time (s)');ylabel('amplitude');title('原始波形','Color','r'