測繪程序設(shè)計(jì)課程實(shí)習(xí)報(bào)告

上傳人：奔*** IP屬地：河北更新時(shí)間：2024-03-01 格式：doc 頁數(shù)：44 大小：317.50KB 人氣指數(shù)：12 舉報(bào) 版權(quán)申訴

已閱讀1頁，還剩43頁未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

1、　　測繪程序設(shè)計(jì)課程實(shí)習(xí)報(bào)告　　一、實(shí)習(xí)目的　　鞏固之前所學(xué)的c++類定義、文件操作等各方面知識，掌握矩陣計(jì)算的基本算法，熟悉并熟練運(yùn)用平差模型，掌握測量平差的基本原理和公式，并處理測量觀測數(shù)據(jù)。　　二、實(shí)習(xí)內(nèi)容<

2、/p>　　該課程的實(shí)習(xí)主要為c語言及c++編程基礎(chǔ)，矩陣操作算法，條件平差，間接平差，以及在水準(zhǔn)網(wǎng)平差中的應(yīng)用。　　三、實(shí)習(xí)過程　　實(shí)習(xí)1（文件的基本操作以及角度、弧度的轉(zhuǎn)化）　　準(zhǔn)備&

3、lt;p>　　C++文件操作函數(shù)包含于頭文件<fstream><ofstream><ifstream>，一般情況下使用<fstream>即可。在計(jì)算過程中需要文件輸出格式控制，函數(shù)包含于<iomainip>。c++的流操作函數(shù)包含于<iostream>，相應(yīng)的流函數(shù)在命名空間“std”。　　角度轉(zhuǎn)弧度的過程中，首先判斷

4、正負(fù)，對符號進(jìn)行不同情況的處理。將度數(shù)的分，秒位分別除以60、3600，然后把得到的結(jié)果與度位絕對值相加，判斷符號。　　弧度轉(zhuǎn)角度的過程中，同樣判斷正負(fù)，取小數(shù)位乘以60取整數(shù)得到相應(yīng)角度的分，再取剩余的小數(shù)乘以60即得到相應(yīng)角度的秒。　　代碼　　#include <i

5、ostream>　　#include <math.h>　　#include <fstream>　　#include <string>　　using namespace std;　　cons

6、t double pi = 3.14159;　　int main()　　{　　fstream file("a.txt"); //定義文件操作對象　　string aa;&

7、lt;b>　　file>>aa;　　cout<<aa<<endl; //輸出文件中的內(nèi)容　　int a = 1;//定義重復(fù)輸入使用的變量　　while(a) //限制循環(huán)，重復(fù)輸入<b&

8、gt;　　{　　float a1,a2,a3;　　cout<<endl;　　double deg;　　cout<<"輸入度"<<endl;

9、;　　cin>>a1;　　cout<<"輸入分"<<endl;　　cin>>a2;　　cout<<"輸入秒"<<endl;

10、　　cin>>a3;　　cout<<"角度為 "<<a1<<"° "<<a2<<"′ "<<a3<<"″"<<endl;<

11、;/p>　　if(a1 >= 0) //判斷輸入度數(shù)的正負(fù)　　deg = a1 + a2/60 + a3/3600;　　else　　deg = a1 - a2/60 - a3/3600;<p&g

12、t;　　double radian;　　radian = deg/180*pi;　　cout<<"轉(zhuǎn)換得弧度為 "<<radian<<endl;　　cout<<endl;　　cout<&l

13、t;"輸入弧度： "<<endl;　　cin>>radian;　　bool tm = 0;　　deg = radian*180/pi;　　if(deg<0)　　　　<b&

14、gt;　　{　　tm = 1;　　deg = -deg;　　}　　a1 = int(deg);　　a2 = int( (deg-a1

15、)*60 );　　a3 = ( (deg-a1-a2/60)*3600 );　　if(tm == 1)　　{　　a1 = -a1;

16、tm = 0;　　}　　cout<<"轉(zhuǎn)換后角度為 "<<a1<<"° "<<a2<<"′ "<<a3<<"″"<

17、<endl;　　　　　　}　　return 0;　　}　　運(yùn)行結(jié)果&l

18、t;b>　　輸入度　　20　　輸入分　　55　　輸入秒

19、;　　32　　角度為 20° 55′ 32″　　轉(zhuǎn)換得弧度為 0.36522　　輸入弧度：　　3.14&l

