基于matlab有噪聲語音信號的處理畢業(yè)論文-在線瀏覽
2024-08-04 15:46本頁面
【正文】 Fs=8000。wp=2*pi*fp/Fs。ws=2*pi*fs/Fs。was=2*Fs*tan(ws/2)。As=100。s39。[z,p,k]=buttap(M)。[bs,as]=lp2lp(bp,ap,wap)。[h,w]=freqz(bz,az)。sound(f1,45000)。%濾波后聲音，因為是低通濾波器，所以原聲音信號要比以前聲音低沉%圖示：%1：時域波形figure(1)。plot(y)。原信號時域波形39。subplot(3,1,2)。title(39。39。subplot(3,1,3)。title(39。)%2頻域波形figure(2)。plot(abs(Y))。原信號頻域波形39。subplot(3,1,2)。title(39。39。subplot(3,1,3)。title(39。)%濾波器波形figure(3)。axis([0,600,0,])。一般來說，濾波器的性能要求往往以頻率響應(yīng)的幅度特性的允許誤差來表征。如圖11所示：在通帶內(nèi)：通帶內(nèi)具有等波動特性得稱為切比雪夫型，阻帶內(nèi)具有等波動振幅特性得稱為切比雪夫型。切比雪夫I型的振幅平方函數(shù)為189。=1/[1+e2C2N(w/wc)]這里，ε是紋波系數(shù)，用來描述波動得大小。一般情況下，波動范圍在1/(1+ε2)～1之間（）。x189。1. CN(x)=cosh[Narccosh(x)], 189。1.首先根據(jù)通帶截頻確定wc wc=wp 。[arccosh(1/e)()]/arccosh(ws/wp)。H（jw）189。H（jw）189。[y,fs,nbits]=wavreadc:\39。Nlength%求出語音信號的長度%sound(y,fs,nbits)。%Noise=ones(N,2)。%載波頻率Noise=modulate(y,fc,fs,39。)。%隨機函數(shù)產(chǎn)生噪聲si=y+Noise。%sound(si,fs,nbits)。Ft=8000。Fs=3500。wp=1。[n13,wn13]=cheb1ord(wp,ws,1,50,39。)。%求模擬的低通濾波器階數(shù)和截止頻率[b13,a13]=cheby1(n13,1,wn13,39。)。%plot(w*21000*,abs(h))。IIR高通濾波器39。legend(39。)。%圖示：%1：時域波形figure(1)。plot(y)。原信號時域波形39。subplot(3,1,2)。title(39。39。subplot(3,1,3)。title(39。)subplot(3,1,1)。title(39。)39。plot(abs(SI))。加躁信號頻域波形)。plot(abs(F1))。濾波信號頻域波形39。%濾波器波形figure(3)。title(39。)。sound(si,fs)。，頻域，波形圖加噪前后時域波形比較加噪前后頻域波形比較切比雪夫I型濾波器切比雪夫濾波器的特點誤差值在規(guī)定的頻段上等幅變化。切比雪夫濾波器的 在通帶范圍內(nèi)是等幅起伏的，所以同樣的通帶衰減，其階數(shù)較巴特沃茲濾波器要小。它的振幅平方函數(shù)為|H(jw)|2=1/[1+ε2R2N(w/wc)]由通帶截頻wp確定wc wc=wp由通帶衰減Ap確定e e=（）由阻帶截頻wc確定k k=wp/ws由阻帶衰減As確定k1 k1=e/()確定濾波器階數(shù)N N=K(k)K(())/K((1k2))K(k1)K(k)=(1k2sin2x) 由得到的參數(shù)N，e，wc,k和k1 確定系統(tǒng)函數(shù)。圖1 N=5時雅可比橢圓函數(shù)的特性曲線clear。[y,fs,nbits]=wavreadc:\39。Nlength%求出語音信號的長度%sound(y,fs,nbits)。%Noise=3*modulate(y,20000,fs,39。)。%隨機函數(shù)產(chǎn)生噪聲si=y+2*Noise。%sound(si,fs,nbits)。%傅里葉變換%Rp=。%Wp=[10%Ws=[8%[n,Wn]=ellipord(Wp,Ws,Rp,Rs,39。)。high39。%f1=filter(b,a,si)。%wp1=*pi。wp2=*pi。%ws1=*pi。%Rp=3。%Wp=[wp1/pi,wp2/pi]。%[N,Wn]=ellipord(Wp,Ws,Rp,Rs)。
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1