(数学建模)小行星的轨迹问题
12.建模作业_行星运行
《数学建模》课程作业题—12第四章 微分方程模型—行星运行1水星到太阳的最远距离为0.68982×1011m ，此时水星绕太阳运行的线速度为3.886×104m/s ，试求: (1）水星到太阳的最近的距离； （2)水星绕太阳运行的周期;（3)画出水星绕太阳运行的轨道曲线；（4）求从远日点开始的第50天(地球天)结束时水星的位置.(1)建立的模型及求解：水星运行的数学模型2212321200000t t t C d r GMdt r r C d dt r r r drr dtθθ===⎧-=-⎪⎪⎪=⎪⎪=⎨⎪⎪=⎪⎪=⎪⎩行星运动轨迹方程 1cos pr e θ=-其中，2110001,,C pe p C r v r MG=-==利用MATLAB 进行计算求解（附录一），得到水星到太阳的最近的距离为104.456510r m =⨯（2）建立的模型及求解：水星绕太阳运行的周期模型为:()2322121p T C eπ=-其中，2110001,,C pe p C r v r MG=-==。
利用MATLAB 进行计算求解(附录二），得到水星绕太阳运行的周期为T = 1.58×107s ，约为183.4天.(3)利用matlab （附录三)画水星绕太阳运行的轨道曲线：（4）建立的模型及求解:水星的位置模型为()2211p d t C ecoc θθθ=-⎰如果要求出1t T =是水星的位置,即求出相应的θ和r ，则意味着先要解方程()211221C T p d p ecoc θθθ=-⎰ 求出θ后,在求出1cos p r e θ=-，进而得到线速度12C d dt rθ=。
最后利用MATLAB 计算求解（附录四），解得从远日点开始的第50天(地球天）结束时水星的位置如图所示：2 (地中海鲨鱼问题) 20世纪20年代中期，意大利生物学家D'Ancona 偶然注意到第一次世界大战期间在原南斯拉夫的里耶卡港，人们捕获的鱼类中，鲨鱼等软骨鱼的百分比大量增加(见下表），而供其捕食的食用鱼的百分比却明显下降. 显然，战争使捕鱼量下降，食用鱼应该增加，鲨鱼等软骨鱼也随之增加，但为何其比例大幅度增加呢？请对软骨鱼及食用鱼的增长情况建立一个数学模型来解释这一现象。
线性代数在数学建模中的应用举例
线性代数在数学建模中的应用举例1 基因间“距离”的表示在ABO 血型的人们中，对各种群体的基因的频率进行了研究。
如果我们把四种等位基因A 1，A 2，B ，O 区别开，有人报道了如下的相对频率，见表1.1。
表1.1基因的相对频率问题 一个群体与另一群体的接近程度如何？换句话说，就是要一个表示基因的“距离”的合宜的量度。
解 有人提出一种利用向量代数的方法。
首先，我们用单位向量来表示每一个群体。
为此目的，我们取每一种频率的平方根，记ki ki f x =.由于对这四种群体的每一种有141=∑=i ki f ，所以我们得到∑==4121i kix .这意味着下列四个向量的每个都是单位向量.记.44434241,34333231,24232221,141312114321⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=x x x x a x x x x a x x x x a x x x x a在四维空间中，这些向量的顶端都位于一个半径为1的球面上. 现在用两个向量间的夹角来表示两个对应的群体间的“距离”似乎是合理的.如果我们把a 1和a 2之间的夹角记为θ，则由于| a 1|=| a 2|=1，再由内只公式，得21cos a a ⋅=θ而.8307.03464.02943.03216.0,8228.01778.00000.05398.021⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=a a 故 9187.0cos 21=⋅=a a θ 得 2.23=θ°. 按同样的方式，我们可以得到表1.2.表1.2基因间的“距离”爱斯基摩人班图人 英国人 朝鲜人 爱斯基摩人 0° 23.2° 16.4° 16.8° 班图人 23.2° 0° 9.8° 20.4° 英国人 16.4° 9.8° 0° 19.6° 朝鲜人16.8°20.4°19.6°0°由表1.2可见，最小的基因“距离”是班图人和英国人之间的“距离”，而爱斯基摩人和班图人之间的基因“距离”最大.2 Euler 的四面体问题问题 如何用四面体的六条棱长去表示它的体积？这个问题是由Euler （欧拉）提出的.解 建立如图2.1所示坐标系，设A ，B ，C 三点的坐标分别为（a 1,b 1,c 1）,( a 2,b 2,c 2)和（a 3,b 3,c 3），并设四面体O-ABC 的六条棱长分别为.,,,,,r q p n m l 由立体几何知道，该四面体的体积V 等于以向量→→→OC OB OA ,,组成右手系时，以它们为棱的平行六面体的体积V 6的16.而)(.3332221116c b a c b a c b a OC OB OA V =⋅⨯= 于是得 .6333222111c b a c b a c b a V = 将上式平方，得.362323233232323231313232322222221212131313121212121212133322211133322211122c b a c c b b a a c c b b a a c c b b a a c b a c c b b a a c c b b a a c c b b a a cb ac b a c b a c b a c b a c b a c b a V ++++++++++++++++++=⋅=根据向量的数量积的坐标表示，有.,,,,232323323232222222313131212121212121c b a OC OC c c b b a a OC OB c b a OB OB c c b b a a OC OA c c b b a a OB OA c b a OA OA ++=⋅++=⋅++=⋅++=⋅++=⋅++=⋅ 于是362OC OC OB OC OB OBOB OBOA OB OA OAV ⋅⋅⋅= （2.1）由余弦定理，可行.2cos 222n q p q p OB OA -+=⋅⋅=⋅θ同理.2,2222222l r q OC OB m r p OC OA -+=⋅-+=⋅将以上各式代入（2.1）式，得.222222362222222222222222222222r l r p m r p l r p p n q p m r p n q p pV -+-+-+-+-+-+=（2.2）这就是Euler 的四面体体积公式.例 一块形状为四面体的花岗岩巨石，量得六条棱长分别为l =10m, m =15m, n =12m, p =14m, q =13m, r =11m.则.952222,462222,5.1102222=-+=-+=-+l r p m r p n q p代入（2.1）式，得.75.13698291219546951695.110465.110196236==V 于是.)195(82639.38050223m V ≈≈即花岗岩巨石的体积约为195m 3.古埃及的金字塔形状为四面体，因而可通过测量其六条棱长去计算金字塔的体积.3 动物数量的按年龄段预测问题问题 某农场饲养的某种动物所能达到的最大年龄为15岁，将其分成三个年龄组：第一组，0～5岁；第二组，6～10岁；第三组，11～15岁.动物从第二年龄组起开始繁殖后代，经过长期统计，第二组和第三组的繁殖率分别为4和3.第一年龄和第二年龄组的动物能顺利进入下一个年龄组的存活率分别为12 和14 .假设农场现有三个年龄段的动物各100头，问15年后农场三个年龄段的动物各有多少头？问题分析与建模 因年龄分组为5岁一段，故将时间周期也取为5年.15年后就经过了3个时间周期.设)(k i x 表示第k 个时间周期的第i 组年龄阶段动物的数量（k =1，2，3；i =1，2，3）.因为某一时间周期第二年龄组和第三年龄组动物的数量是由上一时间周期上一年龄组存活下来动物的数量，所以有).3,2,1(41,21)1(2)(3)1(1)(2===--k x x x x k k k k又因为某一时间周期，第一年龄组动物的数量是由于一时间周期各年龄组出生的动物的数量，所以有).3,2,1(34)1(3)1(2)(1=+=--k x x x k k k于是我们得到递推关系式：⎪⎪⎪⎩⎪⎪⎪⎨⎧==+=----.41,21,34)1(2)(3)1(1213)1(2)(1k k k k k k k x x x x x x x 用矩阵表示).3,2,1(0410021340)1(3)1(2)1(1)(3)(2)(1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---k x x x x x x k k k k k k则).3,2,1()1()(==-k Lx x k k其中.100010001000,04100021340)0(⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=x L 则有),3,2,1()(3)(2)(1)(=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=k x x x x k k k k,250500700010001000100004100021340)0()1(⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡==Lx x,12535002750250500700004100021340)1()2(⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡==Lx x .8751375143751253500275004100021340)2()3(⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡==Lx x 结果分析 15年后，农场饲养的动物总数将达到16625头，其中0～5岁的有14375头，占86.47%，6～10岁的有1375头，占8.27%，11～15岁的有875头，占 5.226%.15年间，动物总增长16625-3000=13625头，总增长率为13625/3000=454.16%.注 要知道很多年以后的情况，可通过研究式)0()1()(x L Lx x k k k ==-中当趋于无穷大时的极限状况得到.关于年龄分布的人口预测模型 我们将人口按相同的年限（比如5年）分成若干年龄组，同时假设各年龄段的田、女人口分布相同，这样就可以通过只考虑女性人口来简化模型.人口发展随时间变化，一个时间周期的幅度使之对应于基本年龄组间距（如先例的5年），令)(k i x 是在时间周期k 时第i 个年龄组的（女性）人口，i =1,2,…,n .用1表示最低年龄组，用n 表示最高年龄组，这意味着不考虑更大年龄组人口的变化.假如排除死亡的情形，则在一个周期内第i 个年龄组的成员将全部转移到i +1个年龄组.但是，实际上必须考虑到死亡率，因此这一转移过程可由一存活系数所衰减. 于是，这一转移过程可由下述议程简单地描述：),1,,2,1()1()(1-==-+n i x b x k ii k i其中i b 是在第i 个年龄组在一个周期的存活率，因子i b 可由统计资料确定.惟一不能由上述议程确定的年龄组是,)(1k x 其中的成员是在后面的周期内出生的，他们是后面的周期内成员的后代，因此这个年龄组的成员取决于后面的周期内各组的出生率及其人数.于是有方程,)1(122)1(11)(1---+++=k n n k k k x a x a x a x （3.1）这里),,2,1(n i a i =是第i 个年龄组的出生率，它是由每时间周期内，第i 个年龄组的每一个成员的女性后代的人数来表示的，通常可由统计资料来确定.于是我们得到了单性别分组的人口模型，用矩阵表示便是,00000000000)1()1(3)1(2)1(11211321)()(3)(2)(1⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡------k n k k k n n n k n k k k x x x x b b b a a a a a x x x x 或者简写成.)1()(-=k k Lx x （3.2）矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--000000000001211321n n n b b b a a a a a L称为Leslie 矩阵.由（3.2）式递推可得)0()1()(x L Lx x k k k ==-这就是Leslie 模型.4 企业投入产生分析模型问题 某地区有三个重要产业，一个煤矿、一个发电厂和一条地方铁路.开采一元钱的煤，煤矿要支付0.25元的电费及0.25元的运输费.生产一元钱的电力，发电厂要支付0.65元的煤费，0.05元的电费及0.05元的运输费.创收一元钱的运输费,铁路要支付0.55元的煤费及0.10元的电费.在某一周内，煤矿接到外地金额为50000元的定货，发电厂接到外地金额为25000元的定货，外界对地方铁路没有需求.问三个企业在这一周内总产值多少才能满足自身及外界的需求？数学模型 设x 1为煤矿本周内的总产值，x 2为电厂本周的总产值，x 3为铁路本周内的总产值，则⎪⎩⎪⎨⎧=⨯++-=++-=++⨯-,0)005.025.0(,25000)10.005.025.0(,50000)55.065.00(321332123211x x x x x x x x x x x x （4.1） 即.02500050000005.025.010.005.025.055.065.00321321⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡x x x x x x 即.025********,005.025.010.005.025.055.065.00,321⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=Y A x x x X 矩阵A 称为直接消耗矩阵，X 称为产出向量，Y 称为需求向量，则方程组（4.1）为,Y AX X =-即Y X A E =-)(， （4.2）其中矩阵E 为单位矩阵，（E-A ）称为列昂杰夫矩阵，列昂杰夫矩阵为非奇异矩阵.投入产出分析表 设,00000,)(3211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=--=-x x x A C E A E B D=（1，1，1）C.矩阵B 称为完全消耗矩阵，它与矩阵A 一起在各个部门之间的投入产生中起平衡作用.矩阵C 可以称为投入产出矩阵，它的元素表示煤矿、电厂、铁路之间的投入产出关系.向量D 称为总投入向量，它的元素是矩阵C 的对应列元素之和，分别表示煤矿、电厂、铁路得到的总投入.由矩阵C ，向量Y ，X 和D ，可得投入产出分析表4.1.表4.1 投入产出分析表 单位：元 煤矿电厂铁路外界需求总产出煤矿 11c 12c 13c 1y 1x电厂 21c 22c 23c 2y 2x 铁路 31c32c33c 3y3x总投入1d 2d 3d计算求解 按（4.2）式解方程组可得产出向量X ，于是可计算矩阵C 和向量D ，计算结果如表4.2.表4.2 投入产出计算结果 单位：元 煤矿 电厂 铁路 外界需求 总产出 煤矿 0 36505.96 15581.51 50000 102087.48 电厂 25521.87 2808.15 2833.00 25000 56163.02 铁路 25521.87 2808.15 0 0 28330.02总投入51043.7442122.2718414.525 交通流量的计算模型问题 图5.1给出了某城市部分单行街道的交通流量（每小时过车数）.假设：（1）全部流入网络的流量等于全部流出网络的流量；（2）全部流入一个节点的流量等于全部流出此节点的流量.试建立数学模型确定该交通网络未知部分的具体流量.建模与计算 由网络流量假设，所给问题满足如下线方程组：234457612157891091083630050020080080010004002006001000x x x x x x x x x x x x x x x x x x x x -+=⎧⎪+=⎪⎪-=⎪+=⎪⎪+=⎪⎨+=⎪⎪=⎪-=⎪⎪=⎪++=⎪⎩ 系数矩阵为11100000000011000000000011000110000000010001000000000001100000000001000000000110000000001010010100A -⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 增广矩阵阶梯形最简形式为1000100000800010010000000010000000200000110000050000000101008000000001100100000000000104000000000001600000000000000000000000B ⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦其对应的齐次方程组为1525345687891000000000x x x x x x x x x x x x x +=⎧⎪-=⎪⎪=⎪+=⎪⎨+=⎪⎪+=⎪=⎪⎪=⎩取（x 5,x 8）为自由取值未知量，分别赋两组值为（1,0），（0,1），得齐次方程组基础解系中两个解向量()11,1,0,1,1,0,0,0,0,0,'η=--()20,0,0,0,0,1,1,1,0,0'η=--其对应的非齐次方程组为1525345687891080002005008001000400600x x x x x x x x x x x x x +=⎧⎪-=⎪⎪=⎪+=⎪⎨+=⎪⎪+=⎪=⎪⎪=⎩赋值给自由未知量（x 5,x 8）为（0,0）得非齐次方程组的特解()800,0,200,500,0,800,1000,0,400,600'x *=于是方程组的通解,*2211x k k x ++=ηη其中k 1，k 2为任意常数，x 的每一个分量即为交通网络未知部分的具体流量，它有无穷多解.6 小行星的轨道模型问题 一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：1.4959787×1011m ）.在5个不同的时间对小行星作了5次观察，测得轨道上5个点的坐标数据如表6.1.表6.1 坐标数据由Kepler （开普勒）第一定律知，小行星轨道为一椭圆.现需要建立椭圆的方程以供研究（注：椭圆的一般方程可表示为012225423221=+++++y a x a y a xy a x a .问题分析与建立模型 天文学家确定小行星运动的轨道时，他的依据是轨道上五个点的坐标数据：（x 1, y 1）, （x 2, y 2）, （x 3, y 3）, （x 4, y 4）, （x 5, y 5）.由Kepler 第一定律知，小行星轨道为一椭圆.而椭圆属于二次曲线，二次曲线的一般方程为012225423221=+++++y a x a y a xy a x a .为了确定方程中的五个待定系数，将五个点的坐标分别代入上面的方程，得2211211314151221222232425222132333343532214244344454221525535455522212221222122212221a x a x y a y a x a y a x a x y a y a x a y a x a x y a y a x a y a x a x y a y a x a y a x a x y a y a x a y ⎧++++=-⎪++++=-⎪⎪++++=-⎨⎪++++=-⎪⎪++++=-⎩这是一个包含五个未知数的线性方程组，写成矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡11111222222222222222543215525552544244424332333232222222211211121a a a a a y x y y x x y x y y x x y x y y x x y x y y x x y x y y x x 求解这一线性方程组，所得的是一个二次曲线方程.为了知道小行星轨道的一些参数,还必须将二次曲线方程化为椭圆的标准方程形式:12222=+bY a X 由于太阳的位置是小行星轨道的一个焦点,这时可以根据椭圆的长半轴a 和短半轴b 计算出小行星的近日点和远日点距离,以及椭圆周长L .根据二次曲线理论,可得椭圆经过旋转和平移两种变换后的方程如下:[]22120D X Y C λλ++=所以,椭圆长半轴:C D a 1λ=;椭圆短半轴: CDb 2λ=;椭圆半焦矩:22b ac -=.计算求解 首先由五个点的坐标数据形成线性方程组的系数矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=7200.69600.142896.112656.509504.550520.53360.143807.62127.363802.516460.35180.133233.36433.246841.454040.25720.124448.11115.155138.39292.1528.114199.04701.72237.33A使用计算机可求得12345(,,,,)(0.6143,0.3440,0.6942, 1.6351,0.2165)a a a a a =---从而⎪⎪⎭⎫⎝⎛--=⎥⎦⎤⎢⎣⎡=6942.03440.03440.06143.03221a a a a C C C ,3081.0=的特征值120.3080, 1.0005λλ==123235450.61430.3440 1.63510.34400.69420.21651 1.63510.21651a a a D a a a a a ---⎡⎤⎡⎤⎢⎥⎢⎥==--⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦.8203.1-=D于是,椭圆长半轴a=19.1834,短半轴b=5.9045,半焦距c=18.2521.小行星近日点距和远日点距为039313,37.4355h a c H a c =-==+=最后,椭圆的周长的准确计算要用到椭圆积分,可以考虑用数值积分解决问题,其近似值为84.7887.7 人口迁移的动态分析问题 对城乡人口流动作年度调查,发现有一个稳定的朝向城镇流动的趋势:每年农村居民的2.5%移居城镇,而城镇居民的1%迁出.现在总人口的60%位于城镇.假如城乡总人口保持不变,并且人口流动的这种趋势继续下去,则一年以后住在城镇人口所占比例是多少两年以后呢十年以后呢最终呢解 设开始时,令乡村人口为,0y 城镇人口为,0z 一年以后有乡村人口,10011000975100y z y =+ 城镇人口 ,10099100025100z z y =+或写成矩阵形式⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡00111009910002510011000975z y z y . 两年以后,有.100991000251001100097510099100025100110009750021122⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y z y . 十年以后,有.100991000251001100097500101010⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y 事实上,它给出了一个差分方程:k k Au u =+1.我们现在来解这个差分方程.首先,1009910002510011000975⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=Ak 年之后的分布(将A 对角化):.75757275100200193115210000⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎪⎭⎫ ⎝⎛⎥⎥⎦⎤⎢⎢⎣⎡-=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y A z y k k k k 这就是我们所要的解,而且容易看出经过很长一个时期以后这个解会达到一个极限状态.7572)(00⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=⎥⎦⎤⎢⎣⎡∞∞z y z y 总人口仍是00z y +,与开始时一样,但在此极限中人口的75在城镇,而72在乡村.无论初始分布是什么样,这总是成立的.值得注意这个稳定状态正是A 的属于特征值1的特征向量.上述例子有一些很好的性质:人口总数保持不变,而且乡村和城镇的人口数决不能为负.前一性质反映在下面事实中:矩阵每一列加起来为1;每个人都被计算在内,而没有人被重复或丢失.后一性质则反映在下面事实中:矩阵没有负元素;同样地0y 和0z 也是非负的,从而1y 和21,y z 和2z 等等也是这样.8 常染色体遗传模型为了揭示生命的奥秘,遗传学的研究已引起了人们的广泛兴趣.动植物在产生下一代的过程中,总是将自己的特征遗传给下一代,从而完成一种“生命的延续”.在常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对.人类眼睛颜色即是通过常染色体控制的,其特征遗传由两个基因A 和a 控制.基因对是AA 和Aa 的人,眼睛是棕色,基因对是aa 的人,眼睛为蓝色.由于AA 和Aa 都表示了同一外部特征,或认为基因A 支配a ,也可认为基因a 对于基因A 来说是隐性的(或称A 为显性基因,a 为隐性基因).下面我们选取一个常染色体遗传——植物后代问题进行讨论.某植物园中植物的基因型为AA ,Aa ,aa .人们计划用AA 型植物与每种基因型植物相结合的方案培育植物后代.经过若干年后,这种植物后代的三种基因型分布将出现什么情形我们假设),2,2,0(,, =n c b a n n n 分别代表第n 代植物中,基因型为AA ,Aa 和aa 的植物占植物总数的百分率,令),,()('=n n n n c b a x为第n 代植物的基因分布, ),,(000)0('=c b a x 表示植物基因型的初始分布,显然,我们有.1000=++c b a (8.1)先考虑第n 代中的AA 型，第1-n 代AA 型与AA 型相结合，后代全部是AA 型；第1-n 代的Aa 型与和与AA 相结合，后代是AA 型的可能性为21；1-n 代的aa 型与AA 型相结合，后代不可能是AA 型。
小行星运动轨迹的Runge-Kutta法模拟
小行星运动的Runge-Kutta 法模拟一、背景介绍由于两个恒星作用下行星运动问题没有解析解，只能用数值方法求解微分方程。
但是在用一阶近似求解微分方程的时候存在严重的误差累积。
当只考虑一个恒星引力影响时的模型如下： (1)当初始值是00001,0,'0,'1x y x y ====时，行星做圆周运动。
此时，微分方程的解是cos()sin()x t y t =⎧⎨=⎩。
在后面的讨论中，用这个初始条件的方程作为测试方程。
如果采用一阶近似，(1)()'(),'(1)'()''()(1)()'(),'(1)'()''()x n x n hx n x n x n hx n y n y n hy n y n y n hy n +=++=+⎧⎨+=++=+⎩，就会有严重的误差累积。
如下图所示当行星偏离理想轨道很小的量以后，之后的偏差就会越来越大，直至脱离恒星的束缚。
在离散化以后，原来临界稳定的系统变得发散了。
二、用高阶系统去求解单恒星问题当用高于一阶的方法近似求解以上方程时，会取得较好一些的近似。
把二阶常微分方程组(1)转化为一阶常微分方程组：3222322200000000''()''(),'',''x x x y y y x y x x x x y y y y -⎧=⎪⎪+⎪-⎪=⎨⎪+⎪==⎪⎪==⎩32223222''()''()x x y y x v x v x y y v y v x y =⎧⎪-⎪=⎪+⎪⎨=⎪⎪-=⎪⎪+⎩,初始条件是00001001x y x y v v =⎧⎪=⎪⎨=⎪⎪=⎩ 一阶常微分方程组00'(,)()x x =⎧⎨=⎩Y F Y Y Y 的经典4阶RK 法的公式是112341213243()6(,)(,)22(,)22(,)n n n n n n n n n n h x h h x h h x x h h +⎧=++++⎪⎪=⎪⎪⎪=++⎨⎪⎪=++⎪⎪=++⎪⎩Y Y K K K K K F Y K F Y K K F Y K K F Y K 当0.01h =时，迭代100000次，模拟行星绕行星100000*0.011592π≈圈的轨迹图如下：从上图中可以看出，当模拟绕中心159圈后，轨道的偏移依然很小。
（数学建模）小行星的轨迹问题
（数学建模）⼩⾏星的轨迹问题问题15 ⼩⾏星的轨迹问题⼀、问题⼀天⽂学家要确定⼀颗⼩⾏星绕太阳运⾏的轨道，他在轨道平⾯内建⽴以太阳为远点的直⾓坐标系，在两坐标轴上取天⽂测量单位（⼀天⽂单位为地球到太阳的平均距离：1.4959787*10^11m ），在5个不同的时间对⼩⾏星作了5次观察，测得轨道上5个点的坐标数据如表2.15.1. 表2.15.1由开普勒第⼀定律知，⼩⾏星轨道为⼀椭圆，现需要建⽴椭圆的⽅程以供研究。
（注：椭圆的⼀般⽅程可表⽰为012225423221=+++++y a x a y a xy a x a ）。
⼆、实验⽬的利⽤5个点确定⼆次曲线的⼀般⽅程，并求出椭圆的重要参数。
三、预备知识线性代数⽅程组理论，椭圆的有关概念及性质。
四、实验内容与要求1.⽤表中5个点的坐标数据分别代⼊椭圆的⼀般⽅程可建⽴5个⽅程的线性代数⽅程组，该⽅程组的系数矩阵为A ，右端项为b ，这⾥，21x 112y x 21y 12x 12y -122x222y x 22y 22x22y -1A= 23x 332y x23y32x 32y b= -124x 442y x24y 42x 42y -1 25x 552y x 25y 52x52y -1试依据题⽬所给的5个点的坐标，⽤计算机计算出矩阵的A 的5*5个数据。
2.利⽤Matalb 指令A\b 求解5元线性代数⽅程组，写出椭圆⽅程012225423221=+++++y a x a y a xy a x a 中的5个待定系数54321,,,,a a a a a 及⼩⾏星多所对应的曲线⽅程。
3.写出曲线表达式中系数所对应的⼆阶矩阵和三阶矩阵：1a2a 1a 2a 3aC= D= 2a 3a 5a2a 3a 4a 5a 1并利⽤Matlab 指令eig （C ）求出矩阵C 的特征值，记录数据=1λ（），=2λ（）利⽤Matlab 指令det （D ）计算⼆阶导数和三阶⾏列式的值；=C （）， =D （）4.利⽤公式计算椭圆的下列参数：长半轴：a=CD1λ=（），短半轴：b=CD2λ=（），办焦距：c=22b a -=（），写出椭圆标准⽅程（），5：利⽤上⾯的椭圆有关数据求出⼩⾏星轨道的参数：⼩⾏星的近⽇点距离：h=a-c=（），⼩⾏星的远⽇点距离：H=2c+h=（），椭圆轨道周长近似值：L ≈π-+ab b a )(23=（）6．*试在Matlab 环境下利⽤参数⽅程： x=acos(t) (t ∈[0, 2π]) y=bsin(t)绘制出以椭圆中⼼为原点的椭圆图形（图2.15.1）图2.15.1五．思考问题你能否利⽤定积分求弧长公式推导出椭圆的周长公式？如果你所得到的是⼀个定积分表达式，利⽤这⼀表达式计算⼩⾏星轨道的椭圆周长。
数学建模-小行星的轨道模型
小行星的轨道模型问题:天文学家要确定一颗小行星绕太阳运行的轨道。
他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（1.4959787e11m ）。
在5个不同的时间对小行星作了5次观察，测的轨道上5个点的坐标数据如表：一、问题分析与建立模型（1） 求解椭圆一般方程的系数由开普勒第一定律知，小行星的轨道为椭圆，现需要建立椭圆方程(221234522210a x a xy a y a x a y +++++=)以供研究。
天文学家确定小行星运动的轨迹时，他的依据是轨道上5个点的坐标数据。
由椭圆轨道属于二次曲线，是一般方程。
为了确定方程中的5个待定系数，需要将上述5个点的坐标代入上面的方程22123452221a x a xy a y a x a y ++++=-,得：22112113141512212222324252221323333435322142443444542215255354555a x 2a x y a y 2a x 2a y 1a x 2a x y a y 2a x 2a y 1a x 2a x y a y 2a x 2a y 1a x 2a x y a y 2a x 2a y 1a x 2a x y a y 2a x 2a y 1⎛++++=- ++++=- ++++=-++++=- ++++=-⎝， 将这一包含5个未知数的线性方程组，写成矩阵的形式221111112222222222333333?2244444422555555x 2x y y 2x 2y x 2x y y 2x 2y x 2x y y 2x 2y x 2x y y 2x 2y x 2x y y 2x 2y ⎛⎫⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭12345a a a a a ⎛⎫⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭= 11111-⎛⎫ ⎪- ⎪⎪- ⎪- ⎪ ⎪-⎝⎭。
求解这一线性方程组，即可得到曲线方程的系数。
数学建模期末大作业-2013年
数学建模期末大作业-2013年期末大作业题目一、小行星的轨道问题一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立了以太阳为原点的直角坐标系，在两坐标轴上取天文观测单位。
在5个不同的时间对（1）建立小行星运行的轨道方程并画出其图形；（2）求出近日点和远日点及轨道的中心（是太阳吗？）；（3）计算轨道的周长。
二、发电机使用计划为了满足每日电力需求（单位：兆瓦），可以选用四种不同类型的发电机。
每日电力需求如下所示：一最小输出功率。
所有发电机都存在一个启动成本，以及工作于最小功率状态时的固定的每小时成本，并且如果功率高于最小功率，则超出部分的功率每兆瓦每小时还存在一个成本，即边际成本。
这些数据均列于下表中。
电机不需要付出任何代价。
我们的问题是：（1）在每个时段应分别使用哪些发电机才能够使每天的总成本最小？（2）如果增加表3中的关闭成本，那么在每个时段应分别使用哪些发电机才能够使每天的总成本最小？（3）如果增加表4中的关闭成本，那么在每个时段应分别使用哪些发电机才能够使每天的总成本最小？三、合理计税问题所以此人一年上税为：245×12+__=__元在实际的执行过程中，每月的岗位津贴和年末一次性奖金实际上是放在一起结算给个人的，而具体每月发放多少岗位津贴和年末一次性发放多少奖金可以由职工本人在年初根据自己的需要进行选择。
显然，不同的选择发放方式所缴纳的税是不同的，这就产生一个合理计税的问题。
假定该事业单位一年中的津贴与奖金之和的上限是__元，试解决下面这个问题：四、光伏电池的选购问题早在1839年，法国科学家贝克雷尔（Becqurel）就发现，光照能使半导体材料的不同部位之间产生电位差。
这种现象后来被称为“光生伏特效应”，简称“光伏效应”。
1954年，美国科学家恰宾和皮尔松在美国贝尔实验室首次制成了实用的单晶硅太阳电池，诞生了将太阳光能转换为电能的实用光伏发电技术。
据预测，太阳能光伏发电在未来会占据世界能源消费的重要席位，不但要替代部分常规能源，而且将成为世界能源供应的主体。
小行星轨道模型(论文)
出椭圆的重要参数。因此将资料锁定在线性代数方程组理论和椭圆的有关概念上。即 用 5 个点的坐标数据分别带入椭圆的一般方程的线性代数方程组������1������2 + 2������2������������ + ������3������2 + 2������4������ + 2������5������ + 1 = 0，该方程组的系数为 A，右端顶为 b，则：
a 椭圆周长的计算：
将托圆周长计算转化为数学模型，已知椭圆长半轴 a ,偏心率 e ，求其周长。
椭圆的方程为：
x asin，y bcos
其微分为
dx acosd，y bsind
弧元为
ds dx2 dy2 a2 cos2 b2 sin2 d a2 (1 sin2 ) b2 sin2 d
a3 a4
1 1
。
x22 2x5 y5 y52
2x5
2y5 a5 1
求解这一线性方程组，即可得到曲线方程的系数。
三、模型的建立与求解
根据上述建立的模型，对数据进行处理，利用 MATLAB 对数据进行处理编程得：
求解结果：
������������ = 4������������ − ������2
������2 = 2(������������������2 + ������������������2 + ������������������������������ − 1) ������ + ������ − √(������ − ������)2 + ������2
小行星轨迹问题
所以代入所观测的5个数据，就会得到5个方程，组成一个方程组，并观察，发现此线性方程组可以用matlab来求解，并且该方程组的矩阵是非奇异的矩阵。下面的是未代入实数的矩阵式子：
= (式子2)
将上列矩阵经过变换得到下面的式子
-1
-1
A= b= -1
-1
-1
(
对于矩阵运算，我们代入测得的5个数据后得到椭圆的方程式，并且可以得到此时的椭圆的在matlab中的图形（图形2）（程序见附录2）
椭圆方程：
0.050 -0.0702xy+0.0381 -0.4530x+0.2642y+1=0;
（模型1）
椭圆图形：
图
为了知道小行星轨道的一些参数比如焦点、长半轴、短半轴，还必须将二次曲线方程化为椭圆的标准方程形式： + =1(式子4);
由于太阳的位置是小行星的一个焦点，这时可以根据椭圆的长半轴a和短半轴b计算出小行星的近日点和远日点距离，以及椭圆周长L.根据二次曲线理论，可得椭圆经过旋转和平移两种变换后的方程如下：
(式子5)
， 、 为C的特征值。
利用matlab计算（见程序3）
C的行列式值为m=7.0347e-004
C的特征值分别为 =0.0088 =0.0801
D的行列式为n=0.0010
又 ， ，所以由matlab（程序4）算得：
长半轴a =12.8320
短半轴b =4.2532
椭圆半焦距： =12.1066
A=[x.^2,2*x.*y,y.^2,2*x,2*y];
>> b=[-1;-1;-1;-1;-1];
>> Y=inv(A)*b%Y是系数
Y =
小行星轨道问题
题目：小行星轨道问题姓名：刘天华班级：车辆工程1106班学号：0121102910819任课老师：陈建业题目：小行星轨道问题摘要本文针对小行星轨道问题提出合理假设，利用开普勒定律和二次曲线理论对模型进行优化，将复杂的方程式用矩阵表示并使用Matlab软件进行求解的出小行星轨道的椭圆标准方程，有利于进一步分析小行星的轨道特征。
关键词：开普勒定律二次曲线理论矩阵 Matlab一、问题重述2013年2月16日，一颗直径大约50米的小行星与地球擦肩而过，小行星撞击地球危险可能再度引起公众的关注。
已知：要确定一颗小行星绕太阳运行的轨道，需要在轨道平面内建立以太阳为原点的空间直角坐标系，然后在不同时刻对小行星进行观测，以确定其轨道。
现在已经在5个不同时刻对某颗小行星进行了5次观测，表1给出了相应的观测数据。
表1：某小行星的5次观测数据（单位：天文单位）其中一个天文单位等于地球到太阳的平均距离，即11101.4959787 米。
要求确定这颗小行星的轨道，如椭圆的半长轴、半短轴、半焦距、近日点、远日点，以及椭圆轨道的周长等。
二、模型的基本假设1：小行星稳定绕太阳运行，不会因为撞击改变轨道。
2：小行星运行符合开普勒第一定律，即为一椭圆。
三、问题的分析及模型的建立由开普勒第一定律知，小行星轨道为一椭圆.现需要建立椭圆的方程，椭圆的一般方程为：012225423221=+++++y a x a y a xy a x a现在已经由上述表格知道轨道上五个点的坐标数据：（x 1, y 1）, （x 2, y 2）, （x 3, y 3）, （x 4, y 4）, （x 5, y 5） 分别对应坐标数据：（5.764，0.648），（6.286 ，1.202），（6.759，1.832），（7.168，2.526），（7.480，3.360）问题就变成了球方程的五个待定系数a1,a2,a3,a4,a5。
为了确定方程中的五个待定系数，将五个点的坐标分别代入上面的方程，得⎪⎪⎪⎩⎪⎪⎪⎨⎧-=++++-=++++-=++++-=++++-=++++.1222122212221222122255542535522514544243442241353423333223125242232222211514213112211y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a由于直接求解需要大量的计算工作，我们可以利用矩阵这一数学工具来优化模型：⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡11111222222222222222543215525552544244424332333232222222211211121a a a a a y x y y x x y x y y x x y x y y x x y x y y x x y x y y x x求解这一线性方程组，所得的是一个二次曲线方程.为了知道小行星轨道的一些参数,还必须将二次曲线方程化为椭圆的标准方程形式:12222=+bY a X 由于太阳的位置是小行星轨道的一个焦点,这时可以根据椭圆的长半轴a 和短半轴b 计算出小行星的近日点和远日点距离,以及椭圆周长L .根据二次曲线理论,可得椭圆经过旋转和平移两种变换后的方程如下:[].02221=++C DY X λλ 所以,椭圆的长半轴:C D a 1λ=;椭圆的短半轴: CDb 2λ=;椭圆的半焦矩:22b ac -=.所以只要求出参数a1,a2,a3,a4,a5，并应用二次曲线理论，即可求出小行星轨道椭圆的半长轴、半短轴、半焦距、近日点、远日点，以及椭圆轨道的周长等数据。
数学实验____方程模型及其求解算法_参考答案
实验2 方程模型及其求解算法一、实验目的及意义[1] 复习求解方程及方程组的基本原理和方法；[2] 掌握迭代算法；[3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)；[4] 通过范例展现求解实际问题的初步建模过程；通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。
这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。
二、实验内容1．方程求解和方程组的各种数值解法练习2．直接使用MATLAB命令对方程和方程组进行求解练习3．针对实际问题，试建立数学模型，并求解。
三、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口；2．根据各种数值解法步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．根据观察到的结果写出实验报告，并浅谈学习心得体会。
四、实验要求与任务基础实验1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。
画出图形程序：x=-10::10;y=x.*sin(x)-1;y1=zeros(size(x));plot(x,y,x,y1)MATLAB运行结果：-10-8-6-4-20246810-8-6-4-22468扩大区间画图程序：x=-50::50;y=x.*sin(x)-1;y1=zeros(size(x));plot(x,y,x,y1)MATLAB 运行结果：由上图可知，该方程有偶数个无数的根。
2．将方程x 5+5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。
（1）画图：x1=-6::6;x2=-3::3;x3=-1::1;x4=::;y1=x1.^5 +5*x1.^3-2*x1+1;y2=x2.^5 +5*x2.^3-2*x2+1;y3=x3.^5 +5*x3.^3-2*x3+1;y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1),title('子图 (1)') ,grid on, subplot(2,2,2),plot(x2,y2),title('子图 (2)'),grid on, subplot(2,2,3),plot(x3,y3),title('子图 (3)'),grid on, subplot(2,2,4),plot(x4,y4),title('子图 (4)') ,grid on,由图可知 x 的初值应在（，）之间。
太阳系内小行星运动规律数学建模与实证
太阳系内小行星运动规律数学建模与实证太阳系是我们所在的宇宙家园，由太阳和围绕其运动的行星、卫星、彗星、小行星等组成。
小行星作为太阳系中的重要组成部分，对于研究太阳系的形成和演化过程以及行星间的相互作用具有重要意义。
研究小行星的运动规律对于探索宇宙的奥秘和保护地球的安全具有重要的科学意义。
小行星运动规律的研究需要建立数学模型，以描述和预测小行星在太阳系中的运动。
数学建模是科学研究中的一种重要方法，通过建立数学模型来描述现实世界的物理过程和规律，从而推导出有关性质和行为的定量预测。
在小行星运动规律的数学建模中，我们需要考虑以下几个关键因素：首先，我们需要考虑太阳的引力对小行星的作用。
根据牛顿的万有引力定律，太阳对于小行星的引力与两者之间的距离和质量有关。
因此，在数学建模中，我们可以利用行星动力学方程来描述小行星的运动轨迹。
该方程即为万有引力定律和牛顿第二定律的结合，能够描述小行星受到的合力与加速度之间的关系。
其次，我们还需考虑其他行星对于小行星的引力影响。
除了太阳的引力外，行星之间也会相互作用，这会对小行星的运动轨迹产生微小的扰动。
这些扰动将导致小行星的轨道发生变化，进而影响其运动规律。
因此，在数学建模中，我们需要考虑多体动力学问题，并将其纳入到行星动力学方程中。
此外，在数学建模中，我们还需考虑相对论效应的影响。
当小行星的速度接近光速时，相对论效应将显著影响其运动规律。
相对论效应主要包括时空弯曲和时间膨胀两个方面。
时空弯曲将导致小行星在太阳系中的运动轨迹发生偏离，而时间膨胀则会影响小行星的时钟速率。
因此，在数学建模中，我们需要考虑相对论效应的修正，并将其纳入到行星动力学方程中。
除了数学建模，我们还需要进行实证研究，以验证数学模型的有效性和准确性。
实证研究主要基于观测数据，通过观测小行星的位置和速度等参数，来评估数学模型的拟合程度和预测精度。
通过与实际观测结果的对比，我们可以检验和修正数学模型，并进一步提高其预测能力。
小行星轨道问题MATLAB数学实验(2)
11.4599 11.9771 12.5177 13.0827 13.6731 14.2902
18/18
思考题与练习题
1. 用线性函数与指数函数两种模型对中国人口数据 进行拟合的结果差异是否很大？哪一种模型的残差 平方和更小？
2.用马尔萨斯模型预测中国人口达到20亿，将会在多 少年以后发生？以中国960万平方公里国土面积计算， 人均占多少平方米？
特征向量 1 [0.12,0.1]T
X=[120;150]; A=[0.9,0.12;0.1,0.88];
Cars=X; for k=1:6
X=A*X; Cars=[Cars,X]; end
figure(1),bar(Cars(1,:)) figure(2),bar(Cars(2,:))
150 100 50
alpha=[0.12;0.1];R=alpha/sum(alpha)
X1=R*270;
taxs=X1;
for k=1:6
X1=A*X1;
taxs=[taxs,X1];
end
figure(3),bar(taxs(1,:))
11/18
离散数据的多项式拟合方法
x
x1 x2 …… xm
f(x)
y1 y2 …… ym
求 n 次多项式 ( n < m )
P(x) = a1xn + a2 xn-1 + …… + an x + an+1
m
使得 S(a1,, an1 ) [ y j P( x j )]2 min j 1
MATLAB求解多项式拟合方法如下:
P =polyfit(x，y，n) Y=polyval(P,x)
分析：根据马尔萨斯人口理论 ,人口增长率与人口 数量N成正比，用微分方程描述为
行星运动的数学模型及轨道分析
行星运动的数学模型及轨道分析引言：行星运动一直以来都是人们关注和研究的话题之一。
通过数学建模和轨道分析，我们可以更深入地了解行星运动的规律和特性。
本文将探讨行星运动的数学模型以及轨道分析的方法。
一、行星运动的数学模型行星运动的数学模型是基于牛顿万有引力定律的。
该定律表明，两个物体之间的引力与它们的质量成正比，与它们之间的距离的平方成反比。
根据这一定律，我们可以用以下公式表示行星的运动：F =G * (m1 * m2) / r^2其中，F表示两个物体之间的引力，G是万有引力常数，m1和m2分别是行星和太阳的质量，r是行星和太阳之间的距离。
数学模型的关键是求解行星的运动轨迹，通常使用数值模拟的方法。
将行星的质量、初始位置和速度输入计算机程序，利用欧拉法或其他数值计算方法，我们可以得到行星在各个时间点上的位置和速度，并通过连续的计算得到整个行星运动轨迹。
二、行星轨道的分析行星轨道分析的目的是了解行星运动的规律和特性。
通过分析行星的轨道，我们可以研究行星的运动周期、轨道形状以及行星之间的相互影响。
1. 运动周期行星的运动周期是指行星绕太阳完成一次运动所需的时间。
根据开普勒第三定律，行星的运动周期与它的轨道半长轴的立方成正比。
因此，通过观测行星的运动轨迹，我们可以计算出它的运动周期，并进一步了解行星运动的规律。
2. 轨道形状行星的轨道形状通常是椭圆。
根据开普勒第一定律，行星绕太阳运动的轨道是一个椭圆，太阳位于椭圆的一个焦点上。
通过测量行星在不同时刻的位置和速度，利用椭圆的方程，我们可以确定行星的轨道参数，如椭圆的长轴、短轴以及离心率等，从而推测行星的轨道形状。
3. 相互影响行星之间的相互影响是行星运动的另一个重要方面。
根据牛顿引力定律，行星之间的引力会影响它们的运动轨道。
特别是在行星相互靠近的时候，它们的轨道可能会发生变化。
通过数值模拟和轨道计算，我们可以研究行星之间的相互作用，进一步了解行星运动的规律。
数学建模典型例题(二)
6 小行星的轨道模型问题 一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：1.4959787×1011m ）.在5个不同的时间对小行星作了5次观察，测得轨道上5个点的坐标数据如表6.1.表6.1 坐标数据由Kepler （开普勒）第一定律知，小行星轨道为一椭圆.现需要建立椭圆的方程以供研究（注：椭圆的一般方程可表示为012225423221=+++++y a x a y a xy a x a .问题分析与建立模型 天文学家确定小行星运动的轨道时，他的依据是轨道上五个点的坐标数据：（x 1, y 1）, （x 2, y 2）, （x 3, y 3）, （x 4, y 4）, （x 5, y 5）.由Kepler 第一定律知，小行星轨道为一椭圆.而椭圆属于二次曲线，二次曲线的一般方程为012225423221=+++++y a x a y a xy a x a .为了确定方程中的五个待定系数，将五个点的坐标分别代入上面的方程，得⎪⎪⎪⎩⎪⎪⎪⎨⎧-=++++-=++++-=++++-=++++-=++++.1222122212221222122255542535522514544243442241353423333223125242232222211514213112211y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a ，y a x a y a y x a x a 这是一个包含五个未知数的线性方程组，写成矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡11111222222222222222543215525552544244424332333232222222211211121a a a a a y x y y x x y x y y x x y x y y x x y x y y x x y x y y x x 求解这一线性方程组，所得的是一个二次曲线方程.为了知道小行星轨道的一些参数,还必须将二次曲线方程化为椭圆的标准方程形式:12222=+b Y a X 由于太阳的位置是小行星轨道的一个焦点,这时可以根据椭圆的长半轴a 和短半轴b 计算出小行星的近日点和远日点距离,以及椭圆周长L .根据二次曲线理论,可得椭圆经过旋转和平移两种变换后的方程如下:[].02221=++C DY X λλ 所以,椭圆长半轴:C D a 1λ=;椭圆短半轴: CDb 2λ=;椭圆半焦矩:22b ac -=.计算求解 首先由五个点的坐标数据形成线性方程组的系数矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=7200.69600.142896.112656.509504.550520.53360.143807.62127.363802.516460.35180.133233.36433.246841.454040.25720.124448.11115.155138.39292.1528.114199.04701.72237.33A使用计算机可求得).2165.0,6351.1,6942.0,3440.0,6143.0(),,,,(54321---=a a a a a从而⎪⎪⎭⎫⎝⎛--=⎥⎦⎤⎢⎣⎡=6942.03440.03440.06143.03221a a a a C C C ,3081.0=的特征值.0005.1,3080.021==λλ.12165.06351.12165.06942.03440.06351.13440.06143.0154532321⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-------=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=a a a a a a a a D .8203.1-=D于是,椭圆长半轴1834.19=a ,短半轴9045.5=b ,半焦距2521.18=c .小行星近日点距和远日点距为.4355.37,039313=+==-=c a H c a h最后,椭圆的周长的准确计算要用到椭圆积分,可以考虑用数值积分解决问题,其近似 值为84.7887.7 人口迁移的动态分析问题 对城乡人口流动作年度调查,发现有一个稳定的朝向城镇流动的趋势:每年农村居民的2.5%移居城镇,而城镇居民的1%迁出.现在总人口的60%位于城镇.假如城乡总人口保持不变,并且人口流动的这种趋势继续下去,那么一年以后住在城镇人口所占比例是多少?两年以后呢?十年以后呢?最终呢?解 设开始时,令乡村人口为,0y 城镇人口为,0z 一年以后有乡村人口,10011000975100y z y =+ 城镇人口,10099100025100z z y =+ 或写成矩阵形式⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡00111009910002510011000975z y z y . 两年以后,有.100991000251001100097510099100025100110009750021122⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y z y . 十年以后,有.100991000251001100097500101010⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y事实上,它给出了一个差分方程:k k Au u =+1.我们现在来解这个差分方程.首先,1009910002510011000975⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=Ak 年之后的分布(将A 对角化):.75757275100200193115210000⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎪⎭⎫ ⎝⎛⎥⎥⎦⎤⎢⎢⎣⎡-=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡z y z y A z y k k k k 这就是我们所要的解,而且容易看出经过很长一个时期以后这个解会达到一个极限状态.7572)(00⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=⎥⎦⎤⎢⎣⎡∞∞z y z y 总人口仍是00z y +,与开始时一样,但在此极限中人口的75在城镇,而72在乡村.无论初始分布是什么样,这总是成立的.值得注意这个稳定状态正是A 的属于特征值1的特征向量.上述例子有一些很好的性质:人口总数保持不变,而且乡村和城镇的人口数决不能为负.前一性质反映在下面事实中:矩阵每一列加起来为1;每个人都被计算在内,而没有人被重复或丢失.后一性质则反映在下面事实中:矩阵没有负元素;同样地0y 和0z 也是非负的,从而1y 和21,y z 和2z 等等也是这样.8 常染色体遗传模型为了揭示生命的奥秘,遗传学的研究已引起了人们的广泛兴趣.动植物在产生下一代的过程中,总是将自己的特征遗传给下一代,从而完成一种“生命的延续”.在常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对.人类眼睛颜色即是通过常染色体控制的,其特征遗传由两个基因A 和a 控制.基因对是AA 和Aa 的人,眼睛是棕色,基因对是aa 的人,眼睛为蓝色.由于AA 和Aa 都表示了同一外部特征,或认为基因A 支配a ,也可认为基因a 对于基因A 来说是隐性的(或称A 为显性基因,a 为隐性基因).下面我们选取一个常染色体遗传——植物后代问题进行讨论.某植物园中植物的基因型为AA ,Aa ,aa .人们计划用AA 型植物与每种基因型植物相结合的方案培育植物后代.经过若干年后,这种植物后代的三种基因型分布将出现什么情形?我们假设),2,2,0(,, =n c b a n n n 分别代表第n 代植物中,基因型为AA ,Aa 和aa 的植物占植物总数的百分率,令),,()('=n n n n c b a x 为第n 代植物的基因分布, ),,(000)0('=c b a x 表示植物基因型的初始分布,显然,我们有.1000=++c b a (8.1)先考虑第n 代中的AA 型，第1-n 代AA 型与AA 型相结合，后代全部是AA 型；第1-n 代的Aa 型与和与AA 相结合，后代是AA 型的可能性为21；1-n 代的aa 型与AA 型相结合，后代不可能是AA 型。
小行星探测任务中的轨迹规划优化算法研究
小行星探测任务中的轨迹规划优化算法研究随着科技的不断进步，人们对于宇宙的探索也变得越来越深入。
在这个过程中，小行星成为了一个备受关注的领域，因为它们不仅可以为我们认识宇宙的大爆炸时期提供线索，还可能带有地外生命的存在。
在小行星勘探任务中，轨迹规划是关键的一步。
因为小行星的运行是不稳定的，其轨迹的复杂性也很高，所以为了使小行星探测任务更加高效和准确，需要对轨迹进行规划和优化。
下面就让我们来探讨小行星探测任务中的轨迹规划优化算法研究。
一、小行星探测任务中的轨迹规划在小行星勘探任务中，轨迹规划是探测任务过程中的重要部分。
根据轨迹规划的特点和小行星运动的规律性，轨迹规划通常可以分为单目标轨迹规划和多目标轨迹规划两种方式。
单目标轨迹规划是为了达到探测任务的一个主要目标而设计的，例如获取小行星表面的样本，或是监测小行星的磁场等。
而多目标轨迹规划则需要考虑多个目标之间的平衡关系，例如同时监测小行星的多个特征或以不同等级搜寻小行星等。
二、轨迹规划优化算法的研究轨迹规划问题本身是一个复杂的优化问题。
因此，众多的优化算法也开始应用到小行星探测任务的轨迹规划中。
以下是一些常用的轨迹规划优化算法：1. 遗传算法（GA）遗传算法起源于进化论的概念，其对生物进化过程的模拟思想被应用到了轨迹规划优化问题中。
遗传算法是通过模拟种群的进化过程，不断从中筛选出最优的轨迹规划方案。
2. 人工神经网络（ANN）人工神经网络是一种计算模型，它能够通过一系列的节点和连接线模拟人类神经元的结构。
在小行星探测任务中，人工神经网络可以用于预测小行星运动模型，进而实现更加准确的轨迹规划。
3. 粒子群优化算法（PSO）粒子群优化算法则是通过随机初始化一群粒子，引导其不断地寻找最优轨迹规划方案。
其核心思想是通过让多个粒子互相作用，达到最终最优目标的效果。
以上三种优化算法仅是小行星探测任务中的部分轨迹规划优化算法。
在实际应用中，我们可以根据任务需求和数据情况及时调整和改进算法。
小行星轨道问题MATLAB 数学实验(2)共21页文档
6、纪律是自由的第一条件。——黑格 尔 7、纪律是集体的面貌，集体的声音， 集体的 动作， 集体的 表情， 集体的 信念。 ——马 卡连柯
8、我们现在必须完全保持党的纪律， 否则一 切都会 陷入污 泥中。 ——马 克思 9、学校没有纪律便如磨坊没有水。— —夸美 纽斯
7/18
例. 简单迁移模型:每年A镇的人口10%迁往B镇;B镇的 人口15%迁往A镇. 假设某年A、B两镇人口各有120人 和80人.问两年后两镇人口数量分布如何? 设两镇总人口不变,人口流动只限于两镇之间.引入变量:
x1(k) 表示 A 镇第 k 年人口数量; x2(k) 表示 B 镇第 k 年人口数量. 由第 k 年到第 k+1 年两镇人口数量变化规律如下
X(0)
120
80
是特征值1的特征向量
9/18
例4.5 出租汽车问题。
出租汽车公司在仅有A城和B城的海岛上,设了A,B两营 业部。如果周一A城有120辆可出租汽车，而B城有150 辆。统计数据表明，平均每天A城营业部汽车的10% 被顾客租用开到B城 ,B城营业部汽车的12%被开到了 A城。假设所有汽车正常，试计算一周后两城的汽车 数量。寻找方案使每天汽车正常流动而A城和B城的汽 车数量不增不减。
a1x12 + 2a2x1y1 + a3 y12 +2a4 x1 + 2a5 y1 = –1 a1x22 + 2a2x2y2 + a3 y22 +2a4 x2 + 2a5 y2 = –1 a1x32 + 2a2x3y3 + a3 y32 +2a4 x3 + 2a5 y3 = –1 a1x42 + 2a2x4y4 + a3 y42 +2a4 x4 + 2a5 y4 = –1 a1x52 + 2a2x5y5 + a3 y52 +2a4 x5 + 2a5 y5 = –1
小行星轨道问题建模
建模作业五A题小行星轨道问题一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：1.4959787 1011m）。
在5个不同的时间对小行星作了5次观察，测得轨道上5个点的坐标数据如表：x 5.764 6.286 6.759 7.168 7.408y 0.648 1.202 1.823 2.526 3.306由开普勒第一定律知，小行星轨道为一椭圆。
请给出分析、合理假设、建立模型，并根据模型求解确定出椭圆的方程以供研究。
分析：由开普勒定律可知，小行星的轨道为椭圆。
由观测得出五个点坐标可建立五个五元一次方程组。
联立方程组便可得出小行星椭圆轨道方程。
假设：1、假设小行星的轨道为标准椭圆。
2、小行星可视为质点3、小行星轨道不会变化4、忽略其他星体对小行星轨道的影响5、小行星不会与其他星体相撞6、观测准确建立模型： 假设椭圆一般方程为椭圆一般方程：01225432221=+++++y a x a xy a y a xax 5.764 6.286 6.759 7.168 7.408 y0.6481.2021.8232.5263.306将表格中的数带入上面椭圆一般方程：1648.0764.52648.0764.52648.0764.55432221=+⨯+⨯+⨯⨯+⨯+⨯a a a a a1202.1286.62202.1286.62202.1286.65432221=+⨯+⨯+⨯⨯+⨯+⨯a a a a a1823.1759.62823.1759.62823.12759.65432221=+⨯+⨯+⨯⨯+⨯+⨯a a a a a1526.2168.72526.2168.72526.2168.75432221=+⨯+⨯+⨯⨯+⨯+⨯a a a a a1306.3408.72306.3408.72306.3408.75432221=+⨯+⨯+⨯⨯+⨯+⨯a a a a a模型求解程序X=[5.764 6.286 6.759 7.168 7.408]; Y=[0.648 1.202 1.823 2.526 3.360]; A=zeros(5);X(1); for i=1:5A(i,1)=X(i)*X(i);A(i,2)=2*X(i)*Y(i);A(i,3)=Y(i)*Y(i);A(i,4)=2*X(i);A(i,5)=2*Y(i); end;B=[-1 -1 -1 -1 -1]';x=A\BB 题 密码的加密、解密与破译1. 请小组的成员，每人自拟一句话（六字以上）以汉语拼音用密钥矩阵编程写成信息传输密文发布。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
问题15 小行星的轨迹问题
一 、问题
一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道
平面内建立以太阳为远点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：1.4959787*10^11m ），在5个不同的时间对小行星作了5次观察，测得轨道上5个点的坐标数据如表2.15.1. 表2.15.1
由开普勒第一定律知，小行星轨道为一椭圆，现需要建立椭圆的方程以供研究。
（注：椭圆的一般方程可表示为
012225423221=+++++y a x a y a xy a x a ）。
二、实验目的
利用5个点确定二次曲线的一般方程，并求出椭圆的重要参数。
三、预备知识
线性代数方程组理论，椭圆的有关概念及性质。
四、实验内容与要求
1.用表中5个点的坐标数据分别代入椭圆的一般方程可建立5
个方程的线性代数方程组，该方程组的系数矩阵为A ，右端项为b ，这里，
21x 112y x 21y 12x 12y -1
22x
222y x 2
2y 22x
22y -1
A= 23
x 332y x
23y
32x 32y b= -1
2
4x 442y x
24y 42x 42y -1 25x 552y x 2
5y 52x
52y -1
试依据题目所给的5个点的坐标，用计算机计算出矩阵的A 的5*5个数据。
2.利用Matalb 指令A\b 求解5元线性代数方程组，写出椭圆
方程
012225423221=+++++y a x a y a xy a x a 中的
5个待定系数
5
4321,,,,a a a a a 及小行星多所对应的曲
线方程。
3.写出曲线表达式中系数所对应的二阶矩阵和三阶矩阵：
1a
2a 1a 2a 3a
C= D= 2a 3a 5a
2a 3a 4a 5a 1
并利用Matlab 指令eig （C ）求出矩阵C 的特征值，记录数据
=1λ（ ） ，
=2λ（ ）
利用Matlab 指令det （D ）计算二阶导数和三阶行列式的值；
=C （
） ， =D （ ）
4.利用公式计算椭圆的下列参数：
长半轴：a=C
D
1λ=（ ）， 短半轴：b=
C
D
2λ=（ ），
办焦距：c=22b a -=（ ）， 写出椭圆标准方程（ ），
5：利用上面的椭圆有关数据求出小行星轨道的参数：
小行星的近日点距离：h=a-c=（ ）， 小行星的远日点距离：H=2c+h=（ ）， 椭圆轨道周长近似值：L ≈π
⎥⎦
⎤
⎢⎣⎡-+ab b a )(23=（ ）
6．*试在Matlab 环境下利用参数方程： x=acos(t) (t ∈[0, 2π]) y=bsin(t)
绘制出以椭圆中心为原点的椭圆图形（图2.15.1）
图2.15.1
五．思考问题
你能否利用定积分求弧长公式推导出椭圆的周长公式？如果你所得到的是一个定积分表达式，利用这一表达式计算小行星轨道的椭圆周长。
实验15 小行星的轨道问题
一． 问题分析与建立模型
天文学家确定小行星运行的轨道时，他的依据是轨道上五点的坐标数据（x 1,y 1）,(x 2,y 2) ,(x 3,y 3) ,(x 4,y 4) ,(x 5,y 5) 由开普勒大一定律知，小行星轨道为一椭圆，而椭圆属于二次曲线，二次曲线的一般方程为f （x j ，y j ）=0，或a 1x 2+2a 2xy+a 3y 2+2a 4x+2a 5y+1=0.为了确定方程中的五个待定系数，将五个点的坐标分别代入上面的方程，得 f （x j ，y j ）=0 ， （j=1,2,3,4,5） 即
1222
1514213112211-=++++y a x a y a y x a x a
122225242
23222221-=++++y a x a y a y x a x a 12223534233332231-=++++y a x a y a y x a x a 12224544243442241-=++++y a x a y a y x a x a 12225554253552251-=++++y a x a y a y x a x a
这是一个包含五个未知数的线性方程组，写成矩阵形式
21x 112y x 21y 12x 12y a 1 -1
22x
222y x 2
2y 22x 22y a 2 -1
23
x 332y x
23y
32x 32y a 3 = -1
2
4x 442y x
24y 42x 42y a 4 -1 25x 552y x 2
5y 52x
52y a 5 -1
求解这一线性方程组，所得的是一个二次曲线方程，为了指导小行星轨道的一些参数，还必须将二次曲线方程化为椭圆的标准方
程形式：
122
22=+b
y
a x 由于太阳的位置是小行星轨道的一个焦点，这是可以根据椭圆的长半轴a 和短半轴
b 计算出小行星的近日点和远日点距离，以及椭圆周长L
根据二次曲线理论，课的椭圆经过旋转和平移两种变换后的
方程如下：
2
221=+
+C
D y x λλ
所以，圆长半轴：C
D
a 1λ= ；椭圆短半轴：b=
C
D
1λ；椭圆的半焦
距：c=22b a -。
二，计算过程
首先由五个点的坐标数据形成线性方程组的系数矩阵，用
Matlab 的指令求解线性代数方程组。
程序2.17.1.
x=[5.764 6.286 6.759 7.168 7.408] ′; y=[0.648 1.202 1.823 2.526 3.360] ′; b=[-1 -1 -1 -1 -1] ′;
A=[x.^2 2*x.*y y.^2 2*x 2*y],pause a=A\b
在Matlab 环境下运行这段程序，计算机将显示出矩阵A 的内容：
A=
7200
.69600
.142896.112656.509504.550520.53360
.143807
.62127
.363802.516460.35180.133233.36433.246841.454040.25720.124448.11115
.155138.39292.1528.114199.04701
.72237
.33
再击任意键，计算机显示a=0.6143 -0.3400 0.6942 -1.6351 -0.2165 这五个数据即是方程式
01222542
3221=+++++y a x a y a xy a x a
中的五个特定系数：5
4321,,,,a a a a a 的
值。
于是椭圆方程为
01)2165.0(2)6351.1(26942.0)3440.0(26143.022=+-+-++-+y x y xy x [a 1 a 2 a 3 a 4 a 5 ]=[ 0.6143 -0.3400 0.6942 -1.6351 -0.2165] 键入：C=[a(1) a(2); a(2) a(3) ]（回车），可得
1a
2a 0.6143 -0.3440
C= =
2a 3a -0.3440 0.6942
键入：eig （C ）（回车）可得
3080.01=λ，0005.12=λ
键入：det （C ）（回车）可得 |C|=0.3081
键入：D=[a(1) a(2) a(4) ; a(2) a(3) a(5) ; a(4) a(5) 1]
（回车） 可得
D=
15
4532421a a a a a a a a = 1
2165.06351.12165.06942.03440.06351.13440.06143.0------ 键入：（D ） |D|=-1.8203
键入：d=abs(det(D)/det(C)) ;d/0.3080(回车)得 椭圆长半轴 a=19.1834
再键入：d/1.0005（回车）得 椭圆短半轴 b=5.9045 键入：sqrt(19.1834^2-5.9045^2)（回车）得 半焦距：
c=18.2521
小行星近日点距和远日点距可类似得出：h=a-c=0.9313，
H=a+c=37.4355
键入：L=pi*（1.5*（19.1834+5.9045）-sqrt （19.1834*5.9045））
（回车）的 椭圆周长近似值：84.7887
最后用参数反复成绘制出以椭圆圆心为原点的椭圆图形 程序2.17.2 t=0:0.1:2pi+0.1; a=19.1834; b=5.9045; x=a*cos(t); y=b*sin(t); plot(x,y)
Matlab环境下运行程序2.17.2可得椭圆图形
三，结果分析
由五个点的坐标确定椭圆方程中的五个系数实际上只用了计算机求解线性代数方程组的一条简单指令A\b由于二次曲线化为标准性的过程较为复繁，加之椭圆的参数较多即使是用计算机也没有十分简捷的方案。
最后椭圆的周长的准确计算要用到椭圆积分，可以考虑用数值积分解决问题。