20、t;p>　　轉(zhuǎn)化后角度為 179° 54′ 32.0387″　　Press any key to continue　　出現(xiàn)的問題　　角度轉(zhuǎn)換過程中需要考慮角度的正負(fù)號問題，角度轉(zhuǎn)化可以添加是否超過59的判斷。<

21、;b>　　實(shí)習(xí)2（概率計(jì)算）　　準(zhǔn)備　　在誤差檢驗(yàn)、粗差探測等平差計(jì)算中經(jīng)常需要計(jì)算某些特定分布函數(shù)的反函數(shù)值，為滿足計(jì)算需要，使用CProbability類對相關(guān)計(jì)算進(jìn)行封裝。　　代碼及運(yùn)行結(jié)果與書上相同。

22、　　實(shí)習(xí)3（矩陣基本運(yùn)算）　　準(zhǔn)備　　矩陣的運(yùn)算有加法、減法、乘法、轉(zhuǎn)置、求逆，封裝在一個(gè)類中。為了方便之后的計(jì)算與使用，類中還添加了從文件中讀入矩陣、輸出矩陣，向控制臺(tái)輸入輸出矩陣的函數(shù)。　　代入運(yùn)算的矩陣均用動(dòng)態(tài)的一維數(shù)組表示。先定義double型指針，在確定矩陣大小后動(dòng)

23、態(tài)申請空間，row行col列的矩陣定義為：double *mat = new double[row*col]; 第i行j列的元素為mat[i*col+j]。　　矩陣的加法、減法只需遍歷兩個(gè)矩陣的每個(gè)元素，相加之后放到結(jié)果矩陣中（如有必要先判斷相加的兩矩陣行數(shù)與列數(shù)是否相等）。　　乘法計(jì)算中設(shè)相成的兩個(gè)矩陣mat1、mat2的行列分別為row1、col,

24、col、col2，則所求結(jié)果矩陣第i行j列的元素用三個(gè)循環(huán)即可計(jì)算矩陣中所有的元素。　　矩陣轉(zhuǎn)置將結(jié)果矩陣中的i行j列元素賦值為需要計(jì)算的矩陣中的j行i列個(gè)元素即可。　　矩陣求逆，由于平差計(jì)算中需要求逆的矩陣一般為實(shí)對稱正定矩陣，為了提高大型矩陣的計(jì)算效率，可以使用對實(shí)對稱正定矩陣求逆的方法。　　設(shè)C為n階

25、對稱正定矩陣，y、x為n維向量，其關(guān)系式為：　　y=C·x （3-1）　　確定了上的一個(gè)映像，如能寫出逆關(guān)系：　　x=B·y （3-2）　　則B為C的逆陣，即B=

26、　　現(xiàn)將式（3-1）寫成　　...... 　　...... （3-3）　　.......　　......<

27、p>　　因C對稱正定，必有0,用除式（3-3）第一個(gè)方程的兩端，解出，把和的位置交換，并將代入其他各式得　　...... 　　...... （3-4）　　......

28、 　　...... &

29、lt;p>　　事實(shí)上式（3-4）可以改寫為　　...... +　　...... +　　........ （3-5）　　......

30、+　　如果對式（3-5）中的變量按如下規(guī)則重新編號　　（3-6）　　經(jīng)n次變換后恢復(fù)原狀，采用變量循環(huán)重新編號法的計(jì)算公式如下：　　對于

31、;　　由于變量循環(huán)重新編號法求逆均在下三角陣（包括對角元素）進(jìn)行，因而運(yùn)行速快，在大型矩陣的計(jì)算中可以很好的適用。　　矩陣運(yùn)算的函數(shù)封裝在類mat中，函數(shù)在mat.h中定義，在mat.cpp中實(shí)現(xiàn)。　　代碼　　文件mat.h：</p

32、>　　#include <fstream>　　#include <iostream>　　#include<math.h>　　using namespace std;　　class mat&l

