测绘程序设计课程实习报告

测绘程序设计课程实习报告本文简介：测绘程序设计课程实习报告一、实习目的巩固之前所学的c++类定义、文件操作等各方面知识，掌握矩阵计算的基本算法，熟悉并熟练运用平差模型，掌握测量平差的基本原理和公式，并处理测量观测数据。二、实习内容该课程的实习主要为c语言及c++编程基础，矩阵操作算法，条件平差，间接平差，以及在水准网平差中的应用。三

测绘程序设计课程实习报告本文内容：

测绘程序设计课程实习报告

一、实习目的

巩固之前所学的c++类定义、文件操作等各方面知识，掌握矩阵计算的基本算法，熟悉并熟练运用平差模型，掌握测量平差的基本原理和公式，并处理测量观测数据。

二、实习内容

该课程的实习主要为c语言及c++编程基础，矩阵操作算法，条件平差，间接平差，以及在水准网平差中的应用。

三、实习过程

（1）

实习1（文件的基本操作以及角度、弧度的转化）

1.

准备

C++文件操作函数包含于头文件，一般情况下使用即可。在计算过程中需要文件输出格式控制，函数包含于。c++的流操作函数包含于，相应的流函数在命名空间“std”。

角度转弧度的过程中，首先判断正负，对符号进行不同情况的处理。将度数的分，秒位分别除以60、3600，然后把得到的结果与度位绝对值相加，判断符号。

弧度转角度的过程中，同样判断正负，取小数位乘以60取整数得到相应角度的分，再取剩余的小数乘以60即得到相应角度的秒。

2.

代码

#include

#include

#include

#include

using

namespace

std;

const

double

pi

=

3.14159;

int

main()

{

fstream

file(“a.txt“);

//定义文件操作对象

string

aa;

file>>aa;

cout>a1;

cout>a2;

cout>a3;

cout=

0)

//判断输入度数的正负

deg

=

a1

+

a2/60

+

a3/3600;

else

deg

=

a1

-

a2/60

-

a3/3600;

double

radian;

radian

=

deg/180*pi;

cout>radian;

bool

tm

=

0;

deg

=

radian*180/pi;

