最新基于matlab的語音信號處理系統(tǒng)設計程序仿真圖畢業(yè)設計(編輯修改稿)
2025-07-16 13:33 本頁面
【文章內容簡介】 %做加噪音信號的頻譜圖plot(w/pi,abs(y2))。title(39。加噪音信號頻譜39。)xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。程序結果：function [db,mag,pha,w]=freqz_m(b,a)。[H,w]=freqz(b,a,1000,39。whole39。)。 %在02*pi之間選取N個點計算頻率響應H=(H(1:501))39。 %頻率響應 w=(w(1:501))39。 %頻率mag=abs(H)。 %響應幅度db=20*log10((mag+eps)/max(mag))。 %增益pha=angle(H)。function [b,a]=afd_butt(Wp,Ws,Rp,As)N=ceil((log10((10^(Rp/10)1)/(10^(As/10)1)))/(2*log10(Wp/Ws)))。 %濾波器階數(shù)N為整數(shù)fprintf(39。\n Butterworth Filter Order=%\n39。,N)OmegaC=Wp/((10^(Rp/10)1)^(1/(2*N))) %求對應于N的3db截止頻率[b,a]=u_buttap(N,OmegaC)。function [b,a]=imp_invr(c,d,T)[R,p,k]=residue(c,d)。 %部分分式展開p=exp(p*T)。 %從模擬到數(shù)字極點對應關系,部分分式系數(shù)相同[b,a]=residuez(R,p,k)。 %將部分分式的形式變換成多項式之比的形式b=real(b39。)。 %求出數(shù)字濾波器系數(shù)a=real(a39。)。程序：wp=*pi。ws=*pi。Rp=1。As=30。T=1。 %定義濾波器參數(shù) OmegaP=wp/T。OmegaS=ws/T。[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As)。[b,a]=imp_invr(cs,ds,T)。[db,mag,pha,w]=freqz_m(b,a)。figure(1)subplot(2,1,1)。plot(w/pi,mag)。 %作出濾波器的圖形title(39。digital filter Magnitude Response39。)axis([0,1,0,])subplot(2,1,2)。plot(w/pi,db)。title(39。digital filter Magnitude in DB39。)axis([0,1,40,5])。x=wavread(39。39。)39。 %打開語音信號 N=length(x)。 %長度n=0:N1。w=2*n*pi/N。z=*cos(10*n)。 %噪音y=x+z。 y1=filter(b,a,y)。 %用濾波器濾波y2=fft(y1)。 %對濾波后信號做FFT變換sound(y1)。 %讀濾波后信號figure(2)subplot(2,1,1)。 plot(n,y1) %做濾波后信號的時域波形圖title(39。濾波后信號時域圖39。)。xlabel(39。時間t39。)。ylabel(39。幅值39。)。subplot(2,1,2)。 %做濾波后信號的頻譜圖plot(w/pi,abs(y2))。title(39。濾波后信號頻譜39。)xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。 GUI的函數(shù)1 開語音global x。global fs。 global nbits。 global w。global N。 H={39。*.wav39。}。 [x,fs,nbits]=wavread(39。39。)。 y=fft(x)。N=length(x)。 n=0:N1。w=2*n*pi/N。subplot(2,1,1)。 plot(n,x) title(39。原始語音信號時域圖39。)。xlabel(39。時間t39。)。ylabel(39。幅值39。)。subplot(2,1,2)。 plot(w/pi,abs(y))。title(39。原始語音信號頻譜39。)xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。點開語音鍵結果：
數(shù)學相關推薦
鄂ICP備17016276號-1