误差理论与测量平差上机指导书.docVIP

误差理论与测量平差上机指导书 辽宁工程技术大学 测绘与地理科学学院测绘工程系  目录 实验1 矩阵加法与乘法运算 1 实验2 矩阵转置与求逆运算 3 实验3 误差椭圆元素计算 11 实验4 水准网间接平差程序设计 12  实验1 矩阵加法与乘法运算 实验名称：矩阵加法与乘法运算。 实验目的和任务：掌握矩阵加法与乘法通用程序的编写。 实验要求： 每人独立编写出矩阵加法与乘法的程序，并上机调试通过； 采用VC++6.0开发平台，C或者C++语言编写程序； 写出矩阵运算的结果。 实验内容： 矩阵加法的示例函数(C语言) void JZjiafa(double a[15][15],double b[15][15],double c[15][15],int m,int n) { for (int i=0;i=m-1;i++) for(int j=0;j=n-1;j++) { c[i][j]=a[i][j]+b[i][j]; } return; } 矩阵乘法的示例程序(C语言) #include stdafx.h void matrixMultiply(double a[14][15],double b[15][13], double c[14][13],long m,long n,long k) { for (long i = 0; i= m-1; i++) { for (long j=0; j=k-1; j++) { c[i][j] =0.0; for (long q=0; q=n-1;q++) { c[i][j] = c[i][j] + a[i][q] * b[q][j]; } } } return; } int main(int argc, char* argv[]) { long n,m,k,i,j; double a[14][15],c[14][13],b[15][13]; FILE *stream; stream = fopen(矩阵输入.txt,r); fscanf(stream,%ld %ld,n,m); for (i=0;in;i++) { for(j=0;jm;j++) { fscanf(stream,%lf,a[i][j]); } } fscanf(stream,%ld %ld,m,k); for(i=0;im;i++) { for(j=0;jk;j++) { fscanf(stream,%lf,b[i][j]); } } fclose(stream); matrixMultiply(a,b,c,4,5,3); stream = fopen(矩阵计算结果.txt,w); for (i=0;i=3;i++) { for(j=0;j=2;j++) fprintf(stream,%16.7e ,c[i][j]); fprintf(stream,\n); } fprintf(stream,\n); fclose(stream); return 0; } 实验2 矩阵转置与求逆运算 实验名称：矩阵转置与求逆运算。 实验目的和任务：掌握矩阵转置的编写，会调用矩阵求逆函数。 实验要求： 每人独立编写出矩阵转置的程序，并上机调试通过； 每人独立完成矩阵求逆函数的调用，并调试通过； 采用VC++6.0开发平台，C或者C++语言编写程序； 写出矩阵运算的结果。 实验内容： 矩阵的转置示例函数(C语言) double JZzhuanzhi(double a[15][15], double b[15][15], int m,int n) { { for(int i=0;im;i++) for(int j=0;jn;j++) b[j][i]=a[i][j]; } return 0.0; } 矩阵求逆的示例函数(C语言) int invGJ(double **a,int n) { int *is,*js,i,j,k,l,u,v; double d,p; is=(int *)malloc(n*sizeof(int)); js=(int *)malloc(n*sizeof(int)); for(k=0;k=n-1;k++) { d=0.0; for(i=k;i=n-1;i++) for(j=k;j=n-1;j++) { l=i*n+j;p=fabs(a[i][j]