【正文】 時(shí)域信號(hào)的 FFT分析 信號(hào)的頻譜分析就是計(jì)算信號(hào)的傅立葉變換。 IIR 數(shù)字濾波器的設(shè)計(jì)方法是利用模擬濾波器成熟的理論及設(shè)計(jì)圖進(jìn)行的，因而保留了一些典型模擬濾波器優(yōu)良的幅度特性。title(‘ 原始信號(hào)波形 ’ )。)。39。)。 %過渡帶寬度 M=ceil(*pi/wdelta)。 title(39。 xlabel(39。 鄭州輕工業(yè)學(xué)院 MATLAB 課程設(shè)計(jì)說明書 9 p=110.^(rp/20)。)。)。39。xlabel(39。)。 grid。在編寫程序的時(shí)候遇到了很多錯(cuò)誤，最后經(jīng)過查找一些資料和自己的理解都一一糾正了。 grid。)。濾波后信號(hào)頻譜 39。)。)。幅度 39。 grid。 wp=2*Fp/Ft。頻率響應(yīng) 39。)。)。時(shí)間 39。ylabel(39。原始信號(hào)波形 39。首先畫出語音信號(hào)的時(shí)域波形，然后對(duì)語音 信號(hào)進(jìn)行頻譜分析。因此， H(z)永遠(yuǎn)穩(wěn)定。 通俗的講采樣頻率是指計(jì)算機(jī)每秒鐘采集多少個(gè)聲音樣本，是描述聲音文件的音質(zhì)、音調(diào)，衡量聲卡、聲音文件的質(zhì)量標(biāo)準(zhǔn)。連續(xù)信號(hào)與系統(tǒng)的傅立葉分析顯然不便于用計(jì)算機(jī)進(jìn)行計(jì)算，使其應(yīng)用受到限制。但設(shè)計(jì)中只考慮了幅度特性，沒考慮相位特性，所設(shè)計(jì)的濾波器一般是某種確定的非線性相位特性。 subplot(2,1,2)。ylabel(39。)。 S=fft(Si)。 %濾波器長度，朝正無窮方向舍入 N=2*M+1。脈沖響應(yīng) 39。頻率 (Hz)39。 s=10.^(rs/20)。 N=length(y)。 axis([0 50 2 2])。)。時(shí)間 39。ylabel(39。 subplot(2,2,4)?？偟膩碚f，通過這次的課程設(shè)計(jì)我對(duì)語音信號(hào)有的更多的認(rèn)識(shí)，對(duì)數(shù)字信號(hào)處理的知識(shí)有了更多的理解，鍛煉了我的理論與實(shí)際相結(jié)合的能力，并且熟練的應(yīng)用 MATLAB 也加深了我對(duì)課程的理解。)。時(shí)間 39。 title(39。g39。濾波后信號(hào)頻譜 39。 ylabel(39。)。 Fs=1200。 title(39。n39。幅度 39。xlabel(39。)。 title(39。當(dāng) N 小于 y的長度時(shí)， fft 函數(shù)計(jì)算 y的前面 N個(gè)元素構(gòu)成的 N點(diǎn) DFT，忽略 y后面的元素。（ ） H（ z）是 z1 的 N1次多項(xiàng)式，它在 z平面上有 N1個(gè)零點(diǎn)，在原點(diǎn) z=0處有一個(gè) N1 重極點(diǎn)。 采樣頻率只能用于周期性采樣的采樣器，對(duì)于非周期性采樣的采樣器沒有規(guī)則限制。而 FFT 是一種時(shí)域和頻域均離散化的變換，適合數(shù)值運(yùn)算，成為用計(jì)算機(jī)分析離散信號(hào)和系統(tǒng)的有力工具。為了得到線性相位特性，對(duì) IIR 濾波器必須另外增加相位校正網(wǎng)絡(luò)。plot(abs(Y))。幅度 39。 N = length (y) 。 %傅里葉變換 subplot(2,1,2)。 %窗口長度 wc=(ws+wp)/2。)。)。 fpts=[wp,ws]。 Noise=*randn(N,1)。 subplot(2,1,2)。 N = length (y) 。)。幅度 39。 plot(z21)。這次設(shè)計(jì)使我更深刻的了解 MATLAB 的使用方法，同時(shí)我相信，進(jìn)一步加強(qiáng)對(duì) MATLAB 的學(xué)習(xí)對(duì)我今后的學(xué)習(xí)將會(huì)起到很大的幫助。幅度 39。xlabel(39。)。 plot(abs(S),39。 title(39。)。FIR 低通濾波器 39。 Fp=1000。)。 xlabel(39。ylabel(39。)。時(shí)間 39。 plot(y)。 sound（ y,fs,nbits） ； 回放語音信號(hào) N=l