【正文】 )。wavwrite(y,fs,39。)。title(39。subplot(2,2,4)。延時(shí)后幅頻39。plot(abs(Z))。)。title(39。subplot(2,2,1:2)。%延遲3000?Z=fft(z,N)。N=length(y)z=[zeros(3000,1)。H:\39。)。ylabel(39。Hz39。原始語(yǔ)音信號(hào)頻譜39。)figure(3)plot(f,abs(y1(1:2048)))。)figure(2)freqz(x)title(39。title(39。原始信號(hào)幅值39。)subplot(3,1,2),plot(magy1)。title(39。angy1=angle(y1)。f=fs*(0:2047)/N。N=length(x)sound(x,fs,bits)。H:\39。)。wavwrite(f1,fs,39。fuzhi39。)。)xlabel(39。title(39。)。ylabel(39。Hz39。濾波前的頻譜39。plot(f,abs(y2(1:512)))。figure(3)y2=fft(x,1024)。F0=fft(f1,1024)。)。title(39。)。figure(2)subplot(2,1,1)plot(t,x)title(39。figure(1)freqz(b,1,512)[h1,w1]=freqz(b,1)。b=fir1(N,wn/pi,hamming(N+1))。t=0:(size(x)1)。wdelta=wswp。wp=*pi。h:\39。% Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fs=22050。h:\39。)。xlabel(39。濾波后的信號(hào)頻譜39。plot(f,abs(Z(1:n/2)))。Z=fft(z)。Hz39。)。title(39。subplot(2,1,1)。)。title(39。subplot(212)。原始信號(hào)的波形39。plot(x)。z=filter(bd,ad,x)。f=fs*(0:(n/21))/n。)。figure(2)[x,fs,bits]=wavread(39。濾波器的性能分析39。plot(w*fs/(2*pi),abs(h))grid。[h,w]=freqz(bd,ad)。[b,a]=lp2lp(Bap,Aap,Wn)。[Z,P,K]=buttap(N)。s39。ws1=2/Ts*tan(ws/2)。Ts=1/Fs。Rs=15。ws=*pi。fs=22050。% Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fb=1000。H:\39。 axis([0 4500 0 300])。fudu39。)。) xlabel(39。title(39。)。ylabel(39。Hz39。原始語(yǔ)音信號(hào)頻譜39。figure(2)subplot(2,1,1),plot(abs(y1))。)。ylabel(39。time n39。加高斯噪聲后語(yǔ)音信號(hào)時(shí)域圖39。)subplot(2,1,2) ,plot(x2)。title(39。sound(x2,fs,bits)。x1=[*sin(2*pi*10000*t)]39。f=fs*(0:511)/N。N=length(x)。H:\39。axis([0 4500 0 400])% Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fs=22050。fudu39。)。)xlabel(39。title(39。頻率響應(yīng)圖39。原始信號(hào)相位39。)subplot(3,1,3),plot(angy1)。title(39。原始信號(hào)波形39。subplot(3,1,1),plot(x)。figure(1)magy1=abs(y1)。y1=fft(x,N)。)。% Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[x,fs,bits]=wavread(39。% Outputs from this function are returned to the mand line.function varargout = GUIDE_OutputF(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT)。% Update handles structureguidata(hObject, handles)。else gui_mainf(gui_State, varargin{:})。 ischar(varargin{1}) = str2func(varargin{1})。if nargin amp。gui_Callback39。gui_LayoutF39。gui_OutputF39。gui_OpeningF39。gui_Singleton39。gui_Name39。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。Value39。Property39。CALLBACK39。我們還能從設(shè)計(jì)中檢驗(yàn)我們所學(xué)的理論知識(shí)到底有多少，鞏固我們已經(jīng)學(xué)會(huì)的，不斷學(xué)習(xí)我們所遺漏的新知識(shí)，把這門課學(xué)的扎實(shí)。 另外，此次課程設(shè)計(jì)讓我更加明白平時(shí)所學(xué)的知識(shí)如果不加以實(shí)踐的話等于紙上談兵。經(jīng)過(guò)此次課程設(shè)計(jì)，使我對(duì)數(shù)字信號(hào)處理的理論知識(shí)有了更深了的了解，更近一步的理解了用MATLAB語(yǔ)言對(duì)信號(hào)進(jìn)行采樣，濾波等處理。拿到課程設(shè)計(jì)任務(wù)的時(shí)候，不知道從哪兒入手，選定題目后，上網(wǎng)查閱了相關(guān)資源，去圖書館搜集資料，確定了設(shè)計(jì)的基本方案，步驟，開始用軟件進(jìn)行實(shí)現(xiàn)，實(shí)現(xiàn)過(guò)程中出現(xiàn)了很多問(wèn)題，上網(wǎng)查閱資料，在專業(yè)的論壇里面去看，下載MATLAB的視頻，慢慢的開始有了進(jìn)一步的理解。在做設(shè)計(jì)的過(guò)程中總會(huì)出現(xiàn)各種問(wèn)題，在這種情況下我們都會(huì)努力尋求最佳路徑解決問(wèn)題，無(wú)形間提高了我們的動(dòng)手，動(dòng)腦能力，并且同學(xué)之間還能相互探討問(wèn)題，研究解決方案，增進(jìn)大家的團(tuán)隊(duì)意識(shí)。 濾除工頻干擾后的時(shí)域圖和頻譜圖帶阻濾波器濾除工頻干擾在28Hz和74Hz的地方最為顯著，帶阻濾波器由于工頻干擾信號(hào)是對(duì)稱的，而且與源信號(hào)的高頻部分在相同的頻段上，因此加噪后高頻部分的頻段上的幅值會(huì)更高，故我們選擇帶阻濾波器來(lái)濾除噪聲干擾的作用。幅值(dB)39。)。xlabel(39。濾除工頻干擾后的頻譜圖39。plot(f1,abs(y5))。y5=fft(y4,1000)。)。ylabel(39。時(shí)間（s）39。)。title(39。subplot(2,1,1)。y4=filter(d,1,x2)。)。d=fir1(1000,[ ],39。（8） 帶阻濾波器的設(shè)計(jì) 設(shè)計(jì)帶阻濾波器，再畫出濾波器的頻率響應(yīng)圖，然后用filter實(shí)現(xiàn)對(duì)工頻干擾信號(hào)的濾波，調(diào)用plot函數(shù)分別畫出濾除工頻干擾后的時(shí)域圖和頻譜圖。 切比雪夫1型數(shù)字濾波器的頻率響應(yīng)圖 濾除白噪聲后的時(shí)域圖和頻域圖由于白噪聲在整個(gè)頻段上都存在，對(duì)與源信號(hào)共存的低頻信號(hào)用選頻濾波器是無(wú)法濾除的，故而采用低通濾波器將噪聲的高頻部分去掉。幅值(dB)39。)。xlabel(39。濾除白噪聲后的頻譜圖39。plot(f1,abs(s3))。f1=100*(0:999)/1000。)。ylabel(39。時(shí)間（s）39。)。title(39。subplot(2,1,1)。y=filter(num,den,s)。[num,den]=bilinear(BT,AT,)。)。[B,A]=cheby1(n1,1,WN,39。s39。WS=WP*(W0^2Ws1^2)/(Ws1*BW)。W00=sqrt(W0)。BW=Omegap2Omegap1。Omegas1=2*fs*tan(Ws1/2)。Omegap1=2*fs*tan(Wp1/2)。Wp2=(f12/fs)*2*pi。Wp1=(f11/fs)*2*pi。f21=5。f11=10。具體過(guò)程如下：figure(7)。 Butterworth型低通數(shù)字濾波器頻率響應(yīng)圖 濾除單頻正弦波后的時(shí)域圖和濾除噪聲后的頻譜圖由于信號(hào)處于頻段的低頻部分，而工頻信號(hào)的頻譜在整個(gè)上是呈對(duì)稱分布的，與源信號(hào)的頻段分布是相似的，采用低通濾波器將噪聲信號(hào)的高頻部分濾掉，濾波后時(shí)域圖的幅值比原來(lái)變小了，濾波后的頻譜圖頻率在0~10Hz和90~100Hz之間的幅值變化較大，而在10~90Hz的區(qū)間，頻譜圖基本趨于直線。幅值(dB)39。)。xlabel(39。濾除白噪聲后的頻譜圖39。plot(f1,abs(y3))。f1=100*(0:999)/1000。)。ylabel(39。時(shí)間（s）39。)。title(39。figure(6)。freqz(bz,az,512,fs)。)。[b,a]=butter(n1,Wn,39。s39。Omegas=2*fs*tan(ws/2)。ws=(f2/fs)*2*pi。f2=10。fs=100。（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)對(duì)工頻干擾的濾波，用plot函數(shù)畫出濾除工頻干擾后的時(shí)域圖和濾除白噪聲后的頻譜圖。)。ylabel(39。幅度39。)。title(39。subplot(2,1,2)。幅值（A）39。)。xlabel(39。加入帶限噪聲后的時(shí)域圖39。plot(k(1:1000))。figure(4)。[hc,wc]=freqz(k,1,1024)。（6） 帶限chirp噪聲 同樣用k表示加入帶限chirp噪聲后的信號(hào)，再用freqz進(jìn)行濾波器的頻率響應(yīng)，調(diào)用plot對(duì)k畫出加入帶限噪聲后的時(shí)域和頻域圖，具體描述如下：p=*chirp(a(:,1),0,a(1000,1),200)。)。ylabel(39。幅度39。)。title(39。subplot(2,1,2)。幅值（A）39。)。xlabel(39。加入工頻干擾后的時(shí)域圖39。plot(t,x)。figure(3)。y2=fft(x2,1000)。x1=a(:,2)。具體步驟如下：x2=[sin(2*pi*50*t)]。 加入白噪聲后的時(shí)域圖和頻譜圖通過(guò)觀察加入白噪聲后的時(shí)域圖和頻域圖，將它與未加入白噪聲進(jìn)行比較，可以發(fā)現(xiàn)頻譜圖在0Hz時(shí)的幅度增加的很大，而且又在沒(méi)有譜線的頻率上竟然出現(xiàn)了頻譜，這是由于白噪聲在所有頻率上都有頻率造成的。Hz39。)。xlabel(39。加入白噪聲后的頻譜圖39。plot(ws/pi,abshs)。)。ylabel(39。xlabel(‘時(shí)間（s）39。加入白噪聲后的