資源描述：
《工程測試與信號處理作業(yè)》由會員上傳分享，免費在線閱讀，更多相關內容在工程資料-天天文庫。
1、每個同學根據(jù)所具備的實驗測量條件，采集一個多頻率成分的信號，要求:1.介紹試驗測試方案，并利用MATLAB對采集的信號進行預處理（掐頭去尾等）（5分）。根據(jù)試驗測量方案、現(xiàn)場照片和采集的信號，對實驗測試方案的設計和實驗過程進行自我評價（5分）。（共10分）2.對信號進行你認為必要的時域分析，并說明分析結果的意義。（共10分）3.對信號做快速傅里葉變換（FFT）。寫出MATLAB程序（5分）；畫出信號的FFT結果圖形（5分）；對FFT結果進行分析，并寫出分析結論（10分）。（共20分）4.對信號做窗口傅里葉變換。寫出MATL
2、AB程序（5分）；畫出分析結果圖形（5分）；對變換結果進行分析，并寫出分析結論（10分）。（共20分）5.對信號做小波變換。寫出MATLAB程序（5分）；畫出分析結果圖形（5分）；對變換結果進行分析，并寫出分析結論（10分）。（共20分）6?設計能濾除信號中你認為不需要的頻率成分的數(shù)字濾波器，并求信號通過該濾波器后的輸出。寫出MATLAB程序（5分）；繪出輸出的時域結果和頻域結果（5分）；對結果進行分析，并寫出分析結論（10分）。（共20分）一采集信號并進行預處理：試驗測量方案：實驗現(xiàn)場圖片如圖1?1所示。首先用一根細繩了
3、將一塊鐵質零件懸掛起來，繩子的另一端固定，打開錄音設備，然后用錘子迅速敲擊零件，然后迅速地收回錘子以免出現(xiàn)連擊現(xiàn)象，關閉錄音設備，并將其保存成WAV格式的文件。經(jīng)過多次實驗并采集聲音信號，篩選出噪音對所采集的聲音信號干擾最小的一個，在MATLAB環(huán)境下對所采集的信號進行分析圖1.1實驗現(xiàn)場因為在信號圖中要顯示其采樣頻率，采用如下MATLAB程序段:[x,fs]=wavread(,F:14.wav,);%調入錄制的聲音信號fs%顯不米樣頻率運行結果為：fs=48000說明采樣頻率為48000.木次試驗，選擇挑選所做實驗屮效
4、果最好的咅頻信號，在MATLAB屮調用如下程序進行分析：[x,fs]=wavread('F:l4.wav')%調用錄制的音頻文件y=x(:,l)sound(y,fs);[m,n]=size(y);t=(O:m-l)/fs;subplot(3,l,l);plot(t,y)titleCH始信號波形J;xlabelCTime(s)');ylabeK'Amplitude^gridon;y=y(40000:m-15000)t=t(40000:m-15000)subplot(3,l,2)plot(t,y)titled截取后的信號波形
5、J;xlabel('Time(s)');ylabel(Amplitude');gi*idon;yft=fft(y,m);Ay=abs(yft);f=(0:m-l)*fs/m;subplot(3,l,3)；plot(f(l:(m-1)/2),2*Ay(1:(m-l)/2))titled幅頻譜圖J;xlabelCPrequency(kHz)');ylabelCAmplitude');gridon;savefFAwp.matVx1);如圖1.2所示，在MATLAB中的commandwindow窗口中輸入m。得到采樣點數(shù)m=141
6、600個。m-141600fx、、圖1.2commandwindow窗口輸入命令通過本次試驗，會發(fā)現(xiàn)由于各種原因都會導致音頻信號有噪音出現(xiàn)，此時需進行一定的掐頭去尾的處理。前邊去掉40000個點，后邊去掉15000個點，得到如下所示初始波形、截取后的波形和頻域分析圖形。Figure1□回岡a)pn七-dlu<原始信號波形截取后的信號波形0.811.21.41.61.822.22.42.62.8Time(s)aprn--duj<4x10振幅頻譜圖00.511.5Frequency(kHz)22.5x10410FileEdit
7、ViewInsertToolsDesktopWindowHelp圖1?3處理后信號2效果圖二時域分析利用MATLAB對信號進行時域分析，得出結論，相應的程序如下:[x,fs,bits]二wavread(F14.wav');%調入錄制的咅頻信號[N,M]=size(y)%顯示音頻信號的尺寸，并保存t=O:N-l;t=t*l/fs;%計算時間序列subplot(2,2,1);plot(t,y);gridon%繪制波形圖title('吋域分析，fs=48000Hz,);xlabel(,time/s');ylabel('Ampl
8、itude,);meany=mean(y)%求、卜均值，信虧的H流分卑subplot(2,2,2);plot(t,meany);gridon%繪制均值圖vary=var(y)%求方差，繞均值波動的程度subplot(2,2,3)；plot(t,vary);gridon%繪制方差圖meany2=mean(
