【正文】 abs(H3)。[N,wp0]=ellipord(wp,ws,rp,rs)。)。title(39。xlabel(39。subplot(2,1,1)。wp=[2*fpl/Fs,2*fpu/Fs]。ylabel(39。subplot(2,1,2)。低通濾波損耗函數(shù)曲線39。%調(diào)用ellip計算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和Ay1t=filter(B,A,st)。低通濾波fp=280。)axis([0,Fs/5,0,])。title(39。grid。fc3=Fs/40。f=k/Tp。調(diào)用filter，用三個濾波器分別對信號st進行濾波，分離出三路不同載波頻率的調(diào)幅信號y1(n)、y2(n)和y3(n)繪圖顯示y1(n)、y2(n)和y3(n)的時域波形和幅頻特性曲線結(jié)束通過觀察，三路信號的載波頻率分別為250，500，1000Hz,所以低通濾波器的截至頻率250Wc500，帶通濾波器的截至頻率250W1500,500W21000。subplot(3,2,6)。title(39。)。%全通濾波器系統(tǒng)函數(shù)分子Az1=[1,zeros(1,R1),a]。y=y(:,1)。plot(angle(Y2))。無限個回聲濾波器時域圖 39。subplot(3,2,2)。%無限個回聲濾波器系統(tǒng)函數(shù)分母y2=filter(Bz,Az,y)。z=[zeros(3000,1)。subplot(3,2,6)。title(39。)。%多重回聲濾波器系統(tǒng)函數(shù)分子Az1=[1,zeros(1,R1),a]。z=[zeros(3000,1)。subplot(3,2,6)。title(39。)。%單回聲濾波器系統(tǒng)函數(shù)分子Az1=(1)。)。)。plot(y1)。)。)。plot(z)。)。plot(abs(Y))。采樣后語音信號的時域波形和頻譜[y,fs,bits]=wavread()。title(39。%播放原始信號Y=fft(y)。（3）無限個回聲濾波器的系統(tǒng)函數(shù): , 1 （4）全通結(jié)構(gòu)的混響器的系統(tǒng)函數(shù):全通濾波器的傳遞函數(shù)公式為H(Z)=y(Z)/X(Z)=(K+Z^(m))/(1K*Z^(m))其中m為回聲延時取樣，k為反饋系數(shù)。（2） 頻率計算。若x為向量且長度小于N，則函數(shù)將x補零至長度N。下面介紹這些函數(shù)。*反饋高頻比（Hi Ratio），即反饋回路上的高頻衰減控制。調(diào)治深度：指上述調(diào)頻電路的調(diào)治深度。高頻滾降：此項參數(shù)用于模擬自然混響當(dāng)中，空氣對高頻的吸收效應(yīng)，以產(chǎn)生較為自然的混響效果。 語音信號的延時與混響一、信號采樣（1）采樣頻率采樣頻率是指計算機每秒鐘采集多少個聲音樣本，采樣頻率越高，即采樣的間隔時間越短，則在單位時間內(nèi)計算機得到的聲音樣本數(shù)據(jù)就越多，對聲音波形的表示也越精確。Hz39。)。f=fs*(0:511)/1024。濾波前的時域波形39。N=ceil(*pi/wdelta)。（4）漢明窗的FIR低通濾波使用窗函數(shù)法,選用海明窗設(shè)計了數(shù)字FIR低通濾波器對加了正弦噪聲的語音信號進行濾波，并繪制了濾波器濾波后的語音信號時域圖和頻譜圖。濾波后的信號頻譜39。)。濾波后信號的波形39。subplot(211)。 pause。[b,a]=lp2lp(Bap,Aap,Wn)。Ts=1/Fs。fb = 1000。fudu39。)。figure(2)subplot(2,1,1)plot(abs(y1))title(39。加正弦噪聲后語音信號時域圖39。 %對信號做n點FFT變換f=fs*(0:511)/n。 axis([0 4500 0 400])（2）wav語音信號加噪聲在MATLAB軟件平臺下，給原始的語音信號疊加上噪聲，噪聲類型分為如下幾種：（1）單頻噪色（正弦干擾）；（2）高斯隨機噪聲。title(39。title(39。 %播放語音信號y1=fft(x,n)。在頻譜特性中分析最大值的位置（可能有幾個），它代表的頻率和時域的采樣時間有關(guān)，相鄰的兩點之間的距離為。,[N1 N2])。脈沖響應(yīng)不變法巴特沃斯濾波器對語音信號去噪流程圖GUI 設(shè)計即可以基本的MATLAB 程序設(shè)計為主,也可以鼠標(biāo)為主利用GUIDE 工具進行設(shè)計。 語音信號去噪處理主要要求：1) 在windows系統(tǒng)下的錄音機錄制一段1s左右的語音信號作為原聲信號，在MATLAB軟件平臺下，利用函數(shù)wavread對語音信號進行采樣，記住采樣頻率和采樣點數(shù)；2) 畫出語音信號的時域波形，對采樣后的語音進行fft變換，得到信號的頻譜特性；對語音信號分別加入正弦噪聲和白噪聲，畫出加噪信號的時域波形和頻譜圖；3) 根據(jù)對加噪語音信號譜分析結(jié)果，確定濾除噪聲濾波器的技術(shù)指標(biāo)，設(shè)計合適的數(shù)字濾波器，并畫出濾波器的頻域響應(yīng)；4) 用所設(shè)計的濾波器對加噪的信號進行濾波，在同一個窗口畫出濾波前后信號的時域圖和頻譜圖，對濾波前后的信號進行對比，分析信號變化；5) 利用sound(x)回放語音信號，驗證設(shè)計效果。用MATLAB提供的GUI界面開發(fā)環(huán)境設(shè)計用戶圖形界面使布局編程簡化。圖形界面的設(shè)計 圖形用戶界面GUI (Graphics User Interface) 是由各種圖形對象,如圖形窗口、圖軸、菜單、按鈕、文本框等構(gòu)建的用戶界面,是人機交流信息的工具和方法。快速傅里葉變換，并且畫頻譜圖設(shè)定濾波器性能指標(biāo)，通帶截止頻率fb=，阻帶截止頻率fc=，通帶波紋Ap=，阻帶波紋As=脈沖響應(yīng)不變法法設(shè)計巴特沃斯濾波器窗函數(shù)法設(shè)計FIR濾波器驗證并進行頻譜分析Blip39。在MATLAB中，可以利用函數(shù)fft對信號進行快速傅里葉變化，得到信號的頻譜特性。 N=length(x) fssound(x,fs,bits)。)subplot(3,1,2),plot(magy1)。)figure(3)plot(f,abs(y1(1:2048)))。)。y1=fft(x,n)。)subplot(2,1,2)plot(x2) %做原始語音信號的時域圖形title(39。y2=fft(x2,n)。fudu39。ylabel(39。用雙線性變換法設(shè)計了巴特沃斯數(shù)字低通IIR濾波器對兩加噪語音信號進行濾波，并繪制了巴特沃斯低通濾波器的幅度圖和兩加噪語音信號濾波前后的時域圖和頻譜圖。Fs=22050。 %創(chuàng)建butterworth模擬濾波器[Bap,Aap]=zp2tf(Z,P,K)。)。z=filter(bd,ad,x)。title(39。原始信號的頻譜39。 title(39。)。wdelta=wswp。figure(2)subplot(2,1,1)plot(t,x)title(39。 %播放濾波后的語音信號F0=fft(f1,1024)。Hz39。)xlabel(39。)。常見參數(shù)有以下幾種：混響時間：能逼真的模擬自然混響的數(shù)碼混響器上都有一套復(fù)雜的程序，其上雖然有很多技術(shù)參數(shù)可調(diào)，然而對這些技術(shù)參數(shù)的調(diào)整都不會比原有的效果更為自然，尤其是混響時間。此項技術(shù)可以有效的消除延時聲陣列的段裂聲，可以增加混響聲的柔和感。*反饋增益（FB Gain），即延時反饋的增益控制。三、離散傅立葉變換在MATLAB的信號處理工具箱中函數(shù)FFT和IFFT用于快速傅立葉變換和逆變換。函數(shù)執(zhí)行N點的FFT。用MATLAB工具箱函數(shù)fft進行頻譜分析時需注意：（1） 函數(shù)fft返回值y的數(shù)據(jù)結(jié)構(gòu)對稱性。（2）多重回聲濾波器為了產(chǎn)生以間隔R個抽樣周期分開的具有指數(shù)衰減振幅的多重回聲，可用一個以下形式傳輸函數(shù)的FIR濾波器： 無限個振幅以指數(shù)衰減間隔為R個抽樣周期的多重回聲可用以下形式傳輸函數(shù)的IIR濾波器生成： IIR多重回聲濾波器的基本頻率FR=Fs/R，通常鎖定在伴音設(shè)備的基頻上，比如基鼓拍子。sound(y,fs,bits)。plot(abs(Y))。)。subplot(2,2,3)。采樣信號相頻39。subplot(2,2,1:2)。延時后幅頻39。%回放延遲信號對采樣后的信號混響[y,fs,bits]=wavread(39。subplot(2,2,1:2)?；祉懙姆l39。%回放混響信號單回聲濾波器程序設(shè)計[y,fs,bits]=wavread(39。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯Bz1=[1,zeros(1,R1),a]。單回聲濾波器幅頻響應(yīng)39。plot(y2)。)。)。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯Bz1=[1,zeros(1,N*R1),a^N]。多重回聲濾波器幅頻響應(yīng)39。plot(y2)。)。])。%無限個回聲濾波器系統(tǒng)函數(shù)分子Az=[1,zeros(1,R1),a]。)。title(39。subplot(3,2,6)。)。%濾波器階數(shù)設(shè)置，其值越高，回聲越明顯Bz1=[a,zeros(1,R1),1]。全通濾波器幅頻響應(yīng)39。plot(yy2)。)。顯示st的時域波形和幅頻特性曲線調(diào)用ellipord和ellip分別設(shè)計三個橢圓濾波器，并繪圖顯示其幅頻響應(yīng)特性曲線。k=0:N1。fm2=fc2/10。 %計算信號st的頻譜subplot(2,1,1)plot(t,st)。axis([0,Tp/8,min(st),max(st)])。(b) s(t)的頻譜39。%三路調(diào)幅信號的載波頻率分別為250Hz、500Hz、1000Hz。%調(diào)用ellipord計算橢圓DF階數(shù)N和通帶截止頻率wp[B,A]=ellip(N,rp,rs,wp0)。title(39。)。)。fsu=900。figure。)。plot(t,y2t)。y2t39。rs=60。[H3,w]=freqz(B,A,1000)。w39。高通濾波后的波形39。部分結(jié)果展示：心電信號的處理開始向心電信號中加入白噪聲，50HZ工頻干擾，帶限chirp噪聲使用Butterworth型低通數(shù)字濾波器對信號過濾使用切比雪夫I型數(shù)字低通濾波器對信號過濾使用帶阻濾波器對信號過濾進行效果比較結(jié)束一、設(shè)計原理（1）、巴特沃斯低通濾波器巴特沃斯低通濾波器的平方幅度響應(yīng)為其中，n為濾波器的階數(shù)，為低通濾波器的截止頻率。 ωs=, a s=40Db（6）模擬低通的技術(shù)指標(biāo)為：歸一化截止角頻率wp=2pi*Fs/Ft。,39。IIR數(shù)字濾波器的特征是：①具有無限持續(xù)時間沖激響應(yīng)；②需要用遞歸模型來實現(xiàn)，這可以從其差分方程得出，也可以從其系數(shù)函數(shù)為：得出。C:\Documents and Settings\hnie\桌面\心電信號數(shù)據(jù)\39。原始波形圖39。y1=fft(a(:,2),1000)。頻率（Hz）39。s=a(:,2)+q。)。title(39。)。f2=100*(0:999)/1000。時間（s）39。加入工頻干擾后的頻譜圖39。 加入工頻干擾后的時域圖和頻譜圖，發(fā)現(xiàn)加噪聲后的幅值有比較微小的變化，而在頻譜圖在38Hz和63Hz附近增值很大，這是由于50Hz的工頻噪聲造成的，而在其他頻率范圍內(nèi)也有比較明顯的變化。title(39。)。)。f1=5。)。subplot(2,1,1)plot(y)。幅值（A）39。)。（7）切比雪夫I型數(shù)字低通濾波器用Wp1,Wp2,Ws1,Ws2表示分別將通帶，阻帶截止頻率的角頻率表示,算出頻帶寬帶，計算階數(shù)n1和截止頻率WN，再設(shè)計切比雪夫I型模擬濾波器，采用雙線性法將模擬濾波器系數(shù)變?yōu)閿?shù)字濾波器系數(shù)，畫出切比雪夫I型數(shù)字濾波器的頻率響應(yīng)，調(diào)用filter實現(xiàn)對白噪聲的濾波，再最后調(diào)用plot函數(shù)畫出濾除白噪聲后的時域圖和頻域圖。Ws1=(f21/fs)*2*pi。W0=Omegap1*Omegap2。s39。plot(y)。幅值（A）39。)。它們在時域圖上差不多是一致的，在濾波后的就明顯變小了。figure(10)。)。title(39。)。初次接觸了MATLAB提供的用戶界面的設(shè)計問題，對其中的一些基本問題有了了解，為以后的實踐打下了一定的基礎(chǔ)，相信在以后的實踐過程中會有更好的提升。,39。, gui_Singleton, ... 39。, [])。% UIWAIT makes GUIDE wait for user response (see UIRESUME)% uiwait()。f=fs*(0:2047)/N。原始信號幅值39。原始語音信號頻譜39。[x,fs,bits]=wavread(39。x2=x+x1。)xlabel(39。title(39。subplot(2,1,2),plot(abs(y2))。)。wp=*pi。[N,Wn]=buttord(wp1,ws1,Rp,Rs,39。figure(1)subplot(111)。n=length(x)。)。plot(f,abs(X(1:n/2)))。subplot(2,1,2)。wavwrite(z,fs,39。ws=*pi。plot(w1*fs/(2*pi),20*log10(abs(h1)))。sound(f1)。)xlabel(39。濾波后的頻譜39。h:\39。figure(1)magy1=abs(y1)。)subplot(3,1,3),plot(angy1)。)xlabel(39。)。延時后時域波形39。plot(angle(Z))。% Executes on button press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)% hObject handle to pus