【正文】 2)。ylabel(39。figure(2)subplot(2,1,1)。幅值A(chǔ)39。)。接下來我們要用設(shè)計的FIR低通濾波器對上一節(jié)中加噪語音信號進行濾波處理。sound(y,fs)n=length(y)y_p=fft(y,n)。時間軸39。xlabel(39。sound(y_z,fs)n=length(y)。xlabel(39。)。Fp=1000。fpts=[wp ws]。FIR低通濾波器39。濾波前信號的頻譜39。plot(y_z)。加噪后的語音信號的時域波形比原始語音信號要模糊得多，頻譜圖則是在頻率5000Hz以后出現(xiàn)了明顯的變化。wp=2*pi*Fp/Ft。s39。grid。 %傅里葉變換z11=filter(num11,den11,s)。濾波前信號的頻譜39。濾波后信號的頻譜39。grid。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。這是很好的性質(zhì)。2) 處理速度較快。沒有老師的幫助也就沒有今天的這篇論文。感謝我的輔導(dǎo)員鄭偉老師，謝謝他在這四年中為我們?nèi)嗨龅囊磺?，他不求回報，無私奉獻的精神很讓我感動，再次向他表示由衷的感謝。3) 占用內(nèi)存空間比較小。FIR幅頻特性精度較之于iir低，但是線性相位，就是不同頻率分量的信號經(jīng)過FIR濾波器后他們的時間差不變。 （1）IIR數(shù)字濾波器的相位特性不好控制，對相位要求較高時，需加相位校準網(wǎng)絡(luò)。plot(z11)。grid。grid。m11=fft(z11)。39。 %求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,)。fp=2*Ft*tan(wp/2)。濾波后信號的頻譜圖也在5000Hz以后開始逐漸接近原始語音信號的頻譜圖。濾波前信號的時域波形39。subplot(2,2,2)。x=fftfilt(b21,y_z)。dev=[p s]。wp=2*Fp/Ft。頻率Hz39。)ylabel(39。f=fs*(0:n/21)/n。)。幅值A(chǔ)39。figure(1)subplot(2,1,1)。函數(shù)fftfilt用的是重疊相加法實現(xiàn)線性卷積的計算。頻率幅值39。plot(f,abs(y_zp(1:n/2)))。title(39。)。title(39。原始語音信號采樣后的時域波形39。 下面一段程序?qū)崿F(xiàn)了利用randn函數(shù)把一段隨機噪音信號加入原始語音信號的信號處理過程：[y,fs,bits]=wavread(39。)。plot(magX)。plot(x)。它的FFT也是N點序列，即X(k)(k=0,1,2,…,N1)，則第k點所對應(yīng)實際頻率值為f=k*f /N. 下面一段程序是語音信號在matlab中的簡單體現(xiàn)，通過該段程序，我們實現(xiàn)了語音信號的讀入與打開，并繪出了語音信號的波形頻譜圖。若向量x的長度大于N，則函數(shù)截短x使之長度為N。其中，x是序列，y是序列的FFT，x可以為一向量或矩陣，若x為一向量，y是x的FFT。利用MATLAB中的wavread命令來讀入（采集）語音信號，將它賦值給某一向量。量化過程中由于存在不同程度的量化誤差，由此會導(dǎo)致濾波器的頻率響應(yīng)出現(xiàn)偏差，嚴重時會使IIR濾波器的極點移到單位圓之外，系統(tǒng)因而失去穩(wěn)定性。遞歸濾波器，也就是IIR數(shù)字濾波器，顧名思義，具有反饋，一般認為具有無限的脈沖響應(yīng)。FIR濾波器基本結(jié)構(gòu)FIR濾波器有以下幾種基本結(jié)構(gòu)： 式的系統(tǒng)的差分方程表達式為 y(n)=∑h(m)x(nm) ( 24）很明顯，這就是線性移不變系統(tǒng)的卷積和公式，也是x (n)的延時鏈的橫向結(jié)構(gòu)，稱為橫截型結(jié)構(gòu)或卷積型結(jié)構(gòu)，也可稱為直接型結(jié)構(gòu)。然后，通過分析經(jīng)典譜減法的原理及其一般改進形式，研究了一種基于噪聲殘差的譜相減改進算法和一種可以不以噪聲是零均值的高斯分布為前提的譜減法改進算法。 去噪的方法數(shù)字信號處理技術(shù)經(jīng)過幾十年的發(fā)展，在國內(nèi)外已經(jīng)取得了很大的成績。采樣頻率與聲音頻率之間有一定的關(guān)系，根據(jù)采樣定理，只有采樣頻率高于聲音信號最高頻率的兩倍時，才能把數(shù)字信號表示的聲音還原成為原來的聲音。2Δt)，...來表示,只要這些采樣點的時間間隔Δt≤1/2F，便可根據(jù)各采樣值完全恢復(fù)原來的信號f(t)。3. 繪制頻譜圖并回放原始語音信號。近年來人工神經(jīng)網(wǎng)絡(luò)(ANN)的研究取得了迅速發(fā)展，語音信號處理的各項課題是促進其發(fā)展的重要動力之一，他的各項成果也體現(xiàn)在語音信號處理的各項技術(shù)之中。s most important, most effective, most popular and most convenient form of exchange of information.. Matlab language is a data analysis and processing functions are very powerful puter application software, sound files which can be transformed into discrete data files, then use its powerful ability to process the data matrix operations, such as digital filtering, Fourier transform, when domain and frequency domain analysis, sound playback and a variety of map rendering, and so on. Its signal processing and analysis toolkit for voice signal analysis provides a very rich feature function, use of these functions can be quick and convenient features plete voice signal processing and analysis and visualization of signals, makes puter interaction more convenient . Matlab Signal Processing is one of the important areas of application. The design of voiceprocessing software for most of the content are numerous, easy to maneuver and so on, using prehensive use GUI interface design, various function calls to voice signals such as frequency, amplitude, Fourier transform and filtering, the program interface concise, simple, has some significance in practice. Finally, the speech signal processing further develop