課程四設(shè)計(jì)_基于iir的有噪聲語音信號處理-文庫吧資料
2025-06-13 00:24本頁面
【正文】 um11,den11,s)。 %隨機(jī)函數(shù)產(chǎn)生噪聲 s=y+noise。 n = length (y) 。speech39。)。 title(39。 %傅 里葉變換 subplot(2,1,2)。)。title(39。 subplot(2,1,1)。 %隨機(jī)函數(shù)產(chǎn)生噪聲 Si=y+Noise。 N = length (y) 。D:\MATLAB\39。在這里，我們選用 Randn（ m,n）函數(shù)。 程序結(jié)果如下圖： 8 圖 5 原始信號波形及頻譜 本文中，利用 MATLAB 中的隨機(jī)函數(shù) (rand 或 randn)產(chǎn)生噪聲加入到語音信號中，模仿語音信號被污染，并對其頻譜分析。幅值 n39。)。 xlabel(39。原始語音信號波形 39。) subplot(2,1,2)。 plot(abs(y(1:512))) %做原始語音信號的 FFT頻譜圖 title(39。 %對信號做 FFT變換 f=fs*(0:511)/1024。 sound(y1,44000)。C:\Program Files\MATLAB\R2021a\39。 在本次設(shè)計(jì)中，我們利用 fft 對語音信號進(jìn)行快速傅里葉變換，就可以 得 到信號的頻譜特性。 MATLAB 提供了快速傅里葉變換算法 FFT計(jì)算 DFT的函數(shù) fft,其調(diào)用格式如下： Xk=fft(xn,N) 參數(shù) xn 為被變換的時(shí)域序列向量， N是 DFT變換區(qū)間長度，當(dāng) N大于 xn 的長度時(shí)，fft 函數(shù)自動在 xn 后面補(bǔ)零。其程序如下： [y,fs,nbits]=wavered (‘ OriSound’ )。 （ 4） y=wavread（ file， [N1,N2]） 功能說明：讀取從 N1 到 N2 點(diǎn)的采樣值放在向 量 y 中。 （ 2） [y,fs,nbits]=wavread(file) 功能說明：采樣值放在向量 y中， fs 表示采樣頻率（ hz）， nbits 表示采樣位數(shù)。下面介紹 Wavread 函數(shù)幾種調(diào)用格式。 第三章 仿真及其結(jié)果分析 語音信號的時(shí)頻分析 利用 MATLAB中的 “ wavread” 命令來讀入（采集）語音信號，將它賦值給某一向量。 函數(shù) filter 的調(diào)用格式： yn=filter(B， )，它是按照直線型結(jié)構(gòu)實(shí)現(xiàn)對 xn 的濾波。調(diào) 用 格式 為 ：y=fftfilter(h,x,M)。 6 圖 5 IIR 高通濾波器 用濾波器對加噪語音信號進(jìn)行濾波： 用自己設(shè)計(jì)的各濾波器分別對加噪的語音信號進(jìn)行濾波，在 Matlab 中， IIR 濾波器利用函數(shù) filter 對信號進(jìn)行濾波。用 cheby1設(shè)計(jì) 39。)。 title(39。 %利用雙線性變換實(shí)現(xiàn) S域到 Z域轉(zhuǎn)換 [h,w]=freqz(num13,den13)。 %求 S域的頻率響應(yīng)的參數(shù) [num,den]=lp2hp(b13,a13,wn13)。s39。)。 [n13,wn13]=cheb1ord(wp,ws,1,50,39。 wp=1。 wp1=tan(pi*Fp/Ft)。 Fp=4000。)。 legend(39。 %根據(jù)參數(shù)求出頻率響應(yīng) plot(w*8000*,abs(h))。%將 S域低通參數(shù)轉(zhuǎn)為帶通的 [num12,den12]=bilinear(num2,den2,)。)。 %求低通濾波器階數(shù)和截止頻率 [b12,a12]=butter(n12,wn12,39。s39。 ws=(wp1*wp2w.^2)/(bw*w)。 bw=wp2wp1。 ws2=tan(pi*Fs2/Ft)。 %帶通到低通濾波器的轉(zhuǎn)換 wp2=tan(pi*Fp2/Ft)。)。 legend(39。 %雙線性變換實(shí)現(xiàn) S域到 Z域的變換 [h,w]=freqz(num11,den11)。)。 %求低通濾波器的階數(shù)和截止頻率 [b11,a11]=butter(n11,wn11,39。s39。 fs=2*Fs*tan(wp/2)。 ws=2*pi*Fs/Ft。下面我們將給出 FIR 和 IIR 數(shù)字濾波器的主要程序。 MATLAB信號處理工具箱函數(shù) cheblap,cheblord和 cheeby1是切比雪夫 I型濾波器設(shè)計(jì)函數(shù)。其中可選的窗函數(shù)有 Rectangular Barlrtt Hamming Hann Blackman窗，其相應(yīng)的都有實(shí)現(xiàn)函數(shù)。 hn=fir1(M， wc， window)，可以指定窗函數(shù)向量 window。 （ 4）加窗得到設(shè)計(jì)結(jié)果。 （ 2）構(gòu)造希望逼近的頻率響應(yīng)函數(shù)。先按照阻帶衰減選擇窗函數(shù)類型。下面我們著重介紹用窗函數(shù)法設(shè)計(jì) FIR濾波
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1