基于matlab的語音信號的采集和處理-資料下載頁
2025-06-18 16:08本頁面
【正文】 夫高通濾波器系數(shù) yh=filter(BZ,AZ,y)。subplot(211)。plot(yh)。xlabel(39。time n39。)。ylabel(39。amplitude39。)。title(39。高通濾波后語音信號時域波形39。)。ah=fftshift(fft(yh,l)/l)。fdh=fs/l。fxh=fdh*(l/2:l/21)。fyh=abs(ah)。 subplot(212)。plot(fxh,fyh)。xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。title(39。高通濾波后信號頻譜39。)。fph=[1000,2000]。 %帶通濾波程序開始 fsh=[800,2500]。 rp=3。 rs=20。 omegaph=fph/(fs/2)。 omegash=fsh/(fs/2)。 [N,wn]=cheb1ord(omegaph,omegash,rp,rs)。 [BZ,AZ]=cheby1(N,rp,wn,39。bandpass39。)。 yb=filter(BZ,AZ,y)。 subplot(211)。plot(yb)。xlabel(39。time n39。)。ylabel(39。amplitude39。)。title(39。帶通濾波后語音信號時域波形39。)。ab=fftshift(fft(yb,l)/l)。 fdb=fs/l。fxb=fdb*(l/2:l/21)。fyb=abs(ab)。subplot(212)。plot(fxb,fyb)。xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。title(39。帶通濾波后信號頻譜39。)。wavplay(y,fs)。 %播放原始語音信號wavplay(yl,fs)。 %播放低通濾波后語音信號wavplay(yh,fs)。 %播放高通濾波后語音信號wavplay(yb,fs)。 %播放帶通濾波后語音信號Au=1。t=0:1/fs:(l1)/fs。d=Au*sin(2*pi*3800*t)。 %構(gòu)造頻率為3800hz的正弦信號，且取其一段d=d39。yy=y+d。 %將所得噪音信號加入原始信號中subplot(211)。plot(yy)。 %繪制加入噪音后波形xlabel(39。time n39。)。ylabel(39。amplitude39。)。title(39。加入噪音后語音信號時域波形39。)。ay=fftshift(fft(yy,l)/l)。fdy=fs/l。fxy=fdy*(l/2:l/21)。fyy=abs(ay)。subplot(211)。plot(yy)。 %繪制加入噪音后頻譜xlabel(39。time n39。)。ylabel(39。amplitude39。)。title(39。加入噪音后語音信號時域波形39。)。fph=3400。 %對加入噪音后信號進(jìn)行低通濾波，濾除噪音 fsh=3600。 rp=2。 rs=60。 omegaph=fph/(fs/2)。 omegash=fsh/(fs/2)。 [N,wn]=cheb1ord(omegaph,omegash,rp,rs)。 [BZ,AZ]=cheby1(N,rp,wn,39。low39。)。 yyy=filter(BZ,AZ,yy)。subplot(211)。plot(yyy)。xlabel(39。time n39。)。ylabel(39。amplitude39。)。title(39。對加噪信號濾波后語音信號時域波形39。)。ayy=fftshift(fft(yyy,l)/l)。fdyy=fs/l。fxyy=fdyy*(l/2:l/21)。fyyy=abs(ayy)。subplot(212)。plot(fxyy,fyyy)。xlabel(39。頻率Hz39。)。ylabel(39。幅度39。)。title(39。加入噪音濾波后信號頻譜39。)。
試題試卷相關(guān)推薦
鄂ICP備17016276號-1