33、t;/b>　　{　　public:　　void add(int row,int col,double m[],double n[],double result[]);//加法　　vo

34、id sub(int row,int col,double m[],double n[],double result[]);//減法　　double* multiply(int row, int col, int col2, double m[],double n[]);//乘法

35、　　void trans(int row, int col, double m[], double result[]);//轉(zhuǎn)置　　bool inverse(double a[], int n);//求逆　　void fileIn(ifstream &in, int size, double *mat);//從文件中讀入&

36、lt;p>　　void fileOut(ofstream out, int row, int col, double *mat);　　//向文件輸出　　void outPut(int row, int col, double *mat);//向控制臺(tái)輸出　　void i

37、nPut(int row, int col, double *mat);//從控制臺(tái)輸入　　int ij(int i, int j)//下三角矩陣元素判斷　　{　　return (i>=j)? i*(i+1)/2+j : j*(j+1)/2+i;

38、　　}　　};　　文件mat.cpp：　　#include "mat.h"　　void mat::add(int row,int col,double m[]

39、,double n[],double result[])　　{　　for(int i=0; i<row*col; i++)　　result[i] = m[i]+n[i];　　}</p&

40、gt;　　void mat::sub(int row,int col,double m[],double n[],double result[])　　{　　for(int i=0; i<row*col; i++)　　result[i] = m[i]

41、-n[i];　　}　　double* mat::multiply(int row, int col, int col2, double m[],double n[])　　{　　int i,j,k;

42、　　double *result = new double[row*col2];　　for(i=0; i<row; i++)　　{　　for(j=0; j<col2; j++)&l

43、t;b>　　{　　result[i*col2+j] = 0;　　for(k=0; k<col; k++)　　result[i*col2+j] += m[i*col+k] * n[k*col2+j]; 　　}</b&g

44、t;　　}　　return result;　　}　　void mat::trans(int row, int col, double m[], double result[])

45、　　{　　int i,j;　　double *temp = new double[row*col];　　for(i=0; i<col; i++)　　{

46、　　for(j=0; j<row; j++)　　{　　result[i*row+j] = m[j*col+i];　　}<b

47、>　　}　　}　　bool mat::inverse(double mat[], int n)　　{　　int i;</p&g

48、t;　　int j;　　int k;　　double *c = new double[n];　　for(k=0; k<n; k++)　　{</b

49、>　　double c11 = mat[0];　　if(c11 + 1.0 == 1.0)　　{　　delete []c;　　return false;<p

50、>　　}　　for(i=1; i<n; i++)　　{　　double ai0 = mat[i*(i+1)/2];　　if(i<=n-k-1)<p&

51、gt;　　c[i] = -ai0/c11;　　else　　c[i] = ai0/c11;　　for(j=1; j<=i; j++)　　{　　mat[

52、(i-1)*i/2 + j -1] = mat[i*(i+1)/2+j] + ai0*c[j];　　}　　}　　for(i=1; i<n; i++)　　{<

53、;/b>　　mat[(n-1)*n/2+i-1] = c[i];　　}　　mat[n*(n+1)/2-1] = 1.0/c11;　　}　　delet

54、e []c;　　return true;　　}　　void mat::fileIn(ifstream &in, int size, double *mat)　　{&

55、lt;p>　　for(int i=0; i<size; i++)　　in>>mat[i];　　}　　void mat::fileOut(ofstream out, int row, int col, double *mat)&l

56、t;p>　　{　　for(int i=0; i<row; i++)　　{　　for(int j=0; j<col; j++)　　out<<mat[i*col+j]&l

57、t;<" "<<endl;　　out<<endl;　　}　　}　　void mat::outPut(int row, int col, doub

58、le *mat)　　{　　int i,j;　　for(i=0; i<row; i++)　　{　　for(j=0

59、; j<col; j++)　　cout<<mat[i*col+j]<<" ";　　cout<<endl;　　}　　}</p&g

60、t;　　void mat::inPut(int row, int col, double *mat)　　{　　int i,j;　　for(i=0; i<row; i++)

61、;　　for(j=0; j<col; j++)　　cin>>mat[i*col+j];　　}　　測試文件main.cpp:　　#include "mat.h"<p&g

