地震波觀測系統(tǒng)的matlab仿真報(bào)告(參考版)
2025-03-26 06:29本頁面
【正文】 在老師和同學(xué)的耐心指導(dǎo)下我發(fā)現(xiàn)了自己在選擇巴特沃斯、切比雪夫?yàn)V波器上的問題，經(jīng)過修改和調(diào)試，終于得到了應(yīng)有的效果，這讓我看到了理論與實(shí)踐相結(jié)合的優(yōu)勢與用處，讓我受益匪淺。讓我們在學(xué)習(xí)理論知識的同時(shí)，明白如何把這些應(yīng)用于實(shí)際。在我們身邊處處都能看到數(shù)字信號處理的相關(guān)知識的應(yīng)用，從語音的識別采集處理到地震波觀測，這直觀的證實(shí)了數(shù)字信號處理這門課程的重要性。 %頻率軸只畫出2Hz頻率之前的部分grid on。)。ylabel(39。幅頻圖39。頻率/Hz39。subplot(2,1,2),plot([0:N1]/N*Fs,abs(Yf))。)。title(39。時(shí)間域39。時(shí)間/s39。)。 %繪制原始值序列title(39。 %對信號進(jìn)行快速Fourier變換(FFT) figure(3)。 %得到序列的長度t=[0:N1]*dt。 %設(shè)定采樣率 單位(Hz)dt=1/Fs。 %讀取數(shù)據(jù)序列Xt=hns3。 %頻率軸只畫出2Hz頻率之前的部分grid on。)。ylabel(39。幅頻圖39。頻率/Hz39。subplot(2,1,2),plot([0:N1]/N*Fs,abs(Yf))。)。title(39。時(shí)間域39。時(shí)間/s39。)。 %繪制原始值序列title(39。 %對信號進(jìn)行快速Fourier變換(FFT)figure(2)。 %得到序列的長度t=[0:N1]*dt。 %設(shè)定采樣率 單位(Hz)dt=1/Fs。 %讀取數(shù)據(jù)序列Xt=hns2。 %頻率軸只畫出2Hz頻率之前的部分grid on。)。ylabel(39。幅頻圖39。頻率/Hz39。subplot(2,1,2),plot([0:N1]/N*Fs,abs(Yf))。)。title(39。時(shí)間域39。時(shí)間/s39。)。 %繪制原始值序列title(39。 %對信號進(jìn)行快速Fourier變換(FFT)figure(1)。 %得到序列的長度t=[0:N1]*dt。 %設(shè)定采樣率 單位(Hz)dt=1/Fs。 %讀取數(shù)據(jù)序列Xt=hns1。grid on。振幅39。)。),title(39。xlabel(39。r39。grid on。振幅39。)。),title(39。xlabel(39。g39。figure(14)。)。title(39。相位/^o39。)。subplot(2,1,2),plot(f,angle(H)) %畫出寬帶濾波器的相頻特性xlabel(39。)。title(39。A(\lambda)/db39。)。%畫出寬帶濾波器的幅頻特性xlabel(39。%按傳遞函數(shù)系數(shù)、數(shù)據(jù)點(diǎn)數(shù)和采樣頻率求得濾波器的頻率特性y3=filtfilt(b,a,Xt)。%按最小階數(shù)、截止頻率、通帶波紋和阻帶衰減設(shè)計(jì)濾波器figure(13)。[Order,Wn]=cheb2ord(wp,ws,Rp,Rs)。Rs=20。wp=[ ]*2/Fs。grid on。振幅39。)。,1)xlabel(39。,39。) %繪制輸入信號legend(39。plot(t,y1,t,real(x1),39。for ii=1:N %得到仿真結(jié)果 if (abs(H1(ii))) XX1(ii)=Yf(ii).*H1(ii)/H(ii)。)。 %在窄帶濾波器上的輸出[H1,f]=freqz(b,a,N,Fs,39。grid on。相位/^o39。)。subplot(2,1,2),plot(f,180/pi*unwrap(angle(H1)))xlabel(39。)。ylabel(39。頻率/Hz39。 %按最小階數(shù)、截止頻率、通帶波紋和阻帶衰減設(shè)計(jì)濾波器figure(11)[H1,f]=freqz(b,a,Nn,Fs)。 %通帶波紋和阻帶衰減以及繪制頻率特性的數(shù)據(jù)點(diǎn)數(shù)[order,Wn]=ellipord(wp,ws,Rp,Rs)。Rs=20。wp=[ ]*2/Fs。grid on。振幅39。)。),title(39。xlabel(39。r39。grid on。振幅39。)。),title(39。xlabel(39。g39。 figure(10)。)。ylabel(39。時(shí)間/s39。)。title(39。 subplot(1,2,2)。)。ylabel(39。時(shí)間/s39。)。title(39。subplot(1,2,1)。for
試題試卷相關(guān)推薦
鄂ICP備17016276號-1