岩性地震学作业(-物探专业)

岩性地震作业1岩性地震学作业题目：根据下列地震道1000-1500ms时窗内各采样点的振幅值（采样间隔为1ms），求取以下10种类型振幅属性值。①均方根振幅，②平均绝对值振幅，③平均波峰振幅，④平均波谷振幅，⑤总绝对值振幅，⑥总振幅，⑦平均能量，⑧总能量，⑨平均振幅，⑩振幅平方差。要求：每人用其中一组数据计算，其所使用数据组序号与学生名单序号相符合（同学间所使用的数据应不相同），采用计算机编程计算，最后提交编程程序、各种振幅值数据、该地震道振幅曲线，并注明所采用的原始数据序号。采用数据编号：53号编程软件：MATLAB程序实现：clc;loaddata53.txt%加载数据s=0;s1=0;s2=0;s3=0;fori=1:501;T=data53(:,1);%取时间数据A=data53(:,2);%取振幅数据s=s+A(i)*A(i);%求振幅平方和，也就是总能量A1=abs(A);%取振幅绝对值s1=sum(A1);%绝对值求和A2=A(A0);%取正值s2=sum(A2);%波峰值求和A3=abs(A(A0));%取负值,取绝对值s3=sum(A3);%求和endA_rms=sqrt(s/501)%求均方根振幅A_aaa=s1/501%求绝对值振幅n1=length(A2);%求波峰长度A_apa=s2./n1%平均波峰振幅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)%平均振幅岩性地震作业2V_a=sum((A-A_ma).^2)/(n1+n2)%振幅平方差plot(T,A,'r')%画图，振幅曲线title('53号数据地震道振幅曲线');ylabel('振幅值');xlabel('时间ms');输出结果：①均方根振幅：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⑨平均振幅：A_ma=-0.0661⑩振幅平方差：V_a=3.2883e+03振幅曲线：过程截图：岩性地震作业3