62、t;　　mat g_mat;　　void main()　　{　　double *mat1;　　double *mat2;　　double *mat3;<p&g

63、t;　　double *mat4;　　double *mat_result;　　mat1 = new double[2*3];　　mat2 = new double[2*3];　　mat3 = new double[3*1];

64、mat4 = new double[2*2];　　cout<<"從文件中讀取"<<endl;　　ifstream in("in.txt");　　g_mat.fileIn(in,2*3,mat1);　　cou

65、t<<"矩陣mat1為："<<endl;　　g_mat.outPut(2,3,mat1);　　cout<<endl;　　g_mat.fileIn(in,2*3,mat2);　　cout<<"矩

66、陣mat2為："<<endl;　　g_mat.outPut(2,3,mat1);　　cout<<endl;　　g_mat.fileIn(in,3*1,mat3);　　cout<<"矩陣mat3為："<

67、<endl;　　g_mat.outPut(3,1,mat1);　　cout<<endl;　　g_mat.fileIn(in,2*2,mat4);　　cout<<"矩陣mat4為："<<endl;</p&g

68、t;　　g_mat.outPut(2,2,mat1);　　cout<<endl;　　in.close();　　cout<<"mat1 + mat2 = "<<endl;　　mat_result

69、= new double[2*3];　　g_mat.add (2, 3, mat1, mat2, mat_result);　　g_mat.fileIn(in, 2*3, mat_result);　　g_mat.outPut(2, 3, mat_result);　　cout

70、<<"mat1 - mat2 = "<<endl;　　g_mat.sub (2, 3, mat1, mat2, mat_result);　　g_mat.fileIn(in, 2*3, mat_result);　　g_mat.outPut(2, 3, mat_result

71、);　　cout<<"mat1 * mat3 = "<<endl;　　delete mat_result;　　mat_result = g_mat.multiply (2, 3, 1, mat1, mat3);　　g_mat.f

72、ileIn(in, 2*3, mat_result);　　g_mat.outPut(2, 1, mat_result);　　cout<<"inv(mat4) = "<<endl;　　g_mat.inverse (mat4, 2);

73、;　　g_mat.outPut (2, 2, mat4);　　cout<<"轉(zhuǎn)置 mat3 為"<<endl;　　delete mat_result;　　mat_result = new double[1*3];　　g_mat.

74、trans (3, 1, mat3, mat_result);　　g_mat.outPut(1, 3, mat_result);　　delete []mat1;　　delete []mat2;　　delete []mat3;　　d

75、elete []mat_result;　　}　　運(yùn)行結(jié)果　　測試文件“in.txt”：　　1 2 3

76、;　　4 5 6　　2 3 4　　5 6 7　　1　　2&

77、lt;b>　　3　　1 3　　4 5　　測試結(jié)果：　　從文件中讀取&l

78、t;p>　　矩陣mat1為：　　1 2 3　　4 5 6　　矩陣mat2為：　　1 2

79、 3　　4 5 6　　矩陣mat3為：　　1　　2&

80、lt;b>　　3　　矩陣mat4為：　　1 2　　3 4　　mat1 + mat2 =

81、　　3 5 7　　9 11 13　　mat1 - mat2 =　　-1 -1 -1　　-1 -1 -1　　mat1 * mat3 =<

82、;/p>　　14　　32　　inv(mat4) =　　-0.8 0.6　　-0.2 5&

83、lt;/p>　　轉(zhuǎn)置 mat3 為　　1 2 3　　Press any key to continue　　出現(xiàn)的問題　　矩陣的

84、輸入與輸出都要與一維數(shù)組相對應(yīng)，在用動(dòng)態(tài)指針操作內(nèi)存容易出現(xiàn)數(shù)組越界的問題，直接造成程序崩潰，可以在關(guān)鍵的代碼上使用try……catch排錯(cuò)。定義的類可以將函數(shù)體放在類中形成內(nèi)聯(lián)函數(shù)以增加程序的效率。函數(shù)除了使用類外還可以考慮使用一個(gè)命名空間對函數(shù)進(jìn)行封裝，或是在類中封裝矩陣的其他元素，重載操作符元素（在簡單的運(yùn)算中比較方便，在較為復(fù)雜的運(yùn)算中不太適用，比如用下三角存放的矩陣）。

85、;　　實(shí)習(xí)3（間接平差）　　準(zhǔn)備　　間接平差法（參數(shù)平差法）是通過選定t個(gè)與觀測值有一定關(guān)系的獨(dú)立未知量作為參數(shù)，將每個(gè)觀測值都分別表達(dá)成這t個(gè)參數(shù)的函數(shù)，建立函數(shù)模型，按最小二乘原理，用求自由極值的方法解出參數(shù)的最或然值，從而求得各觀測值的平差值。　　設(shè)

