普通物理建模：4物理过程模型 -2014.pptVIP

物理过程模型 单摆 刚体旋转 陀螺进动 单摆 上机作业 单摆程序 普通物理建模 * 理想情况 摆动周期 m 重力分力： 空气阻力： 驱动力： 运动方程： 无量纲化 常微方程求解 一阶常微分方程的初值问题 要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 初值问题的数值解法的基本特点是：求解过程是沿节点排列的顺序依次向前推进，即按递推方法由已知的y0, y1, … ,yn求出yn+1. 初值问题的数值解法的基本特点是：求解过程是沿节点排列的顺序依次向前推进，即按递推方法由已知的y0, y1, … ,yn求出yn+1. 计算 函数f(x)在区间[a,b]上连续，且原函数为F(x)，则可用牛顿―莱布尼兹公式来求得定积分。 有些原函数是无法用初等函数来表示的?? 有些原函数表达式复杂，不便使用?? 有些只用图形或一张数据表来表示 数值方法 求近似解！ 设想：积分的值由被积函数的值直接决定 回顾积分中值定理： 难以准确算出f(ξ)的值 左矩形公式 中矩形公式 右矩形公式 y x 欧拉公式 左矩形公式 ? 欧拉公式几何意义 x0 x1 (x0,y0) f(x0,y0) (x1,y1) f(x1,y1) y(x) 梯形公式 中矩形公式 Step 1: 先用显式欧拉公式作预测，算出 ) , ( 1 i i i i y x f h y y + = + Step 2: 再将 代入隐式梯形公式的右边作校正，得到 1 + i y )] , ( ) , ( [ 2 1 1 1 + + + + + = i i i i i i y x f y x f h y y 改进欧拉公式 预测-校正系统 欧拉公式 梯形公式 一阶方程组与高阶方程 ? 一阶常微分方程组 一般形式为： = ? = ? )) ( , ... ), ( , ( ) ( ... ... ... )) ( , ... ), ( , ( ) ( 1 1 1 1 x y x y x f x y x y x y x f x y m m m m 初值 0 0 0 2 0 2 0 1 0 1 ) ( , ... , ) ( , ) ( m m y x y y x y y x y = = = 将问题记作向量形式，令： 前述所有公式皆适用于向量形式。 ? 高阶微分方程 = = ? = ? = - - - 1 0 ) 1 ( 1 0 0 0 ) 1 ( ) ( ) ( , ... , ) ( , ) ( ) , ... , , , ( n n n n a x y a x y a x y y y y x f y 化作一阶微分方程组求解。 引入新变量 初值条件为： 二阶方程初值问题 引入新变量z = y’ 改进的欧拉法: （1）无阻尼无驱动的情况 周期与初始角度的关系？ 相图 （2）有阻尼无驱动的情况 相图 （3）有阻尼有驱动的情况 共振现象： q=0.5 b=0.1 Periodic 相图 普通物理建模 *