岩性地震学作业 答案.doc

岩性地震学作业导航目录：一、振幅曲线和振幅属性值2二、程序源代码之振幅属性值3三、程序源代码之振幅曲线91一、振幅曲线和振幅属性值8二、程序源代码之振幅属性值#include #include #include #include #define N 0struct Amplitudedataint nnn;float aaa;struct Amplitudedata *xxx;#define LENGTH sizeof(struct Amplitudedata)struct Amplitudedata *yyy(void) struct Amplitudedata *bbb,*ccc,*ddd;int i,j=0,number1502;float number502;FILE *fp;if(fp=fopen(63.txt,r)=N)printf( ERROR n);elsefor(i=0;iaaa=number0; ccc-nnn=number10;int n=0; for(i=1;ixxx=ccc; ddd=ccc;ccc=(struct Amplitudedata *)malloc(LENGTH);ccc-nnn=number1i;ccc-aaa=numberi; ccc-xxx=N;ddd-xxx=ccc;return bbb;float ggg(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,*ddd,ggg;bbb=eee;ggg.aaa=0;ddd=(struct Amplitudedata *)malloc(LENGTH);for(ccc=bbb;ccc!=N;ccc=ccc-xxx)if(ggg.aaaaaa)ggg.aaa=ccc-aaa;return ggg.aaa;void jjj(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc;ccc=bbb=eee; float total(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,*ddd;bbb=eee;float fff=0;ddd=(struct Amplitudedata *)malloc(LENGTH);for(ccc=bbb;ccc!=N;ccc=ccc-xxx)fff=fff+ccc-aaa;return fff;float jdzf(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,ddd;float fff=0;bbb=eee;int i=0;for(ccc=bbb;ccc!=N;ccc=ccc-xxx)if(ccc-aaaaaa;elseddd.aaa=ccc-aaa;fff=fff+ddd.aaa;return fff;float pjjdzf(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,ddd;int i=0;bbb=eee;float fff=0;for(ccc=bbb;ccc!=N;ccc=ccc-xxx)if(ccc-aaaaaa;elseddd.aaa=ccc-aaa;fff=fff+ddd.aaa;i+;fff=fff/(float)i;return fff;float Averaging(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,ddd;int i=0;bbb=eee;float fff=0;for(ccc=bbb;ccc!=N;ccc=ccc-xxx)ddd.aaa=ccc-aaa;fff=fff+ddd.aaa;i+;fff=fff/(float)i;return fff;void main()int number1(struct Amplitudedata *eee);int number2(struct Amplitudedata *eee);float iii,fff,jfff,averagejfff,hhh,variance,average,averagepeak,averagebottom;struct Amplitudedata *bbb;int i,j;bbb=yyy();jjj(bbb);iii=ggg(bbb);fff=total(bbb);jfff=jdzf(bbb);i=number1(bbb);j=number2(bbb);averagejfff=jfff/(float)i;average=fff/(float)j;printf(|总振幅：%.3f |n| |n|总绝对值振幅：%f |n| |n|平均绝对值振幅：%.3f |n| |n|平均振幅：%f |n| |n,fff,jfff,averagejfff,average);float peak(struct Amplitudedata *eee);float bottom(struct Amplitudedata *eee);averagepeak=peak(bbb);if(averagepeak!=0)printf(|平均波峰振幅：%f |n| |n,averagepeak);averagebottom=bottom(bbb);if(averagebottom!=0)printf(|平均波谷振幅：%f |n| |n,averagebottom);float variety(struct Amplitudedata *eee);variance=variety(bbb);printf(|振幅平方差：%f |n| |n,variance);float kkk(struct Amplitudedata *eee);hhh=kkk(bbb);printf(|均方根振幅：%f |n| |n,sqrt(hhh/(float)i);printf(|平均能量：%f |n| |n,hhh/(float)i);printf(|总能量：%f |n|-|n,hhh);getchar();float peak(struct Amplitudedata *eee) struct Amplitudedata *lll(struct Amplitudedata *bbb1,struct Amplitudedata *mmm,int j);float averagepeakbottom(struct Amplitudedata *bbb3);struct Amplitudedata *bbb,*ccc;int k;bbb=eee;float fff=0;struct Amplitudedata *bbb2=N;for(ccc=bbb,k=1;ccc!=N;ccc=ccc-xxx)if(ccc-aaa0)bbb2=lll(bbb2,ccc,k);k+;if(k=1)return fff=0;elsefff=averagepeakbottom(bbb2); return fff;float bottom(struct Amplitudedata *eee) struct Amplitudedata *lll(struct Amplitudedata *bbb1,struct Amplitudedata *mmm,int j);float averagepeakbottom(struct Amplitudedata *bbb3);struct Amplitudedata *bbb,*ccc;int k;bbb=eee;float fff=0;struct Amplitudedata *bbb3=N;for(ccc=bbb,k=1;ccc!=N;ccc=ccc-xxx)if(ccc-aaaaaa;ccc1=(struct Amplitudedata *)malloc(LENGTH);if(i=1) ccc1-aaa=eee; ccc1-xxx=N; bbb1=ccc1; else ccc1-aaa=eee; ccc1-xxx=N; ddd1=bbb1; while(ddd1-xxx!=N) ddd1=ddd1-xxx; ddd1-xxx=ccc1; return bbb1;float averagepeakbottom(struct Amplitudedata *bbb3)struct Amplitudedata *ccc2;int i2=0;float eee=0;ccc2=bbb3;while(ccc2!=N) eee=eee+ccc2-aaa;i2=i2+1;ccc2=ccc2-xxx;eee=eee/(float)i2;return eee;float variety(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,ddd;int i=0;bbb=eee;float average,fff=0;float Averaging(struct Amplitudedata *eee);average=Averaging(bbb);for(ccc=bbb;ccc!=N;ccc=ccc-xxx)fff=fff+(ccc-aaa-average)*(ccc-aaa-average);if(ccc-aaa!=0)i+;fff=fff/(float)i;return fff;int number1(struct Amplitudedata *eee)struct Amplitudedata *ccc;int i=0;for(ccc=eee;ccc!=N;ccc=ccc-xxx)i+;return i;int number2(struct Amplitudedata *eee)struct Amplitudedata *ccc;int i=0;for(ccc=eee;ccc!=N;ccc=ccc-xxx)if(ccc-aaa!=0)i+;return i;float kkk(struct Amplitudedata *eee)struct Amplitudedata *bbb,*ccc,ddd;int i=0;bbb=eee;float fff=0;for(ccc=bbb;ccc!=N;ccc=ccc-xxx)fff=fff+ccc-aaa*ccc-aaa;i+;return fff;三、程序源代码之振幅曲线/ imageView.cpp : implementation of the CImageView class/#include stdafx.h#include image.h#include imageDoc.h#include imageView.h#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CImageViewIMPLEMENT_DYNCREATE(CImageView, CScrollView)BEGIN_MESSAGE_MAP(CImageView, CScrollView)/AFX_MSG_MAP(CImageView)/ NOTE - the ClassWizard will add and remove mapping macros here./ DO NOT EDIT what you see in these blocks of generated code!/AFX_MSG_MAP/ Standard printing commandsON_COMMAND(ID_FILE_PRINT, CScrollView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT, CScrollView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_PREVIEW, CScrollView:OnFilePrintPreview)END_MESSAGE_MAP()/ CImageView construction/destructionCImageView:CImageView()/ TODO: add construction code hereCImageView:CImageView()BOOL CImageView:PreCreateWindow(CREATESTRUCT& cs)/ TODO: Modify the Window class or styles here by modifying/ the CREATESTRUCT csreturn CScrollView:PreCreateWindow(cs);/ CImageView drawingvoid CImageView:OnDraw(CDC* pDC)CImageDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data hereint i;double k501=66,50.957,29,8.057,-11,-29.972,-48,-63.538,-76,-85.031,-90,-90.998,-88,-81.716,-73,-62.835,-52,-40.863,-30,-19.676,-11,-4.253,0,2.467,3,2.265,0,-3.548,-9,-16.034,-24,-31.635,-38,-42.502,-45,-45.404,-43,-37.6,-30,-21.535,-13,-4.787,3,10.187,16,19.899,22,22.955,23,22.472,22,21.907,22,21.854,22,22.604,24,25.658,28,30.815,35,40.554,48,56.788,67,77.743,88,96.392,102,104.145,102,95.219,83,65.882,45,21.859,-3,-29.234,-57,-84.086,-105,-115.846,-119,-119.635,-119,-117.835,-119,-122.164,-119,-103.552,-80,-55.694,-33,-11.009,10,29.67,47,60.901,70,74.642,76,75.633,73,67.742,61,54.28,48,41.807,36,30.347,25,19.641,15,11.217,9,8.192,9,10.997,14,17.508,21,23.795,25,23.996,20,13.146,4,-6.419,-18,-30.218,-42,-51.886,-59,-62.962,-64,-62.541,-59,-53.82,-48,-42.425,-38,-34.881,-33,-31.7,-31,-30.358,-29,-26.09,-22,-17.205,-12,-6.307,0,6.74,13,18.046,22,25.392,28,29.61,30,29.231,27,23.312,18,11.367,4,-3.513,-11,-18.358,-25,-30.183,-33,-33.143,-31,-27.424,-23,-18.278,-14,-10.806,-9,-8.551,-9,-9.796,-10,-9.033,-6,-0.614,8,19.61,34,49.747,66,81.217,95,106.26,114,117.661,119,119.512,119,116.346,111,102.868,92,78.497,62,42.444,21,-1.249,-24,-46.981,-68,-85.14,-97,-103.475,-104,-98.98,-89,-75.456,-60,-44.286,-30,-17.951,-9,-3.031,-1,-2.339,-6,-10.119,-15,-20.363,-25,-27.463,-28,-27.203,-25,-21.254,-16,-9.811,-4,0.715,5,9.61,13,13.871,12,8.371,3,-3.8,-12,-20.959,-30,-37.949,-44,-47.618,-49,-48.55,-47,-44.74,-42,-39.044,-37,-36.505,-37,-37.314,-37,-35.647,-32,-25.481,-17,-7.514,4,18.026,33,46.547,58,67.03,73,75.591,76,75.32,74,71.963,70,68.399,67,65.087,63,60.46,57,51.689,45,37.199,29,20.581,13,6.787,3,1.683,3,6.033,10,13.929,18,21.992,25,26.005,25,22.603,19,14.651,10,5.409,0,-6.562,-14,-21.405,-29,-36.819,-45,-53.039,-61,-68.318,-75,-80.464,-85,-88.307,-91,-92.996,-95,-96.853,-99,-100.865,-102,-101.411,-99,-94.333,-87,-76.497,-63,-47.098,-30,-12.579,5,22.607,39,52.977,64,72.279,78,81.515,83,82.672,81,78.354,75,70.802,66,60.359,54,46.639,39,31.252,24,17.207,12,8.536,7,6.679,8,10.683,15,20.048,25,28.907,33,37.867,43,47.166,51,55,59,62.309,65,66.859,67,64.657,60,53.371,44,31.784,18,4.216,-10,-24.968,-40,-53.89,-67,-79.322,-90,-98.184,-105,-111.176,-116,-118.188,-118,-115.683,-111,-103.353,-93,-80.23,-66,-50.787,-35,-18.847,-4,8.639,19,27.926,35,40.072,43,44.251,44,42.758,41,39.102,37,34.551,32,29.301,26,21.548,16,9.834,4,-1.099,-6,-11.401,-17,-22.039,-26,-28.717,-30,-30.155,-30,-29.916,-29,-26.481,-23,-19.461,-16,-12.259,-8,-3.477,0,1.73,3,5.195,8,10.44,13,16.127,20,24.349,30,36.974,45,53.016,61,68.477,75,79.778,83,84.51,84,81.028,76,69.447,62,53.831,45,35.338,25,14.292,4,-5.598,-15,-24.72,-34,-42.019,-49,-55.178,-60,-62.77,-64,-64.216,-64,-63.362,-63; /画坐标pDC-MoveTo(50,50);/纵轴pDC-LineTo(50,750);pDC-MoveTo(50,50);/箭头pDC-LineTo(40,60);pDC-MoveTo(50,50);pDC-LineTo(60,60);pDC-MoveTo(50,400);/横轴pDC-LineTo(1100,400);pDC-MoveTo(1100,400);/箭头pDC-LineTo(1090,390);pDC-MoveTo(1100,400);pDC-LineTo(1090,410);for(i=0;iMoveTo(50