活动断层区域地壳垂直运动监测数据分析之速度面多项式模型建立程序代码.docVIP

附录一： 程序代码如下： #includestdio.h #includestdlib.h #includemath.h #define N 8//全局变量 double *mul(double *a,double *b,double *c,int m,int p,int n)//矩阵乘法 { int i,j,k; double temp; for(i=0;im;i++) { for(j=0;jn;j++) { temp=0.0; for(k=0;kp;k++) temp+=a[i*p+k]*b[k*n+j]; c[i*n+j]=temp; } } return c; } double *inv(double *m1, int n)//矩阵的逆矩阵 { int *is,*js; int i=0,j,k,l,u,v; double temp,max_v; is=(int *)malloc(n*sizeof(int)); js=(int *)malloc(n*sizeof(int)); if(is==NULL||js==NULL) { printf(out of memory!\n); return m1; } for(k=0;kn;k++) { max_v=0.0; for(int i=k;in;i++) { for(j=k;jn;j++) { temp=fabs(m1[i*n+j]); if(tempmax_v) { max_v=temp; is[k]=i; js[k]=j; } } } if(max_v==0.0) { free(is); free(js); printf(invers is not availble!\n); return m1; } if(is[k]!=k) { for(j=0;jn;j++) { u=k*n+j; v=is[k]*n+j; temp=m1[u]; m1[u]=m1[v]; m1[v]=temp; } } if(js[k]!=k) for(int i=0;in;i++) { u=i*n+k; v=i*n+js[k]; temp=m1[u]; m1[u]=m1[v]; m1[v]=temp; } l=k*n+k; m1[l]=1.0/m1[l]; for(j=0;jn;j++) { if(j!=k) { u=k*n+j; m1[u]*=m1[l]; } } for(i=0;in;i++) { if(i!=k) { for(j=0;jn;j++) { if(j!=k) { u=i*n+j; m1[u]-=m1[i*n+k]*m1[k*n+j]; } } } } for(i=0;in;i++) { if(i!=k) { u=i*n+k; m1[u]*=-m1[l]; } } } for(k=n-1;k=0;k--) { if(js[k]!=k) for(j=0;jn;j++) { u=k*n+j; v=js[k]*n+j; temp=m1[u]; m1[u]=m1[v]; m1[v]=temp; } if(is[k]!=k) for(int i=0;in;i++) { u=i*n+k; v=i*n+is[k]; temp=m1[u]; m1[u]=m1[v]; m1[v]=temp; } } free(is); free(js); return m1; // int i,j,k; // double p,q,*c; // c=(double *)malloc(n*sizeof(double)); // for(k=0;kn;k++) // 循环求解 // { // p=mat[0]; // p=1.0/p; // for(i=1;in-k;i++) // { // q=mat[i*n]; // c[i]=-q*p; // for(j=1;ji+1;j++) //