基于matlab的語(yǔ)音信號(hào)的時(shí)、頻域分析課程設(shè)計(jì)-文庫(kù)吧資料
2025-01-24 14:16本頁(yè)面
【正文】 。 plot(w/pi,180/pi*unwrap(angle(hf)))。)。title(39。幅度(dB)39。)。xlabel(39。subplot(3,1,2)。脈沖響應(yīng)39。)。 ylabel(39。n39。 stem(n,hi)。%得到脈沖響應(yīng)[hf w]=freqz(b,1,512)。n=0:1:N。 %截止頻率win=hamming(N)。 %濾波器長(zhǎng)度，朝正無(wú)窮方向舍入N=2*M+1。 %性能指標(biāo)wdelta=wswp。程序結(jié)果如下圖：4 數(shù)字濾波器設(shè)計(jì) 用窗函數(shù)法設(shè)計(jì)IIR帶通濾波器程序如下：wp=*pi。幅度39。)。 xlabel(39。加噪語(yǔ)音信號(hào)的頻域波形39。plot(abs(S))。S=fft(Si)。幅度39。)。xlabel(39。加噪語(yǔ)音信號(hào)的時(shí)域波形39。plot(Si)。 %語(yǔ)音信號(hào)加入噪聲 sound(Si)。 %求出語(yǔ)音信號(hào)的長(zhǎng)度noise=*randn(N,1)。)。程序結(jié)果如下圖：3構(gòu)造受干擾信號(hào)并對(duì)其進(jìn)行FFT頻譜分析程序如下：[y,fs,nbits]=wavread (39。幅度39。)。xlabel(39。原始信號(hào)頻譜39。plot(abs(Y))。)。ylabel(39。時(shí)間39。)。title(39。subplot(2,1,1)。N=length(y)。)。title(‘原始信號(hào)頻譜’)；程序如下：[y,fs,nbits]=wavread(39。subplot(2,1,2)。plot(y)。在MATLAB中利用fft對(duì)新號(hào)進(jìn)行快速傅立葉變換，得到信號(hào)的頻譜特性。當(dāng)N小于y的長(zhǎng)度時(shí)，fft函數(shù)計(jì)算y的前面N個(gè)元素構(gòu)成的N點(diǎn)DFT，忽略y后面的元素。[y,fs,nbits]=wavread(‘’)；把語(yǔ)音信號(hào)進(jìn)行載MATLAB仿真軟件平臺(tái)中，采樣值放在向量y中，fs表示采樣頻率（hz），nbits表示采樣位數(shù)。在MATLAB軟件平臺(tái)下，利用wavread函數(shù)對(duì)語(yǔ)音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。 其中,N為DFT變換的區(qū)間長(zhǎng)度。FIR濾波器在保證幅度特性滿足技術(shù)要求的同時(shí)，很容易做到有嚴(yán)格的線性相位特性。但設(shè)計(jì)中只考慮了幅度特性，沒(méi)考慮相位特性，所設(shè)計(jì)的濾波器一般是某種確定的非線性相位特性。IIR數(shù)字濾波器最大的優(yōu)點(diǎn)是給定一組指標(biāo)時(shí)，它的階數(shù)要比相同組的FIR濾波器低
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1