【文章內(nèi)容簡介】 后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——10——1（6） IIR帶通濾波程序：fp1=1200。fp2=3000。fs1=1000。fs2=3200。ft=10000。As=100。Rp=1。wp1=2*pi*fp1/ft。wp2=2*pi*fp2/ft。ws1=2*pi*fs1/ft。ws2=2*pi*fs2/ft。Omgp1=2*ft*tan(wp1/2)。Omgp2=2*ft*tan(wp2/2)。Omgp=[Omgp1,Omgp2]。Omgs1=2*ft*tan(ws1/2)。Omgs2=2*ft*tan(ws2/2)。Omgs=[Omgs1,Omgs2]。bw=Omgp2Omgp1。w0=sqrt(Omgp1*Omgp2)。%通帶寬和中心頻率[n,Omgn]=cheb1ord(Omgp,Omgs,Rp,As,39。s39。)。[z,p,k]=cheb1ap(n,Rp)。B=k*real(poly(z))。A=real(poly(p))。[b,a]=lp2bp(B,A,w0,bw)。[ba,aa]=bilinear(b,a,ft)。[Ha,w]=freqz(ba,aa)。[y,fs,nbits]=wavread(39。39。)。%加載語音信號f=filter(ba,aa,y)。wavplay(f)。f1=fft(f)。subplot(2,1,1)plot(f)。title(39。IIR帶通濾波后的波形39。)。subplot(2,1,2)plot(abs(f1))。title(39。IIR帶通濾波后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——11——15. 異性語音采集分析 運(yùn)用相關(guān)MATLAB函數(shù)和所學(xué)的相關(guān)知識，采集同樣內(nèi)容的異性的語音信號，對它的時(shí)頻域進(jìn)行分析處理，比較兩種相同內(nèi)容而不同性別語音信號的特點(diǎn)。步驟如下：（1） 異性語音錄制運(yùn)用MATLAB中的wavrecord()、wavwrite()、wavplay()等函數(shù)進(jìn)行語音的錄制、播放、存儲。把錄制好的語音文件存為“C:\Users\PKL\Desktop\”，具體程序如下： fs=16000。%設(shè)計(jì)抽樣頻率 channel=1。%設(shè)計(jì)錄音通道 t=5。 fprintf(39。按任意鍵開始%d秒錄音:39。,t)。 pause。%暫停命令 fprintf(39。正在錄音39。)。 y=wavrecord(t*fs,fs,channel,39。double39。)。%錄制語音 fprintf(39。錄音結(jié)束\n39。)。 wavwrite(y,fs,39。C:\Users\PKL\Desktop\39。)。 fprintf(39。按任意鍵回放語音39。)。 pause。 wavplay(y,fs)。（2） 異性語音的采樣分析 運(yùn)用MATLAB中的wavread()對錄制好的異性語音信號進(jìn)行加載采樣，然后用MATLAB的相關(guān)函數(shù)對采集的語音信號進(jìn)行時(shí)頻域畫圖分析，作好的圖為圖——12。具體程序如下： [y,fs,nbits]=wavread(39。39。)。%加載語音信號 wavplay(y,fs)。%回放語音信號 N=length(y)。%語音信號長度 Y=fft(y,N)。%用傅里葉變換求頻域 subplot(2,1,1)。 plot(y)。 title(39。異性語音原始信號時(shí)域波形圖39。)。 subplot(2,1,2)。 plot(abs(Y))。 title(39。異性語音原始信號頻譜圖39。)。6. 語音信號加噪、濾噪處理分析運(yùn)用MATLAB相關(guān)函數(shù)錄制一段相同長度的背景噪聲，再把此噪聲疊加到原始語音信號中，然后再選擇合適的濾波器對加噪語音信號進(jìn)行濾波。步驟如下：（1） 噪聲錄制與疊加運(yùn)用MATLAB隨即錄制一段與語音信號相同長度的噪聲，再疊加到原始語音信號上，把加噪后的語音信號存為“C:\Users\PKL\Desktop\”。實(shí)現(xiàn)程序如下：[y1,fs,nbits]=wavread(39。39。)。Y1=fft(y1)。n=length(y1)。noise=*randn(n,1)。%隨機(jī)加同樣長度的噪聲y2=y1+noise。Y2=fft(y2)。sound(y2)wavwrite(y2,fs,39。C:\Users\PKL\Desktop\39。)。subplot(2,2,1)。plot(y1)。title(39。疊加前語音信號波形39。)。subplot(2,2,3)。plot(abs(Y1))。title(39。疊加前語音信號頻譜圖39。)。subplot(2,2,2)。plot(y2)。title(39。疊加后語音信號波形39。)。subplot(2,2,4)。plot(abs(Y2))。title(39。疊加后語音信號頻譜圖39。)。程序所得圖形為附件：圖——13（2） 濾除噪聲通過（1）中的加噪前后時(shí)頻域的分析，濾除此噪聲應(yīng)用低通濾波器進(jìn)行濾波，實(shí)現(xiàn)程序如下：ft=5500。fpts=[1000 1200]。dev=[ ]。%通帶與阻帶波動mag=[1 0]。[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,ft)。b21=fir1(n21,wn21,Kaiser(n21+1,beta))。[h,w]=freqz(b21,1)。[y,fs,nbits]=wavread(39。39。)。%加載語音信號Y=fft(y)。f=filter(b21,1,y)。F=fft(f)。subplot(2,2,1)plot(y)。title(39。低通噪聲濾波前的波形39。)。subplot(2,2,3)plot(abs(Y))。title(39。低通噪聲濾波前的頻譜圖39。)subplot(2,2,2)plot(f)。title(39。低通噪聲濾波后的波形39。)。subplot(2,2,4)plot(abs(F))。title(39。低通噪聲濾波后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——147. 信號處理系統(tǒng)界面設(shè)計(jì)建立GUI文件：FileNewGUI:如圖——15，后選擇Blan GUI(Default),得到如圖——16的GUI開發(fā)環(huán)境。布置好各控件，為控件編寫程序來實(shí)現(xiàn)兩數(shù)相加的功能?？丶O(shè)計(jì)如圖——17,各按鈕設(shè)計(jì)程序如下：GUI設(shè)計(jì)程序：function varargout = untitled1(varargin)gui_Singleton = 1。 gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @untitled1_OpeningF, ... 39。gui_OutputF39。, @untitled1_OutputF, ... 39。gui_LayoutF39。, [] , ...