《误差理论课程设计指导书》DOC

1、误差理论与测量平差基础课程设计编写：张建霞闽江学院地理科学系2015 年 3 月适用专业：测绘工程（6、7）学时数：1周一、课程设计的性质、目的和任务误差理论与测量平差课程设计是完成测量平差基础课程教学后进行的 综合应用该课程基本知识和技能的一个教学环节，通过课程设计培养学生解决生 产实际问题的能力和所学基本知识的综合应用能力。二、课程设计的基本要求1. 巩固和加深课堂所学理论知识，培养学生理论联系实际、实际动手能力；2. 熟练掌握平差编程语言的使用；3. 掌握利用条件平差和间接平差进行平差实例计算。三、课程设计主要内容与安排利用VB可视化编程软件并结合Matlab软件实现平差中条件平差和间接

2、平差 的程序编写，并通过水准测量实例验证所编写程序的正确性。基本内容与要求如下：1. 条件平差原理与计算步骤；2间接平差原理与计算步骤；3.条件平差的算法流程；4间接平差的算法流程；5. 条件平差与间接平差编程实现6. 水准平差+GPSF差实例验证程序的正确性。四、进度安排序号设计主要内容学时分配1VB语言熟悉1天2Matlab编程熟悉1天3条件平差与间接平差理论与算法学习1天4水准网平差编程实现（条件平差与间接平差）+GPS平差1天5设计报告编写与成果制作1天五、课程设计所需的仪器设备1、具备容纳 60 人以上的实验机房2、每台电脑上需配备 VB6.0 软件以及 Matlab 软件六、课程设

3、计考核方式成绩按百分制记载。评定学生成绩主要依据以下几项：1. 平时成绩。主要包括：出勤率、态度、是否遵守学校及教师所规定的各 项纪律等。该项总分为 30 分。2. 报告成绩。主要包括： 报告的编写格式和内容是否符合要求， 分析问题、 解决问题的能力及有无独特见解， 是否工整清晰， 数据计算是否正确等。 该项总分为 40 分。3. 考核成绩。使用软件的编程能力。主要包括：对理论知识的掌握程度， 将所学理论用计算机编程语言加以编程实现。该项总分为 30 分。七、参考资料误差理论与测量平差基础 ，武汉大学出版社， 2003附录一：平差课程设计主要操作内容、vb编程学习1、启动Visual Basi

4、c6.0创建一个“标准EXE”类型的应用程序，要求:在屏幕上显示“欢迎学习Visual BasiC，并在“请输入你的姓名”标签后的文 本框Textl中输入姓名；单击“你输入的姓名是”按钮，在Label3标签显示在文本框Textl输入的姓名。程序运行效果如图示。程序以 Sy1-1.frm和 Sy1-1.vbp为文件名保存在外存上。属性设置：Label1的字体：隶书，字号：二号。Label2的字体：楷体，有下 划线。Label3有边框线。其它参见图。2、创建一个工程，由3个窗体组成。Form1用于输入用户名和口令， 当输入正确时单击“确定”按钮显示 Form2，当输入出错时则显示 Form3， 单

5、击“退出”按钮结束程序；Form2中用标签显示系统日期和时间（用 now 函数），单击“返回”按钮回到Form1； Form3为退出窗体，单击窗体则结束 程序。3、设计一个窗体，输入两个数求它们的和、差、积、商1】上机验证并分析结果： 1、设 a=2, b=3, c=4, d=5，表达式 ab and cc 的值是o2、设 a=2, b=3, c=4, d=5，表达式 32*b or a=c and bc or cd 的值是3、设 a=2, b=3, c=4, d=5，表达式 not a=c or 4*c=bA2 and ba+c 的值是4、执行语句Print chr$(65)输出结果是 5、

