确定高精度参数问题.pptxVIP

确定高精度参数问题  (2006研究生数学建模竞赛); 包括“神舟六号”载人航天宇宙飞船、人造地球卫星等航天器围绕地球在轨运行的过程中，要受到很多力的作用，其中主要的是地球万有引力和航天器发动机作用力。 一:考虑航天器在仅受到地球万有引力、航天器自身发动机作用力的作用下作平面运动，将地球和航天器视为质点，试建立航天器运动的数学模型（只要列出模型，不要求解）。 ; 显然这样的数学模型在精度上是远远不能满足实际需要的，在其他要求精确制导等有关高科技的实际问题中，我们都面临着类似的问题：我们必须建立高精度的数学模型，必须高精度地估计模型中的大批参数，因为只有这样的数学模型才能解决实际问题，而不会出现差之毫厘，结果却失之千里的情况。 这时所建立数学模型的精度就成了数学模型的生命线。 ;例如上述问题中的航天器还要受到地球质量分布不均匀所引起的摄动力，大气阻力，日、月及其它星球的摄动引力的影响，以及航天器发动机为调整航天器自身姿态运作时作用力的影响。这样不但数学模型十分复杂，而且在这些数学模型中还要涉及到许多重要的参数，如地球的引力场模型就有许多待定参数。不仅如此，在对航天器进行测量时，还涉及到观测站的地理位置以及设备的系统误差等参数。 为此人们要设法利用长期积累的丰富的观测资料，高精度确定这些重要的参数。由于航天器的问题太复杂，下面本题仅考虑较简单的确定高精度参数问题。; 本题选择了众所周知的Lotka-Volterra生态系统模型.; 通过对此生态系统的观测，可以得到相关的观测数据。观测数据的格式依次为： 观测时刻 、A生物数目 、B生物数目 ;3) 在观测资料有误差（时间变量不含有误差）的情况下，请分别利用观测数据DATA2.TXT和DATA3.TXT，确定参数 在某种意义下的最优解，并与仿真结果比较，进而改进你们的数学模型。 4) 假设连观测资料的时间变量也含有误差，试利用数据DATA4.TXT，建立数学模型，确定参数 在某种意义下的最优解。;二、模型假设;三、符号说明;四、??题分析; 第三个问题，讨论在观测数据有误差的情况下，模型参数的求解。这是一个典型的非线性模型的参数辩识问题，可以采用经典的高斯-牛顿算法进行求解；需要解决两个问题： （1） 的初值确定问题，由于所给的数据量比较大，而且数据观测时间间隔比较短，可以采用灰色估计的方法获得； （2）各参数偏导数的计算，由于无法得到模型的显式表达式，它们无法通过解析的方法获得，考虑通过原微分方程导出各参数偏导数关于时间的一组微分方程，从而利用四阶龙格-库塔求得其数值解。;模型的改进: （1）为了减少算法对初值的依赖性，可以采用改进的高斯-牛顿迭代法进行求解，其思想是通过适当的选取步长，避免迭代过程的波动性； （2）为了避免由于系数矩阵的条件数过大，而导致算法不稳定，可以考虑通过引入阻尼因子，使用阻尼最小二乘法进行迭代求解。;五、模型的建立与求解;;C为任意常数，如果给定了初始条件; 表示初始时刻，令 ，从观测数据中可直接得到 。由于 已知，所以共有四个未知数，即只要四组数据便可求。 也可在数据量充足的情况下, 通过最小二乘法得到一组最优值。（类似多元线性回归）; 利用data1.txt中提供的6组观测数据，可以得到式中 的系数矩阵;ezplot( -2*log(y)+0.2*y-12*log(x)+x+13.82 ,[0,30,0,60]);5.2 问题二的求解——①奇异值分解模型 在该问题中， 未知，共有 和C七个未知数，除去 是初始时刻的值，可用观测数据直接估算，对于剩下的五个未知数，根据线性方程组有解的必要条件，至少需要五组数据才能确定五个参数，这样总共需要六组数据(包括初始数据).;设观测数据共有n组，将数据带入可以得到如下方程：;由线性代数知识可知，当A为非奇异矩阵时，X只有零解；只有当A为奇异矩阵时，X才有非零通解；任意代入五组数据，经过简单计算可知A为近似奇异的，可利用其最小的右奇异向量[1]来确定 和C的值。;5.2 问题二的求解——②迭代模型;　由已知数据data1, 对方程(9) 进行线性方程组的最小二乘拟合可得到:;5.2 问题二的求解——③牛顿迭代法;可利用mat