【正文】 形，通過多波形的分析，使學(xué)生理解時(shí)域信號(hào)好頻域信號(hào)的特性，從而進(jìn)一步加深學(xué)生對(duì)理論知識(shí)的理解?！拔⒄n”在闡述問題的方式上除了教師講解外，還結(jié)合動(dòng)畫、視頻等各種信息技術(shù)，使得對(duì)知識(shí)的講解過程更為直觀。[Bz,Az]=bilinear(b,a,Fs)。[n,Wc]=buttord(OmegaP,OmegaS,Rp,Rs,39。Wp1=(fp1/Fs)*2*pi。)。Ws=[Ws1,Ws2]。% 通帶截止頻率 fs2=265。high39。Wp=(fp/Fs)*2*pi。[b,a]=butter(n,Wc,39。s39。Wp=(fp/Fs)*2*pi。Ws=[Ws1,Ws2]。% 通帶截止頻率 fs2=60。Ws1=(fs1/Fs)*2*pi。注：要好好看實(shí)驗(yàn)中關(guān)于低通，高通，帶通，帶阻的設(shè)計(jì)代碼。% 頻率預(yù)畸 OmegaS=2*Fs*tan(Ws/2)。2）雙線性法濾除第三個(gè)信號(hào)： Fs=256。OmegaS=Ws*Fs。1）脈沖響應(yīng)不變法濾除第三個(gè)信號(hào)： Fs=256。最后可得（級(jí)聯(lián)型）：(s)=(s+)(s2++)五、脈沖響應(yīng)不變法（P177 ）156158頁脈沖響應(yīng)不變法的優(yōu)點(diǎn)：l 時(shí)域逼近。(2)阻帶截止頻率：Ωs=；阻帶最小衰減：As=16dB。式中x(n)和y(n)分別表示系統(tǒng)的輸入和輸出信號(hào)。x(n)=x(n)RN(n)=x((n))NRN(n)x(n)=IDFT[X(k)]=X(k)W 229。有限長(zhǎng)序列的離散傅里葉變換（DFT）：N1236。))set(hObject,39。)returns contents of edit12 as text % str2double(get(hObject,39。)。stem(n1,yn,39。)。stem(n1,yn,39。ylabel(39。n1=[0:length(yn)1]/length(yn)*2。w/π39。yn=20*log10(yn)。elseBt=wpws。String39。)。)。stem(n1,yn,39。)。stem(n1,yn,39。ylabel(39。n1=[0:length(yn)1]/length(yn)*2。w/π39。yn=20*log10(yn)。elseBt=wpws。String39。)。)。stem(n1,yn,39。)。stem(n1,yn,39。ylabel(39。n1=[0:length(yn)1]/length(yn)*2。w/π39。yn=20*log10(yn)。wp=2*fp/fre*pi。String39。相位39。)。|X(ejw)|39。)。)。title(39。ylabel(39。axes()。Bt=wswp。fs=str2num(get(,39。defaultUicontrolBackgroundColor39。endfunction edit11_Callback(hObject, eventdata, handles)function edit11_CreateF(hObject, eventdata, handles)if ispcset(hObject,39。,39。)。stem(n1,yn,39。)。stem(n1,yn,39。)。)。title(39。%求階數(shù)，截止頻率%[B,A]=butter(N,wc,39。rs=str2num(get(,39。String39。))。BackgroundColor39。white39。elseset(hObject,39。,get(0,39。)。title(39。)。%坐標(biāo)軸單位控制% title(39。,xn)。))。f(3)=str2num(get(,39。))。BackgroundColor39。white39。elseset(hObject,39。,get(0,39。))。BackgroundColor39。amp。, df_OpeningF,...39。每次課程設(shè)計(jì)都收獲頗多，而且每次都對(duì)自己的學(xué)習(xí)態(tài)度做一次調(diào)整，這個(gè)的作用確實(shí)是很大的，希望以后更加珍惜的課程設(shè)計(jì)的時(shí)間。對(duì)數(shù)字信號(hào)本身內(nèi)容的理解來說，全書的內(nèi)容確實(shí)是可以融合在這兩個(gè)課程設(shè)計(jì)題目中的，一個(gè)是DFT運(yùn)用，另一個(gè)是濾波器的設(shè)計(jì)和利用。幾天下來收獲是很多的，對(duì)數(shù)字信號(hào)處理有了一次全面的回顧，而且也看到了matlab的神秘面容。但是在一些函數(shù)操作行會(huì)受一些限制，比如函數(shù)參數(shù)的傳遞，目前來說還是采用全局變量才能完成。,wn)。%wdel為窗函數(shù)的過渡帶寬，h對(duì)應(yīng)不同窗函數(shù)的值 wn=boxcar(N+1)。濾波器的設(shè)計(jì)：根據(jù)以上技術(shù)指標(biāo)（通帶截止頻率、通帶最大衰減、阻帶截止頻率、阻帶最小衰減），設(shè)計(jì)數(shù)字濾波器，生成相應(yīng)的濾波器系數(shù)，并畫出對(duì)應(yīng)的濾波器幅頻特性與相頻特性。% 頻率1f2=50。在輸入信號(hào)確定后，要對(duì)其進(jìn)行頻譜分析，從而確定濾波器的參數(shù)和類型（低通、高通、帶通、帶阻），此模塊作用也就完成，將數(shù)據(jù)分別用全局變量傳遞給下面的模塊。特別在實(shí)現(xiàn)窗函數(shù)濾波器時(shí)，各個(gè)函數(shù)的主要區(qū)別是不同的頻率采樣，可以通過選擇結(jié)構(gòu)實(shí)現(xiàn)，簡(jiǎn)化程序。動(dòng)態(tài)輸入，即輸入由一系列頻率的正弦信號(hào)相加的組成的信號(hào)，需要經(jīng)過采樣的，注意在設(shè)置采樣頻率時(shí)一定要符合奈奎斯特準(zhǔn)則，提高采樣點(diǎn)數(shù)，增加頻譜分辨率。第二章 設(shè)計(jì)分析此題目的實(shí)現(xiàn)可分為三個(gè)某塊的設(shè)計(jì)實(shí)現(xiàn)：輸入信號(hào)模塊，設(shè)計(jì)濾波器模塊，濾波模塊。⑶ 濾波器設(shè)計(jì)—根據(jù)IIR/FIR 數(shù)字濾波器技術(shù)指標(biāo)設(shè)計(jì)濾波器，生成相應(yīng)的濾波器系數(shù)，并展示對(duì)應(yīng)的濾波器幅頻（衰減）、相頻特性。主要分為兩種：有限脈沖響應(yīng)FIR和無限脈沖響應(yīng)IIR。IIR濾波器的設(shè)計(jì)主要有兩類，一是借助于模擬濾波器設(shè)計(jì)進(jìn)行，二是直接在頻域或時(shí)域中進(jìn)行設(shè)計(jì)。① IIR DF：要求通過差分方程迭代實(shí)現(xiàn)濾波，未知初值置零處理； ② FIR DF：要求通過快速卷積實(shí)現(xiàn)濾波。先設(shè)計(jì)輸入模塊。第一步需要總結(jié)設(shè)計(jì)濾波器需要哪些參數(shù)，通過學(xué)習(xí)可以總結(jié)，所有濾波器的參數(shù)有四個(gè)：通帶截止頻率、阻帶截止頻率、通帶最大衰減、阻帶最小衰減。如果性能符合要求，則將系數(shù)傳遞給濾波模塊。頻率f3=200。% 采樣頻率Ts=1/fs。[N,Wc]=cheb1ord(wp,ws,rp,rs)。wn=hamming(N+1)。需要注意的是，窗函數(shù)對(duì)濾波參數(shù)的使用，需要通過運(yùn)算得到一窗函數(shù)階數(shù)，對(duì)通帶最大衰減無使用。其余和直接用函數(shù)創(chuàng)建GUI界面無太大區(qū)別。這次課程設(shè)計(jì)也使我認(rèn)識(shí)到了matlab的強(qiáng)大，或者進(jìn)一步說工具軟件有著你意想不到的功效，能節(jié)省你的時(shí)間，同時(shí)又能夠讓你從實(shí)踐上更深層次的認(rèn)識(shí)到所學(xué)知識(shí)的奇妙，而且可以同時(shí)明白了學(xué)習(xí)與實(shí)踐的相輔相成。對(duì)于課程設(shè)計(jì)中出現(xiàn)的問題，解決的方式有兩種，一是自己解決。gui_Name39。gui_LayoutF39。elsegui_mainf(gui_State, varargin{:})。)。,39。endfunction edit3_Callback(hObject, eventdata, handles)function edit3_CreateF(hObject, eventdata, handles)if ispcset(hObject,39。defaultUicontrolBackgroundColor39。BackgroundColor39。)。String39。f(4)=str2num(get(,39。))。%橫坐標(biāo)% axes()。xlabel(39。axes()。xlabel(39。,39。endfunction edit7_Callback(hObject, eventdata, handles)function edit7_CreateF(hObject, eventdata, handles)if ispcset(hObject,39。defaultUicontrolBackgroundColor39。BackgroundColor39。)。)。String39。wp=2*fp/fre。%巴特沃茲模擬低通濾波器系數(shù)%[Bz,Az]=bilinear(B,A,1)。xlabel(39。title(39。)。title(39。)。title(39。)。elseset(hObject,39。white39。userdata39。fre=str2num(get(,39。wc=(wp+ws)/2/pi。矩形窗的損耗函數(shù)39。axis([0,1,min(yn),max(yn)])。xlabel(39。%重疊相加法求濾波后序列% yn=abs(fft(yn2))。)。yn=angle(fft(yn2))。)。pushbutton1_Callback(hObject, eventdata, handles)。String39。N0=ceil(*pi/Bt)。plot(n1,yn)。20log|Hg(w)|39。三角窗的相頻特性39。axis([0,1,min(yn),max(yn)])。title(39。)。title(39。)。))。ws=2*fs/fre*pi。wc=(wp+ws)/2/pi。plot(n1,yn)。20log|Hg(w)|39。漢寧窗的相頻特性39。axis([0,1,min(yn),max(yn)])。title(39。)。title(39。)。))。ws=2*fs/fre*pi。wc=(wp+ws)/2/pi。plot(n1,yn)。20log|Hg(w)|39。哈明窗的相頻特性39。axis([0,1,min(yn),max(yn)])。title(39。)。title(39。)。amp。white39。x(n)WNx(n)=x((n))N236。Nk=0238。y(n)=(a+b)y(n－1)－aby(n－2)+x(n－2)+(a+b)x(n－1)+ab 式中, |a|解：(1)直接型結(jié)構(gòu)。 N=234。模擬頻率Ω和數(shù)字頻率ω之間呈線性關(guān)系ω=ΩT。% 阻帶截止頻率 Rp=1。)。% 阻帶截止頻率 Rp=1。)。% 阻帶截止頻率 fp2=60。Wp=[Wp1,Wp2]。Wp1=(fp1/Fs)*2*pi。p70+p180。OmegaS=Ws*Fs。（2）雙線性變換法OmegaP=2*Fs*tan(Wp/2)。[Bz,Az]=bilinear(b,a,Fs)。% 頻率預(yù)畸 OmegaS=2*Fs*tan(Ws/2)。)。Wp1=(fp1/Fs)*2*pi。[n,Wc]=buttord(OmegaP,OmegaS,Rp,Rs,39。% 通帶截止頻率 fs1=240。Ws2=(fs2/Fs)*2*pi。[b,a]=butter(2*n,Wc,39。針對(duì)該問題，本文提出了一種基于微課的教學(xué)方式，將課程內(nèi)容制作成微課，分享給學(xué)生，以提高學(xué)生的學(xué)習(xí)興趣。 教學(xué)方法方面?zhèn)鹘y(tǒng)教學(xué)多采用“灌輸式”的方法進(jìn)行授課，教師作為課堂主體，單向地向?qū)W生灌輸知識(shí)，教學(xué)的方式過于單調(diào)，不足以吸引學(xué)生的注意力，和學(xué)習(xí)興趣的提高。 雙向互動(dòng)式研究型學(xué)習(xí)針對(duì)課程部分重點(diǎn)難點(diǎn)教學(xué)內(nèi)容，提取出合適的討論課題，進(jìn)行定期或不定期的專題討論?！緟⒖嘉墨I(xiàn)】[1]曾偉梁，劉穎，鮑曼，徐宏艷， MATLAB 實(shí)踐數(shù)字信號(hào)處理課程的輔助教學(xué)[J].黑龍江生態(tài)工程職業(yè)學(xué)院學(xué)報(bào)，2014（1）.[2]高靜，王鳳文，舒冬梅.《數(shù)字信號(hào)處理》課程教學(xué)實(shí)踐與探索[J].科技創(chuàng)新導(dǎo)報(bào)，2011（4）：153154.[3]曾偉梁，劉穎，鮑曼，徐宏艷， 實(shí)踐數(shù)字信號(hào)處理課程的輔助教學(xué)[J].黑龍江生態(tài)工程職業(yè)學(xué)院學(xué)報(bào)，2014（1）.[4]劉伯紅，閻英，[J].實(shí)驗(yàn)室研究與探索，2012，31（12）：121123，150.[責(zé)任編輯：楊玉潔]第四篇：數(shù)字信號(hào)處理課程設(shè)計(jì)..課程設(shè)計(jì)報(bào)告課程名稱： 數(shù)字信號(hào)處理 課題名稱： 語音信號(hào)的處理與濾波姓 名： 學(xué) 號(hào)： 院 系： 專業(yè)班級(jí)： 指導(dǎo)教師： 完成日期： 2013年7月2日目錄第1部分 課程設(shè)計(jì)報(bào)告………………………………………3 一．設(shè)計(jì)目的……………………………………………3 二．設(shè)計(jì)內(nèi)容……………………………………………3 三．設(shè)計(jì)原理……………………………………………3 四．具體實(shí)現(xiàn)……………………………………………5 …………………………………5 ………………………8 …………………13 ………………………………19 ……………………19 ……………………………22 五． 結(jié)果分析……………………………………………27 第2部分 課程設(shè)計(jì)總結(jié)………………………………28 一． 參考文獻(xiàn)……………………………………………28第1部分 課程設(shè)計(jì)報(bào)告一．設(shè)計(jì)目的綜合運(yùn)用本課程的理論知識(shí)進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為工具進(jìn)行實(shí)現(xiàn)，從而復(fù)習(xí)鞏固課堂所學(xué)的理論知識(shí)，提高對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)對(duì)數(shù)字信號(hào)的處理。n=165。xa(t)DT(jW)171。如果：236。2Wh 或 fs179。%暫停命令 fprintf(39。)。)。%錄制信號(hào)的時(shí)域波形 title(39。%快速傅立葉變換 figure(2)。)。%加網(wǎng)格 原始信號(hào)的時(shí)域與頻譜圖針對(duì)實(shí)際發(fā)出聲音落后錄制動(dòng)作半拍的現(xiàn)象，如何拔除對(duì)無效點(diǎn)的采樣的問題: 出現(xiàn)這種現(xiàn)象的原因主要是錄音開始時(shí)，人的反應(yīng)慢了半拍，導(dǎo)致出現(xiàn)了一些無效點(diǎn)，而后而出現(xiàn)的無效的點(diǎn)，主要是已經(jīng)沒有聲音的動(dòng)作，先讀取聲音出來，將原始語音信號(hào)時(shí)域波形圖畫出來，根據(jù)己得到的信號(hào)，可以在第二次讀取聲音的后面設(shè)定采樣點(diǎn)，取好有效點(diǎn)，畫出濾除無效點(diǎn)后的語音信號(hào)時(shí)域波形圖，對(duì)比可以看出。)。%畫網(wǎng)格 去除無效點(diǎn)MATLAB程序如下。Fs1=1。)。xlabel(39。%高通濾波fp=4800。Fs1=1。,39。巴特沃斯高通濾波器39。grid on。ws=2*pi*fs/Fs。[B,A]=butter(N,wc,39。巴特沃斯帶通濾波器39。grid on。f=[(n11)*(2*pi/N)]/pi。幅值/A39。)。原始信號(hào)的頻譜圖39。wn=(ws+wc)/2。%將濾波后的語音信號(hào)進(jìn)行快速傅立葉變換 figure(4)。)。頻率w/pi39。%標(biāo)題 grid。n=ceil(*pi/width)。%信號(hào)的時(shí)域波形 xlabel(39。)。%繪制模值 axis([0 1 0 50])。高通濾波器濾波后信號(hào)的頻譜圖39。ws1=2*pi*1200/fs。band39