86、誤差方程 V = Bx-l 其中V為觀測值改正數(shù)，A為系數(shù)矩陣，x為參數(shù)向量，L為觀測值向量，則x的最小二乘解　　其中P為觀測值權(quán)陣　　單位權(quán)中誤差　　在間接平差的計(jì)算函數(shù)中，需要輸入已知參數(shù)：

87、　　n（觀測值數(shù)）　　t（參數(shù)個(gè)數(shù)）　　B（誤差方程系數(shù)矩陣，大小n*t）　　L（觀測值向量，大小n*1）　　P（觀測值權(quán)陣簡化為對角數(shù)組）<b&g

88、t;　　需要求解的參數(shù)：　　X（參數(shù)平差值向量，大小t*1）　　V（觀測值平差值向量，大小n*1）　　u（單位權(quán)中誤差）　　計(jì)算過程：

89、　　通過B、P、L計(jì)算、　　求的逆矩陣　　計(jì)算平差值 　　將X代入計(jì)算 V = Bx-l 　　代入計(jì)算 　　計(jì)

90、算函數(shù)void adjust(int n, int t, double V[], double B[], double X[], double L[], double P[],double &u)　　代碼　　文件“mat.cpp” “mat.h”同上；　　文件“mai

91、n.cpp”:　　#include"mat.h"　　#include<iomanip>　　mat mat_c;　　/**　　V = BX-L;

92、　　V: n*1　　B: n*t　　X: t*1　　L: n*1<p&

93、gt;　　P: n*1　　已知：B、L、P　　求：X、V、　　*/　　void adjust(int n, int t,

94、 double V[], double B[], double X[], double L[], double P[],double &u)　　{　　int i,j,k;　　double *BTPB = new double[(t+1)*t/2];

95、　　double *BTPB2 = new double[(t+1)*t/2];　　for(i=0; i<(t+1)*t/2; i++)　　BTPB2[i] = 0;　　cout<<"BTPB2"<<endl;<

96、p>　　for(i=0; i<t; i++)　　{　　for(j=0; j<=i; j++)　　{　　for(k=0; k<n; k++)<p

97、>　　BTPB2[i*(i+1)/2+j] += B[k*t+i]*P[k]*B[k*t+j];　　cout<<BTPB2[i*(i+1)/2+j];　　}　　cout<<endl;<

98、;b>　　}　　mat_c.inverse(BTPB2,t);　　double *BTPL = new double[t*1];　　for(i=0; i<t; i++)　　{&l

99、t;p>　　BTPL[i] = 0;　　for(k=0; k<n; k++)　　{　　BTPL[i] += B[k*t+i]*P[k]*L[k];　　}<

100、p>　　}　　for(i=0; i<t; i++)　　{　　X[i] = 0;　　for(j=0; j<t; j++)

101、　　X[i] += BTPB2[mat_c.ij(i,j)]*BTPL[j];　　}　　//X = mat_c.multiply(t,t,1,BTPB2,BTPL); 下三角矩陣不能用一般乘法　　double *BX;

102、　　BX = mat_c.multiply(n,t,1,B,X);　　mat_c.sub(n,1,BX,L,V);// V = BX-L 　　u = 0;　　for(i=0; i<n; i++)<b&g

103、t;　　{　　u += V[i]*P[i]*V[i];　　}　　u = sqrt(u/(n-t));　　}　　void mai

104、n()　　{　　double* B = new double[3*6];　　double* L = new double[6*1];　　double* P = new double[6];

105、double* X = new double[3];　　double* V = new double[6];　　double u = 0;　　ifstream file("filein.txt");　　mat_c.fileIn(file,3*6,B);

106、　　mat_c.fileIn(file,6,L);　　mat_c.fileIn(file,6,P);　　int n,t;　　file>>n;&

107、lt;b>　　file>>t;　　cout<<setiosflags(ios::fixed) << setprecision(6) ;　　adjust(n,t,V,B,X,L,P,u);　　cout<<u<<endl;</

108、p>　　mat_c.outPut(6,1,V);　　double H[] = {0.023,1.114,1.142,0.079,0.099,1.21};　　for(int i=0; i<6; i++)　　{<p&

109、gt;　　cout<<"改正后H"<<i<<":"<<H[i]+V[i]<<endl;　　H[i] = H[i]+V[i];　　}　　cout<<"A:

110、"<<10+H[3]<<endl;　　cout<<"B: "<<10+H[4]<<endl;　　cout<<"D: "<<10+H[5]<<endl;

111、　　}　　測試文件數(shù)據(jù)：　　-1 1 0 　　0 -1 1　　-1 0 1　　1 0 0　　0

