打靶法(含Matlab程序).docVIP

西京学数学软件实验任务书 课程名称 数学软件实验 班级 数 0901 学号 07 姓名 李亚强 实验课题 微分方程组边值问题数值算法 （打靶法，有限差分法） 熟悉微分方程组边值问题数值算法（打靶法，有限差 实验目的 分法） 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中一 实验要求 种语言完成 实验内容 微分方程组边值问题数值算法 （打靶法，有限差分法） 成绩 教师 动方向控制减速的推力，主要的控制量只有一个减速推力，减速还会消耗燃料让登月器的质量减小。所以在极坐标下系统的状态 就是 x‘=[质量 m ，角度 theta ，高度 r，角速度 omega，径速度 v] 这五个量，输入就是减速力 F。先列微分方程， dx/dt=f(x)+B*F，其 中 x 是 5*1 的列向量，质量 dm/dt=-F/2940 ，剩下几个翻下极坐标的手册。把这个动力学模型放到 matlab 里就能求解了，微分方程数值解用 ode45。第一问 F=0，让你求椭圆轨道非常容易。注意附件 1 里说 15 公里的时候速度是 s。算完以后验证一下对不对，对的话就是他了，不对的话说明这个椭圆轨道有进动，到时再说。 算出轨道就能计算减速力了。这时候你随便给个常数减速力到方程里飞船八成都能降落，但不是最优解。想想整个过程，开始降落之前飞船总机械能就那么多，你需要对飞船做负功让机械能 减到 0。题目里写发动机喷出翔的相对速度是一定的， 直觉告诉我飞船速度快的时候多喷一些速度慢的时候少喷一些，可以提高做 负功的效率。 但是多喷也不能超过上限 7500N，所以这就是一个带约束优化问题， matlab 里边有专用的优化函数，用 fmincon 就好。找出最优解以后把过程画出来， 看看 F 可不可以是那 5 个状态量的线性组合，如果是的话就非常 happy，不是的话再说。 三四阶段你可以扯点图像识别，什么二维复利叶分解找平坦区域，怎么一边下降一边根据自身状态调整路径之类的。五六阶段还真不知道说什么。一二阶段肯定是重点啦 误差分析其实还挺难的。可能的误差来源是地球的引力，月亮 绕地球向心加速度，太阳的引力 (可能会很小 )，对自身速度、角度的 测量 误差 ( 比如你测出 自身 当前 速度 100m/s 但实际上是 105m/s) ，控制的时候 F 大小以及角度的误差 (比如你想朝正前方向 喷 2000N 但实际上偏了 2 度而且 F=2010N之类 )。上一问已经求出了最优控制策略和飞船路线，把这些扰动加进去以后算出新的路 线减掉理想路线求偏差，然后随便用个卡尔曼滤波器把误差给校 正 All for Joy 2014/9/13 11:14:38 老师的思路，求大神解答给我一份呀 实验二十七实验报告 实验名称： 微分方程组边值问题数值算法（打靶法，有 限差分法）。 实验目的：进一步熟悉微分方程组边值问题数值算法 （打 靶法，有限差分法）。 实验要求： 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 实验原理 ： 打靶法： 对于线性边值问题 y p( x) y q( x) y f (x)x [a, b] y(a) , y(b) (1) 假设 L 是一个微分算子使： Ly y p( x) y q( x) y 则可得到两个微分方程： Ly1 f ( x) ， y1 (a) ， y1 (a) 0 y1 p(x) y1 q( x) y1 f ( x) ， y1( a) ， y1 ( a) 0 (2) Ly 2 0 ， y2 (a) 0 ， y2 ( a) 1 y2 p(x) y2 q( x) y2 0 ， y2 (a) 0 ， y2 (a) 1 (3) 方程（2），（3）是两个二阶初值问题 .假设 y1 是问题（2） 的解， y2 是问题（ 3）的解，且 y2 (b) 0 ，则线性边值问 题（ 1）的解为 : y(x) y1 (x) y1 (b) y2 (x) 。 y2 (b) 有限差分法： 基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商 来近似， 积分用积分和来近似， 于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程 组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。实验内容： %线性打靶法 function [k,X,Y,wucha,P]=xxdb(dydx1,dydx2,a,b,alpha,beta,h) n=fix((b-a)/