【正文】 As,39。)。B=k*real(poly(z))。[b,a]=lp2bp(B,A,w0,bw)。[Ha,w]=freqz(ba,aa)。title(39。)。 (2） 用窗函數(shù)法設(shè)計濾波器 窗函數(shù)設(shè)計方法主要有矩形窗、巴特列特窗、漢寧窗、海明窗、布拉克曼窗凱澤窗等。FIR低通濾波器設(shè)計程序：ft=5500。dev=[ ]。[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,ft)。[h,w]=freqz(b21,1)。title(39。)。fpts=[4800 5000]。mag=[1 0]。b22=fir1(n22,wn22,39。,Kaiser(n22+1,beta))。plot(w/pi,abs(h))。圖——7 FI高通濾波器39。頻率響應(yīng)圖為：圖——7FIR帶通濾波器的設(shè)計程序：fp1=1200。fs1=1000。ft=10000。ws2=fs2/(ft/2)*pi。wp2=fp2/(ft/2)*pi。ws=[ws1,ws2]。N0=ceil(11*pi/deltaw)。windows=(blackman(N))39。wc2=(ws2+wp2)/2。b=hd.*windows。plot(w/pi,abs(h))。圖——8 FIR帶通濾波器39。 頻率響應(yīng)圖為：圖——84. 對“”語音信號進(jìn)行濾波 分別采用設(shè)計好的低通，高通，帶通濾波器對“”語音信號進(jìn)行濾波，用MATLAB畫出濾波后的時頻域圖，并回放錄播后的語音信號。fpts=[1000 1200]。mag=[1 0]。b21=fir1(n21,wn21,Kaiser(n21+1,beta))。[y,fs,nbits]=wavread(‘39。%加載語音信號f=filter(b21,1,y)。f1=fft(f)。title(39。)。title(39。)運行結(jié)果圖為附件：圖——9（2） FIR高通濾波程序：ft=25000。dev=[ ]。[n22,wn22,beta,ftype]=kaiserord(fpts,mag,dev,ft)。high39。[h,w]=freqz(b22,1)。39。%加載語音信號f=filter(b22,1,y)。f1=fft(f)。title(39。)。title(39。)運行結(jié)果圖為附件：圖——10（3） FIR帶通濾波程序：fp1=1200。fs1=1000。ft=10000。ws2=fs2/(ft/2)*pi。wp2=fp2/(ft/2)*pi。ws=[ws1,ws2]。N0=ceil(11*pi/deltaw)。windows=(blackman(N))39。wc2=(ws2+wp2)/2。b=hd.*windows。 [y,fs,nbits]=wavread(39。)。 wavplay(f)。 subplot(2,1,1) plot(f)。帶通濾波后的波形39。 subplot(2,1,2) plot(abs(f1))。帶通濾波后的頻譜圖39。 ft=5500。 wp=2*pi*fp/ft。%截止頻率 Omgp=2*ft*tan(wp/2)。 [n,Omgc]=buttord(Omgp,Omgs,1,100,39。)。 B=k*real(poly(z))。 [b,a]=lp2lp(B,A,Omgc)。[Ha,w]=freqz(ba,aa)。39。%加載語音信號f=filter(ba,aa,y)。f1=fft(f)。title(39。)。title(39。)運行結(jié)果圖為附件：圖——9——1（5） IIR高通濾波程序： fp=5000。%取抽樣頻率 fs=4800。As=100。%通帶頻率 ws=2*pi*fs/ft。 Omgs=2*ft*tan(ws/2)。s39。 %計算階數(shù)n和截止頻率 [z,p,k]=ellipap(n,Rp,As)。 A=real(poly(p))。 [ba,aa]=bilinear(b,a,ft)。[y,fs,nbits]=wavread(39。)。wavplay(f)。subplot(2,1,1)plot(f)。IIR高通濾波后的波形39。subplot(2,1,2)plot(abs(f1))。IIR高通濾波后的頻譜圖39。fp2=3000。fs2=3200。As=100。wp1=2*pi*fp1/ft。ws1=2*pi*fs1/ft。Omgp1=2*ft*tan(wp1/2)。Omgp=[Omgp1,Omgp2]。Omgs2=2*ft*tan(ws2/2)。bw=Omgp2Omgp1。%通帶寬和中心頻率[n,Omgn]=cheb1ord(Omgp,Omgs,Rp,As,39。)。B=k*real(poly(z))。[b,a]=lp2bp(B,A,w0,bw)。[Ha,w]=freqz(ba,aa)。39。%加載語音信號f=filter(ba,aa,y)。f1=fft(f)。title(39。)。title(39。)運行結(jié)果圖為附件：圖——11——15. 異性語音采集分析 運用相關(guān)MATLAB函數(shù)和所學(xué)的相關(guān)知識，采集同樣內(nèi)容的異性的語音信號，對它的時頻域進(jìn)行分析處理，比較兩種相同內(nèi)容而不同性別語音信號的特點。把錄制好的語音文件存為“C:\Users\PKL\Desktop\”，具體程序如下： fs=16000。%設(shè)計錄音通道 t=5。按任意鍵開始%d秒錄音:39。 pause。正在錄音39。 y=wavrecord(t*fs,fs,channel,39。)。錄音結(jié)束\n39。 wavwrite(y,fs,39。)。按任意鍵回放語音39。 pause。（2） 異性語音的采樣分析 運用MATLAB中的wavread()對錄制好的異性語音信號進(jìn)行加載采樣，然后用MATLAB的相關(guān)函數(shù)對采集的語音信號進(jìn)行時頻域畫圖分析，作好的圖為圖——12。39。%加載語音信號 wavplay(y,fs)。%語音信號長度 Y=fft(y,N)。 plot(y)。異性語音原始信號時域波形圖39。 subplot(2,1,2)。 title(39。)。步驟如下：（1） 噪聲錄制與疊加運用MATLAB隨即錄制一段與語音信號相同長度的噪聲，再疊加到原始語音信號上，把加噪后的語音信號存為“C:\Users\PKL\Desktop\”。39。Y1=fft(y1)。noise=*randn(n,1)。Y2=fft(y2)。C:\Users\PKL\Desktop\39。subplot(2,2,1)。title(39。)。plot(abs(Y1))。疊加前語音信號頻譜圖39。subplot(2,2,2)。title(39。)。plot(abs(Y2))。疊加后語音信號頻譜圖39。程序所得圖形為附件：圖——13（2） 濾除噪聲通過（1）中的加噪前后時頻域的分析，濾除此噪聲應(yīng)用低通濾波器進(jìn)行濾波，實現(xiàn)程序如下：ft=5500。dev=[ ]。[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,ft)。[h,w]=freqz(b21,1)。39。%加載語音