【正文】 ite(x,Fs,39。%返回采樣點數(shù) df=fs/N。subplot(2,1,1)。)。%顯示縱坐標的表示意義 grid。plot(f,abs(y0(n1)))。%加標題 xlabel(39。幅值 39。)。x=wavread(39。%畫出截取后的時域圖形 title(39。頻率39。)。fs=1200。wp=2*pi*fp/Fs。was=2*tan(ws/2)。[B,A]=butter(N,wc,39。figure(1)。巴特沃斯低通濾波器39。)。gridon。Fs=22050。ws=2*pi*fs/Fs。was=2*tan(ws/2)。[B,A]=butter(N,wc,39。)。plot(w,abs(h))。xlabel(39。耗損（dB）39。%帶通濾波fp=[1200,3000]。rs=100。Fs1=1。s39。)。plot(w,abs(h))。xlabel(39。耗損（dB）39。C:39。%返回采樣點數(shù) df=fs/N。subplot(4,2,1)。)。%顯示縱坐標的表示意義 grid。%原始信號的頻譜圖 title(39。頻率w/pi39。)。%加標題 grid。ws=2*pi*1200/fs。h=fir1(n,wn/pi,39。ys=fftfilt(h,x)。%在四行兩列的第二個窗口顯示圖形 plot(ys)。%加標題 xlabel(39。幅值/A39。%在四行兩列的第四個窗口顯示圖形 plot(f, abs(fftwave(n1)))。ylabel(39。低通濾波器濾波后信號的頻譜圖39。%播放濾波后信號3．2高通濾波器濾波 fs=8000。width=wswc。high39。%將信號送入高通濾波器濾波 subplot(4,2,5)。)。%顯示標表示的意義 title(39。幅值/A39。%將濾波后的語音信號進行快速傅立葉變換 subplot(4,2,7)。頻率w/pi39。)。%標題 grid。wc1=2*pi*1000/fs。width=ws1wc1。n=ceil(*pi)。[h1,w]=freqz(h)。%在四行兩列的第六個窗口顯示圖形 plot(ys1)。%加標題 xlabel(39。幅值/A39。%對濾波后的信號進行快速傅立葉變換 subplot(4,2,8)。頻率w/pi39。)。%加標題 grid。%在Matlab中，函數(shù)sound可以對聲音進行回放，其調(diào)用格式: sound(xhigh, Fs,bits)。%抽樣頻率 time=3。%文字提示 pause。x=wavrecord(time*Fs,Fs,39。錄音結(jié)束39。)。C:39。)。%采樣間隔 n1=1:N/2。plot(x)。%加標題 ylabel(39。%加網(wǎng)格y0=fft(x)。原始女生信號的頻譜圖39。)。%顯示縱坐標表示的意義 grid。% 對語音信號進行采樣wavplay(y,fs)。f=[(n11)*(2*pi/N)]/pi。原始男生信號的時域波形39。)。%在四行兩列的第三個窗口顯示圖形 plot(f,abs(y0(n1)))。%加標題 xlabel(39。幅值 39。%抽樣頻率 time=3。%文字提示 pause。x=wavrecord(time*Fs,Fs,39。錄音結(jié)束39。)。C:39。)。%采樣間隔n1=1:N/2。plot(x1)。grid。y0=fft(x1)。幅值39。)。)。%對噪聲信號進行采樣yy=x+x1。%采樣間隔 n1=1:N/2。plot(yy,39。疊加信號的時域波形39。)。grid。%疊加信號的頻譜圖 title(39。頻率w/pi39。)。%通帶截至頻率 Ws=2*2000/fs。%buttord函數(shù)（n為階數(shù)，Wn為截至頻率）[num,den]=butter(N,Wn)。%將濾波后的語音信號進行快速傅立葉變換 figure(5)。低通濾波后信號的時域波形39。)。%繪制模值 title(39。頻率w/pi39。)。圖形如下：噪聲信號的時域波形1100噪聲信號的頻譜圖幅值/ 10疊加信號的時域波形10101時間/t2200幅值/db34幅值/，常會出現(xiàn)實際發(fā)出聲音落后錄制動作半拍，可在[x,fs,bits]=wavread(39。[100 10000])，窗的長度可根據(jù)需要定義。：時域圖的振幅大小與性別無關(guān)，只與說話人音量大小有關(guān)，音量越大，振幅越大。第2部分 課程設(shè)計總結(jié)通過本次課程設(shè)計，使我們對數(shù)字信號處理相關(guān)知識有了更深刻的理解，尤其是對各種濾波器的設(shè)計。整個過程中學到了很多新的知識，特別是對Matlab的使用終于有些了解。八．參考文獻[1]數(shù)字信號處理教程(第三版)程佩青 清華大學出版社 [2]MATLAB信號處理 劉波 文忠 電子工業(yè)出版社 [3] 王宏 清華大學出版社[4]MATLAB基礎(chǔ)與編程入門 張威 西安電子科技大學出版社[5] 數(shù)字信號處理及其MATLAB實驗 趙紅怡 張常 化學工業(yè)出版社[6]MATLAB信號處理詳解 陳亞勇等 人民郵電出版社 [7] 數(shù)字信號處理錢同惠 機械工業(yè)出版社 29第三篇：數(shù)字信號處理課程設(shè)計參考題目數(shù)字信號處理課程設(shè)計資料使用MATLAB（或其他開發(fā)工具）編程實現(xiàn)下述內(nèi)容并寫出課程設(shè)計報告。(1)采集數(shù)據(jù)x(n)長度取N=16點，編寫程序計算出x(n)的16點DFTX(k),并畫出相應(yīng)的幅頻圖。補零DFT能否提高信號的頻譜分辨率，說明提高頻譜物理分辨率與頻譜頻域分辨率的措施各是什么？設(shè)計二用窗函數(shù)法設(shè)計FIR數(shù)字低通濾波器 主要要求：,在MATLAB命令窗下瀏覽各種窗函數(shù)，繪出（或打?。└鞣N窗函數(shù)圖。，將結(jié)論寫在報告中。二、課程設(shè)計撰寫具體要求 ；（各種理論推導和計算）、系統(tǒng)功能結(jié)構(gòu)及功能說明，并列出相應(yīng)重要的MATLAB程序； ,并做出說明；，歸納得出結(jié)論，并分析設(shè)計中遇到的問題及解決思路和方法； ； ；（放入課程設(shè)計報告冊附錄中）。一、課程設(shè)計目的：綜合運用本課程的理論知識進行頻譜分析以及濾波器設(shè)計，通過理論推導得出相應(yīng)結(jié)論，并利用MATLAB作為工具進行實現(xiàn)，從而復習鞏固課堂所學的理論知識，提高對所學知識的綜合應(yīng)用能力，并從實踐上初步實現(xiàn)對數(shù)字信號的處理。4．針對電話信道（最高3500Hz），設(shè)計一個FIR或IIR濾波器進行濾波，把抽樣率轉(zhuǎn)變?yōu)?000Hz/s，并進行頻譜分析，得到幅度和相位譜。四、相關(guān)原理知識錄音原理采用Matlab中的兩個函數(shù)可以用于實現(xiàn)錄音功能，一個是waverecord函數(shù)調(diào)用windows音頻設(shè)備進行錄音，主要用于實時語音采集。另外語音采集的頻率和精度不能改變。濾波器的設(shè)計原理及設(shè)計方法FIR濾波器的設(shè)計問題在于尋求一系統(tǒng)函數(shù)H(z)，使其頻率響應(yīng)H(ejw)逼近濾波器要求的理想頻率響應(yīng)Hd(ejw)，其對應(yīng)的單位脈沖響應(yīng)hd(n)。（3）漢寧(Hanning)窗，又稱升余弦窗12npw(n)=[1cos()]RN(n)2N1其頻率響應(yīng)和幅度響應(yīng)分別為：j(2p2pjwW(e)={(w)+[WR(w)+WR(w+)]}eN1N1=W(w)ejwaN1)w2W(w)=(w)+[WR(w2p2p)+WR(w+)]N1N1（4）漢明(Hamming)窗，又稱改進的升余弦窗2npw(n)=[()]RN(n)N1其幅度響應(yīng)為：W(w)=(w)+[WR(w2p2p)+WR(w+)] N1N1（5）布萊克曼(Blankman)窗，又稱二階升余弦窗2np4npw(n)=[()+()]RN(n)N1N12p2pW(w)=(w)+[WR(w)+WR(w+)]N1N1其幅度響應(yīng)為：4p4p+[WR(w)+WR(w+)]N1N1（6）凱澤(Kaiser)窗I0(b1[12n/(N1)]2)w(n)=,0163。)是第一類修正零階貝塞爾函數(shù)。%設(shè)定錄音時間為3秒鐘fprintf(39。%暫停命令fprintf(39。double39。)。%存儲錄音信號的文件夾fprintf(39。%作用是暫停wavplay(x,fs)。%讀取原有聲音文件及取點范圍plot(y)。%標題 xlabel(39。幅值39。%．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．取8000點頻譜分析y=wavread(39。%對語音信號做8000點的FFT變換 subplot(2,2,1)。)。%Y軸幅值title(39。%添加網(wǎng)格 subplot(2,2,2)。)。%Y軸angle(d)title(39。%添加網(wǎng)格觀察圖可以看出取8000點和取16000點所得的幅度譜和相位譜的圖形基本上是相同的；不同之處只有8000點圖形線比16000點的圖形線梳松。%歸一化阻帶數(shù)字截止頻率ws Rp=。%傳輸分子和分母的系數(shù) [h,w]=freqz(num,den)。%添加網(wǎng)格 title(39。omega/pi39。)。%添加網(wǎng)格 title(39。omega/pi39。)。)。%讀取濾波后語音信號 yx=fft(y)。濾波前的相對幅度譜39。title(39。不過課程設(shè)計的過程中也遇到了各種困難，但我還是堅持下來了。在本次課程設(shè)計中主要以數(shù)字信號處理來解決問題。數(shù)字信號處理技術(shù)及設(shè)備具有靈活、精確、抗干擾強、設(shè)備尺寸小、造價低、速度快等突出優(yōu)點，這些都是模擬信號處理技術(shù)與設(shè)備所無法比擬的。它在數(shù)學類科技應(yīng)用軟件中在數(shù)值計算方面首屈一指。二、課程設(shè)計題目：（一）：1）生成信號發(fā)生器：能產(chǎn)生頻率（或基頻）為10Hz的周期性正弦波、三角波和方波信號。分析說明選擇理由。將此已調(diào)信號通過信道傳輸，描述該信道的差分方程為y(n)(n1)+(n2)(n3)+(n4)=(n)(n1)+(n2)(n3)+x(n4)得到接收信號y(n)=s(n)*h(n)1）分析Hanning信號p(n)的時域與頻域特性 2）分析已調(diào)信號s(n)的時域與頻域特性 3）分析系統(tǒng)的單位脈沖響應(yīng)h(n)4）分析接收信號y(n)的頻譜5）設(shè)計帶通濾波器從接收信號y(n)中還原出三個已調(diào)信號。3）截取數(shù)據(jù)長度（即信號長度）N=T*fs=2 4）頻率為10Hz的周期性正弦波的幅頻特性和相頻特性：頻率為10Hz的周期性三角波的幅頻特性和相頻特性：頻率為10Hz的周期性方波的幅頻特性和相頻特性：5）頻譜混疊：一個信號的最高頻率比如是fmax,那它的頻譜就是在fmax~fmax之間有值。如果fs頻譜泄露：對于頻率為fs的正弦序列，它的頻譜應(yīng)該只是在fs處有離散譜。N取100時：N取200時：（二）：1）Hanning信號p(n)的時域與頻域特性：2）已調(diào)信號s(n)的時域與頻域特性：3）系統(tǒng)的單位脈沖響應(yīng)h(n)：4）接收信號y(n)的頻譜：5）（設(shè)計帶通濾波器）從接收信號y(n)中還原出三個已調(diào)信號：（三）：1）讀入一幅彩色圖像：2）將彩色圖像進行三原色分解，分解出R、G、B分量，并用圖像顯示出來：3）將彩色圖像灰度化，轉(zhuǎn)換為灰度圖像并顯示：4）對灰度圖像用幾種典型的邊緣檢測算子進行邊緣檢測，顯示檢測出的邊緣：四、源程序：（一）： 1）：周期性正弦波： fs=10。plot(t,y)。正弦信號39。plot(t,y)。三角波信號39。grid on。)4）：周期性正弦波： t=0::。title(39。T39。stem(t*1000,abs(b)/1000,39。Hz39。title(39。)。axis([0 ])grid on。a=sawtooth(10*pi*t,)。title(39。fill39。title(39。)。相頻特性39。grid on。)b=fft(a)。axis([0 50 0 4])。)subplot(313)stem(t*100,angle(b)/100,39。grid on。T=1。f=1/dt。subplot(2,1,1)plot(t,x)title(39。)ylabel(39。Frequency39。)（二）： 1）： N=100。subplot(211)。)。length(H)。2）：N=100。sn=2*pn.*(cos(pi*n/4)+4*pn.*cos(pi*n/2)+6*pn.*cos(3*pi*n/4))。已調(diào)信號39。length(w)。plot(w/pi,magHdB)。hn=impz(b,a,n)。f39。h(n)39。)。4）：a=[ ]。s=p.*(cos(pi*k/4)+2*cos(pi*k/2)+3*cos(3*pi*k/4))。axis([30 30 0 120])。wp1=[,]。wp3=[,]。N=100。sn=2*pn.*(cos(pi*n/4)+4*pn.*cos(pi*n/2)+6*pn.*cos(3*pi*n/4))。[n,wc]=cheb1ord(wp,ws,ap,as)。x1=filter(bz1,az1,yn)。x2=filter(bz2,az2,yn)。x3=filter(bz3,az3,yn)。)。2）：clc。subplot(3,2,1),imshow(a),title(39。ag=a(:,:,2)。)。subplot(3,2,6),imshow(ab),title(39。a=imread(39。原彩色圖像39?；叶葓D像39。)。)。b3=edge(b,39。log39。)?；叶葓D像39。title(39。imshow(b2)。subplot(3,2,4)。)。LoG邊緣檢測39。title(39。我進一步了解到凡事都需要耐心，細心仔細是成功的重要保證之一。還有就是有些不定參數(shù)存在時，可先取定值，用于調(diào)試，這樣可以節(jié)約調(diào)試時間，從而提高效率