数学实验---Mathematica.ppt

数学建模 网络模型与统筹模型 数学建模 数学建模 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）当 时 ，平衡点是稳定的，此时由于甲、乙两国的军事力量趋于 和 ,所以在这一时期内两国会和平相处； （2）当 时 ，平衡点是不稳定的，此时由于一方的军事扩张，可能形成战争的威胁. 例如，1909年到1914年，德奥匈联盟与法俄联盟的军备竞争，双方扩军与抑制的情形是 k=l=0.9, α=β=0.2, αβ –k l=0.04-0.81<0 平衡点(0,0)是不稳定的，双方不能和平相处，两个军事集团最终爆发了第一次世界大战! 甲乙双方军事力量的变化 The End ! Thanks ! 数学实验 微分方程实验 数 学 实 验 —用计算机做数学 科学实验 根据一定目的,运用仪器、设备等物质手段，在人为控制的条件下模拟自然现象，以认识自然界事物的本质和规律为目的和任务。 数学实验 以数据、图形等为“思想材料”，以计算机为手段，以数学软件为实验平台，通过对数学问题和实际问题的探索，得到相应问题的解，并进行计算机模拟。　 【主要内容】首先介绍如何利用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}] 微分方程的积分曲线 求一阶微分方程组 的通解 DSo