112、 1 0　　0 0 1　　0.003　　0.003　　0.011　　0</

113、b>　　0　　0　　2　　2　　2

114、　　1　　1　　1　　6 3　　運(yùn)行結(jié)果

115、　　BTPB2　　5.000000　　-2.0000005.000000　　-2.000000-2.0000005.000000　　0.

116、004243　　0.001000　　0.001000　　-0.003000　　-0.004000&l

117、t;/p>　　-0.000000　　0.004000　　改正后H0:0.024000　　改正后H1:1.115000　　改正后H2:1.139000

118、　　改正后H3:0.075000　　改正后H4:0.099000　　改正后H5:1.214000　　A: 10.075000　　B: 10.099000　　D: 11.214000&l

119、t;p>　　Press any key to continue　　出現(xiàn)的問題　　由于double類型計(jì)算精度問題，輸出時(shí)容易出現(xiàn)很趨近于0的小數(shù)，需要控制輸出流格式，包含頭文件<iomainip>，使用　　cout<<setiosflags(io

120、s::fixed) << setprecision(6) ;　　最后的數(shù)字為控制小數(shù)輸出的位數(shù)。　　為了提高計(jì)算效率，使用下三角矩陣存儲(chǔ)，使用向量存儲(chǔ)權(quán)陣，使用這種方法時(shí)，之前所用的矩陣類中乘法不能使用，需要直接用for循環(huán)進(jìn)行計(jì)算。　　實(shí)習(xí)4（條件平差）<

121、/p>　　準(zhǔn)備　　條件平差是根據(jù)各觀測元素間所構(gòu)成的幾何條件以及起始數(shù)據(jù)間的強(qiáng)制條件，按最小二乘法的原理求得各觀測值的最或然值，以消除由于多次觀測產(chǎn)生的矛盾的平差方法。　　條件平差模型為其中V為觀測值改正數(shù)，A為系數(shù)矩陣，W為條件方程自由項(xiàng)向量，則由最小二乘原理：

122、　　其中P為觀測值權(quán)陣，K為拉格朗日算法中聯(lián)系數(shù)向量，解為：　　改正數(shù)向量V的解為：　　單位權(quán)中誤差　　在條件平差的計(jì)算函數(shù)中，需要輸入已知參數(shù)：　　n（觀測值數(shù)）</p

123、>　　r（條件方程個(gè)數(shù)）　　A（條件方程系數(shù)矩陣，大小r*n）　　W（條件方程自由項(xiàng)向量，大小r*1）　　P（觀測值權(quán)陣簡化為對角數(shù)組）　　需要求解的參數(shù)：</p

124、>　　V（觀測值平差值向量，大小n*1）　　u（單位權(quán)中誤差）　　計(jì)算過程：　　通過A、P計(jì)算　　求的逆矩陣&l

125、t;/b>　　計(jì)算平聯(lián)系數(shù)向量 　　計(jì)算單位權(quán)中誤差 　　代碼　　文件“mat.cpp” “mat.h”同上；

126、　　文件“main.cpp”:　　#include <iostream>　　#include <cmath>　　#include <iomanip>　　#include "mat.h"<

127、p>　　using namespace std;　　mat mat_cul;　　void condition(int n, int r, double A[], double W[],　　double P[], double V[],double &u)&l

128、t;b>　　{　　int i,j,k;　　double *AQAT = new double[r*(r+1)/2];　　//法方程系數(shù)陣　　for(i=0; i<r; i++)&

129、lt;p>　　{　　for(j=0; j<=i; j++)　　{　　AQAT[mat_cul.ij(i,j)] = 0;　　for(k=0; k<n; k++)

130、;　　AQAT[mat_cul.ij(i,j)] += A[i*n+k]*A[j*n+k]/P[k]; 　　}　　}　　//系數(shù)陣求逆

131、　　if(mat_cul.inverse(AQAT,r) == false)　　return ;　　double *K = new double [r];　　for(i=0; i<r; i++)<b&g

132、t;　　{　　K[i] = 0.0;　　for(j=0; j<r; j++)　　K[i] -= AQAT[mat_cul.ij(i,j)]*W[j];　　}<b

133、>　　//計(jì)算V　　double pvv = 0;　　for(i=0; i<n; i++)　　{　　double vi =0.0;　　for(int j=0;

134、j<r; j++)　　vi += A[j*n+i]*K[j];　　vi /= P[i];　　V[i] = vi;　　pvv+= vi*vi*P[i];　　}<

