【正文】 發(fā)現(xiàn)頻譜圖在0Hz時的幅度增加的很大，而且又在沒有譜線的頻率上竟然出現(xiàn)了頻譜，這是由于白噪聲在所有頻率上都有頻率造成的。（5） 工頻干擾（50Hz）用x表示加入工頻干擾后的信號，再對工頻信號的1000個頻率點(diǎn)進(jìn)行FFT運(yùn)算，在進(jìn)行相關(guān)的運(yùn)算后，通過調(diào)用plot函數(shù)直接繪出加入工頻干擾后的時域圖和頻譜圖。具體步驟如下：x2=[sin(2*pi*50*t)]。 t=0::*(10001)。x1=a(:,2)。x=x1+x2。y2=fft(x2,1000)。f2=100*(0:999)/1000。figure(3)。subplot(2,1,1)。plot(t,x)。title(39。加入工頻干擾后的時域圖39。)。xlabel(39。時間（s）39。)。ylabel(39。幅值（A）39。)。subplot(2,1,2)。plot(f2,abs(y2))。title(39。加入工頻干擾后的頻譜圖39。)。xlabel(39。幅度39。)。ylabel(39。Hz39。)。 加入工頻干擾后的時域圖和頻譜圖，發(fā)現(xiàn)加噪聲后的幅值有比較微小的變化，而在頻譜圖在38Hz和63Hz附近增值很大，這是由于50Hz的工頻噪聲造成的，而在其他頻率范圍內(nèi)也有比較明顯的變化。（6） 帶限chirp噪聲 同樣用k表示加入帶限chirp噪聲后的信號，再用freqz進(jìn)行濾波器的頻率響應(yīng)，調(diào)用plot對k畫出加入帶限噪聲后的時域和頻域圖，具體描述如下：p=*chirp(a(:,1),0,a(1000,1),200)。k=x1+p。[hc,wc]=freqz(k,1,1024)。abshc=abs(hc)。figure(4)。subplot(2,1,1)。plot(k(1:1000))。title(39。加入帶限噪聲后的時域圖39。)。xlabel(39。時間（s）39。)。ylabel(39。幅值（A）39。)。subplot(2,1,2)。plot(wc/pi,abshc)。title(39。加入帶限噪聲后的頻譜圖39。)。xlabel(39。幅度39。)。ylabel(39。Hz39。)。 加入帶限噪聲后的時域圖和頻譜圖加入帶限chirp噪聲干擾后，發(fā)現(xiàn)加噪之前時域圖是均勻分布的，而加噪后，則變?yōu)榍懊嫦∈韬竺婷芗那闆r了，對于頻譜圖而言，加噪前是中間凹兩邊凸，加噪后是前凸后平，在個別點(diǎn)上幅值增加很多，這是由于加入chirp噪聲的結(jié)果。（6）濾波器的設(shè)計(jì)Butterworth型低通數(shù)字濾波器用wp和ws表示分別將通帶，阻帶截止頻率的角頻率表示，在分別計(jì)算階數(shù)n1和截止頻率Wn,再設(shè)計(jì)低通Butterworth型模擬濾波器，然后采用雙線性法將模擬濾波器系數(shù)變?yōu)閿?shù)字濾波器系數(shù)，畫出濾波器頻譜圖，調(diào)用filter實(shí)現(xiàn)對工頻干擾的濾波，用plot函數(shù)畫出濾除工頻干擾后的時域圖和濾除白噪聲后的頻譜圖。具體操作如下：figure(5)。fs=100。f1=5。f2=10。wp=(f1/fs)*2*pi。ws=(f2/fs)*2*pi。Omegap=2*fs*tan(wp/2)。Omegas=2*fs*tan(ws/2)。[n1,Wn]=buttord(Omegap,Omegas,1,50,39。s39。)。[b,a]=butter(n1,Wn,39。s39。)。[bz,az]=bilinear(b,a,fs)。freqz(bz,az,512,fs)。y=filter(bz,az,x2)。figure(6)。subplot(2,1,1)plot(y)。title(39。濾除工頻干擾后的時域圖39。)。xlabel(39。時間（s）39。)。ylabel(39。幅值（A）39。)。y3=fft(y,1000)。f1=100*(0:999)/1000。subplot(2,1,2)。plot(f1,abs(y3))。title(39。濾除白噪聲后的頻譜圖39。)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值(dB)39。)。 Butterworth型低通數(shù)字濾波器頻率響應(yīng)圖 濾除單頻正弦波后的時域圖和濾除噪聲后的頻譜圖由于信號處于頻段的低頻部分，而工頻信號的頻譜在整個上是呈對稱分布的，與源信號的頻段分布是相似的，采用低通濾波器將噪聲信號的高頻部分濾掉，濾波后時域圖的幅值比原來變小了，濾波后的頻譜圖頻率在0~10Hz和90~100Hz之間的幅值變化較大，而在10~90Hz的區(qū)間，頻譜圖基本趨于直線。（7）切比雪夫I型數(shù)字低通濾波器用Wp1,Wp2,Ws1,Ws2表示分別將通帶，阻帶截止頻率的角頻率表示,算出頻帶寬帶，計(jì)算階數(shù)n1和截止頻率WN，再設(shè)計(jì)切比雪夫I型模擬濾波器，采用雙線性法將模擬濾波器系數(shù)變?yōu)閿?shù)字濾波器系數(shù)，畫出切比雪夫I型數(shù)字濾波器的頻率響應(yīng)，調(diào)用filter實(shí)現(xiàn)對白噪聲的濾波，再最后調(diào)用plot函數(shù)畫出濾除白噪聲后的時域圖和頻域圖。具體過程如下：figure(7)。fs=100。f11=10。f12=25。f21=5。f22=30。Wp1=(f11/fs)*2*pi。Ws1=(f21/fs)*2*pi。Wp2=(f12/fs)*2*pi。Ws2=(f22/fs)*2*pi。Omegap1=2*fs*tan(Wp1/2)。Omegap2=2*fs*tan(Wp2/2)。Omegas1=2*fs*tan(Ws1/2)。Omegas2=2*fs*tan(Wp2/2)。BW=Omegap2Omegap1。W0=Omegap1*Omegap2。W00=sqrt(W0)。WP=1。WS=WP*(W0^2Ws1^2)/(Ws1*BW)。[n1,WN]=buttord(WP,WS,1,50,39。s39。)。[B,A]=cheby1(n1,1,WN,39。s39。)。[BT,AT]=lp2bp(B,A,W00,BW)。[num,den]=bilinear(BT,AT,)。freqz(num,den,64)。y=filter(num,den,s)。figure(8)。subplot(2,1,1)。plot(y)。title(39。濾除白噪聲后的時域圖39。)。xlabel(39。時間（s）39。)。ylabel(39。幅值（A）39。)。s3=fft(y,1000)。f1=100*(0:999)/1000。subplot(2,1,2)。plot(f1,abs(s3))。title(39。濾除白噪聲后的頻譜圖39。)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值(dB)39。)。 切比雪夫1型數(shù)字濾波器的頻率響應(yīng)圖 濾除白噪聲后的時域圖和頻域圖由于白噪聲在整個頻段上都存在，對與源信號共存的低頻信號用選頻濾波器是無法濾除的，故而采用低通濾波器將噪聲的高頻部分去掉。，它們在時域圖上差不多是一致的，在濾波后的就明顯變小了。（8） 帶阻濾波器的設(shè)計(jì) 設(shè)計(jì)帶阻濾波器，再畫出濾波器的頻率響應(yīng)圖，然后用filter實(shí)現(xiàn)對工頻干擾信號的濾波，調(diào)用plot函數(shù)分別畫出濾除工頻干擾后的時域圖和頻譜圖。具體操作如下：figure(9)。d=fir1(1000,[ ],39。stop39。)。freqz(d,512)。y4=filter(d,1,x2)。figure(10)。subplot(2,1,1)。plot(y4)。title(39。濾除工頻干擾后的時域圖39。)。xlabel(39。時間（s）39。)。ylabel(39。幅值（A）39。)。subplot(2,1,2)。y5=fft(y4,1000)。f1=100*(0:999)/1000。plot(f1,abs(y5))。title(39。濾除工頻干擾后的頻譜圖39。)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值(dB)39。)。 濾除工頻干擾后的時域圖和頻譜圖帶阻濾波器濾除工頻干擾在28Hz和74Hz的地方最為顯著，帶阻濾波器由于工頻干擾信號是對稱的，而且與源信號的高頻部分在相同的頻段上，因此加噪后高頻部分的頻段上的幅值會更高，故我們選擇帶阻濾波器來濾除噪聲干擾的作用。設(shè)計(jì)心得 課程設(shè)計(jì)是培養(yǎng)我們綜合運(yùn)用所學(xué)知識,發(fā)現(xiàn),提出,分析和解決實(shí)際問題,鍛煉實(shí)踐能力的重要環(huán)節(jié),是對我們實(shí)際工作能力的具體訓(xùn)練和考察過程，通過課程設(shè)計(jì)我們能夠比較系統(tǒng)的了解理論知識，把理論和實(shí)踐相結(jié)合，并且用到生活當(dāng)中。在做設(shè)計(jì)的過程中總會出現(xiàn)各種問題，在這種情況下我們都會努力尋求最佳路徑解決問題，無形間提高了我們的動手，動腦能力，并且同學(xué)之間還能相互探討問題，研究解決方案，增進(jìn)大家的團(tuán)隊(duì)意識。 此次課程設(shè)計(jì)綜合了數(shù)字信號處理的理論知識和實(shí)踐經(jīng)驗(yàn)，經(jīng)過這段時間查閱的相關(guān)資料并結(jié)合以前學(xué)過的數(shù)字信號處理課程和實(shí)驗(yàn)，我如期完成了設(shè)計(jì)任務(wù)要求，并加深了對知識的理解，提高了自己的實(shí)踐動手能力。拿到課程設(shè)計(jì)任務(wù)的時候，不知道從哪兒入手，選定題目后，上網(wǎng)查閱了相關(guān)資源，去圖書館搜集資料，確定了設(shè)計(jì)的基本方案，步驟，開始用軟件進(jìn)行實(shí)現(xiàn)，實(shí)現(xiàn)過程中出現(xiàn)了很多問題，上網(wǎng)查閱資料，在專業(yè)的論壇里面去看，下載MATLAB的視頻，慢慢的開始有了進(jìn)一步的理解。在以前的實(shí)驗(yàn)中沒有接觸過GUI界面的問題，根本不知道從哪入手，在查閱了資料，看視頻，有了一定的了解，但是由于初次接觸，出現(xiàn)問題的時候還是很難解決，經(jīng)過和同學(xué)討論，終于解決了這些問題。經(jīng)過此次課程設(shè)計(jì)，使我對數(shù)字信號處理的理論知識有了更深了的了解，更近一步的理解了用MATLAB語言對信號進(jìn)行采樣，濾波等處理。初次接觸了MATLAB提供的用戶界面的設(shè)計(jì)問題，對其中的一些基本問題有了了解，為以后的實(shí)踐打下了一定的基礎(chǔ)，相信在以后的實(shí)踐過程中會有更好的提升。 另外，此次課程設(shè)計(jì)讓我更加明白平時所學(xué)的知識如果不加以實(shí)踐的話等于紙上談兵。課程設(shè)計(jì)主要是我們理論知識的延伸，它的目的主要是要在設(shè)計(jì)中發(fā)現(xiàn)問題，并且自己要能找到解決問題的方案，形成一種獨(dú)立的意識。我們還能從設(shè)計(jì)中檢驗(yàn)我們所學(xué)的理論知識到底有多少，鞏固我們已經(jīng)學(xué)會的，不斷學(xué)習(xí)我們所遺漏的新知識，把這門課學(xué)的扎實(shí)。附錄function varargout = GUIDE(varargin)% GUIDE Mfile for % GUIDE, by itself, creates a new GUIDE or raises the existing% singleton*.%% H = GUIDE returns the handle to a new GUIDE or the handle to% the existing singleton*.%% GUIDE(39。CALLBACK39。,hObject,eventData,handles,...) calls the local% function named CALLBACK in with the given input arguments.%% GUIDE(39。Property39。,39。Value39。,...) creates a new GUIDE or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before GUIDE_OpeningFunction gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to GUIDE_OpeningF via varargin.%% *See GUI Options on GUIDE39。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 GUIDE% Last Modified by GUIDE 12Sep2013 09:13:19% Begin initialization code DO NOT EDITgui_Singleton = 1。gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @GUIDE_OpeningF, ... 39。gui_OutputF39。, @GUIDE_OutputF, ... 39。gui_LayoutF39。, [] , ... 39。gui_Callback39。, [])。if nargin amp。amp。 ischar(varargin{1}) = str2func(varargin{1})。endif nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。else gui_mainf(gui_State, varargin{:})。end% End initialization code DO NOT EDIT% Execut