【正文】 號 y=fft(y1)。 figure(1) subplot(2,1,1)。原始語音信號 FFT頻譜 39。%做原始語音信號的時域圖形 plot(y1) title(39。)。時間 n39。 ylabel(39。)。 Randn 函數(shù)有兩種基本調用格式： Randn(n)和 Randn（ m,n） ,前者產生 n n服從標準高斯分布的隨機數(shù)矩陣，后者產生 m n 的隨機數(shù)矩陣。 加噪程序如下： [y,fs,nbits]=wavread (39。)。 %求 出語音信號 的長度 Noise=*randn(N,2)。 %語音信號加入噪聲 sound(Si)。 plot(Si)。加噪語音信號的時域波形 39。 S=fft(Si)。 plot(abs(S))。加噪語音信號的頻域波形 39。 程序結果如下圖： 9 圖 6 加噪后的波形及頻譜分析 比較濾波前后語音信號的波形及頻譜 % ======================雙線性變換法 ======================= %*************************低通濾波器 ************************ [y,fs,nbits]=wavread (39。)。 %求出語音信號的長度 noise=*randn(n,2)。 %語音信號加入噪聲 S=fft(s)。 sound(z11)。 %求濾波后的信號 subplot(2,2,1)。g39。 title(39。)。 subplot(2,2,2)。r39。 title(39。)。 subplot(2,2,3)。 title(39。)。 subplot(2,2,4)。 title(39。)。speech39。 n = length (y) 。 %隨機函數(shù)產生噪聲 s=y+noise。 %傅里葉變換 z12=filter(num12,den12,s)。 m12=fft(z12)。 plot(abs(S),39。)。濾波前信號的頻譜 39。 subplot(2,2,2)。r39。 title(39。)。 plot(s)。濾波前信號的波形 39。 subplot(2,2,4)。 title(39。)。speech39。 n = length (y) 。 %隨機函數(shù)產生噪聲 s=y+noise。 %傅里葉變換 z13=filter(num13,den13,s)。 m13=fft(z13)。 plot(abs(S),39。)。濾波前信號的頻譜 39。 subplot(2,2,2)。r39。 title(39。)。 plot(s)。濾波前信號的波形 39。 subplot(2,2,4)。 title(39。)。speech39。 n = length (y) 。 %隨機函數(shù)產生噪聲 s=y+noise。 %傅里葉變換 z21=fftfilt(b21,s)。 m21=fft(z21)。 plot(abs(S),39。)。濾波前信號的頻譜 39。 subplot(2,2,2)。r39。 title(39。)。 plot(s)。濾波前信號的波形 39。 subplot(2,2,4)。 title(39。)。speech39。 n = length (y) 。 %隨機函數(shù)產生噪聲 s=y+noise。 %傅里葉變換 z22=fftfilt(b22,s)。 m22=fft(z22)。 plot(abs(S),39。)。濾波前信號的頻譜 39。 subplot(2,2,2)。r39。 title(39。)。 plot(s)。濾波前信號的波形 39。 subplot(2,2,4)。 title(39。)。speech39。 n = length (y) 。 %隨機函數(shù)產生噪聲 s=y+noise。 %傅里葉變換 z23=fftfilt(b23,s)。 m23=fft(z23)。 plot(abs(S),39。)。 plot(abs(m23),39。)。濾波后信號的頻譜 39。 subplot(2,2,3)。 title(39。)。 plot(z23)。濾波后的信號波形 39。 15 圖 14 窗函數(shù)法高通濾波