6、执行下列语句, a$=”Good”输出结果是 b$=”Morning print a$+b$ print a$&b$6、执行下列语句,写出结果printprint“abcd” & 1234, ”abcd”print“12” & 34, ”12” +34print 12 & 34,12+34print“ab” & true,print“12” & false,print 12 & false,12 +falseabcd” &“a34d”，+” 1234” 输出结果是+1234输出结果是输出结果是 。输出结果是。7、执行下列语句, s$=” ABCDEFGHIJ”Kprint left$(s$,4

7、)print right$(s$,4)ab” +true输出结果是 。”12” + false 输出结果是 输出结果是 。输出结果是8、执行下列语句, s$=” ABCDEFGHIJ”K输出结果是print mid$(s$,3,4)print len(s$)9、执行下列语句, s$=” ABCDEFGHIJ”Kprint instr(s$,” efg ”)print lcase$(s$)输 出 结 果是 10、执行下列语句， print day(now)print month(now)print year(now)print weekday(now) 输出结果是 11、执行语句 print i

8、nt(12345. 6789*100+0.5)/100 的输出结果是 。12、执行语句 print int(-12345.6789*100+0.5)/100 的输出结果是 。13、执行语句 print fix(12345 . 6789*100+0.5)/100的输出结果是 。14、执行语句 print fix(-12345 . 6789*100+0.5)/100 的输出结果是 。15、执行下列语句，写出输出结果：Print Abs(-3.9),Fix(-3.9),Int(-3.9),Round(-3.9)输出结果： 、Print Exp(5),Log(5)输出结果： Print Sqr(256

9、)输出结果： 。Print Asc( ” d” ),Chr(100),Asc( “ ASD”) 输出结果： Print Ucase( “ aBcDef ” ),Lcase( “ aBcDef ” )输出结果： 、Print Str(12345),Val(“ 12.3AB4-51”23).,4S5t)r,(Val( “ asd123.5 ”)输出结果： 、Print Len（“ 12学创！” ）输出结果： 。Print Replace（我是教师你是学生”教师”学生”）输出结果：。Print String（9, 我喜”欢 VB” ）输出结果： 。【2】编写简单应用程序1 、每单击窗体一次，就能产生

10、 1 个 100 到 200之间的随机整数，并在窗体上输出。2、从键盘上输入 4 个数，编写程序，计算并输出这 4 个数的和及平均值。通过 inputbox 函数输入数据，在窗体上显示和及平均值。3、通过键盘任意输入两个数，对调后输出。【3】数组与矩阵1、任意输入 10 个数，按从大到小排列。2、任意建立一个二维数组，求数组中最大元素值和所在下标。3、任意建立二维数组，求每行元素的和。4、任意建立二维数组，求各行最大值之和。5、输出 4*4 阶矩阵的主对角线元素之和。6、计算 4*4 阶矩阵的周边元素值之和。7、录入矩阵 A(2*2) 和矩阵 B(2*3 )，求 C=A*B二、平差编程1、求矩

11、阵A(3X 4)与矩阵 B (4X6) 的乘积。2、求矩阵 A( 5X 5) 的转置。3、 以矩阵A(5X5)的下三角元素为基础，将矩阵 A变为对称阵。4、编程实现教材(P78)例题5- 2的条件平差结果，要求如下：(1) 通过inputbox函数输入已知点A B的高程值。(2) 通过文本框(text )接收系数矩阵A和协因素阵Q通过inputbox函 数录入常数项矩阵 W。(3) 对方阵Na求逆通过编程实现。( 4)界面设计如下：矩阵求逆程序示例：Dim n%, a()Private Sub Form_Load()Me.Top = 800: Me.L eft = 800End SubPriv

12、ate Sub Label1_Click()c = 1For i = 1 To nForj = 1 To na(i, j) = Val(Text1(j + (i - 1) * n)Next j, iFor i = 1 To nFor j = n + 1 To 2 * nIf i + n = j The na(i, j) = 1Elsea(i, j) = 0End IfNext j, iCall txjz(a(), n, 2 * n)Call zjjz(a(), n, 2 * n)For i = 1 To nFor j = n + 1 To 2 * nText1(j + (i - 1) * n

13、- n) = Format(a(i, j), 0.#) Next j, iEnd SubPrivate Sub Label2_Click()For Each x In Text1If x.Index 0 Then Unload xNext xsta:n = Val(InputBox( 请你输入方阵的介数 ): GoTo staIf n 10.5 Or n 0.8 Then MsgBox (请你输入不大于十的正整数Scale (0, 0)-(n + 2, n + 2)ReDim a(1 To n, 1 To 2 * n)For i = 1 To nFor j = 1 To nLoad Text1

14、(j + (i - 1) * n)Text1(j + (i - 1) * n).Visible = TrueText1(j + (i - 1) * n).Top = iText1(j + (i - 1) * n).Left = jNext j, iEnd SubSub txjz(a(), m%, n%) 将矩阵化为阶梯型矩阵i% = 1: j% = 1DoFor k = i To mIf a(k, j) 0 ThenFor l = j To nb = a(k, l)a(k, l) = a(i, l)a(i, l) = bNext lFor l = j + 1 To nFor h = i + 1

15、 To ma(h, l) = a(h, l) - a(i, l) * a(h, j) / a(i, j)Next h, lFor l = i + 1 To ma(l, j) = 0Next li = i + 1Exit For End If Next kj = j + 1Loop Until i = m Or j = n + 1End SubSub zjjz(a(), m%, n%) 将梯形矩阵化为行最简型矩阵i% = 1: j% = 1DoFor k = m To i Step -1If a(k, j) 0 Thenb = a(k, j)For l = j To na(k, l) = a(k

16、, l) / bNext lFor l = j + 1 To nFor h = 1 To k - 1a(h, l) = a(h, l) - a(h, j) * a(k, l)Next h, lFor l = 1 To k - 1 a(l, j) = 0Next lEnd IfNext kj = j + 1Loop Until i = m + 1 Or j = n + 1End Sub循环 , 统计变量Sub main()Dim ci As Byte, cj As Byte, i As Byte, j As Byte, n As ByteDim a(2, 2) As Single矩阵数组Dim

17、b(2, 2) As Single逆阵数组Dim t(3) As Single 余子式数组Dim rValue As Single 矩阵值变量 初始矩阵a(0, 0) = 1: a(0, 1) = -1: a(0, 2) = 3a(1, 0) = 2: a(1, 1) = -1: a(1, 2) = 4a(2, 0) = -1: a(2, 1) = 2: a(2, 2) = -4Debug.Print 原始矩阵为 : For i = 0 To 2For j = 0 To 2Debug.Print a(i, j),Next jDebug.PrintNext i 求伴随阵For ci = 0 To

18、 2For cj = 0 To 2 求余子式n = 0For i = 0 To 2For j = 0 To 2If i ci And j cj Then t(n) = a(i, j): n = n + 1Next jNext ib(cj, ci) = (-1) A (ci + cj) * (t(0) * t(3) - t * t(2)Next cjNext ci 求矩阵的值rValue = 0For i = 0 To 2rValue = rValue + a(i, 0) * b(0, i)Next iIf rValue Then 逆阵存在 求矩阵的逆阵Debug.Print 逆阵为 : rVa

19、lue = 1 / rValueFor i = 0 To 2For j = 0 To 2 b(i, j) = b(i, j) * rValue Debug.Print b(i, j), Next jDebug.PrintNext iElse 逆阵不存在Debug.Print 逆阵不存在 ! End IfEnd Sub三、使用 Matlab 间接平差1、使用 Matlab 实现教材 P77 例题 5-1 的平差计算2、使用 Matlab 实现教材 P157 例题 7-14 的平差计算附录二：课程设计报告编写提纲一、目的 通过课程设计，掌握平差的基本原理、主要平差方法，并学会利用所学理论与方法解决实际测量中的平差问题。二、设计主要内容+GPS平差实例验证所编利用VB可视化编程软件并结合Matlab软件实现平差中条件平差和间接平差的程序