大地主题解算-正算-C#
大地主题解算-正算-程序
using System;
using System.Collections.Generic;
using https://www.doczj.com/doc/4816095986.html,ponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
namespace WindowsFormsApplication2
{
public partial class Form1 : Form
{
public Form1()
{
InitializeComponent();
}
private void button1_Click(object sender, EventArgs e) {
int n,du,fen,miao;
double B1, L1, A1, S, dS, dB1, dL1, M1, N1, W1;
double E = 0.006694384999588, a = 6378140, C;
B1 =( Convert.ToDouble(textBox1.Text) + Convert.ToDouble(textBox9.Text)/60 + Convert.ToDouble(textBox10.Text)/3600 )* Math.PI / 180;
L1 = (Convert.ToDouble(textBox2.Text) + Convert.ToDouble(textBox11.Text)/60 + Convert.ToDouble(textBox12.Text)/3600) * Math.PI / 180; ;
A1 = (Convert.ToDouble(textBox3.Text) + Convert.ToDouble(textBox13.Text)/60 + Convert.ToDouble(textBox14.Text)/3600) * Math.PI / 180; ;
S = Convert.ToDouble(textBox5.Text);
n = Convert.ToInt16(textBox4.Text);
double[] L = new double[100000];
double[] B = new double[100000];
double[] dL = new double[100000];
double[] dB = new double[100000];
double[] A = new double[100000];
double[] W = new double[100000];
double[] N = new double[100000];
double[] M = new double[100000];
W1 = Math.Pow(1 - E * Math.Pow(Math.Sin(B1), 2), 0.5);
N1 = a / W1;
M1 = a * (1 - E) / (W1 * W1 * W1);
C = N1 * Math.Sin(A1) * Math.Cos(B1);
dS = S / n;
dB1 = dS * Math.Cos(A1) / M1;
dL1 = dS * Math.Sin(A1) / N1 / Math.Cos(B1);
B[0] = B1;
L[0] = L1;
dB[0] = dB1;
dL[0] = dL1;
A[0] = A1;
M[0] = M1;
N[0] = N1;
W[0] = W1;
int i;
for (i = 1; i < n; i++)
{
dB[i - 1] = dS * Math.Cos(A[i - 1]) / M[i - 1];
B[i] = B[i - 1] + dB[i - 1];
W[i] = Math.Pow(1 - E * Math.Pow(Math.Sin(B[i]), 2), 0.5);
N[i] = a / W[i];
M[i] = a * (1 - E) / (W[i] * W[i] * W[i]);
dL[i - 1] = dS * Math.Sin(A[i - 1]) / (N[i - 1] * Math.Cos(B[i - 1])); L[i] = L[i - 1] + dL[i - 1];
A[i] = Math.Asin(C / (Math.Cos(B[i]) * N[i]));
}
du = Convert.ToInt16(Math.Floor(B[i - 1] * 180 / Math.PI));//du
fen = Convert.ToInt16(Math.Floor((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16(((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60 - Math.Floor((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60)) * 60));//秒?
textBox6.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
du = Convert.ToInt16(Math.Floor(L[i - 1] * 180 / Math.PI));//du
fen = Convert.ToInt16(Math.Floor((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16(((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60 - Math.Floor((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60)) * 60));//秒?
textBox7.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
du = Convert.ToInt16(Math.Floor((180 + A[i - 1] * 180 / Math.PI)));//du
fen = Convert.ToInt16(Math.Floor(((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16((((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60 - Math.Floor(((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60)) * 60));//秒?
textBox8.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
}
private void Form1_Load(object sender, EventArgs e)
{
} }
}
第2章矩阵与线性代数方程组 一般的线性代数方程组，A非奇异可根据Cramer法则求解方程唯一解但是它的计算量很大。 高斯消元法的算法时间复杂度是O(n3)，可以解一系列的线性方程；所占数据空间符合原地工作的原则。但是算法对数值计算不稳定（当分母为0或很小时）。可以用在计算机中来解决数千条等式及未知数。不过，如果有过百万条等式时，这个算法会十分费时。 解决高斯法中的不稳定性，在每次归一化前增加选主元（列选主元、全选主元）过程。但是列选主元法仍不稳定，不适求解大规模线性代数方程组。全选主元的高斯消去法，则在复杂度降低的同时能够避免舍入误差，保证数值稳定性。 高斯-约当消去法算法产生出来的矩阵是一个简化行梯阵式，而不是高斯消元法中的行梯阵式。相比起高斯消元法，此算法的效率比较低，却可把方程组的解用矩阵一次过表示出来。线性代数方程组的迭代解法 简单迭代法：迭代格式发散但迭代值序列不一定发散，但收敛格式收敛，迭代值序列收敛于方程组的准确解与选取迭代初值无关。 雅可比迭代法: 计算公式简单，且计算过程中原始矩阵A始终不变，比较容易并行计算。但是收敛速度较慢，而且占据的存储空间较大，所以工程中一般不直接用雅克比迭代法，而用其改进方法。 高斯-赛德尔迭代法：较上面的迭代复杂，但是矩阵的条件相对宽松。 松弛法：需要根据经验去调整，收敛速度依赖松弛参数的选择，收敛条件的要求更宽松。共轭梯度法：是介于最速下降法与牛顿法之间的一个方法，它仅需利用一阶导数信息，但克服了最速下降法收敛慢的缺点，又避免了牛顿法需要存储和计算Hesse矩阵并求逆的缺点，在各种优化算法中，共轭梯度法是非常重要的一种。其优点是所需存储量小，具有步收敛性，稳定性高，而且不需要任何外来参数。共轭梯度法不仅是解决大型线性方程组最有用的方法之一，也是解大型非线性最优化最有效的算法之一。 第3章矩阵特征值 乘幂法计算绝对值最大的特征值：其收敛速度受限于最大与次大特征值比值绝对值的大小，实际应用中采用加速技术。 求对称特征值的雅克比方法96：每进行一次选装变换钱都需要在飞对角线的元素中选取绝对值最大的元素，很费时间，雅克比过关法对此做了改进。 QR方法求一般实矩阵的全部特征值98下100下：重复多次进行QR分解费时，计算工作量很大。一般先进行相似变换然后进行QR分解。但是这样仍然收敛速度慢，一般是线性收敛。实际应用中使用双重步QR变换将带原点的QR算法中相邻两步合并一步，加速收敛避免复数运算。 第4章非线性方程与方程组 二分法：每次运算后，区间长度减少一半，是线形收敛。优点是简单，但是不能计算复根和重根。 简单迭代法：直接的方法从原方程中隐含的求出x，从而确定迭代函数 (x)，这种迭代法收敛速度较慢，迭代次数多。 埃特金迭代法113中：对简单迭代进行改进，使在其不满足收敛条件下迭代过程也收敛，在其收敛时加快收敛速度，减少迭代次数降低时间复杂度。 牛顿迭代法：其最大优点是在方程f(x) = 0的单根附近具有平方收敛，收敛速度快。而且该法还可以用来求方程的重根、复根。缺点：初值的选择会影响收敛结果。 牛顿下山法：保证函数值稳定下降，且有牛顿法的收敛速度。
实验一多功能计算器 一、问题描述 设计一个多功能计算器，可以完成基本的计算。 设计要求： 1. 具备整型数据、浮点型数据的算术（加、减、乘、除）运算功能。依次输入第一个运算数、运算符（+，- ，* ，/ ）、第二个运算数，然后输出结果。结果可以作为下一个运算的第一运算数。按’C清 屏，按’R'返回菜单。 例如：输入：2 + 5 输出：7 2. 实现单运算符表达式计算的功能。输入的操作数可以包含整数或浮点数。输入表达式如下： 例如：输入：2+5 输出：7 二、算法说明 1．数据结构说明（可以图示说明，也可以文字说明） 本程序主要根据选择菜单编写了六个自定义函数，用于在main（）函数中调用，在main（）中，用一个字符变量num1 来记录下菜单选项的标号，根据num1 的值来决定调用哪个函数。 程序要完成的功能及所要用到的函数如下：
使用各菜单分别调用的函数 '≡ua∩ Ξ huyun≡u an O ■> SUanShUyUnSLLan2 panduan() > bi aodashiyunsuan 0 > QirLgPing () > fanhuicaidan() 菜 I Y l - ▽ 算术运算 a 卑运篡符表达 b 清屏 C 返回菜单 E 退岀程序 t F 面是整个程序的流程图:
2?算法说明(即函数说明) Void SUa nshuyu nsua n() //做算术运算时调用的函数 Void SUanshuyunsuan2() //选择继续做算术运算调用的函数，将上次运算的结果做为下次算术运算的第一个操作数〃判断算术运算是否继续 void Pan dua n() 〃判断算术运算是否继续 void biaodashiyu nsua n() 〃单运算符表达式实现函数 void qin gpi ng() 〃清除屏幕 void fan huicaida n() // 显示菜单 三、测试结果(这部分需文字与图示结合) 1 ?第一组测试用例 (1) 测试输入： (≡ΓF;I学习\匚语言程序设计实习'实验一Ilkb II丛匚语言实验一一 测试目的：测试算术运算的功能结果输出:
￥ 数值分析思考题1 1、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。 4、取 ，计算 ，下列方法中哪种最好为什么（1）(3 3-，（2）(2 7-，（3） ()3 1 3+ ，（4） ()6 1 1 ，（5）99- , 数值实验 数值实验综述：线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式，则为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。对正定对称矩阵，采用平方根方法无需选主元。方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进行求解。 数值计算方法上机题目1 1、实验1. 病态问题 实验目的： 算法有“优”与“劣”之分，问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。希望读者通过本实验对此有一个初步的体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 $ r e x x e x x ** * ** - == 141 . ≈)61
课程设计名称：C语言课程设计课程设计题目：简单计算器程序
目录 第1章需求分析 (1) 1.1设计要求 (1) 1.2任务 (1) 第2章总体设计 (2) 2.1设计简介及设计方案论述 (2) 2.2功能模块层次图 (2) 第3章详细设计 (3) 3.3由（后缀）逆波兰表达式计算中缀表达式原理 (8) 3.3.1算法描述 (8) 第4章调试分析 (10) 4.1程序设计中所遇到的错误及犯错的原因 (10) 4.2错误的解决方法 (10) 第5章用户手册 (11) 总结 (15) 参考文献 (16) 附录（程序清单） (17)
第1章需求分析 1.1 设计要求 （1）用 C 语言数据结构实现程序设计； （2）利用结构体、栈、进行相关信息处理； （2）系统的各个功能模块要求用函数的形式实现； （4）界面简单，可操作性高。 1.2任务 （1）定义一个结构体类型数组，输入0~9 及+、--、*等符号的信息，将其信息存储起来； （2）输入简单的加减乘除算术计算式，并在屏幕上显示逆波兰（后缀式）表达式和计算结果； （3）编写代码； （4）程序分析与调试。 说明： 本课程设计将实现一个简单计算器。在功能上尽量模仿windows 的计算器。系统界面不做牵制要求。该程序能实现标准型中+、-、*、/、（、）、.、的混合运算表达式（一般意义上的中缀表达式），将其转换成逆序波兰表达式（后缀表达式）并计算输出结果。在进行运算后可以选择继续运算或者结束当前运算。即时准确地获得需要的计算的结果，充分降低了数字计算的难度和节约了时间，对人们的生活有一定的帮助。
第2章 总体设计 2.1设计简介及设计方案论述 逆波兰表达式又叫做后缀表达式。在通常的表达式中，二元运算符总是置于与之相 关的两个运算对象之间，所以，这种表示法也称为中缀表达式。波兰逻辑学家 J.Lukasiewicz 于 1929 年提出了另一种表示表达式的方法。按此方法，每一运算符都置 于其运算对象之后，故称为后缀表达式。 后缀表达式的优点是显而易见的， 编译器在处理时候按照从左至右的顺序读取逆波 兰表达式，遇到运算对象直接压入堆栈，遇到运算符就从堆栈提取后进的两个对象进行计算，这个过程正好符合了计算机计算的原理。后缀表达式比前缀表达式更加易于转换，并且它的最左面一定为数字，这一点在实 际编程的时候就会体会到它的好处了。 逆波兰表达式有一个更大的优点，就是拆括号，根据运算符的级别将中缀表达式转 换成逆波兰表达式后，运算顺序就已经替代了运算符的级别，这样也避免了括号提高运 算级别的特殊处理。 2.2功能模块层次图 将算术表达式转化为逆波兰表达式 计算逆波兰表达式的值 简单计算器 表 达 式 格 式 转 换 系统 求 值 计 算 系 统
白塞尔大地主题解算 方向： 学号： 姓名：
一．基本思路： 基本思想：将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后在将球面上的计算结果换算到椭球面上。其关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式，同时解决在球面上进行大地主题解算的方法。 正算流程： 1.计算起点的归化纬度 2.计算辅助函数值，解球面三角形 3.按公式计算相关系数A,B,C 以及α,β 4.计算球面长度 5.计算纬度差改正数 6.计算终点大地坐标及大地方位角 011122S B C A {sin (cos )} σσσ=-+10101022222sin ()sin sin cos cos σσσσσσ+=+101010 22222cos ()cos cos sin sin σσσσσσ+=-001101 522B C A [cos ()]sin ()σσσσσσ=+ +++010122L A sin [(sin ()sin )] λδασβσσσ-==++-2111u u u A sin sin cos cos cos sin σσ = +222 222 22 22 2222221111u B u B W W u B u B B arctan e u sin sin cos cos tan sin tan cos ?==????=????==--1 111 A arctan u u A sin sin [] cos cos sin sin cos σλσσ=-21L L λδ =+ -
反算流程： 1.辅助计算 2.用逐次趋近法同时计算起点大地方位角、球面长度及经差，第一次趋近时，取δ＝０。 计算下式,重复上述计算过程2. 3.计算大地线长度S 4.计算反方位角 二．已知数据 L λδ =+211212u p A u u u u q sin cos tan cos sin sin cos cos λλ= = -2121p p u q b b A arctan q sin cos cos λλ==-=11 p A q A sin cos tan cos σσ += 11 p A q A sin sin cos σ=+12a a cos cos σλ =+arctan sin cos σσσ?? = ??? 011 A u A sin cos sin =111u A tan tan sec σ=21+σσσ =02122L A sin [(sin sin )] λδασβσσ-==+-L λδ =+11222222S A B C B C sin (cos )sin (cos ) σσσσσ=++-+
淮阴工学院 《数值分析》考试 ──基于Matlab的方法综合应用报告 班级：金融1121 姓名：姚婷婷 学号：1124104129 成绩： 数理学院 2014年6月7日
《数值分析》课程综述报告 前言： 数值分析也称计算方法，它与计算工具的发展密切相关。数值分析是一门为科学计算提供必需的理论基础和有效、实用方法的数学课程，它的任务是研究求解各类数学问题的数值方法和有关的理论。 正文： 第一章 近似计算与误差分析 1、产生误差的原因：①模型误差；②观测误差；③截断误差；④舍入误差。 2、四则运算的误差： ①加减法运算 ()()()****x y x y δδδ±=+ ②乘法运算 ()()() ****** *** ******xy x y xy xy xy x y x y y y x x x y x y y x δδδ-=-+-≤-+-?=+ ③ 除法运算： ()()() () () ***** ******* * * ** * * ** * *2 ** x x xy x y y y yy xy x y x y x y yy x x y y y x yy x y y x x y y δδ δ--=-+-=-+-= +?? ?≈ ??? 3、科学表示法、有效数字、近似值的精度 任何一个实数都可以表示成如下的形式： 其中：是正整数，是整数， 如果是数的近似值 并且 则称该近似值具有位有效数字（significant digit ）。
此时，该近似值的相对误差为 另一方面，若已知 ()() *111 1021n r x a δ-≤ + 那么， ()()***1112110.10 211 102 r m n n m n x x x x a a a a δ----≤?=+≤ 即：*x 至少有n 位有效数字。 例： 3.141592653589793...π= 取其近似值如下： x*=3.14 x * =3.14159 x*=3.1415 x*=3.141 **213 100.314 110.0016...0.005101022 x x π--=?-=<=?=? **516 100.314159 110.0000026...0.00000510102 2 x x π--=?-=<=?=? **314 100.31415 110.000092...0.0001101022 x x π--=?-=<
目录 一.课程设计目的 (1) 二.设计环境 (1) 三.设计内容 (1) 四.设计说明 (2) 五.设计程序流程图 (2) 六.调试 (4) (1)错误原因分析一 (4) (2)语法错误 (5) (3)逻辑错误 (5) 七. 调试结果图 (6) 八. 结论与心得体会 (7) 九.附录 (8) 具体代码实现 (8) 十．参考文献 (18)
一.课程设计目的 1.通过一个学期的学习，我认为要学号C语言程序这门课程，不仅要认真阅读课本知识，更重要的是要通过上机实践来巩固我 们的知识，特别是学计算机专业的，我们更应该注重这一环节， 只有这样我们才能成为一个合格的计算机人才。通过这一个课程 设计，进一步来巩固所学的语句，如：循环，和分支结构的运用。还要熟悉四则运算和函数的算法。 2.通过这次课程设计扩展自己的知识面，课本上的东西是远 远不够的，可以通过上网或去图书馆查资料等方式得到一些新的 知识， 3.通过课程设计，加深对课程化设计思想的理解，能进行一 个系统功能分析，并设计一个合理的模块化结构，提高程序开发 能力。 二.设计环境 1.硬件：一台完整的电脑，包括键盘、鼠标，最小硬盘空间1GHz 2.软件：安装有Microsoft visual c++6.0 三.设计内容 以简易计算器为例，通过对简单应用软件计算器的设计，编制、调试，实现
简单的加，减，乘，除等运算，以学习应用MFC库类编写对话框的原理，加深对C++类的学习及应用。 （1）定义一个结构体类型数组，输入0~9及+、--、*等符号的信息，将其信息存入文件中； （2）输入简单的加减乘除算术计算式，并在屏幕上显示计算结果； （3）画出部分模块的流程图； （4）编写代码； （5）程序分析与调试。 四.设计说明 1)包含的功能有：加、减、乘、除运算，开方、平方等功能。 (2)计算器上数字0—9为一个控件数组，加、减、乘、除为一个控件数组，其余为单一的控件。 (3)输入的原始数据、运算中间数据和结果都显示在窗口顶部的同一个标签中。 (4)计算功能基本上是用系统内部函数。 (5)程序可以能自动判断输入数据的正确性，保证不出现多于一个小数点、以0开头等不正常现象。 (6)“CE”按钮可以清除所有已输入的数据从头计算 五.设计程序流程图
大地主题结算实验报告 姓名：高成 学号：20104214 测绘2班 一、实验目的： 1.提高运用计算机语言编程开发的能力； 2.加深对大地主题解算计算公式及辅助参数的理解并掌握计算步骤； 3.通过编程语言实现大地主题解算。 二、工具： Windows XP Mode 环境下的Microsoft Visual C++ 6.0 三、注意事项： 1.计算所需变量多，容易混淆； 2.正反算函数的编写； 3.函数调用； 4.弧度与角度之间的转化。 四、实验要求： 1.提交报告，实验总结，编写代码； 2.独立编程，调试运行； 3.上交成果：编写思想，编写过程，问题分析，源代码，计算结果；
实验源代码： #include
毕业论文文献综述 信息与计算科学 导数的数值计算方法 一、 前言部分 导数概念的产生有着直觉的起源，与曲线的切线和运动质点的速度有密切的关系．导数用于描述函数变化率，刻画函数的因变量随自变量变化的快慢程度．比如说，物理上考虑功随时间的变化率（称为功率），化学上考虑反应物的量对时间的变化率（称为反应速度），经济学上考虑生产某种产品的成本随产量的变化率（称为边际成本）等等，这些变化率在数学上都可用导数表示． 导数由于其应用的广泛性，为我们解决所学过的有关函数问题提供了一般性的方法，导数是研究函数的切线、单调性、极值与最值等问题的有力工具；运用它可以简捷地解决一些实际问题，导数的概念是用来研究函数在一点及其附近的局部性质的精确工具，而对于函数在某点附近的性质还可以应用另一种方法来研究，就是通过最为简单的线性函数来逼近，这就是微分的方法．微分学是数学分析的重要组成部分，微分中值定理作为微分学的核心，是沟通导数和函数值之间的桥梁， Rolle 中值定理， Lagrange 中值定理， Cauchy 中值定理， Taylor 公式是微分学的基本定理， 统称为微分学的中值定理，这四个定理作为微分学的基本定理，是研究函数形态的有力工具 ] 1[．在微分学中，函数的导数是通过极限定义的，但 当函数用表格给出时，就不可用定义来求其导数，只能用近似方法求数值导数] 2[．最简单 的数值微分公式是用差商近似地代替微商，常见的有 [3] ． ()()() 'f x h f x f x h +-≈ ， ()()() 'f x f x h f x h --≈， ()()() '2f x h f x h f x h +--≈ ． 需要注意的是微分是非常敏感的问题，数据的微小扰动会使结果产生很大的变化] 4[．
#include"stdio.h" #include"math.h" #include"stdlib.h" struct optrstyle { int top; char stack[20]; } struct opndstyle { int top; int stack[20]; } void main() { char ch; while(1) { printf("-----------------------------------------\n"); printf("|-----------欢迎来到速算24--------------|\n"); printf("|----------是否开始游戏?????------------|\n"); printf("|--------- 开始游戏请按'Y'-------------|\n"); printf("|--------- 退出游戏请按其他键-------------|\n"); printf("-----------------------------------------\n"); scanf("%c",&ch); if(ch=='Y' or ch=='y') { } else break; } } int EvaluateExpression() { // 算术表达式求值的算符优先算法。设OPTR 和OPND 分别为运算符栈和运算数栈，OP 为运算符集合。 InitStack(OPTR); Push(OPTR,'#'); InitStack(OPND); c = getchar(); while(c!='#' || GetTop(OPTR)!='#') { if(!In(c,OP))
地球科学与环境工程学院实验报告书 课程名： 学号： 姓名： 指导老师： 日期：
目录 一、目的与要求 (1) 二、实验内容 (1) 三、计算公式整理 (1) 四、程序代码 (4) 五、计算结果 (15) 六、实验体会 (16)
一、目的与要求 参考椭球面是大地测量计算的基准面。大地坐标是椭球面上的基本坐标系，根据大地测量的观测成果（如距离与方向），从大地原点出发，逐点计算在椭球面上的大地坐标；或根据两点的大地坐标，计算它们之间的大地线长度和大地方位角，这类计算称为大地问题解算（或称为大地主题解算）。大地问题解算的用途是多方面的，随着现代空间技术和航空航天、航海等领域的发展，大地问题解算（尤其是大地反算）有着更为重要的作用，因此需要熟练掌握其计算。 二、实验内容 在《大地测量学基础》教材中，介绍了高斯平均引数法与白塞尔方法的计算过程、步骤。鉴于此，需要熟练掌握高斯平均引数法与白塞尔方法解大地主题问题的基本方法与原理。采用所熟悉的计算机语言编程计算。计算时采用克拉索夫椭球参数，至少完成其中一种方法正反算，按照数据序号选取不同的已知数据，在计算结果中注明所选取的数据序号，选取其它数据作为无效数据处理。 三、计算公式整理 、高斯平均引数正算计算公式（S< 200 km） 、高斯平均引数正算计算公式（S< 200 km）
四、程序代码 、角度转换类的头文件： #pragma once const double Pi=; class AngleTrans { public: AngleTrans(void); ~AngleTrans(void); double D,F,M,DFM,Rad,Ten; double trans1(double DFM), f"),B2); (_T("%.7f"),L2), (_T("%.7f"),A21); UpdateData(false); } 、反算的计算按钮代码： void C大地主题高斯引数正反算Dlg::OnBnClickedButton1() { UpdateData(true); ZhengFanSuan _ZhengFanSuan; double B1=_wtof(fB1), L1=_wtof(fL1), B2=_wtof(fB2), L2=_wtof(fL2), S=(B1,L1,B2,L2), A12=(B1,L1,B2,L2), A21=(B1,L1,B2,L2); (_T("%.7f"),S); (_T("%.7f"),A12); (_T("%.7f"),A21); UpdateData(false); } 、清零按钮代码： void C大地主题高斯引数正反算Dlg::OnBnClickedButton3() { UpdateData(true);
课程设计（论文） 题目: 三次样条插值问题 学院: ___ 理学院 _ 专业: __ _ 数学与应用数学 班级：数学08-2班 学生姓名: 魏建波 学生学号: 080524010219 指导教师：李文宇 2010年12月17日
课程设计任务书
目录 摘要……………………………………………………………………… 一、前言………………………………………………………………… （一）Lagrange插值的起源和发展过程……………………………………… （二）本文所要达到的目的……………………………………………………… 二、插值函数…………………………………………………………… （一）函数插值的基本思想…………………………………………………… （二）Lagrange插值的构造方法……………………………………………… 三、MATLAB程序………………………………………………………… （一）Lagrange程序…………………………………………………………… （二）龙格程序………………………………………………………………… 四、理论证明…………………………………………………………… 五、综述……………………………………………………………………谢辞………………………………………………………………………参考文献…………………………………………………………………
摘要
前言 要求：500字以上，宋体小四，行距20磅，主要内容写该算法的产生及发展、应用领域等。 题目 整体要求：报告页数，正文在8页以上 字体：宋体小四（行距20磅） 内容：1、理论依据 2、问题描述 3、问题分析 4、求解计算（程序） 5、结论 注：（1）页码编号从正文页开始 （2）标题可根据情况自己适当改动 示例见下： 2判别…………………… 2.1 判……………… 2.1.1 判别……………… 所谓的判别分析，………………………………………………方法[3]。 2.1.2 判………………………… 常用的有四种判别方法：…………………………………………………步判别法[6]。 1. 马氏………………
超详细 一、因为计算器设计的控件太多，不便使用控制台应用程序完成，所以这里 使用Windows窗体应用程序，并命名为Calc,如下图所示： 二、向窗体中拖入需要的控件，如下图所示： （完成效果图） 结果显示区（作者博客左边的文本框）是TextBox控件，并修改其name为txtShow，按键0~9为Button控件，并将其name分别修改为btn_0、btn_1、btn_2、btn_3、btn_4、btn_5、btn_6、btn_7、btn_8、btn_9；按键【负数】的name值修改为btn_sign，按键【.】的name修改为btn_dot，按键【+-*/】的name值分别修改为btn_add、btn_sub、btn_mul、btn_div，按键【=】的name值修改为btn_equ，按键【倒数】的name值修改为btn_rev，按键【平方】的name值修改为btn_sqr，按键【开方】的name值修改为btn_sqrt。右边的计算器图片空间是PictureBox，作者博客控件是LinkLabel，可以不添加，以上所有控件均可按照需求添加，只保留自己需要的按钮控件和textbox控件即可。 三、代码部分（含解释），采用switch多分支语句编写 using System; using System.Drawing; using System.Collections; using https://www.doczj.com/doc/4816095986.html,ponentModel; using using System.Data; namespace Calc { ///
白塞尔大地主题解算。 方向： 学号： 姓名： \
一．基本思路： ； 基本思想：将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后在将球面上的计算结果换算到椭球面上。其关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式，同时解决在球面上进行大地主题解算的方法。 正算流程： 1.计算起点的归化纬度 2.计算辅助函数值，解球面三角形 3.按公式计算相关系数A,B,C 以及α,β 4.计算球面长度 # 5.计算纬度差改正数 6.计算终点大地坐标及大地方位角 、 011122S B C A {sin (cos )} σσσ=-+101010 22222sin ()sin sin cos cos σσσσσσ+=+101010 22222cos ()cos cos sin sin σσσσσσ+=-001101 522B C A [cos ()]sin () σσσσσσ=++++010122L A sin [(sin ()sin )] λδασβσσσ-==++-2111u u u A sin sin cos cos cos sin σσ = +222 222 2222 2222221111u B u B W u B u B B arctan e u sin cos cos tan sin tan cos ?==????=????==--1 111 A arctan u u A sin sin [] cos cos sin sin cos σλσσ=-21L L λδ =+ -
反算流程： 1.辅助计算 2.用逐次趋近法同时计算起点大地方位角、球面长度及经差，第一次趋近时，取δ＝０。 ! 计算下式,重复上述计算过程2. 3.计算大地线长度S ; 4. 计算反方位角 二．已知数据 L λδ=+211212u p A u u u u q sin cos tan cos sin sin cos cos λλ= = -2121p p u q b b A arctan q sin cos cos λλ==-=11 p A q A sin cos tan cos σσ += 11 p A q A sin sin cos σ=+12a a cos cos σλ =+arctan sin cos σσσ?? = ??? 011 A u A sin cos sin =111u A tan tan sec σ=21+σσσ =02122L A sin [(sin sin )] λδασβσσ-==+-L λδ =+11222222S A B C B C sin (cos )sin (cos ) σσσσσ=++-+
第四章 问题一 一、问题综述 在离地球表面高度为y处的重力加速度如下： 计算高度y=55000m处的重力加速度值。 二、问题分析 以高度y作为自变量，重力加速度的值为因变量。得到以下信息： f(0)=9.8100; f(30000)=9.7487; f(60000)=9.6879; f(90000)=9.6278; f(120000)=9.5682; 本题要求的就是f(55000)的值。 以下将采用课堂中学到的Lagrange插值多项式法、Newton插值多项式法、分段低次插值法和样条插值法求解该问题。 三、问题解决 1. lagrange插值多项式法 对某个多项式函数，已知有给定的k+ 1个取值点： 其中对应着自变量的位置，而对应着函数在这个位置的取值。 假设任意两个不同的x j都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：
其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为： 拉格朗日基本多项式的特点是在上取值为1，在其它的点上取值为0。 源程序lagrange.m function [c,f]=lagrange(x,y,a) % 输入：x是自变量的矩阵；y是因变量的矩阵；a是要计算的值的自变量； % 输出：c是插值多项式系数矩阵；f是所求自变量对应的因变量； m=length(x); l=zeros(m,m); % l是权矩阵 f=0; for i=1:m v=1; for j=1:m if i~=j v=conv(v,poly(x(j)))/(x(i)-x(j)); % v是l_i(x)的系数矩阵 end end l(i,:)=v; % l矩阵的每一行都是x从高次到低次的系数矩阵 end c=vpa(y*l,10); % 对应阶次的系数相加，乘以y，显示10位有效数字 for k=1:m f=f+c(k)*a^(m-k); end 输入矩阵 x=[0 30000 60000 90000 120000] y=[9.81 9.7487 9.6879 9.6278 9.5682] a=55000 再运行源函数，可得： c = [ -2.057613169e-23, 4.938271605e-18, -3.703703702e-14, -0.000002046111111, 9.81] f = 9.6979851723251649906109417384537
C语言简单计算器的编 写课程设计 IMB standardization office【IMB 5AB- IMBK 08- IMB 2C】
2013-2014学年第二学期 《高级语言程序设计》 课程设计报告 题目：简单计算器 专业：数字媒体技术 班级：13级1班 姓名：邢洪波 指导教 师：王辉 成绩： 计算机与信息工程系 二0一四年二月二十日 1设计内容及要求 利用本学期所学的《C语言程序设计》课程，运用相关知识，查阅相关资料，编写C语言程序，设计一个简单计算器，要求编写的简单计算器能够模拟windows系统的计算器，用户能够用键盘输入
相关数据，能够进行简单的加、减、乘、除运算，并且在程序运行过程中能够正常的退出程序。 在计算器程序运行中，输入数据时如果遇到输入错误的情况，能够能过键盘上的退格键进行删除，并且重新输入正确的数据。在数据输入完成后，如果需要放弃本次计算操作，可以利用程序中设置好的按键进行清零，并为下一次运算作准备。 运用相关知识，查阅相关资料，编写一个简单的计算器，能够实现简单的基本的加减乘除运算，在计算器工作时用户能利用键盘或鼠标进行相应操作。程序基本功能要求实现完整，并有简单的验证。程序还 须有符合标准的程序设计报告。
2概要设计 系统流程图 系统流程图主要功能表
3设计过程或程序代码 #include<> #include<> #include<>/*malloc的头文件代表在内存中开辟存储空间*/ doublejisuan(chara[]) { inti=1,j,k,m,cnt=0,t1=0,t2=0,t3=0;/*定义整型变量 i,j,k,m,n,cnt,t1,t2,t3*/ charnibo[50],zhan2[50];/*定义字符型数组*/ doublex,n,l,z=0,zhan3[50];/*定义双精度变量x,n,l,z,zhan3[50]*/ typedefstruct/*表示可以用con定义结构体变量*/ { doubled1; intd2; }dd; typedefstruct { dddata[50]; inttop; }zhan1; zhan1*shu;/*定义指针变量*/ shu=(zhan1*)malloc(sizeof(zhan1)); shu->top=0; while(a[i]!='\0')/*外循环是用while语句*/ {
简易四则运算计算器 计算机工具地历史而言，中国古代最早采用地一种计算工具叫筹策，又叫做算筹，这种算筹多用竹子制成，也有用木头，兽骨充当材料的，大约270枚一束，放在布袋里可随身携带。直到今天仍在使用的住算盘，是中国古代计算工具领域中的另一项发明明代时的住算盘已经与现代的住算盘几乎相同。 17世纪初，西方国家的计算工具有了较大的发展。英国数学家纳尔发现了“纳皮尔算筹”；英国牧师奥却德发明了圆柱形对数计算尺，这种计算尺不仅能做加减乘涂，乘方和开方运算，甚至可以计算三角函数，指数函数和对数函数。这些计算工具不仅带动了计算的发展，也为现代计算器发展尊定了良好的基础，成为现代社会应用广泛的计算工具，如手机，操作系统上附带的多功能计算器。 项目设计目的 掌握模块的划分； 掌握结构体定义； 掌握栈的操作； 掌握函数的定义； 掌握栈的基本操作； 掌握VC的基本绘图库； 掌握鼠标的响应操作； 目设计内容 设计一个含有界面的计算器，输入一个包含括号的表达式，使用栈数据类型实现整数的四则运算操作，开方运算。 项目设计要求 根据分析，系统要求实现以下几点基本功能： 可通过按钮输入数字、运算符； 能通过按钮实现退格、清除功能； 实现加、减、乘、除、开方运算功能； 实现括号运算； 显示运算结果； 良好的交互界面。
项目详细设计 功能模块图 根据系统功能要求，主要含有四大模块，分别是计算器界面设计模块、计算器按键控制模块、四则混合运算模块、计算器记忆处理模块。 计算器界面设计模块：调用VC图形系统函数和字符函数画出计算器的界面，包括24个按钮和一个文本输入框。 计算机按键控制模块：计算器通过键盘按键值的判断，执行相应的操作，如接收数字输入等。 计算器计算处理模块。计算处理模块主要完成可以包括括号的表达式运算，运算包括加，减，乘，除、开方。计算处理模块在按键控制模块中被调用执行。处理模块是通过栈来实现的，分别定义了两个顺序栈，一个用来存放操作符栈，一个用来存放操作数栈。 计算器记忆处理模块。记忆处理模块使计算其具有记忆功能。记忆当前输入的表达式，给出一个最后的运算结果。可以清楚部分或全部记忆的值。 图简易计算器模块图 数据结构 该系统主要的数据有操作符集，用一维数组存放；存放操作符的优先级关系集，操作符的优先级关系用一个二维数组来体现；另外定义两个栈，一个用来存放操作符，一个用来存放操作数。 char ops[7]={'+','-','*','/','(',')','#'}; int cmp[7][7]={ {2,2,1,1,1,2,2}, {2,2,1,1,1,2,2}, {2,2,2,2,1,2,2}, {2,2,2,2,1,2,2}, {1,1,1,1,1,3,0}, {2,2,2,2,0,2,2}, {1,1,1,1,1,0,3}}; typedef struct { char elem[Stack_Size]; int top; }SeqStack; typedef struct { int elem[Stack_Size]; int top; }SeqStackOperand;
大地测量编程实习报告 班号：XXXX 学号：XXXXXXXXXX 姓名：XXX 大地主题解算（分段累加法正算） 结果截图： 主要代码： privatevoid button1_Click(object sender, EventArgs e) { double bx, by, bz, B1, lx, ly, lz, L1, ax, ay, az, A1, S; double dB, dL; double e2 = 0.006693421622966, a = 6378245.0000000000; double M, N, W, C; double B2 = 0, L2 = 0, A2 = 0; int Bx, By, Lx, Ly, Ax, Ay; double Bz, Lz, Az; bx = Convert.ToDouble(du1.Text); by = Convert.ToDouble(fen1.Text); bz = Convert.ToDouble(miao1.Text); lx = Convert.ToDouble(du2.Text); ly = Convert.ToDouble(fen2.Text); lz = Convert.ToDouble(miao2.Text); ax = Convert.ToDouble(du3.Text); ay = Convert.ToDouble(fen3.Text); az = Convert.ToDouble(miao3.Text); S = Convert.ToDouble(m.Text); double PI = Math.PI; B1 = (bx + by / 60 + bz / 3600) * PI / 180; L1 = (lx + ly / 60 + lz / 3600) * PI / 180; A1 = (ax + ay / 60 + az / 3600) * PI / 180; W = Math.Sqrt(1 - e2 * (Math.Sin(B1) * Math.Sin(B1)));