c语言误差椭圆,2021误差理论与测量平差上机指导书.docx

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

辽宁工程技术大学

测绘与地理科学学院测绘工程系

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

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

目录

TOC \o "1-5" \h \z \o "Current Document" 实验1矩阵加法与乘法运算 3

\o "Current Document" 实验2矩阵转置与求逆运算 5

\o "Current Document" 实验3误差椭圆元素计算13

\o "Current Document" 实验4水准网间接平差程序设计 14

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

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

{ PAGE

{

PAGE #

实验1矩阵加法与乘法运算

一、实验名称：矩阵加法与乘法运算。

二、实验目的和任务：掌握矩阵加法与乘法通用程序的编写。

三、实验要求：

1每人独立编写出矩阵加法与乘法的程序，并上机调试通过；

2米用VC++6.0开发平台，C或者C++语言编写程序；

3写出矩阵运算的结果。

四、实验内容：

1矩阵加法的示例函数(C语言)

void JZjiafa(double a[15][15],double b[15][15],doublec[15][15],i nt

m,i nt n)

{

for (i nt i=0;i<=m-1;i++)

for(i nt j=0;j<=n-1;j++)

{

c[i][j]=a[i][j]+b[i][j];

}

return;

}

2矩阵乘法的示例程序(C语言)

#in clude "stdafx.h"

void matrixMultiply(double a[14][15],doubleb[15][13], double

c[14][13],lo ng m,lo ng n,lo ng k)

{

for (l ong i = 0; i<= m-1; i++)

for (l ong j=0; j<=k-1; j++)

{

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

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

PAGE

PAGE #

c[i][j] =0.0;

for (long q=0; qv=n-1;q++)

{

c[i][j] = c[i][j] + a[i][q] * b[q][j];

}

}

}

return;

}

int main (i nt 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");

fsca nf(stream,"%ld %ld",&n,&m);

for (i=0;i

{

for(j=0;j

{

fscan f(stream,"%lf",&a[i][j]);

}

}

fsca nf(stream,"%ld %ld",&m,&k);

for(i=0;i

{

for(j=0;j

{

fscan f(stream,"%lf",&b[i][j]);

}

}

fclose(stream);

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

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

{ PAGE

{

PAGE #

matrixMultiply(a,b,c,4,5,3);

stream = fopen(”矩阵计算结果.txt","w"); for (i=0;i<=3;i++)

{

for(j=0;j<=2;j++)

fprin tf(stream,"%16.7e",c[i][j]);

fprin tf(stream,"\n");

}

fprin tf(stream,"\n");

fclose(stream);

return 0;

}

实验2矩阵转置与求逆运算

一、实验名称：矩阵转置与求逆运算。

二、实验目的和任务：掌握矩阵转置的编写，会调用矩阵求逆函数。

三、实验要求：

1每人独立编写出矩阵转置的程序，并上机调试通过；

2每人独立完成矩阵求逆函数的调用，并调试通过；

3米用VC++6.0开发平台，C或者C++语言编写程序；

4写出矩阵运算的结果。

四、实验内容：

矩阵的转置示例函数(C语言)

double JZzhua nzhi(double a[15][15], double b[15][15], i nt m,i nt n)

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

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

{ PAGE

{

PAGE #

for(int i=O;i

for(i nt j=O;j

b[j][i]=a[i][j];

}

return 0.0;