巖性地震學(xué)作業(yè)(-物探專業(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上巖性地震學(xué)作業(yè)題目：根據(jù)下列地震道1000-1500ms時(shí)窗內(nèi)各采樣點(diǎn)的振幅值（采樣間隔為1ms），求取以下10種類型振幅屬性值。均方根振幅，平均絕對值振幅，平均波峰振幅，平均波谷振幅，總絕對值振幅，總振幅，平均能量，總能量，平均振幅，振幅平方差。要求：每人用其中一組數(shù)據(jù)計(jì)算，其所使用數(shù)據(jù)組序號(hào)與學(xué)生名單序號(hào)相符合（同學(xué)間所使用的數(shù)據(jù)應(yīng)不相同），采用計(jì)算機(jī)編程計(jì)算，最后提交編程程序、各種振幅值數(shù)據(jù)、該地震道振幅曲線，并注明所采用的原始數(shù)據(jù)序號(hào)。采用數(shù)據(jù)編號(hào)： 53號(hào)編程軟件： MATLAB程序?qū)崿F(xiàn)：clc;load data53.txt %加載數(shù)據(jù)s=0;s1=0;s

2、2=0;s3=0;for i=1:501; T=data53(:,1); %取時(shí)間數(shù)據(jù) A=data53(:,2); %取振幅數(shù)據(jù) s=s+A(i)*A(i); %求振幅平方和，也就是總能量 A1=abs(A); %取振幅絕對值 s1=sum(A1); %絕對值求和 A2=A(A>0); %取正值 s2=sum(A2); %波峰值求和 A3=abs(A(A<0); %取負(fù)值,取絕對值 s3=sum(A3); %求和endA_rms=sqrt(s/501) %求均方根振幅A_aaa=s1/501 %求絕對值振幅n1=length(A2); %求波峰長度A_apa=s2./n1 %平均

3、波峰振幅n2=length(A3); %求波谷長度A_ata=s3/n2 %平均波谷振幅A_taa=s1 %總絕對值振幅A_ta=sum(A) %總振幅E_a=s/501 %平均能量E_t=s %總能量A_ma=A_ta/(n1+n2) %平均振幅V_a=sum(A-A_ma).2)/(n1+n2) %振幅平方差 plot(T,A,'r') %畫圖，振幅曲線title('53號(hào)數(shù)據(jù)地震道振幅曲線');ylabel('振幅值');xlabel('時(shí)間ms');輸出結(jié)果：均方根振幅：A_rms = 57.1145平均絕對值振幅：A_aaa = 45.1646平均波峰振幅：A_apa = 44.8307平均波谷振幅：A_ata = 46.2455總絕對值振幅：A_taa = 2.2627e+04總振幅：A_ta = -32.8290平均能量：E_a = 3.2621e+03總能量：E_t = 1.6343e+06平均振幅