【正文】 信號在matlab中的簡單體現(xiàn)，通過該段程序，我們實現(xiàn)了語音信號的讀入與打開，并繪出了語音信號的波形頻譜圖。plot(magX)。 下面一段程序?qū)崿F(xiàn)了利用randn函數(shù)把一段隨機噪音信號加入原始語音信號的信號處理過程：[y,fs,bits]=wavread(39。title(39。title(39。頻率幅值39。figure(1)subplot(2,1,1)。)。)ylabel(39。wp=2*Fp/Ft。x=fftfilt(b21,y_z)。濾波前信號的時域波形39。fp=2*Ft*tan(wp/2)。39。grid。plot(z11)。FIR幅頻特性精度較之于iir低，但是線性相位，就是不同頻率分量的信號經(jīng)過FIR濾波器后他們的時間差不變。感謝我的輔導(dǎo)員鄭偉老師，謝謝他在這四年中為我們?nèi)嗨龅囊磺校磺蠡貓?，無私奉獻的精神很讓我感動，再次向他表示由衷的感謝。2) 處理速度較快。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準(zhǔn)型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。濾波后信號的頻譜39。 %傅里葉變換z11=filter(num11,den11,s)。s39。加噪后的語音信號的時域波形比原始語音信號要模糊得多，頻譜圖則是在頻率5000Hz以后出現(xiàn)了明顯的變化。濾波前信號的頻譜39。fpts=[wp ws]。)。sound(y_z,fs)n=length(y)。時間軸39。接下來我們要用設(shè)計的FIR低通濾波器對上一節(jié)中加噪語音信號進行濾波處理。幅值A(chǔ)39。ylabel(39。plot(y)。title(39。angX=angle(X)。函數(shù)執(zhí)行N點的FFT。 語音文件在MATLAB平臺上的錄入與打開單擊自己的電腦開始程序，選擇所有程序，接著選擇附件，再選擇娛樂，最后選擇錄音。但是這種結(jié)構(gòu)所需要的系數(shù)B2k（I = 0，1，2，k，= 1，2，．．．，[N/2]）比卷積型的系數(shù)h (n)要多，因而所需的乘法次數(shù)也比卷積型的要多。在分析了語音增強算法理論的基礎(chǔ)上，本文首先研究了語音激活檢測算法。 通俗的講采樣頻率是指計算機每秒鐘采集多少個聲音樣本，是描述聲音文件的音質(zhì)、音調(diào)，衡量聲卡、聲音文件的質(zhì)量標(biāo)準(zhǔn)。為什么把采樣頻率設(shè)為8kHz?在數(shù)字通信中，根據(jù)采樣定理, 最小采樣頻率為語音信號最高頻率的2倍，頻帶為F的連續(xù)信號f(t)可用一系列離散的采樣值f(t1),f(t1177。 國內(nèi)外同行研究現(xiàn)狀 20世紀(jì)60年代中期形成的一系列數(shù)字信號處理的理論和算法，如數(shù)字濾波器、快速傅立葉變換（FFT）等是語音信號數(shù)字處理的理論和技術(shù)基礎(chǔ)。宿州學(xué)院畢業(yè)論文設(shè)計基于matlab聲音信號的濾波去噪處理 摘要濾波器設(shè)計在數(shù)字信號處理中占有極其重要的地位。其強大的數(shù)據(jù)處理能力可以極大程度上削弱噪聲影響還原出真實的語音信號相符度在90%以上。 1924年奈奎斯特(Nyquist)就推導(dǎo)出在理想低通信道的最高大碼元傳輸速率的公式: 理想低通信道的最高大碼元傳輸速率=2W*log2 N (其中W是理想低通信道的帶寬,N是電平強度)。 采樣頻率的常用的表示符號是 fs。譜相減語音增強算法的核心是噪聲檢測和譜減規(guī)則。 這種結(jié)構(gòu)的每一節(jié)控制一對零點，因而再需要控制傳輸零點時，可以采用它。在MATLAB下設(shè)計不同類型IIR濾波器均有與之對應(yīng)的函數(shù)用于階數(shù)的選擇。式中，x，y意義同前，N為正整數(shù)。magX=abs(X)。plot(angX)。figure(1)subplot(2,1,1)。)。)ylabel(39。 去噪及仿真（1）FIR濾波器法去噪通過對上一節(jié)中加噪語音信號和原始語音信號頻譜圖對比可以知道，噪音大部分是Hz大于5000的部分，故設(shè)計低通濾波器進行濾波處理。xlabel(39。y_z=y+noise。加噪語音信號頻譜圖39。s=10.^(rs/20)。title(39。)sound(x,fs,bits)而后得到圖像如下：分析：從以上四圖可以很明顯和直觀的看出原始語音信號和加噪語音信號時域波形和頻譜圖的區(qū)別。 %求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,39。 %語音信號加入噪聲S=fft(s)。title(39。 結(jié)合去噪后的頻譜圖對比兩種方式濾波的優(yōu)缺點IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。在進行語音信號處理的操作界面中，菜單按鍵明了，每個功能只對應(yīng)一個按鍵，省去了大型軟件的操作復(fù)雜步驟。在這四年的學(xué)期中結(jié)識的各位生活和學(xué)習(xí)上的摯友讓我得到了人生最大的一筆財富。這是很好的性質(zhì)。 title(39。subplot(2,2,2)。)。fs=2*Fs*tan(wp/2)。)subplot(2,2,4)。X=fft(x,n)。ws=2*Fs/Ft。幅值A(chǔ)39。ylabel(39。plot(y)。)。加噪語音信號時域波形39。原始語音信號采樣后的頻譜圖39。39。title(39。 [x,fs,bits]=wavread(39。且和x相同長度。為了獲得最優(yōu)的濾波器系數(shù)，采用以下步驟進行量化。將轉(zhuǎn)置定理用于可得到的轉(zhuǎn)置直接型結(jié)構(gòu)。到目前為止，已經(jīng)比較成熟的去噪方法比較典型的有切比雪夫去噪法、雙線性變換去噪法、窗函數(shù)去噪法、譜相減去噪法、巴特沃茲去噪法等有名的去噪方法。（22） 采樣頻率采樣頻率，也稱為采樣速度或者采樣率，定義了每秒從連續(xù)信號中提取并組成離散信號的采樣個數(shù)，它用赫茲（Hz）來表示。4. 利用MATLAB編程加入一段正弦波噪音，設(shè)計濾波器去噪。在信號傳輸過程中，由于實驗條件或各種其他主觀或客觀條件的原因，語音處理系統(tǒng)都不可避免地要受到各種噪聲的干擾。利用MATLAB信號處理工具箱可以快速有效地設(shè)計各種數(shù)字濾波器。本論文針對以上問題，運用數(shù)字信號學(xué)基本原理實現(xiàn)語音信號的處理，幅頻變換以及傅里葉變換、濾波等技術(shù)來進行語音信號處理。 時域采樣定理的另一種表述方式是：當(dāng)時間信號函數(shù)f(t)的最高頻率分量為fM時,f(t)的值可由一系列采樣間隔小于或等于1/2fM的采樣值來確定,即采樣點的重復(fù)頻率f≥2fM。采樣位數(shù)和采樣率對于音頻接口來說是最為重要的兩個指標(biāo)，也是選擇音頻接口的兩個重要標(biāo)準(zhǔn)。 由于譜相減算法去噪比較復(fù)雜，我們在這里就不多做討論，并且本論文也不采用此種方法。需要將MATLAB設(shè)計出的IIR數(shù)字濾波器進一步分解和量化，從而獲得可用FPGA實現(xiàn)的濾波器系數(shù)。利用Sound函數(shù)可以清晰的聽到“畢業(yè)設(shè)計聲音文件”的語音。經(jīng)函數(shù)fft求得的序列y一般是復(fù)序列，通常要求其幅值和相位。原始信號波形39。產(chǎn)生均值為0，方差 σ^2 = 1，標(biāo)準(zhǔn)差σ = 1的正態(tài)分布的隨機數(shù)或矩陣的函數(shù)。xlabel(39。y_z=y+noise。加噪語音信號頻譜圖39。其中，h是系統(tǒng)單位沖擊響應(yīng)向量；x是輸入序列向量；y是系統(tǒng)的輸出序列向量；M是有用戶選擇的輸入序列的分段長度，缺省時，默認(rèn)的輸入向量的重長度M=512。plot(f,abs(y_p(1:n/2)))。plot(y_z)。ylabel(39。b21=fir1(n21,wn21,Kaiser(n21+1,bet