基于matlab的語音信號的時、頻域分析課程設計(編輯修改稿)
2024-07-09 08:41 本頁面
【文章內(nèi)容簡介】 最好的，但通帶和阻帶內(nèi)均呈現(xiàn)等波紋幅頻特性，線性特性的非線性稍重。 IIR 數(shù)字濾波器最大的優(yōu)點是給定一組指標時，它 的階數(shù)要比相同組的FIR 濾波器低得多。 IIR 數(shù)字濾波器的設計方法是利用模擬濾波器成熟的理論及設計圖進行的，因而保留了一些典型模擬濾波器優(yōu)良的幅度特性。但設計中只考慮了幅度特性，沒考慮相位特性，所設計的濾波器一般是某種確定的非線性相位特性。為了得到線性相位特性，對 IIR 濾波器必須另外增加相位校正網(wǎng)絡。 FIR 濾波器在保證幅度特性滿足技術(shù)要求的同時，很容鄭州輕工業(yè)學院 MATLAB 課程設計說明書 3 易做到有嚴格的線性相位特性。穩(wěn)定和線性相位是 FIR濾波器最突出的優(yōu)點。 離散傅立葉變換 ????? 10 )()]([)( Nn knNWnxnxD F TkX 其 中 NjeW ?2?? ,N 為 DFT 變換的區(qū)間長度。 2 信號采集及讀取 用錄音工具錄一段自己的聲音，轉(zhuǎn)化為 wav 格式的文件，然后保存到MATLAB 的文件夾下的 work文件夾中。在 MATLAB 軟件平臺下，利用 wavread函數(shù)對語音信號進行采樣，記住采樣頻率和采樣點數(shù)。 Wavread 函數(shù)調(diào)用格式： y=wavread(file),讀取 file所規(guī)定的 wav文件，返回采樣值放在向量 y中。 [y,fs,nbits]=wavread(‘ ’ )； 把語音信號進行載 MATLAB 仿真軟件平臺中，采樣值放在向量 y中， fs 表示采樣頻率（ hz）， nbits 表示采樣位數(shù)。 sound（ y,fs,nbits） ； 回放語音信號 N=length(y)； 求出語音信號的長度 Y=fft(y,N)； 傅立葉變換 調(diào)用函數(shù) y為被變換的時域系列向量， N是 DFT 變換區(qū)間長度，當 N大于 y的長度時， fft 函數(shù)自動在 y后面補零，函數(shù)返回 y 的 N 點 DFT 變換結(jié)果向量 Y。當 N 小于 y的長度時， fft 函數(shù)計算 y的前面 N個元素構(gòu)成的 N點 DFT，忽略 y后面的元素。首先畫出語音信號的時域波形，然后對語音 信號進行頻譜分析。在MATLAB 中利用 fft 對新號進行快速傅立葉變換，得到信號的頻譜特性。 subplot(2,1,1)。plot(y)。title(‘ 原始信號波形 ’ )。 subplot(2,1,2)。plot(abs(Y))。title(‘ 原始信號頻譜 ’ )； 鄭州輕工業(yè)學院 MATLAB 課程設計說明書 4 程序如下： [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。)。 程序結(jié)果如下圖： 鄭州輕工業(yè)學院 MATLAB 課程設計說明書 5 程序如下： [y,fs,nbits]=wavread (39。39。)。 N = length (y) 。 %求出語音信號的長度 noise=*randn(N,1)。 %隨機函數(shù)產(chǎn)生噪聲 Si=y+noise。 %語音信號加入噪聲 sound(Si)。 subplot(2,1,1)。 plot(Si)。
畢業(yè)設計相關(guān)推薦
鄂ICP備17016276號-1