数学实验[—用计算机做数学].pptVIP

数学实验 用计算机做数学 科学实验根据一定目的 运用仪器 设备等物质手段 在人为控制的条件下模拟自然现象 以认识自然界事物的本质和规律为目的和任务 数学实验以数据 图形等为 思想材料 以计算机为手段 以数学软件为实验平台 通过对数学问题和实际问题的探索 得到相应问题的解 并进行计算机模拟 主要内容 首先介绍如何利用mathematica求解微分方程 然后讨论与微分方程相关的几个问题 微分方程的线素场 欧拉方法 最陡下降法 微分方程的稳定性 主要目的 本实验将借助mathematica研究与一阶常微分方程相关的几个问题 让我们了解几类微分方程模型 利用mathematica求解微分方程的方法以及通过几何直观和数值计算的方法了解微分方程解的性态 第二讲微分方程实验 牛顿最早研究天体的运动规律时 假定行星是各自独立的对象 只受到太阳引力的作用 因此牛顿的分析只涉及两个对象 这类问题统称为 二体问题 牛顿采用数学方法研究二体问题时 需要求解的运动方程就是微分方程 1864年英国数学家勒维耶由引力定律建立了一个微分方程模型 经严密分析后断定了海王星的存在性 并确定了海王星在天空中的位置 利用mathematica解微分方程 dsolve y x f x y x y x x 求一阶微分方程的特解 dsolve y x f x y x y x0 y0 y x x 求一阶微分方程的通解 例1利用mathematica求微分方程的通解及在初始条件下的特解 dsolve y x 1 x 2 y x y x x dsolve y x 1 x 2 y x y 0 1 y x x s dsolve y x 1 x 2 y x y 0 1 y x x y x y x s 1 1 plot y x x 1 2 微分方程的积分曲线 求一阶微分方程组的通解 dsolve x t f t x t y t y t g t x t y t x t y t t 求一阶微分方程组的特解 dsolve x t f t x t y t y t g t x t y t x t0 x0 y t0 y0 x t y t t 例2利用mathematica求微分方程组的通解和在初始条件下的特解 dsolve x t 3x t 2y t y t 2x t y t x t y t t dsolve x t 3x t 2y t y t 2x t y t x 0 1 y 0 0 x t y t t 轨线 解曲线 s dsolve x t 3x t 2y t y t 2x t y t x 0 1 y 0 0 x t y t t x t x t s 1 1 y t y t s 1 2 f1 parametricplot3d 0 x t y t t 0 1 f2 parametricplot3d t x t y t t 0 1 show f1 f2 ndsolve y x f x y x y x0 y0 y x x x0 x1 求一阶微分方程的数值解 例3求黎卡提 riccati 方程在初始条件下的数值解 ndsolve y x x 2 y x 2 y 0 0 5 y x x 0 1 实验一 磁场的分布 微分方程的线素场 磁力线演示 磁力线是用来形象地描述磁场状态的一种工具 以力线上某一点的切线方向表示该点的磁场强度的方向 以力线的疏密程度表示磁场的强度 斜率 方程的积分曲线与线素 积分曲线 由mathematica描绘方程的线素场 方向场 clearall f h directionalfield f x xmin xmax y ymin ymax m n h opts module s s1 x y a b k i j lines s solve y b k x a x a 2 y b 2 h 2 x y simplify s1 flatten s x a b k x s1 4 y a b k y s1 3 h1 xmax xmin m h2 ymax ymin n a xmin ih1 b ymin jh2 k f x a y b lines table line a b x a b k y a b k i 0 m j 0 n show graphics lines opts axes automatic aspectratio automatic 这是一个通用程序 只要给定右端函数 线素场的区域和选取点的个数 就可自动生成该方程对应的线素场 f x y x 2 y 2 1directionalfield f x y x 2 2 y 2 2 20 20 0 1 在平面和处分别放置两个正 负单位电荷 则它们在平面上产生一磁场 磁场的分布图 实验二 炮弹飞行的轨迹 欧拉方法 考虑初值问题 取一列点 记则有 欧拉 euler 方法 用欧拉方法和mathematica求微分方程近似解 clearall f euler f x x0 x1 h y y0 opt module points y 0 y0 x n x0 nh y n y n 1 hf x x n 1 y y n 1 n m floor x1 x0 h points table x n y n n 0 m listplot points opt 线素场和欧拉方法近似解的比较 比较方程的解析解与近似解 欧拉方法的改进 实验三 黄土高原上的沟壑 最陡下降法 也许你去过西北的黄土高原 也许从各种媒体上了解那里的景象 到处是沟壑交错 这主要是由于长年雨水的冲刷而形成的 你也许会想到 这些沟壑的分布走向与地形有关 如果你手中有一张当地的地形图 你会发现这些沟壑往往是垂直穿过所倚山峦的等值线 西点军校附近地区地形等高线图 例设有一表面光滑的橄榄球 它的表面形状为由长半轴为6 短半轴为3的椭圆绕其长轴旋转一周所得的旋转椭球面 在无风的细雨天 将球放在室外草坪上 使长轴水平放置 求雨水从椭球面流下的路线方程 建立方程 解方程 其中 quickestdescent a b c x0 y0 h1 h2 opt module f fx fy g ellipsoid s t r l p d descentroad projection generator y f x y csqrt 1 x 2 a 2 y 2 b 2 fx x y d f x y x fy x y d f x y y g x y fy x y fx x y ellipsoid parametricplot3d asin u cos v bsin u sin v ccos u u 0 pi 2 v 0 2pi displayfunction identity s dsolve y x g x y x y x0 y0 y x x y x y x s 1 1 print y x y x t solve f x y 0 y y x s 1 1 x y for i 1 i0 l r p table x y x s 1 1 0 x x0 l h1 projection graphics3d rgbcolor 1 0 0 thickness 0 015 line p d table x y x s 1 1 f x y x s 1 1 x x0 l h1 descentroad graphics3d rgbcolor 1 0 0 thickness 0 015 line d generator table graphics3d rgbcolor 0 1 0 line x y x s 1 1 0 x y x s 1 1 f x y x s 1 1 x x0 l h2 startpoint graphics3d pointsize 0 03 rgbcolor 1 0 0 point x0 y0 f x0 y0 show projection ellipsoid startpoint descentroad generator plotrange a a 0 b 0 c opt x0 0 1 y0 1 5 a 3 b 6 c 3 f1 quickestdescent a b c x0 y0 0 001 0 2 viewpoint 1 2 1 displayfunction identity boxed false f2 quickestdescent a b c x0 y0 0 001 0 2 viewpoint 1 2 1 displayfunction identity boxed false f3 quickestdescent a b c x0 y0 0 001 0 2 viewpoint 1 2 0 7 displayfunction identity boxed false show graphicsarray f1 f2 f3 实验四 战争能爆发吗 微分方程稳定性 两国势均力敌的军事力量互相制约 能保证相互之间的和平共处 一旦某一国的军事力量无限扩张 便会对另一国构成威胁 爆发战争的机会大大增加 影响军事力量增加的速度有三个方面 1 对方国家的军事力量可以刺激本国军事力量的发展 2 本国已有的军事力量会对本国继续扩军产生抑制作用 3 本国人民 军队对对方国人民 军队产生仇视的程度会增加本国扩军 将甲乙两国的军事力量量化为和 解方程得平衡点 1 当时 平衡点是稳定的 此