語(yǔ)音數(shù)字信號(hào)處理與分析及matlab實(shí)現(xiàn)_數(shù)字信號(hào)課程設(shè)計(jì)-資料下載頁(yè)
【導(dǎo)讀】用,為下一步學(xué)習(xí)和工作打下良好的基礎(chǔ)。語(yǔ)音信號(hào)處理的發(fā)展依賴于這些學(xué)科的發(fā)展，而語(yǔ)。音信號(hào)處理技術(shù)的進(jìn)步也會(huì)促進(jìn)這些領(lǐng)域的進(jìn)步?？梢詫?duì)數(shù)字化的語(yǔ)音信號(hào)進(jìn)行處理。由于MATLAB是一種面向科學(xué)。的庫(kù)函數(shù)，所以編程簡(jiǎn)單、編程效率高、易學(xué)易懂。信號(hào)進(jìn)行加噪和去噪、濾波、截取語(yǔ)音等。設(shè)計(jì)合適的濾波器濾除噪聲，恢復(fù)原信號(hào)。法，各種不同類型濾波器的性能比較），闡明設(shè)計(jì)原理。采集語(yǔ)音信號(hào)，把語(yǔ)音信號(hào)保存為.wav文件，長(zhǎng)度小于3秒，錄，并對(duì)其進(jìn)行FFT頻譜分析，畫出信號(hào)的時(shí)域波形圖和頻譜圖。行分析，闡明原因，得出結(jié)論。波器的幅頻和相頻特性。設(shè)計(jì)處理系統(tǒng)的用戶界面，在所設(shè)。N=128點(diǎn)DFT的幅頻圖和N=1024點(diǎn)DFT幅頻圖。由于采樣頻率f=100Hz，故Nyquist頻率為50Hz，、或可見(jiàn)，整個(gè)頻譜圖是以Nyquist頻率為軸對(duì)稱的。半）范圍的幅頻特性。小與fft選用點(diǎn)數(shù)N有關(guān)，但只要點(diǎn)數(shù)N足夠不影響研究結(jié)果。從上圖幅頻譜可見(jiàn)，信號(hào)中包括15Hz和40Hz的正弦分量。
【正文】 :N/2),magy(1:N/2))。 xlabel(39。頻率 (Hz)39。)。ylabel(39。幅值 39。)。 II title(39。N=128(b)39。)。grid fs=100。N=1024。 y=fft(x,N)。 magy=abs(y)。 f=(0:length(y)1)39。*fs/length(y)。 subplot(223)。plot(f,magy)。 xlabel(39。頻率 (Hz)39。)。ylabel(39。幅值 39。)。 title(39。N=1024(c)39。)。grid subplot(224)。plot(f(1:N/2),magy(1:N/2))。 xlabel(39。頻率 (Hz)39。)。ylabel(39。幅值 39。)。 title(39。N=1024(d)39。)。grid [x,fs,bits]=wavread(39。39。)。 N=256。 f=0:fs/N:fs/21/N。 X=fft(x)。 X=abs(X)。 subplot(211) plot(f(45:60),X(45:60))。grid xlabel(39。Hz39。),ylabel(39。|H(ejw)|39。) N=N*4。 f=0:fs/N:fs/21/N。 X=fft(x)。 X=abs(X)。 III subplot(212) plot(f(45*4:4*60),X(4*45:4*60))。grid xlabel(39。Hz39。),ylabel(39。|H(ejw)|39。) x=wavread(39。39。)。 sound(x)。 N=5。wc=。 [b,a]=butter(N,wc,39。high39。)。 X=fft(x)。 subplot(321)。plot(x)。title(39。濾波前信號(hào)的波形 39。)。 subplot(322)。plot(X)。title(39。濾波前信號(hào)的頻譜 39。)。 y=filter(b,a,x)。 Y=fft(y)。 subplot(323)。plot(y)。title(39。IIR 濾波后信號(hào)的波形 39。)。 subplot(324)。plot(Y)。title(39。IIR 濾波后信號(hào)的頻譜 39。)。 z=fftfilt(b,x)。 Z=fft(z)。 subplot(325)。plot(z)。title(39。FIR 濾波后信號(hào)的波形 39。)。 subplot(326)。plot(Z)。title(39。FIR 濾波后信號(hào)的頻譜 39。)。 x=wavread(39。39。)。 sound(x)。 N=5。wc=[,]。 [b,a]=butter(N,wc)。 X=fft(x)。 IV subplot(321)。plot(x)。title(39。濾波前信號(hào)的波形 39。)。 subplot(322)。plot(X)。title(39。濾波前信號(hào)的頻譜 39。)。 y=filter(b,a,x)。 Y=fft(y)。 subplot(323)。plot(y)。title(39。IIR 濾波后信號(hào)的波形 39。)。 subplot(324)。plot(Y)。title(39。IIR 濾波后信號(hào)的頻譜 39。)。 z=fftfilt(b,x)。 Z=fft(z)。 subplot(325)。plot(z)。title(39。FIR 濾波后信號(hào)的波形 39。)。 subplot(326)。plot(Z)。title(39。FIR 濾波后信號(hào)的頻譜 39。)。 x=wavread(39。39。)。 sound(x)。 N=5。wc=[,]。 [b,a]=butter(N,wc,39。stop39。)。 X=fft(x)。 subplot(321)。plot(x)。title(39。濾波前信號(hào)的波形 39。)。 subplot(322)。plot(X)。title(39。濾波前信號(hào)的頻譜 39。)。 y=filter(b,a,x)。 Y=fft(y)。 subplot(323)。plot(y)。title(39。IIR 濾波后信號(hào)的波形 39。)。 subplot(324)。plot(Y)。title(39。IIR 濾波后信號(hào)的頻譜 39。)。 z=fftfilt(b,x)。 V Z=fft(z)。 subplot(325)。plot(z)。title(39。FIR 濾波后信號(hào)的波形 39。)。 subplot(326)。plot(Z)。title(39。FIR 濾波后信號(hào)的頻譜 39。)。 x=wavread(39。39。)。 sound(x)。 b=100。a=5。 y=filter(b,a,x)。 X=fft(x,4096)。 subplot(221)。plot(x)。title(39。濾波前信號(hào)的波形 39。)。 subplot(222)。plot(abs(X))。title(39。濾波前信號(hào)的幅值 39。)。 Y=fft(y， 4096)。 subplot(223)。plot(y)。title(39。濾波后信號(hào)的波形 39。)。 subplot(224)。plot(abs(Y))。title(39。濾波后信號(hào)的幅值 39。)。 sound(y)。 figure, subplot(211)。plot(angle(X))。title(39。濾波前信號(hào)相位 39。)。 subplot(212)。plot(angle(Y))。title(39。濾波后信號(hào)相位 39。)。
鄂ICP備17016276號(hào)-1