參與評(píng)論

評(píng)論

登錄后參與評(píng)論

全部評(píng)論(2條)

• 

  汪樂(lè)安 2009-06-25 00:00:00

  太多了，其實(shí)網(wǎng)上這方面現(xiàn)成的程序很多，大概如下： 1：使用wavrecord函數(shù) 2：使用awgn等函數(shù) 3：使用fdatool可以設(shè)計(jì)濾波器 5：wavplay 6: guide

  贊(13)

  回復(fù)(0)

  評(píng)論

  評(píng)論

  登錄后參與評(píng)論

• 

  love_xue_58 2017-12-16 13:42:01

  你自己整合吧，我沒(méi)時(shí)間幫你整合，我給你提供一些程序： 正確的代碼：程序1： fs=22050; %語(yǔ)音信號(hào)采樣頻率為22050 x1=wavread('Windows Critical Stop.wav'); %讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1 sound(x1，22050); %播放語(yǔ)音信號(hào) y1=fft(x1，1024); %對(duì)信號(hào)做1024點(diǎn)FFT變換 f=fs*(0:511)/1024; figure(1) plot(x1) %做原始語(yǔ)音信號(hào)的時(shí)域圖形 title('原始語(yǔ)音信號(hào)'); xlabel('time n'); ylabel('fuzhi n'); figure(2) freqz(x1) %繪制原始語(yǔ)音信號(hào)的頻率響應(yīng)圖 title('頻率響應(yīng)圖') figure(3) subplot(2，1，1); plot(abs(y1(1:512))) %做原始語(yǔ)音信號(hào)的FFT頻譜圖 title('原始語(yǔ)音信號(hào)FFT頻譜') subplot(2，1，2); plot(f，abs(y1(1:512))); title('原始語(yǔ)音信號(hào)頻譜') xlabel('Hz'); ylabel('fuzhi'); 程序2： fs=22050; %語(yǔ)音信號(hào)采樣頻率為22050 x1=wavread('Windows Critical Stop.wav'); %讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1 t=0:1/22050:(size(x1)-1)/22050; y1=fft(x1，1024); %對(duì)信號(hào)做1024點(diǎn)FFT變換 f=fs*(0:511)/1024; x2=randn(1，length(x1)); %產(chǎn)生一與x長(zhǎng)度一致的隨機(jī)信號(hào) sound(x2，22050); figure(1) plot(x2) %做原始語(yǔ)音信號(hào)的時(shí)域圖形 title('高斯隨機(jī)噪聲'); xlabel('time n'); ylabel('fuzhi n'); randn('state'，0); m=randn(size(x1)); x2=0.1*m+x1; sound(x2，22050);%播放加噪聲后的語(yǔ)音信號(hào) y2=fft(x2，1024); figure(2) plot(t，x2) title('加噪后的語(yǔ)音信號(hào)'); xlabel('time n'); ylabel('fuzhi n'); figure(3) subplot(2，1，1); plot(f，abs(y2(1:512))); title('原始語(yǔ)音信號(hào)頻譜'); xlabel('Hz'); ylabel('fuzhi'); subplot(2，1，2); plot(f，abs(y2(1:512))); title('加噪后的語(yǔ)音信號(hào)頻譜'); xlabel('Hz'); ylabel('fuzhi'); 根據(jù)以上代碼，你可以修改下面有錯(cuò)誤的代碼 程序3:雙線性變換法設(shè)計(jì)Butterworth濾波器 fs=22050; x1=wavread('h:\課程設(shè)計(jì)2\shuzi.wav'); t=0:1/22050:(size(x1)-1)/22050; Au=0.03; d=[Au*cos(2*pi*5000*t)]'; x2=x1+d; wp=0.25*pi; ws=0.3*pi; Rp=1; Rs=15; Fs=22050; Ts=1/Fs; wp1=2/Ts*tan(wp/2); %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo) ws1=2/Ts*tan(ws/2); [N，Wn]=buttord(wp1，ws1，Rp，Rs，'s'); %選擇濾波器的Z小階數(shù) [Z，P，K]=buttap(N); %創(chuàng)建butterworth模擬濾波器 [Bap，Aap]=zp2tf(Z，P，K); [b，a]=lp2lp(Bap，Aap，Wn); [bz，az]=bilinear(b，a，F(xiàn)s); %用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換 [H，W]=freqz(bz，az); %繪制頻率響應(yīng)曲線 figure(1) plot(W*Fs/(2*pi)，abs(H)) grid xlabel('頻率／Hz') ylabel('頻率響應(yīng)幅度') title('Butterworth') f1=filter(bz，az，x2); figure(2) subplot(2，1，1) plot(t，x2) %畫(huà)出濾波前的時(shí)域圖 title('濾波前的時(shí)域波形'); subplot(2，1，2) plot(t，f1); %畫(huà)出濾波后的時(shí)域圖 title('濾波后的時(shí)域波形'); sound(f1，22050); %播放濾波后的信號(hào) F0=fft(f1，1024); f=fs*(0:511)/1024; figure(3) y2=fft(x2，1024); subplot(2，1，1); plot(f，abs(y2(1:512))); %畫(huà)出濾波前的頻譜圖 title('濾波前的頻譜') xlabel('Hz'); ylabel('fuzhi'); subplot(2，1，2) F1=plot(f，abs(F0(1:512))); %畫(huà)出濾波后的頻譜圖 title('濾波后的頻譜') xlabel('Hz'); ylabel('fuzhi'); 程序4:窗函數(shù)法設(shè)計(jì)濾波器： fs=22050; x1=wavread('h:\課程設(shè)計(jì)2\shuzi.wav'); t=0:1/22050:(size(x1)-1)/22050; Au=0.03; d=[Au*cos(2*pi*5000*t)]'; x2=x1+d; wp=0.25*pi; ws=0.3*pi; wdelta=ws-wp; N=ceil(6.6*pi/wdelta); %取整 wn=(0.2+0.3)*pi/2; b=fir1(N，wn/pi，hamming(N+1)); %選擇窗函數(shù)，并歸一化截止頻率 figure(1) freqz(b，1，512) f2=filter(bz，az，x2) figure(2) subplot(2，1，1) plot(t，x2) title('濾波前的時(shí)域波形'); subplot(2，1，2) plot(t，f2); title('濾波后的時(shí)域波形'); sound(f2，22050); %播放濾波后的語(yǔ)音信號(hào) F0=fft(f2，1024); f=fs*(0:511)/1024; figure(3) y2=fft(x2，1024); subplot(2，1，1); plot(f，abs(y2(1:512))); title('濾波前的頻譜') xlabel('Hz'); ylabel('fuzhi'); subplot(2，1，2) F2=plot(f，abs(F0(1:512))); title('濾波后的頻譜') xlabel('Hz'); ylabel('fuzhi');

  贊(19)

  回復(fù)(0)

  評(píng)論

  評(píng)論

  登錄后參與評(píng)論