基于matlab的語音信號采集與處理課程設(shè)計(jì)報(bào)告-在線瀏覽
2025-01-13 03:34本頁面
【正文】 號 加上一個(gè)高頻噪聲 在 Matlab中人為設(shè)計(jì)一 個(gè) 固定頻率 5500Hz的噪聲干擾信號。,給出的干擾信號為一個(gè) 正 弦信號，針對上面的語音信號 ,采集了其中一段。 (3). 設(shè)計(jì) 一個(gè) 濾波器，濾除高頻噪聲 由模擬濾波器變換為數(shù)字濾波器時(shí) ,采用的是雙線性變換法 , 它保留的是從模擬到數(shù)字域的系統(tǒng)函數(shù)表示。如果 FIR濾波器情況下 ,最優(yōu)的設(shè)計(jì)將是橢圓濾波器。的 3 本設(shè)計(jì)是用雙線性變換法設(shè)計(jì) BW帶阻型濾波器。用設(shè)計(jì)好的帶阻濾波器對含噪語音信號進(jìn) 行濾波 ,在 Matlab 中 FIR 濾波器利用函數(shù) fftfilt對信號進(jìn)行濾波 ,FIR濾波器利用函數(shù) filter對信號進(jìn)行濾波。從圖中可以看出 ,5500Hz 看到的高峰消失了 ,語音信號與開始的一樣 ,濾波器成功的濾除了干擾信號。 三、實(shí)踐步驟 1． 用 MATLAB對原始 語音 信 號進(jìn)行分析，畫出它的時(shí)域波形和頻譜 程序： fs=22050。C:\Documents and Settings\Administrator\桌面 \聲音文件 .wav39。 %讀取語音信號的數(shù)據(jù)，賦給變量 x1， sound(x1,22050)。 %對信號做 1024點(diǎn) FFT變換 f=fs*(0:511)/1024。原始語音信號 39。 xlabel(39。)。fuzhi n39。 %y軸的名字是“ fuzhi n” figure(2) freqz(x1) %繪制原始語音信號的頻率響應(yīng)圖 title(39。) figure(3) subplot(2,1,1)。原始語音信號 FFT頻譜 39。 plot(f,abs(y1(1:512)))。原始語音信號頻譜 39。Hz39。 ylabel(39。)。 對 加噪后的 語音進(jìn)行分析，并畫出其 信號時(shí)域和頻譜圖 。 x1=wavread(39。)。 %將 0到 511，步長為 1的序列的值與 fs相乘并除以 1024的值，賦值給 f t=0:1/fs:(length(x1)1)/fs。 %噪聲幅值 d=[Au*sin(2*pi*5500*t)]39。 %將正弦信號噪聲加在語音信號上 sound(x2,22050)。 %對信號 y1做 1024點(diǎn) FFT變換 y2=fft(x2,1024)。 %創(chuàng)建圖形窗 1 plot(t,x2)。加噪后的信號 39。 xlabel(39。)。fuzhi n39。 %y軸的名字是“ fuzhi n” figure(2) %創(chuàng)建圖形窗 2 subplot(2,1,1)。 %做原始語音信號的頻譜圖 title(39。)。Hz39。 %x軸的名字是“ time n” ylabel(39。)。 %創(chuàng)建兩行一列繪圖區(qū)間的第 2個(gè)繪圖區(qū)間 plot(f,abs(y2(1:512)))。加噪后的信號頻譜 39。 xlabel(39。)。fuzhi39。 % y軸的名字是“ fuzhi” 時(shí)域波形和頻譜 ： 圖 4 加噪后的信號時(shí)域圖 圖 5 原始信號與加噪后信號頻譜對比 8 與原始信號對 比，區(qū)別 ： 先原始信號沒加噪音之前 0到 2020有幅值，在 4000到 6000之間沒有幅值，但是在加了噪音之后 4000到 6000 之間出現(xiàn)最大幅值 12，超出正常值。 3. 設(shè)計(jì)合適的濾波器 ，濾除高頻噪聲，繪出濾波后的信號頻域和時(shí)域波形： 程序： %用雙線性變換法設(shè)計(jì)巴特沃思低
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1