matlab上机作业绘制椭球面

第七章 椭球面上的基本计算 2006版 控制测量学讲稿

第七章 椭球面上的基本计算

§1 地球椭球的基本知识

离心力 地球的自然表面——不规则；不能在上面进行计算； 大地水准面——平均海水面延伸得到的封闭曲面，最接近大地自然 地心引力 重力G 表面； ∵大地水准面具有性质：大地水准面上任一点处的垂线（重力方向） 与该点处切面正交；

又：重力是离心力与地心引力的合力（离心力与地心引力之比约1： 300），而大地水准面上各点处引力不等，造成各点处垂线方向各异。

∴各点处切面组成的曲面——大地水准面亦不规则，有微小起伏，是一个具有物理性质的曲面。 实践和理论均可证明：1）在各水准面（与大地水准面的不平行性不很明显）上测得的水平角，因归化到大地水准面上改正极微小，完全可以看成大地水准面上的角值；2）各高程面上测得之边长也可化算到大地水准面上；3）地面点的高程亦从大地水准面起算。

结论：大地水准面是测量外业的基准面；但它是物理曲面而非数学曲面，所以不能作为测量计算的基准面。 大地体—

第七章 椭球面上的测量计算

1．试写出椭球的基本元素及其基本关系式。

2．在控制测量的椭球解算中，常引用下列符号：t、?、W、V，试问它们之间函数关系的一个基本共同特点是什么？

3．我国解放后主要采用哪两种参考椭球?其主要参数是什么? 4．绘图并说明表示椭球面上点位的三种常用坐标系统。

5．在报纸上经常看到X X号轮船在东经XXX度，北纬X X度遇险一类的报导，试问这是指的什么坐标系，为什么?

6．写出参考椭球体的五个基本元素及相互间的关系。 7．什么叫子午圈？什么叫平行圈？

8．参考椭球体扁率的变化，椭球体的形状发生怎样的变形？

9．简要说明并图示地面某一点的大地高、正常高以及大地水准面差距的几何意义。

10．什么是大地测量的基本坐标系？有何优点？

11．画图表示地心纬度坐标系和归化纬度坐标系，这两种坐标系在大地测量中有何意义？ 12．用公式表示空间直角坐标系和大地坐标系之间的关系。 13．何为大地纬度、归化纬度、地心纬度？三者间有何关系？ 14．简要叙述M、N、R三种曲率半径之间的关系。

15．大地坐标系和天文坐标系各以什么作基准面和基准线？ 16．试推证卯酉圈、子午圈曲率

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后

数学建模上机作业（3月30日与4月2日1,2班用）

提交作业邮箱

gdshumo2015@163.com;

作业将两题M文件压缩在一个自己中文名+学院+年级+学号命名的压缩文件里面提交

1．课堂练习题（课堂完成，和课后作业一起在周日交）

编写一个函数，实现以下功能： （1）提示“请输入矩阵X=“（使用input）

（2）任意输入一个矩阵X，求矩阵中各行与各列元素的最大值Mrow与Mlist，并求整个矩阵的最大值与最小值M与m；

（3）求该矩阵的每行元素的乘积和全部元素的乘积S1，S2；

（4）对矩阵元素按照行和按照列进行排序，输出结果XR与XL和排序前元素位置矩阵R和L；

（5）将R和L合并成一个新矩阵U1，U1元素中为R和L相同位置对应元素较大者;同理将XR与XL合并成U2，U2为两个矩阵对应位置元素较小者；

（6）计算U1每列的算术平均和U2每行的中位数，将结果分别乘以M,m，再除以S2，得到U11和U22；

（7）将向量Mrow，Mlist，S1，U11，U22全部作列向量（提示：回顾一下上面输出的哪些是行向量）合并成一个新矩阵XX，计算XX的协方差阵CXX； （8）计算CXX所有元素的和SS2和每行元素的和SS1； （9）输出列向量

Matlab上机作业4

说明：1、应用Matlab软件，完成下列各题

2、交实验报告纸文档:实验报告格式 （1）标题：经济模型上机实验4 （2）年级专业学号姓名

（3）实验内容：问题（序号）、源程序（或函数命令）、运行结果

一、概率计算

1、设X~N(350,3502),求：概率P(X?250)；

2、设X~P(4)柏松分布,求：X0为何值时，P(X?X0)达到0.5；

3、某人定点投篮投中率为0.3，求： （1）投篮10次，命中5次的概率；

（2）投篮几次，命中达到或超过5次的概率达到0.5；

二、描述统计

4、在北京奥运场馆某次比赛中收集到观众的调查数据为（文件夹s08s4：data6.m）：

矩阵A，各列数据意义如下：

性别（男1、女2）、年龄（20岁以下1、20—30岁2、30—50岁3、50岁以上4）、坐公交出行（南北方向）、坐公交出行（东西方向）、坐出租出行、开私车出行、坐地铁出行（东向）、坐地铁出行（西向）、中餐馆午餐、西餐馆午餐、商场内餐饮午餐、非餐饮消费额

求：（1）男女各为多少人

（2）非餐饮消费额：最高、最低、平均、标准差 （1）：

（2）

三、回归分析

5、某地区12个气象观测站近8年来各观测站测得的周降水量为data6

《MATLAB语言与应用》上机实验作业

第一部分

1、

在matlab窗口下点击’help’弹出如下窗口，窗口左侧的列表框可以选择各种不同组合的演示内容。

2、

>> A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1] A =

1 2 3 4 4 3 2 1 2 3 4 1 3 2 4 1 >>

B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j] B =

1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i

4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i

2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.000