135、;p>　　delete []AQAT;　　delete []K;　　u = sqrt(pvv/r);　　}　　void main()　　{</

136、p>　　double V[6];　　double A[] = {1,0,0,1,-1,0,　　0,0,-1,-1,0,1,　　0,1,0,0,1,-1};　　double W[] = {0.003,-0.011,0.003};

137、;　　double P[] = {2,2,2,1,1,1};　　double H[] = {0.023,1.114,1.142,0.079,0.099,1.21};　　double u;　　condition(6,3,A,W,P,V,u);</p&

138、gt;　　cout<<"中誤差： "<<u<<endl;　　cout<<setiosflags(ios::fixed) << setprecision(4) ;　　int i;<p&g

139、t;　　for(i=0; i<6; i++)　　cout<<V[i]<<endl;　　for(i=0; i<6; i++)　　{　　cout<<"改正后H"<&

140、lt;i<<":"<<H[i]+V[i]<<endl;　　H[i] = H[i]+V[i];　　}　　cout<<"A: "<<10+H[3]<<endl;<

141、/p>　　cout<<"B: "<<10+H[4]<<endl;　　cout<<"D: "<<10+H[5]<<endl;　　}<

142、b>　　運(yùn)行結(jié)果　　中誤差： 0.00424264　　0.0010　　0.0010　　-0.0030<

143、;p>　　-0.0040　　-0.0000　　0.0040　　改正后H0:0.0240　　改正后H1:1.1150<p&g

144、t;　　改正后H2:1.1390　　改正后H3:0.0750　　改正后H4:0.0990　　改正后H5:1.2140　　A: 10.0750　　B: 10.0990　　D: 11.2140

145、　　Press any key to continue　　出現(xiàn)的問題　　和之前間接平差的問題類似，主要為輸出流控制與矩陣計(jì)算方面的問題。　　實(shí)習(xí)5（水準(zhǔn)網(wǎng)平差）　　準(zhǔn)備<

146、;/b>　　利用前面的間接平差函數(shù)，通過已知的數(shù)據(jù)求出平差計(jì)算所需要的　　n（觀測值數(shù)）　　t（參數(shù)個(gè)數(shù)）　　B（誤差方程系數(shù)矩陣，大小n*t）

147、L（觀測值向量，大小n*1）　　P（觀測值權(quán)陣簡化為對角數(shù)組）　　水準(zhǔn)網(wǎng)平差用類Clevel進(jìn)行封裝，類中主要包含計(jì)算所需的矩陣和數(shù)據(jù)，以及進(jìn)行計(jì)算的函數(shù)。通過類的封裝在建立對象時(shí)輸入水準(zhǔn)網(wǎng)數(shù)據(jù)所在文件名，并計(jì)算后輸出相應(yīng)的結(jié)果。　　類中的函數(shù)成員如下：　　構(gòu)造函數(shù)Cle

148、vel(char *fileName) 輸入文件的名稱，并調(diào)用其它的計(jì)算函數(shù)，輸出平差結(jié)果。　　void inputdata(char *fileName);//從文件中輸入讀入計(jì)算所需的數(shù)據(jù)，存到相應(yīng)的數(shù)組和變量中。　　首先從文件中讀取觀測值n、總點(diǎn)數(shù)pNum、已知點(diǎn)數(shù)p_k、精度m。　　從中得到t = pN

眾賞文庫> 全部分類> 畢業(yè)設(shè)計(jì)

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲(chǔ)空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

測繪程序設(shè)計(jì)課程實(shí)習(xí)報(bào)告

文檔簡介

溫馨提示

最新文檔

評論

測繪程序設(shè)計(jì)課程實(shí)習(xí)報(bào)告

文檔簡介

溫馨提示

最新文檔

評論

免費(fèi)下載