飞行管理与交调的频率设计详解.ppt

数学建模的一些方法 1提问题法 1）这个问题与什么问题相类似？ 2）改变问题的某些条件会怎样？ 3）将问题分解成若干部分再考虑会怎样？ 4）重新组合又会怎样？ 5）我们还可以做些什么工作？ 6) 还有没有需要进一步完善的内容？ 7）可否换一种数学工具来解决问题？ 抓住问题的关键词，然后不受约束的浮想联翩； 把联想到的内容记在卡片上。 通过对卡片上的内容琢磨，激发产生新的主意。 通过对这些卡片的整理，形成解决问题的初步思路。 3 整体把握问题法 将问题分为“三要素”，既初态、目标态和过程。 初态：察觉到的现在状态。 “有什么”：如已知的条件、数据等。 目标态：察觉到的希望目标。 希望什么、避免什么。 过程：在初态和目标态之间发生作用的行动。 “做什么” 形成研究对象的一个整体结构。 3.模型建立与求解 确定输入输出关系 配置频率 确定信噪比 数据拟合 理论分析、计算机计算 系数的确定 + + + + + + + + + x y y=f(x) (xi,yi) ?i 数据拟合的一般提法为: 已知某函数y=f(x)的一组测量数据 (xi,yi)(i=1,2,…n),要寻求一个函数g(x),使g(x)与上述测量数据在某种准则下最为接近。 即g(xi)≈yi,则可以用g(x)近似代替f(x)。 1 线性最小二乘拟合 问题归结为，求 a1,a2, …am 使 s 最小。 因为a1,a2, …am 都以线性形式出现,故称此方法为线性最小二乘拟合。权系数wi根据需要可灵活选择。 上方程称为正规方程组. 对于给定的测量数据,只要函数系 选得合适,就可以从正规方程组中解出a1,a2, …am . 2.多项式插值拟合： 3.样条函数拟合(略) g(x)=a1g1(x)+ …+amgm(x)中函数{g1(x), …gm(x)}的选取 1. 通过机理分析建立数学模型来确定 g(x)； + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + g=a1+a2x g=a1+a2x+a3x2 g=a1+a2x+a3x2 g=a1+a2/x g=aebx g=ae-bx 2. 将数据 (xi,yi) i=1, …n 作图，通过直观判断确定 g(x)： 3. 取g(x)为x 的多项式函数 (1).输入输出关系的确立 根据题目给出的数据条件,首先要确定输入 输出的函数关系.这是一个数据拟合问题。 由于交调是因为输入u(t)的乘方产生的.故此处用多项式拟合输入输出关系是恰当的。 98.50 87.85 75.10 56.40 35.70 20.15 6.80 2.25 0 输出y 80 60 50 40 30 20 10 5 0 输入u 用什么函数拟合？ 那么,我们试用不同次数的多项式进行拟合来比较,结果发现: 用≥4次的多项式进行拟合时,拟合出的多项式,其次数≥4的项的系数非常小(≤10－５),以致不会对结果产生影响。 故用三次多项式进行拟合已达到精度了。 uk(t)可能产生≤k阶类型的交调; 题目要求考虑二阶和三阶类型的交调; 故最高次数必定≥3。 拟合多项式的最高次数是多少? 注: 拟合函数一般只对内插有效, 外推一般不可靠。 设拟合多项式为: 在所给的数据中有u=0时,y=0.故选取a0=0较好。于是拟合多项式化为: (3.1) (3.2) 用线性最小二乘拟合，对y(t)进行三元回归确定系数。记 (3.3) 则(3.2)式变为 令 求 使S为最小。 由 ，得正规方程组: 解得： 高阶误差为2.6716e-007，故三阶拟合是比较精确的。得输入输出关系式为： 拟合曲线如图： (3.4) 也可以用Mathematica求出拟合函数： 函数形式 功能 . Fit[date，funs，vars] 用变量为vars,函数类为funs,按 最小二乘法拟合一组数据date 程序： date = {{0, 0}, {5, 2.25},