基于matlab的語音信號的時、頻域分析課程設計(編輯修改稿)
2025-02-14 14:16 本頁面
【文章內容簡介】 點DFT，忽略y后面的元素。首先畫出語音信號的時域波形，然后對語音信號進行頻譜分析。在MATLAB中利用fft對新號進行快速傅立葉變換，得到信號的頻譜特性。subplot(2,1,1)。plot(y)。title(‘原始信號波形’)。subplot(2,1,2)。plot(abs(Y))。title(‘原始信號頻譜’)；程序如下：[y,fs,nbits]=wavread(39。39。)。sound(y,fs,nbits)。N=length(y)。Y=fft(y,N)。subplot(2,1,1)。plot(y)。title(39。原始信號波形39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。subplot(2,1,2)。plot(abs(Y))。title(39。原始信號頻譜39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。程序結果如下圖：3構造受干擾信號并對其進行FFT頻譜分析程序如下：[y,fs,nbits]=wavread (39。39。)。N = length (y) 。 %求出語音信號的長度noise=*randn(N,1)。 %隨機函數產生噪聲Si=y+noise。 %語音信號加入噪聲 sound(Si)。subplot(2,1,1)。plot(Si)。title(39。加噪語音信號的時域波形39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。S=fft(Si)。 %傅里葉變換subplot(2,1,2)。plot(abs(S))。 title(39。加噪語音信號的頻域波形39。)。 xlabel(39。時間39。)。ylabel(39。幅度39。)。程序結果如下圖：4 數字濾波器設計 用窗函數法設計IIR帶通濾波器程序如下：wp=*pi。 ws=*pi。 %性能指標wdelta=wswp。 %過渡帶寬度M=ceil(*pi/wdelta)。 %濾波器長度，朝正無窮方向舍入N=2*M+1。 %窗口長度wc=(ws+wp)/2。 %截止頻率win=hamming(N)。 %因為衰減不小于50dB，所以選擇海明窗，這里得到海明窗的時域響應b=fir1(N1,w
環(huán)評公示相關推薦
鄂ICP備17016276號-1