if(deg

#include

#include

using

namespace

std;

class

mat

{

public:

void

add(int

row,int

col,double

m[],double

n[],double

result[]);//加法

void

sub(int

row,int

col,double

m[],double

n[],double

result[]);//减法

double*

multiply(int

row,int

col,int

col2,double

m[],double

n[]);//乘法

void

trans(int

row,int

col,double

m[],double

result[]);//转置

bool

inverse(double

a[],int

n);//求逆

void

fileIn(ifstream

//从文件中读入

void

fileOut(ofstream

out,int

row,int

col,doublemat);

//向文件输出

void

outPut(int

row,int

col,doublemat);//向控制台输出

void

inPut(int

row,int

col,doublemat);//从控制台输入

int

ij(int

i,int

j)//下三角矩阵元素判断

{

return

(i>=j)?

i*(i+1)/2+j

:

j*(j+1)/2+i;

}

};

文件mat.cpp：

#include

“mat.h“void

mat::add(int

row,int

col,double

m[],double

n[],double

result[])

{

for(int

i=0;

i>mat[i];

}

void

mat::fileOut(ofstream

out,int

row,int

col,doublemat)

{

for(int

i=0;

i>mat[i*col+j];

}

测试文件main.cpp:

#include

“mat.h“mat

g_mat;

void

main()

{

doublemat1;

doublemat2;

doublemat3;

doublemat4;

doublemat_result;

mat1

=

new

double[2*3];

mat2

=

new

double[2*3];

mat3

=

new

double[3*1];

mat4

=

new

double[2*2];

cout

mat

mat_c;

/**

V

=

BX-L;

V:

n*1

B:

n*t

X:

t*1

L:

n*1

P:

n*1

已知：B、L、P

求：X、V、/

void

adjust(int

n,int

t,double

V[],double

B[],double

X[],double

L[],double

P[],double

doubleBTPB

=

new

double[(t+1)*t/2];

doubleBTPB2

=

new

double[(t+1)*t/2];

for(i=0;

i>n;

file>>t;

cout，使用

cout

#include

#include

#include

“mat.h“using

namespace

std;

mat

mat_cul;

void

condition(int

n,int

r,double

A[],double

W[],double

P[],double

V[],double

doubleAQAT

=

new

double[r*(r+1)/2];

//法方程系数阵

for(i=0;

i

#include

#include

#include

#include

#include

“mat.h“using

namespace

std;

class

Clevel

{

public:

int

n;

//总观测值数

int

p_k;

//已知点数

int

t;

//未知点数

int

pNum;

//高程点总数

double

m;

//单位权

stringPname;

//点名地址数组

intstartP;

//起点号

intendP;

//终点号

doubleh;

//观测值

doubleH;

//高程值

doubleS;

//距离长度

doubleP;

//权

doubleB;

doubleL;

doubleX;

doubleV;

double

u;

//中误差

Clevel(charfileName);

void

inputdata(charfileName);

int

GetPNum(string

name);

void

adjust();

void

outputResult();

void

calculataHBPL();

mat

mat_c;//矩阵计算类

private:

void

calH();

void

calB();

void

calP();

void

calL();

};

文件Clevel.cpp:

#include

“Clevel.h“Clevel::Clevel(char*

fileName)

{

inputdata(fileName);

calculataHBPL();

adjust();

outputResult();

}

/**

V

=

BX-L;

V:

n*1

B:

n*t

X:

t*1

L:

n*1

P:

n*1

已知：B、L、P

求：X、V、/

void

Clevel::adjust()

{

int

i,j,k;

doubleBTPB

=

new

double[(t+1)*t/2];

doubleBTPB2

=

new

double[(t+1)*t/2];

//

下三角求BTPB2

for(i=0;

i>n>>pNum>>p_k>>m;

t

=

pNum-p_k;

startP

=

new

int[n];//起点号

endP

=

new

int[n];//终点号

L

=

new

double[n];//观测值

H

=

new

double[pNum];//高程值

S

=

new

double[n];//距离长度

P

=

new

double[n*n];//权

B

=

new

double[t*n];

h

=

new

double[n];

Pname

=

new

string[pNum];

X

=

new

double[t];

V

=

new

double[n];

int

i;

for(i=0;i>buffer;

int

bb

=

GetPNum(buffer);

in>>H[bb];

//cout>buffer;//读取起点名

startP[i]

=

GetPNum(buffer);//获取起点对应的点号

in>>buffer;//读取终点名

endP[i]

=

GetPNum(buffer);

in>>h[i]>>S[i];//读取观测值高差及对应的距离

}

in.close();

//cout-9999.0

Hcount++;

}

}

if(Hcount

==

t)

break;

}

}

void

Clevel::calB()

{

int

i,j;

for(i=0;

i=p_k

}

}

//输出矩阵B

cout<

for(

i=0;i

{

for(int

j=0;j

{

cout<

}

cout<

}

}

void

Clevel::calP()

{

int

i;

for(i=0;

i

i++)

{

P[i]=m/S[i];

}

//输出矩阵P

cout<

for(i=0;

i

i++)

{

cout<

}

}

void

Clevel::calL()

{

int

i;

for(i=0;i

{

int

k1=startP[i];

int

k2=endP[i];

L[i]=h[i]-(H[k2]-H[k1]);

}

//输出矩阵L

cout<

for(i=0;

i

i++)

{

cout<

<<

setprecision(4)<

}

}

void

Clevel::outputResult()

{

cout<<“中误差：

“<

cout<<“观测值改正数：“<

mat_c.outPut(n,1,V);

int

i;

for(i=0;

i

i++)

{

cout<<“改正后观测值h“<h[i]

=

H[i]+V[i];

}

cout<

for(i=p_k;

i

i++)

{

cout<<“改正后高程值“<

}

}

void

Clevel::calculataHBPL()

{

calH();

calB();

calP();

calL();

}

主函数main.cpp:

#include

“Clevel.h“void

main()

{

Clevel

l(“Dat.txt“);

//建立对象之后又输入文件名的数据计算结果

}

3.

运行结果

B:

1

0

0

0

1

0

1

0

0

0

1

0

-1

1

0

-1

0

1

0

0

-1

P:

0.909091

0.588235

0.434783

0.37037

0.416667

0.714286

0.384615

L:

0.0000

-0.0000

0.0040

0.0030

0.0070

0.0000

0.0020

中误差：

0.0022

观测值改正数：

-0.0002

0.0029

-0.0042

-0.0001

-0.0039

-0.0006

-0.0011

改正后观测值h0:1.3588

改正后观测值h1:2.0119

改正后观测值h2:0.3588

改正后观测值h3:1.0119

改正后观测值h4:0.6531

改正后观测值h5:0.2374

改正后观测值h6:-0.5961

改正后高程值P1:6.