資源描述：
《數(shù)字信號的處理實驗報告材料》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、實用標(biāo)準(zhǔn)文案長春理工大學(xué)電工電子實驗教學(xué)中心學(xué)生實驗報告2014——2015學(xué)年第一學(xué)期實驗課程數(shù)字信號處理實驗實驗地點東1教學(xué)樓414實驗室學(xué)院電子信息工程學(xué)院專業(yè)學(xué)號姓名精彩文檔實用標(biāo)準(zhǔn)文案實驗項目音頻信號采集、分析及處理實驗時間實驗臺號操作成績報告成績一、實驗?zāi)康?.以音頻信號為例，熟悉模擬信號數(shù)字處理過程，進一步理解數(shù)字信號處理概念。2.掌握運用Matlab實現(xiàn)對音頻信號的時頻分析方法；3.初步掌握數(shù)字音頻信號合成的方法。4.掌握運用Matlab設(shè)計IIR和FIR濾波系統(tǒng)的方法；5.掌握運用Matlab實
2、現(xiàn)對加噪音頻信號進行去噪濾波的方法。鍛煉學(xué)生運用所學(xué)知識獨立分析問題解決問題的能力，培養(yǎng)學(xué)生創(chuàng)新能力。二、實驗儀器計算機、Matlab軟件三、實驗性質(zhì)綜合分析、設(shè)計性實驗四、實驗任務(wù)?實驗內(nèi)容二：任意音頻信號的時域和頻域分析及數(shù)字濾波器設(shè)計1.音頻信號的采樣音頻信號的采集可以通過Windows自帶的錄音機也可以用專用的錄制軟件錄制一段音頻信號（盡量保證無噪音、干擾?。?，也可以直接復(fù)制一段音頻信號（時間為1s），但必須保證音頻信號保存為.wav的文件。2.音頻信號時域采樣對選取的音頻信號分別按照欠采樣率、過采樣率及
3、合理采樣率進行時域采樣，對采樣之后的音頻信號進行回放，對比回放結(jié)果，進行討論分析。3.音頻信號的時域、頻域分析運用Matlab軟件實現(xiàn)對音頻信號的打開操作以及時域分析和頻域分析，并畫出相應(yīng)的圖形，打印在實驗報告中。4.音頻信號的頻域分析對音頻信號按照過采樣、欠采樣及合理采樣進行頻域采樣并進行回放，對比結(jié)果，進行討論和分析。5.引入干擾信號基本要求：在原有的音頻信號上，疊加一個頻率為100KHz的正弦波干擾信號（幅度自定，可根據(jù)音頻信號幅度情況而定）。發(fā)揮：在原有音頻信號上，疊加其他復(fù)雜噪聲。6.數(shù)字濾波器設(shè)計精彩
4、文檔實用標(biāo)準(zhǔn)文案基本要求：運用Matlab實現(xiàn)IIR或FIR濾波系統(tǒng)，要求加入干擾的音頻信號經(jīng)過濾波系統(tǒng)后，能夠濾除100KHz的干擾信號，同時保留原有的音頻信號，要求繪制出濾波系統(tǒng)的頻譜圖，同時繪制出濾波前后音頻信號的時域、頻域?qū)Ρ葓D。發(fā)揮：對加有復(fù)雜噪聲的音頻信號進行濾波按照自己所加噪聲特點，確定濾波器類型及濾波器參數(shù)并設(shè)計該濾波器，利用自己所設(shè)計的濾波器對加有復(fù)雜噪聲的音頻信號進行濾波，回放濾波前后音頻信號，對比結(jié)果進行分析。%音頻信號的采集及其頻譜分析clearall;closeall;[y,fs,bit
5、s]=wavread('C:UsersAdministratorDesktopasd3Windowsll.wav');%讀取音頻文件sound(y,fs)%播放音頻文件n=length(y);%選取變換的點數(shù)取ny_p=fft(y,n);%對n進行傅里葉變換到頻域f=fs*(0:n/2-1)/n;%對應(yīng)點subplot(2,1,1);%第一組第一張圖片的頻率figure(1)%創(chuàng)建第一個繪圖窗口plot(y)%繪制音頻信號時域波形title('音頻信號時域波形');%圖片標(biāo)題xlabel('時間(t)'
6、);%x軸為時間軸ylabel('幅值(V)');%y軸為幅值subplot(2,1,2);%第一組第二張圖片plot(f,abs(y_p(1:n/2)));%繪制音頻信號頻譜圖title('音頻信號頻譜圖');%圖片標(biāo)題xlabel('頻率(Hz)');%x軸為頻率ylabel('頻率幅值(V)');%頻率幅值pause(1)%暫停一秒%噪聲合成處理t=0:1/fs:(n-1)/fs;%產(chǎn)生序列noise=0.05*cos(2*pi*100000*t);%噪聲信號Noise=[noise;noise];%創(chuàng)建矩
7、陣Noise=Noise';%對該矩陣進行轉(zhuǎn)置y_z=y+Noise;%加噪聲處理，兩個信號合成sound(y_z,fs);%播放合成后的音頻信號n=length(y);%選取長度y_zp=fft(y_z,n);%對n進行傅里葉變換到頻域f=fs*(0:n/2-1)/n;%對應(yīng)點的頻率精彩文檔實用標(biāo)準(zhǔn)文案figure(2)%創(chuàng)建第二個繪圖窗口subplot(2,1,1);%第二組第一張圖片plot(y_z);%繪制加噪聲處理音頻信號時域波形title('加噪聲處理音頻信號時域波形');%圖片標(biāo)題xlabel('時
8、間(t)');%x軸為時間軸ylabel('幅值(V)');%y軸為幅值subplot(2,1,2);%第二組第二張圖片plot(f,abs(y_zp(1:n/2)));%繪制加噪聲處理后音頻信號頻譜圖title('加噪聲處理音頻信號頻譜圖');%圖片標(biāo)題xlabel('頻率(Hz)');%x軸為頻率ylabel('頻率幅值(V)');%頻率幅值pause(1)%暫停一秒
