【正文】 以后有這樣的機(jī)會(huì)一定會(huì)更加的很好利用，它不僅可以提高學(xué)習(xí)的針對性而且可以很好的鍛煉動(dòng)手能力以及自己的邏輯設(shè)計(jì)能力和處理問題的能力，希望在以后的學(xué)習(xí)生活中會(huì)有更多的機(jī)會(huì)來加強(qiáng)這方面的能力。這次課程設(shè)計(jì)對我各方面的能力有了很大的提高，對我以后的工作、實(shí)踐都有很大的幫助。)。subplot(3,2,6)。imshow(b1)。title(39。Prewitt邊緣檢測39。)。subplot(3,2,2)。imshow(b)。b5=edge(b,39。)。prewitt39。b1=edge(b,39。4）：a=imread(39。b=rgb2gray(a)。)。)。G39。subplot(3,2,4),imshow(ar),title(39。)。F:39。原彩色圖像39。a=imread(39。subplot(313)bz3=bz。subplot(312)bz2=bz。subplot(311)bz1=bz。b=[ 1]。Rn=[ones(1,N1)zeros(1,101N)]。ap=1。wp2=[,]。5）：wp=200。z1=fft(y,256)。k=0::100。單位脈沖響應(yīng)39。xlabel(39。ylabel(39。k39。b=[ 1]。magHdB=20*log10(magH)。[H,w]=freqz(sn,1,200)。stem(n,sn)。Rn=[ones(1,N1)zeros(1,101N)]。subplot(212)。magH=abs(H)。title(39。Rn=[ones(1,N1)zeros(1, 101N)]。ylabel(39。)axis([0,1,1,1])。)xlabel(39。F=X(1:N/2+1)。x=sin(2*pi*10*t)。相頻特性39。)。title(39。fill39。title(39。a=square(2*pi*10*t,50)。grid on。)subplot(313)stem(t*100,angle(b)/100,39。axis([0 50 0 ])。)b=fft(a)。grid on。相頻特性39。Hz39。)subplot(313)stem(t*1000,angle(b)/1000,39。axis([10 50 0 ])。)。b=fft(a)。)。plot(t,a)。title(39。y=square(2*pi*10*t,50)。axis([ ])。周期性三角波： t=::。axis([0 22 2])。t=0::2。整周泣截?。号c周期函數(shù)一樣，周期信號是每隔時(shí)間T以后，信號重復(fù)出現(xiàn)。有一個(gè)定理：對信號進(jìn)行時(shí)域fs的采樣，信號的頻譜就會(huì)在頻域以fs為周期重復(fù)。三、調(diào)試及結(jié)果分析（截圖）：（一）：1）頻率為10Hz的周期性正弦波：頻率為10Hz的周期性三角波：頻率為10Hz的周期性方波：2）采樣頻率不能過低，必須fs=2fm，即采樣頻率必須大于或等于最高截止頻率的二倍（對采樣頻率的要求，即采樣頻率要足夠大，采樣的值要足夠多，才能不失真的恢復(fù)原信號）。（二）：已知三個(gè)信號aip(n)，經(jīng)調(diào)制產(chǎn)生信號s(n)=229。說明選擇理由。一、課程設(shè)計(jì)目的：1．全面復(fù)習(xí)課程所學(xué)理論知識(shí)，鞏固所學(xué)知識(shí)重點(diǎn)和難點(diǎn)，將理論與實(shí)踐很好地結(jié)合起來。而使數(shù)字信號處理從理論走向?qū)嵱玫氖强焖俑盗⑷~變換(FFT)，F(xiàn)FT的出現(xiàn)大大減少了DFT的運(yùn)算量，使實(shí)時(shí)的數(shù)字信號處理成為可能、極大促進(jìn)了該學(xué)科的發(fā)展。而數(shù)字信號處理的輸出經(jīng)常也要變換到模擬域，這是通過數(shù)模轉(zhuǎn)換器實(shí)現(xiàn)的。成功正在盡頭。)。subplot(2,1,2)。subplot(2,1,1)plot(20*log10(abs(yx)))。F:39。y=wavread(39。%加X軸說明（/后表示轉(zhuǎn)義字符）ylabel(39。)。%窗口分成2行1列繪圖區(qū),第2個(gè)區(qū)為當(dāng)前活動(dòng)區(qū) plot(w/pi,20*log10(abs(h)))。%加X軸說明（/后表示轉(zhuǎn)義字符）ylabel(39。)。%窗口分成2行1列繪圖區(qū),第1個(gè)區(qū)為當(dāng)前活動(dòng)區(qū) plot(w/pi,abs(h))。%最小阻帶衰減Rs [N,Wn]=buttord(wp,ws,Rp,Rs)。3．濾波器的設(shè)計(jì)：針對電話信道我們設(shè)計(jì)一個(gè)IIR低通濾波器對我們錄制的語音進(jìn)行濾波,具體程序如下: %低通濾波器設(shè)計(jì)wp=2*1200/7000。)。 angle(d)39。%繪出語音信號的相位譜 xlabel(39。)。幅值39。%繪出語音信號的幅度譜 xlabel(39。)。%y軸坐標(biāo) grid。)。原始聲音信號時(shí)域圖形39。F:數(shù)字信號處理數(shù)字信號課程設(shè)計(jì) 39。)。F:數(shù)字信號處理數(shù)字信號課程設(shè)計(jì) 39。%調(diào)用錄制自己的聲音信號函數(shù)fprintf(39。)。,t)。%抽樣頻率 channel=1。N1I0(b)其中：β是一個(gè)可選參數(shù)，用來選擇主瓣寬度和旁瓣衰減之間的交換關(guān)系，一般說來，β越大，過渡帶越寬，阻帶越小衰減也越大。設(shè)理想濾波器Hd(ejw)的單位脈沖響應(yīng)為hd(n)。但在本設(shè)計(jì)中主要是對已有的信號進(jìn)行處理，所以選用wavread函數(shù)。不過wavread函數(shù)只能讀取計(jì)算機(jī)中已經(jīng)有的“.wav” 文件，因此必須事先使用windows的錄音機(jī)等設(shè)備把要用的語音錄制并存為“.wav” 文件格式。三、基本步驟：1．語音信號的采集熟悉并掌握MATLAB中有關(guān)聲音（wave）錄制、播放、存儲(chǔ)和讀取的函數(shù)，在MATLAB環(huán)境中，使用聲音相關(guān)函數(shù)錄制3秒左右自己的聲音。2．在MATLAB環(huán)境中，使用聲音相關(guān)函數(shù)錄制2秒左右自己的聲音，抽樣率是8000Hz/s。：（35個(gè)）一設(shè)計(jì)目的與要求………………………………………頁碼 二總體設(shè)計(jì)方案…………………………………………頁碼 三設(shè)計(jì)原理、結(jié)果與仿真分析…………………………頁碼 四結(jié)論……………………………………………………頁碼 五心得體會(huì)………………………………………………頁碼 參考文獻(xiàn)…………………………………………………頁碼 附錄………………………………………………………頁碼特別注意：，所有的表也要有編號和表名； ，公式編號右對齊。設(shè)計(jì)四IIR數(shù)字濾波器設(shè)計(jì) 主要要求：、帶通、帶阻和高通四種數(shù)字濾波器（FIR數(shù)字濾波器的指標(biāo)自行選擇）；，并分別闡述上述四類濾波器的設(shè)計(jì)原理、設(shè)計(jì)步驟，并給出所編寫的相應(yīng)的m程序；（頻率區(qū)間[O，π]上的幅頻響應(yīng)特性曲線），并分析各個(gè)濾波器的特點(diǎn)，將結(jié)論寫在報(bào)告中。3根據(jù)指標(biāo)（低通FIR濾波器的指標(biāo)自行選擇）要求選擇窗函數(shù)的形狀與長度N。(3)采集數(shù)據(jù)x(n)長度取為M點(diǎn)（注意不是補(bǔ)零至M），編寫程序計(jì)算出M點(diǎn)采集數(shù)據(jù)x(n)的的頻譜并畫出相應(yīng)的幅頻圖。(2)將(1)中x(n)補(bǔ)零加長至M點(diǎn)，長度M自己選（,為了比較補(bǔ)零長短的影響，M可以取兩次值，一次取較小的整數(shù)，一次取較大的整數(shù)），編寫程序計(jì)算x(n)的M點(diǎn)DFT, 畫出時(shí)域序列圖和兩次補(bǔ)零后相應(yīng)的DFT幅頻圖。在這次的課程設(shè)計(jì)中讓我體會(huì)最深的是：知識(shí)來不得半點(diǎn)的馬虎。通過查閱資料和請教同學(xué)收獲了很多以前不懂的理論知識(shí)。，噪聲與原信號明顯區(qū)分,但通過低通濾波器后，噪聲沒有濾除，信號產(chǎn)生失真。當(dāng)通過高通濾波器后，把低頻濾除，聲音變得比較就尖銳。)加 窗[x,fs,bits]=wavread(39。%加網(wǎng)格wavplay(ys,8000)。ylabel(39。)。%網(wǎng)格 subplot(3,2,6)。%加標(biāo)題 ylabel(39。plot(ys)。%DTFT變換ys=filter(num,den,yy)。%最大衰減 Rs=100。 設(shè)計(jì)一個(gè)合適的濾波器將噪聲濾除 fs=18000。ylabel(39。)。%快速傅立葉變換 subplot(3,2,4)。幅值/A39。xlabel(39。,2)。%頻帶寬度 figure(5)。%播放疊加后信號 N=length(yy)。 C: 39。[x,fs,bits]=wavread(39。title(39。plot(f,abs(y0(n1)))。幅值/A39。噪聲信號的時(shí)域波形39。%頻帶寬度 figure(5)。%播放噪聲信號 N=length(x1)。%存儲(chǔ)語音信號 對噪聲進(jìn)行頻譜的分析[x1,fs,bits]=wavread(39。%暫停命令 wavplay(x,Fs)。%文字提示 fprintf(39。)。錄音中......39。按Enter鍵錄音%ds39。%顯示縱坐標(biāo)表示的意義 grid。)。原始男生信號的頻譜圖39。%加網(wǎng)格y0=fft(y)。%加標(biāo)題 ylabel(39。plot(y)。%返回采樣點(diǎn)數(shù) df=fs/N。 C: 39。幅值 39。%加標(biāo)題 xlabel(39。plot(f,abs(y0(n1)))。)。原始女生信號的時(shí)域波形39。%頻帶寬度 figure(1)。%播放原始信號 N=length(x)。%存儲(chǔ)語音信號 分析男女聲音的頻譜x=wavread(39。%暫停命令 wavplay(x,Fs)。%文字提示 fprintf(39。)。錄音中......39。按Enter鍵錄音%ds39。sound(xdaitong, Fs,bits)。%播放濾波后信號 圖形如下：原始信號的時(shí)域波形幅值/A01012x 10原始信號的頻譜圖34幅值/分析：三個(gè)濾波器濾波后的聲音與原來的聲音都發(fā)生了變化。帶通濾波器濾波后信號的頻譜圖39。ylabel(39。%繪制模值 axis([0 1 0 50])。%顯示縱坐標(biāo)表示的意義 grid。)。帶通濾波后信號的時(shí)域波形39。%將信號送入高通濾波器濾波 figure(4)。band39。wn2=(ws2+wc2)/2。ws1=2*pi*1200/fs。%播放濾波后信號3．3帶通濾波器 fs=8000。高通濾波器濾波后信號的頻譜圖39。ylabel(39。%繪制模值 axis([0 1 0 50])。%顯示縱坐標(biāo)的表示意義 grid。)。幅值/A39。%信號的時(shí)域波形 xlabel(39。[h1,w]=freqz(h)。n=ceil(*pi/width)。ws=2*5000/fs。%標(biāo)題 grid。)。頻率w/pi39。%顯示標(biāo)表示的意義 grid。)。低通濾波后信號的時(shí)域波形39。%將濾波后的語音信號進(jìn)行快速傅立葉變換 figure(4)。,kaiser(n+1,beta))。wn=(ws+wc)/2。beta=。原始信號的頻譜圖39。%顯示橫坐標(biāo)表示的意義 ylabel(39。)。%快速傅立葉變換 subplot(4,2,3)。幅值/A39。%錄制信號的時(shí)域波形title(39。f=[(n11)*(2*pi/N)]/pi。%播放原始信號wavplay(x,fs)。grid on。)。巴特沃斯帶通濾波器39。figure(4)。[B,A]=butter(N,wc,39。was=2*tan(ws/2)。ws=2*pi*fs/Fs。Fs=8000。grid on。)。巴特沃斯高通濾波器39。figure(1)。,39。s39。Fs1=1。rs=100。%高通濾波fp=4800。耗損（dB）39。xlabel(39。plot(w,abs(h))。)。s39。Fs1=1。rp=1。%畫網(wǎng)格 去除無效點(diǎn)MATLAB程序如下。ylabel(39。)。, 7[4000,24000])。%加網(wǎng)格 原始信號的時(shí)域與頻譜圖針對實(shí)際發(fā)出聲音落后錄制動(dòng)作半拍的現(xiàn)象，如何拔除對無效點(diǎn)的采樣的問題: 出現(xiàn)這種現(xiàn)象的原因主要是錄音開始時(shí)，人的反應(yīng)慢了半拍，導(dǎo)致出現(xiàn)了一些無效點(diǎn)，而后而出現(xiàn)的無效的點(diǎn)，主要是已經(jīng)沒有聲音的動(dòng)作，先讀取聲音出來，將原始語音信號時(shí)域波形圖畫出來，根據(jù)己得到的信號，可以在第二次讀取聲音的后面設(shè)定采樣點(diǎn)，取好有效點(diǎn)，畫出濾除無效點(diǎn)后的語音信號時(shí)域波形圖，對比可以看出。%顯示縱坐標(biāo)表示的意義 title(39。)。原始信號的頻譜圖39。%快速傅立葉變換 figure(2)。幅值/A39。%錄制信號的時(shí)域波形 title(39。f=[(n11)*(2*pi/N)]/pi。)。pause。)。double39。%暫停命令 fprintf(39。%錄音時(shí)間 fprintf(39。2Wh 或 fs179。238。如果：236。Xa(jWjk)=229。xa(t)DT(jW)171。229。n=165。，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進(jìn)行濾波。第二篇：數(shù)字信號處理課程設(shè)計(jì)..課程設(shè)計(jì)報(bào)告課程名稱： 數(shù)字信號處理 課題名稱： 語音信號的處理與濾波姓 名： 學(xué) 號： 院 系： 專業(yè)班級： 指導(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í)，提高對所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)對數(shù)字信號的處理。另外，導(dǎo)師嚴(yán)謹(jǐn)求實(shí)的治學(xué)態(tài)度、一絲不茍的工作作風(fēng)和高尚的人格魅力，都給了學(xué)生很大感觸，使學(xué)生終生受益。因此，本次課程設(shè)計(jì)讓我受益匪淺。聲音放送可以加入數(shù)字回聲，