語音信號(hào)處理與濾波課程設(shè)計(jì)-畢業(yè)設(shè)計(jì)-全文預(yù)覽
【正文】 bel(39。 subplot(2,1,1)。） 。 [b,a]=firl(N,wn/pi,hamming(N+1))。 Rs=100。)。)。頻率 /Hz39。plot(f,abs(y1(1:512)))。 ylabel(39。FIR 高通濾波器濾波前的頻譜 39。 y2=fft(x1,1024)。FIR高通濾波器 39。high39。 wdelta=wpws。 wp=2*pi*5000/fs。 fs=22050。)。title(39。幅值 39。)； xlabel(39。 subplot(2,1,1)。)。)。 N=ceil(8*i/wdelta)。 Rs=100。wp2=2*pi*3000/fs。 x1=wavread(39。 else if temp==3 IIR_BP。 if temp==1 IIR_LP。IIR 帶通濾波器 39。IIR 低通濾波器 39。 else main。 else if temp==2 FIR_HP。 main。)。,39。,39。)。,39。,39。北京：電子工業(yè)出版社， 2021 [3]王宏 。 一年之計(jì)在于春，而我們一學(xué)期的學(xué)習(xí)就在于這次課程設(shè)計(jì)了，我覺得這才是真正學(xué)到知識(shí)的，也鍛煉了自己的動(dòng)手實(shí)踐能力，耐心和毅力。 選用凱澤窗設(shè)計(jì)方法，濾波器的幅度和相位響應(yīng)滿足設(shè)計(jì)的指標(biāo)但濾波器長(zhǎng) 度（ N=708）太長(zhǎng)，實(shí)現(xiàn)起來很困難，主要原因是濾波器的指標(biāo)太苛刻，因此，一向不用窗函數(shù)法設(shè)計(jì)這種類型的濾波器。 題目要求最大衰減 Ap=1dB， fp=5000Hz，由圖可看出，基本符合，而最小衰減 As=100dB， fc=4800Hz，也基本符合。. 經(jīng)過調(diào)試分析，觀察后發(fā)現(xiàn)是書寫錯(cuò)誤，將語句 wp=[2*fp1/fs,2*fp2/s]。出錯(cuò)，將其改為 N=ceil(()/(*(wcwp)/2))+1。 （ 9）比較濾波前后語音信號(hào)的波形及頻譜，在同一個(gè)窗體中畫出濾波前后的波形圖 : FIR 低通濾波器濾波前后的波形圖： FIR 高通濾波器濾波前后的波形圖： FIR 帶通濾波器濾波前后的波形圖： IIR 低通濾波器濾波前后的波形圖 : IIR 高通濾波器濾波前后的波形圖： IIR 帶通濾波器濾波前后的波形圖： （ 10）回放語音信號(hào) 在 MATLAB 中，函數(shù) sound 可以對(duì)聲音進(jìn)行回放。 [b,a]=ellip(n,Ap,As,wn)。 fs=22050。 fc1=1000。 [b,a]=ellip(n,Ap,As,wn)。 fs=22050。 程序運(yùn)行結(jié)果如下圖： （ 5）高通濾波器的設(shè)計(jì)過程 程序代碼： fp=5000。 wp=2*fp/fs。 As=100。 b=fir1(N,wc,Win)。 wp=2*fp/fs。 As=100。 ）設(shè)計(jì)上面要求的三種濾波器。幅值 39。 xlabel(39。)。頻率 /Hz39。 figure(1) subplot(2,1,1)。首先畫出語音信號(hào)的時(shí)域波形，然后對(duì)語音信號(hào)進(jìn)行頻譜分析。 subplot(2,1,1)。 y=wavread(file,[N1， N2])，讀取從 N1點(diǎn)到 N2點(diǎn)的采樣值放在向量 y中。通過使用 wavread 函數(shù)，理解采樣頻率、采樣位數(shù)等概念。 曾勝： IIR 帶通濾波器的設(shè)計(jì)和用此濾波器對(duì)語音信號(hào)進(jìn)行濾波。 分組情況： 組長(zhǎng)：張鳳 組員：張鳳、張志廣、李云、黃生涯、曾勝、沈呈潔 分工情況： 張鳳：設(shè)計(jì)全過程的監(jiān)督及 協(xié)助、部分源程序代碼的編寫和整個(gè)源程序代碼的整理。 （ 2）掌握序列快速傅里葉變換 FFT 方法。 二、課程設(shè)計(jì)基本要求及分組 對(duì)知識(shí)點(diǎn)的掌握要求 ： （ 1）熟悉離散信號(hào)和系統(tǒng)的時(shí)域特性。 （ 5）掌握 MATLAB 設(shè)計(jì) FIR 和 IIR 數(shù)字濾波器的方法。 黃生涯： IIR 低通濾波器、 IIR 高通濾波器的設(shè)計(jì)。 詳細(xì)操作步驟和部分運(yùn)行結(jié)果 （ 1） 采集語音信號(hào)： 利用 Windows 下的錄音機(jī)，錄制一段自己的話音，時(shí)間控制在 1s 左右；然后在MATLAB 軟件平臺(tái)下，利用函數(shù) wavread 對(duì)語音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。 y=wavread(file,N)，讀取前 N點(diǎn)的采樣值放在向量 y中。 y=fft(x1,1024)。 原始語音采樣后的的時(shí)域信號(hào)波形圖： （ 2）對(duì)語音信號(hào)進(jìn)行頻譜分析： 在 MATLAB 中，利用函數(shù) fft 對(duì)信號(hào)進(jìn)行快速傅里葉變換，得到信號(hào)的頻譜特性。 f=fs*(0:511)/1024
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1