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