微分方程数值解--上机2.docVIP

微分方程练习 练习1 求解范德堡（vander pol）方程 练习2 单摆运动 图4.3中一根长的细线，一端固定，另一端悬挂质量为 的小球，在重力作用下，小球处于竖直的平衡位置. 现使 小球偏离平衡位置一个小的角度，然后使其自由运动，在不 考虑空气阻力情形下，小球将沿弧线作周期一定的简谐运动. 为平衡位置，在小球摆动过程中，当与平衡位置夹 角为时，小球所受重力在其动运轨迹的分量为 （负号表示力的方向使减少），由牛顿第二定律可得微分方 程 （4.12） 设小球初始偏离角度为，且初速为0，式（4.12）的初 始条件为 （4.13） 当不大时，，式（4.12）化为线性常系数微 分方程 图4.3 （4.14） 解得 （4.15） 简谐运动的周期为. 现在的问题是：当较大时，仍用近似，误差太大，式（4.12）又无解析解，试用数值方法在两种情况下求解，画出的图形，与近似解（4.15）比较，这里设. 练习3 捕食与被捕食 当鲨鱼捕食小鱼，简记为乙捕食甲，在时刻，小鱼的数量