【正文】 sound(y,fs,bits)。title(39。figure。fudu39。頻率響應(yīng)圖39。原始信號波形39。)。ylabel(39。fuzhi39。subplot(2,1,1)。subplot(2,1,2),plot(t,f1)。wn=(+)*pi/2。[x,fs,bits]=wavread(39。)。xlabel(39。濾波后信號的波形39。subplot(211)。)。[Bap,Aap]=zp2tf(Z,P,K)。Fs=22050。)。Hz39。)。fudu39。原語音信號時域圖39。y1=fft(x,N)。ylabel(39。)figure(2)freqz(x)title(39。title(39。H:\39。endif nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。, GUIDE_OutputF, ... 39。gui_State = struct(39。附錄function varargout = GUIDE(varargin)% GUIDE Mfile for % GUIDE, by itself, creates a new GUIDE or raises the existing% singleton*.%% H = GUIDE returns the handle to a new GUIDE or the handle to% the existing singleton*.%% GUIDE(39。 此次課程設(shè)計綜合了數(shù)字信號處理的理論知識和實踐經(jīng)驗，經(jīng)過這段時間查閱的相關(guān)資料并結(jié)合以前學(xué)過的數(shù)字信號處理課程和實驗，我如期完成了設(shè)計任務(wù)要求，并加深了對知識的理解，提高了自己的實踐動手能力。頻率(Hz)39。subplot(2,1,2)。濾除工頻干擾后的時域圖39。stop39。ylabel(39。subplot(2,1,2)。xlabel(39。freqz(num,den,64)。[n1,WN]=buttord(WP,WS,1,50,39。Omegap2=2*fs*tan(Wp2/2)。f12=25。ylabel(39。subplot(2,1,2)。xlabel(39。[bz,az]=bilinear(b,a,fs)。Omegap=2*fs*tan(wp/2)。 加入帶限噪聲后的時域圖和頻譜圖加入帶限chirp噪聲干擾后，發(fā)現(xiàn)加噪之前時域圖是均勻分布的，而加噪后，則變?yōu)榍懊嫦∈韬竺婷芗那闆r了，對于頻譜圖而言，加噪前是中間凹兩邊凸，加噪后是前凸后平，在個別點上幅值增加很多，這是由于加入chirp噪聲的結(jié)果。加入帶限噪聲后的頻譜圖39。時間（s）39。abshc=abs(hc)。)。)。title(39。 t=0::*(10001)。幅度39。幅值（A）39。subplot(2,1,1)。)。title(39。)。t=a(1:1000,1)。2)間接法，目前所用的方法主要是：沖激響應(yīng)不變法、階躍響應(yīng)不變法和雙線性法。s39。通帶最大衰減為：Rp=1。圖1 巴特沃斯低通濾波器的幅頻特性（2） 切比雪夫I型數(shù)字低通濾波器（1）確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶截止頻率ωp、通帶衰減ap、阻帶截止頻率ωs、阻帶衰減as切比雪夫濾波器的振幅平方特性如圖2所示：圖2 切比雪夫濾波器的振幅平方特性(2)將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo)。)。)。title(39。)。fs=900。xlabel(39。ylabel(39。plot(w/pi,20*log(m/max(m)))。rs=60。帶通濾波Fs=10000。低通濾波后的波形39。w39。[H1,w]=freqz(B,A,1000)。ws=2*fs/Fs。)。.39。)。xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t)。k=0:N1。Fs=10000。全通結(jié)構(gòu)相頻 39。subplot(3,2,5)。title(39。YY2=fft(yy2,6001)。y=[y。)。plot(abs(Y2))。無限個回聲濾波器相頻響應(yīng)39。subplot(3,2,1)。zeros(3000,1)]。多重回聲濾波器相頻 39。subplot(3,2,5)。title(39。[h,w]=freqz(Bz1,Az1)。zeros(3000,1)]。單回聲濾波器相頻 39。subplot(3,2,5)。title(39。[h,w]=freqz(Bz1,Az1)。y=[y。title(39。)。y=[y。title(39。)。)。)。plot(y)。subplot(2,2,4)。title(39。原始語音信號時域波形和頻譜[y,fs,bits]=wavread(39。其中向量b和a分別表示系統(tǒng)函數(shù)的分子、分母多項式的系數(shù)，若a＝1，此時表示FIR濾波器，否則就是IIR濾波器。經(jīng)函數(shù)fft求得的序列y一般是復(fù)序列，通常要求其幅值和相位。且和x相同長度。*高頻增益（HF），指高頻均衡控制。當(dāng)延遲時間在3~35ms之間時人耳感覺不到滯后音的存在，并且他與原音源疊加后，會因其相位干涉而產(chǎn)生梳狀濾波效應(yīng)，這就是鑲邊效果。擴(kuò)散度：此項參數(shù)可調(diào)整混響聲陣密度的增長速度，其可調(diào)范圍為0~10，其值較高時，混響效果比較豐厚、溫暖；其值較低時，混響效果則較空曠、冷僻。聲卡的位客觀地反映了數(shù)字聲音信號對輸入聲音信號描述的準(zhǔn)確程度。fuzhi39。)。plot(f,abs(y2(1:512)))。title(39。b=fir1(N,wn/pi,hamming(N+1))。)。Hz39。)。figure(3)sound(z,fs,bits)。原始信號的波形39。)。 figure(1)subplot(111)。 [N,Wn]=buttord(wp1,ws1,Rp,Rs,39。wp=*pi。wavwrite(x2,fs,39。)。xlabel(39。)。x2=x+x1。[x,fs,bits]=wavread(39。Hz39。title(39。angy1=angle(y1)。由于沒有考慮相位和其他的頻譜分量，所以波形和原來的波形相差甚大，但大體的頻率是沒有錯的。sound(y)?？焖俑道锶~變換，并且畫頻譜圖在GUIDE 設(shè)計環(huán)境中,需要用到的工具有屬性編輯器、控件布置編輯器、菜單編輯器、對象瀏覽器、網(wǎng)格標(biāo)尺設(shè)置編輯器以及GUI 應(yīng)用屬性設(shè)置編輯器等。 心電信號的處理主要要求：1) 在MATLAB軟件平臺下，給原始的心電信號疊加上噪聲或干擾，干擾類型分為如下幾種：白噪聲、工頻干擾（50Hz）、諧波干擾（二次、三次諧波為主，分別為100Hz、150Hz）繪出疊加噪聲后的心電信號時域和頻譜圖，在視覺上與原始心電信號圖形對比，分析頻域基本特征變化。DSP課程設(shè)計作者：9月9日10月9日中南大學(xué)通信班目錄緒論1GUI圖形界面設(shè)計4語音信號去噪處理4語音信號的延時與混響4數(shù)字濾波器的設(shè)計與實現(xiàn)4心電信號的處理4設(shè)計小結(jié)4附錄4緒論 《數(shù)字信號處理》課程是通信專業(yè)的一門重要專業(yè)基礎(chǔ)課，是信息的數(shù)字化處理、存儲和應(yīng)用的基礎(chǔ)。2) 給定濾波器的規(guī)一化性能指標(biāo)（參考指標(biāo)，實際中依據(jù)每個同學(xué)所疊加噪聲情況而定），例如：通帶截止頻率wp=*pi, 阻通帶截止頻率ws=*pi。GUI 設(shè)計面板是GUI 設(shè)計工具應(yīng)用的平臺,面板上部提供了菜單和常用工具按鈕,左邊提供了多種GUI 控件,如按鈕、單選按鈕、復(fù)選框、文本框等。設(shè)定濾波器性能指標(biāo)，通帶截止頻率fb=1100，阻帶截止頻率fc=1200，通帶波紋Ap=1，阻帶波紋As=20 用于對聲音的回放。 [x,fs,bits]=wavread(39。subplot(3,1,1),plot(x)。原始信號相位39。)。 H:\39。sound(x2,fs,bits)。ylabel(39。Hz39。xlabel(39。 H:\39。ws=*pi。s39。plot(w*fs/(2*pi),abs(h))grid。n=length(x)。)。subplot(211)。Z=fft(z)。)。n=length(x)。 %選擇窗函數(shù)，并歸一化截止頻率f1=fftfilt(b,x)。濾波后的時域波形39。 %畫出濾波前的頻譜圖title(39。subplot(2,1,2)F1=plot(f,abs(F0(1:512)))。)。聲卡的主要的作用之一是對聲音信息進(jìn)行錄制與回放，在這個過程中采樣的位數(shù)和采樣的頻率決定了聲音采集的質(zhì)量。預(yù)延時：自然混響聲陣的建立都會延遲一段時間，預(yù)延時即為模擬次效應(yīng)而設(shè)置。如果延遲時間在50ms以上時，其延遲音就清晰可辨，此時的處理效果才是回音。*預(yù)延時（Ini Dly），指主延時電路預(yù)延時時間調(diào)整。若x為一矩陣，則y是對矩陣的每一列向量進(jìn)行FFT。MATLAB提供求復(fù)數(shù)的幅值和相位函數(shù)：abs，angle，這些函數(shù)一般和FFT同時使用。該函數(shù)是利用給出的向量b和a，對x中的數(shù)據(jù)進(jìn)行濾波，結(jié)果放入向量y。39。原始信號時域波形39。plot(angle(Y))。title(39。subplot(2,2,4)。z=[zeros(3000,1)。subplot(2,2,3)。延時后相頻39。zeros(3000,1)]。subplot(2,2,3)?；祉懙南囝l39。zeros(3000,1)]。%求設(shè)計的濾波器頻譜subplot(3,2,1)。單回聲濾波器相頻響應(yīng)39。plot(abs(Y2))。)。a=。subplot(3,2,1)。多重回聲濾波器相頻響應(yīng)39。plot(abs(Y2))。)。a=。plot(abs(h))。)。title(39。sound(y2,fs,bits)。zeros(3000,1)]。subplot(3,2,1)。全通濾波器相頻響應(yīng)39。plot(abs(YY2))。)。T=1/Fs。fc1=Fs/10。xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t)。ylabel(39。)。ylabel(39。rp=。%低通的頻率響應(yīng)m=abs(H1)。)。)。fpl=440。[N,wp0]=ellipord(wp,ws,rp,rs)。grid on。H239。t39。wp=2*fp/Fs。y3t=filter(B,A,st)。高通濾波損耗函數(shù)曲線39。subplot(2,1,2)。ylabel(39。如果采用脈沖響不變法，邊界頻率的轉(zhuǎn)換關(guān)系為:如果采用雙線性變換法，邊界頻率的轉(zhuǎn)換關(guān)系為(3) 按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計模擬低通濾波器。阻帶最大衰減為：As=15。)。它們都是借助于 已經(jīng)成熟的現(xiàn)有低通濾波器原型進(jìn)行設(shè)計，即對數(shù)字低通數(shù)字濾波器，先將數(shù)字低通濾波器的技術(shù)指標(biāo)按希望的設(shè)計方法轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標(biāo)，再按指定的模擬低通濾波器的類型設(shè)計模擬濾波器H(s)，然后，將模擬濾波器的系統(tǒng)函數(shù)H(s)從s平面轉(zhuǎn)換到z平面，得到數(shù)字低通濾波器的系統(tǒng)函數(shù)H(z)；如所設(shè)計的數(shù)字濾波器為高通、帶通或帶阻濾波器，則可借助模擬濾波器的頻帶變換轉(zhuǎn)換為低通模擬濾波器。b=a(1:1000,2)。ylabel(39。原始頻譜圖39。圖通過導(dǎo)入的心電信號數(shù)據(jù)發(fā)現(xiàn)在其頻譜圖上的0~20H z和80~100Hz之間的幅值比較大，而在30~70Hz之間的幅值相對較小。plot(s(1:1000))。)。)。x1=a(:,2)。加入工頻干擾后的時域圖39。subplot(2,1,2)。ylabel(39。figure(4)。)。)。（6）濾波器的設(shè)計Butterworth型低通數(shù)字濾波器用wp和ws表示分別將通帶，阻帶截止頻率的角頻率表示，在分別計算階數(shù)n1和截止頻率Wn,再設(shè)計低通Butterworth型模擬濾波器，然后采用雙線性法將模擬濾波器系數(shù)變?yōu)閿?shù)字濾波器系數(shù)，畫出濾波器頻譜圖，調(diào)用filter實現(xiàn)對工頻干擾的濾波，用plot函數(shù)畫出濾除工頻干擾后的時域圖和濾除白噪聲后的頻譜圖。Omegas=2*fs*tan(ws/2)。freqz(bz,az,512,fs)。時間（s）39。plot(f1,abs(y3))。幅值(dB)39。f21=5。Omegas1=2*fs*tan(Ws1/2)。s39。y=filter(num,den,s)。時間（s）39。plot(f1,abs(s3))。幅值(dB)39。)。)。y5=fft(y4,1000)。)。拿到課程設(shè)計任務(wù)的時候，不知道從哪兒入手，選定題目后，上網(wǎng)查閱了相關(guān)資源，去圖書館搜集資料，確定了設(shè)計的基本方案，步驟，開始用軟件進(jìn)行實現(xiàn)，實現(xiàn)過程中出現(xiàn)了很多問題，上網(wǎng)查閱資料，在專業(yè)的論壇里面去看，下載MATLAB的視頻，慢慢的開始有了進(jìn)一步的理解。CALLBACK39。gui_Name39。gui_LayoutF39。else gui_mainf(gui_State, varargin{:})。)。原始信號波形39。頻率響應(yīng)圖39。fudu39。f=fs*(0:511)/N。)subplot(2,1,2) ,plot(x2)。)。ylabel(39。)。% Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved to be de