测量平差基础课程设计.docx
《测量平差基础课程设计.docx》由会员分享,可在线阅读,更多相关《测量平差基础课程设计.docx(14页珍藏版)》请在冰点文库上搜索。
测量平差基础课程设计
测量平差基础课程设计
题目:
导线平差
班级:
测绘二班
姓名:
李文义
学号:
080615380
完成日期:
测量平差课程设计
测量平差简介
测量平差的主要任务是对一系列带有观测误差的观测值,运用概率统计的方法来消除它们之间的不符值,求出未知量的最可靠值,并且评定成果的精度。
测量平差与其他学科一样,是由于生产的需要而产生的,并在生产实践过程中,随着科学技术的进步而发展。
18世纪末,在测量学、天文测量学等实践中提出如何消除由于观测误差引起的观测值之间矛盾的问题,即如何从带有误差的观测值中找到观测值的最优值。
1794年,年仅17岁的高斯(C.FGauss)首先提除了这个问题的解决方法——最小而乘法。
他是根据偶然误差的四个特性,并以算术平均值为待求量的最或然值出发,导出了偶然误差的概率分布,给出了在最小二乘原理下求待定量最或然值的计算方法。
当时高斯没有正式发表。
19世纪初(1801年),天文学家对刚发现的谷神星运行轨道的一段弧长作了一系列观测,后来因故中止了。
这就需要根据这些带有误差的观测结果求出该星运行的实际轨道。
高斯用自己提出的最小二乘法解决了这个当时很大的难题,对谷神星运行轨道进行了预报,使天文学家及时的又找到了这颗彗星。
1809年高斯才在《天体运动的理论》一书中正式发表了他的方法。
在此之前,1806年,勒戎德尔(A.M.Legendre)发表了《决定彗星轨道的新方法》一文,从代数观点上也独立的提出了最小二乘法,并定名为最小二乘法。
所以,后人称它为高斯-勒戎德尔方法。
自19世纪初到20世纪五六十年代的一百多年来,测量平差学者在基于最小二乘原理的平差方法上做了许多研究,提出了一系列解决各类测量问题的平差方法,针对这一时期的计算工具的情况,提出了许多分组解算线性方程组的方法,达到了简化计算的目的。
自20世纪70年代开始,随着计算机技术的进步和生产实践中的高精度要求,测量平差得到了很大发展,主要表现在:
1.从单纯的偶然误差到包含系统误差和粗差。
2.1947年,铁斯特拉(T.M.Tienstra)提出了相关观测值的平差理论。
3.1969年,克拉鲁(T..Krarup)提出最小二乘虑波,也称为拟和推估法。
4.1962年,迈赛尔(P.meissl)提出针对非满秩平差问题的内制约束平差原理。
5.对系统误差的特性、传播、检验、分析的理论的发展
6..测量中粗差理论的研究。
要求:
系统(或顶层文件)结构描述,各个模块(或子程序)的功能描述;
1)用原理图实现的,报告需包含以下内容:
(1)系统原理图
(2)主要模块或符号说明;
1)数据录入:
可以用TXT文本编好过后直接进行数据传输,也可一在EXCEL中直接编写进行计算
2)数据处理:
这一模块住要在EXCEL中做出,对所有的公式进行编写一达到我们所要的计算目的。
3)成图:
是在VB图形控件中显示,即与CAD结合
2)用VB语言实现,报告包含以下内容:
(1)源代码;
DimxlappAsExcel.Application
DimxlbookAsExcel.Workbook
DimxlsheetAsExcel.Worksheet
DimaAsInteger
DimbAsInteger
DimbAsInteger
DimdAsInteger
DimeAsInteger
DimfAsInteger
a=Val(Text1.Text)
b=Val(Text5.Text)
c=Val(Text6.Text)
d=Val(Text7.Text)
e=Val(Text8.Text)
f=Val(Text9.Text)
PrivateSubCommand3_Click()
IfDir("d:
\temp\excel.bz")=""Then
Setxlapp=CreateObject("excel.application")
xlapp.Visible=True
Setxlbook=xlapp.Workbooks.open("d:
\temp\bb.xls")
Setxlsheet=xlbook.Worksheets
(1)
xlsheet.Activate
xlsheet.Cells(2,3)="abc"
xlbook.RunAutoMacros(xlAutoOpen)
Else
MsgBox("excel已经打开!
")
EndIf
EndSub
PrivateSubCommand2_Click()
IfDir("d;\temp\excel.bz")<>""Then
xlbook.RunAutoMacros(xlAutoClose)
xlbook.Close(True)
xlapp.Quit
EndIf
Setxlapp=Nothing
EndSub
PrivateSubCommand4_Click()
IfDir("d:
\temp\excel.bz")=""Then
Setxlapp=CreateObject("excel.application")
xlapp.Visible=True
Setxlbook=xlapp.Workbooks.open("d:
\temp\bb.xls")
Setxlsheet=xlbook.Worksheets
(2)
xlsheet.Activate
xlsheet.Cells(1,1)="abc"
xlbook.RunAutoMacros(xlAutoOpen)
Else
MsgBox("excel已经打开!
")
EndIf
EndSub
程序代码:
OptionBase1
PrivateSubForm_Click()
DimmAsInteger,nAsInteger
DimiAsInteger,jAsInteger
ReDimb(3,3)AsInteger
ReDimc(3,3)AsInteger
Fori=1To3
Forj=1To3
b(i,j)=CInt(Rnd*3)+2
Nextj
Nexti
Print"矩阵b:
"
Print
Fori=1To3
Forj=1To3
PrintTab(4*j);b(i,j);
Nextj
Nexti
Fori=1To3
Forj=1To3
c(i,j)=b(j,i)
Nextj
Nexti
Print
Print"b矩阵的转置:
";‘显示转置后的数组以确认代码的正确性
Print
Fori=1To3
Forj=1To3
PrintTab(4*j);c(i,j);
Nextj
Nexti
EndSub
b)矩阵乘法模块
模块代码:
OptionBase1
PrivateSubForm_Click()
DimmAsInteger
DimnAsInteger
DimiAsInteger
Dima(3,3)AsInteger
Dimb(3,3)AsInteger
Dimc(3,3)AsInteger
Form=1To3
Forn=1To3
a(m,n)=Rnd*20+8
b(m,n)=Rnd*12+7
Nextn
Nextm
Print"矩阵a:
"
Form=1To3
Forn=1To3
PrintTab(5*n);a(m,n);
Nextn
Nextm
Print
Print"矩阵b:
"
Form=1To3
Forn=1To3
PrintTab(5*n);b(m,n);
Nextn
Nextm
Form=1To3
Forn=1To3
c(m,n)=0
Fori=1To3
c(m,n)=c(m,n)+a(m,i)*b(i,n)
Nexti
Nextn
Nextm
Print
Print"矩阵c:
"
Form=1To3
Forn=1To3
PrintTab(5*n);c(m,n);
Nextn
Nextm
EndSub
c)矩阵求逆模块
程序采用了求矩阵的行列式和矩阵伴随矩阵,最后求商得到矩阵的逆的方法,由于采用了大量的循环语句和递归调用的方式,程序在执行大数据量的时候显得有些慢,例如在解算校内导线的时候,是一个8×8矩阵,并且数据量较大,因而在解算时花了好几秒的时间。
因此该方法还有待改进。
程序代码:
OptionBase1
'建立一个函数求下一个矩阵
Functionnext_a(ByRefa()AsInteger,ByValnAsInteger,ByVali
AsInteger)AsInteger
Dimk,jAsInteger
DimcountAsInteger
count=0
Fork=2Ton+1
Forj=1Ton+1
Ifj<>iThen
a(Fix(count/n)+1,(countModn)+1)=a(k,j)
count=count+1
EndIf
Nextj
Nextk
EndFunction
'建立一个函数求取矩阵的行列式值
Functiondvalue(ByRefa()AsInteger,ByValnAsInteger)As
Integer
Dimb()AsInteger
Dimo,pAsInteger
Dimi,k,jAsInteger
DimsignAsInteger:
sign=1
DimsumAsInteger:
sum=0
ReDimb(n,n)AsInteger
Ifn=1Then
sum=a(1,1)
Else:
Fori=1Ton
Forj=1Ton
Fork=1Ton
b(j,k)=a(j,k)
Nextk
Nextj
next_ab(),n-1,i'调用nexta函数
sum=sum+sign*a(1,i)*dvalue(b(),n-1)'递归调用
dvalue函数
sign=(-1)*sign
Nexti
EndIf
dvalue=sum
EndFunction
'建立一个函数求取矩阵的伴随矩阵
Functionaccom_a(ByRefa()AsInteger,ByValnAsInteger,ByVali
AsInteger,ByValjAsInteger)AsInteger
Dimt()AsInteger
Dime,bAsInteger
DimcountAsInteger:
count=0
DimsAsInteger:
s=0
ReDimt(n,n)AsInteger
Fore=1Ton
Forb=1Ton
Ife<>iAndb<>jThen
t(Fix(count/(n-1))+1,(countMod(n-1))+1)=a(e,b)
count=count+1
EndIf
Nextb
Nexte
s=dvalue(t(),n-1)
If((i+j)Mod2)<>0Then
s=(-1)*s
Else:
s=s
EndIf
accom_a=s
EndFunction
abc
闭合导线平差
度
分
秒
X坐标
Y坐标
距离
所测角度
方位角
已知方位角
第1点
第2点
第3点
第4点
第5点
第6点
第7点
第8点
第9点
第10点
第11点
第12点
第13点
第14点
x
y
S=
相对精度=
4.课程设计过程中遇到的问题以及解决方法
在vb中实现导线的平差太难了,而且编制过成也太长,我们在接到通知要交的时候还六下3天时间,所以我们主要编制导线网的平差,在导线那一部分我们使用vb与EXCEL连接计算,本来我们还要把acd连接进入直接进行成图,但是还没找到相关资料,不过我们知道以后会在这一方面努力的
首先在页面设计是我要考虑到我们要的起誓数据是什么,一开始我们并没有在意,可是后来我们在编制程序的时候发现了问题。
所以又耽误了我们好长时间重新来过,包括EXCEL的编制,主要使用起内部函数。
第二在实现编程的时候我们对于vb的好多命令有点生疏了,使得进度更加慢了我们知道在学习编程的时候我们不可以停下脚步,我们要不段的饿学习,不断的更新自己的知识
在最后编制过车工中出现了错误,晕人的挖,不知道是哪里错了,可是用的vb是下载的没有msdn不好查找,只好出现一个在网上搜一个,只怪自己技术不到家。
其实我想如果我们在做课程设计的时候有个老师在边上指导我想,那效果会更好的
5.总结
经过一个星期的上机实践学习,使我对vb语言有了更进一步的认识和了解,要想学好它要重在实践,要通过不断的上机操作才能更好地学习它,通过实践,我也发现我的好多不足之处,首先是自己在指法上还不行,经常按错字母,通过学习也有所改进;再有对vb语言的一些标准库函数不太了解,还有对函数调用的正确使用不够熟悉,还有对vb语言中经常出现的错误也不了解,通过实践,使我在这几个方面的认识有所提高。
通过实践的学习,我认到学好计算机要重视实践操作,不仅仅是学习vb语言,还是其它的语言,以及其它的计算机方面的知识都要重在实践,所以后在学习过程中,我会更加注视实践操作,使自己便好地学好计算机。
6.心得体会
通过该课程设计,对该课程的一些基本概念、词法分析、语法分析、属性文法和语法制导翻译有了深刻的了解并比较掌握了什么是平差,平差工作的基本过程及其各阶段的基本任务,熟悉了编译程序总流程框图,了解了编译程序的生成过程、构造工具及其相关的技术对课本上的知识有了更深的理解,课本上的知识机械的,表面的。
通过把该算法的内容,算法的执行顺序在计算机上实现,把书本知识更好的融入自己的头脑,对实验原理有更深的理解。
通过该课程设计,全面系统的理解了编译原理程序构造的一般原理和基本实现方法。
把学过的计算机平差原理以及相关课程的知识强化,能够把课堂上学的知识通过自己设计的程序表示出来,加深了对理论知识的理解和对自己动手能力的提高。
以前对与计算机平差系统概念上的认识是模糊的,现在通过自己动手做实验,从实践上认识了操作系统是如何处理命令的,如何协调计算机内部各个部件运行记忆是怎么运行的,对计算机平差原理的认识更加深刻。
在这次课程设计中,通过自己动手,加深理解文件系统的内部功能及内部实现,培养实践动手能力和程序开发能力的目的。
把死板的课本知识变得生动有趣,激发了学习的积极性。
参考文献:
互联网