【精品】课件---数学实验第四次讲稿.ppt

观察结果 1、该曲线包含两个“圆盘”，每一个都是由螺线形轨道构成。某些轨道几乎是垂直地离开圆盘中一个而进入另一个。 2、随着t的增加，x(t)先绕一个圆盘几圈，然后“跳”到另一个圆盘中。绕第二个圆盘几圈，又跳回原来的圆盘。 并以这样的方式继续下去，在每个圆盘上绕的圈数是随机的。 思考：该空间曲线与初始点x0的选择有关吗？ 1）x0=[0 0.1 0.1]‘;[t0,tf]=[0,30];解向量y 2）x00=[0.01 0.11 0.11]‘;[t0,tf]=[0,30];解向量x y – x = (y1-x1,y2-x2,y3-x3).注：这是这两个向量必须是同维数，并且所取的时间节点应该相同。解决的办法是将[t0,tf]直接定义为如下的形式，比如0：0.1：30。 三，微分方程建模实例 [1]追击路线问题：一艘缉私舰雷达发现距c km处有一艘走私船正以匀速 a 沿直线行驶。缉私舰立即以最大的速度 b 追赶，若用雷达进行跟踪，保持船的瞬时速度方向始终指向走私船，试求缉私舰追逐路线和追上的时间。 0 (c, 0) x y D R 假设 1、建立坐标系； 2、缉私船在(c,0)处发现走私船在(0,0)处； 3、走私船逃跑方向为y轴方向； 4、t时刻，走私船到达R(0,at) 缉私舰到达D(x,y)； [1] 追击路线问题 0 (c, 0) x y D R y x (0, at) α 化简： 对 t求导 化简： （1） [1] 追击路线问题 (2) 结合(1)、(2)得到如下微分方程： 是否存在解析解？ [1] 追击路线问题 1）当a < b时，k=a/b<1, 走私船所用的时间为：t = y/a. 缉私船一定能追上走私船。 2）当a > b时，显然缉私船不可能追上走私船。 讨论 当a = b时呢？ 思考 考虑数值求解法 [1] 追击路线问题 敌艇逃跑速度：a = 60公里/小时； 缉私舰追击速度：b = 80公里/小时 缉私舰发现敌艇时的相距：c = 500公里； 计算得：k = (a/b) = 0.75; 编程计算，并画出追击路线和追击所花费的时间。 假设 [1] 追击路线问题 转化成一阶微分方程组 2. Zhui.m [x,y]=ode23('zhuiji',[500,1],[0,0]); y1=y(:,1); plot(x,y1) 1. Zhuiji.m function f=zhuiji(x,y) f=[y(2);0.75*sqrt(1-y(2)^2)/x]; MATLAB程序 [1] 追击路线问题 262.7 时间： 其中a=60(公里/小时) (小时) [1] 追击路线问题 缉私舰追击路线 敌艇逃跑路线 思考 1、在本问题的求解过程中，假定了走私艇的逃跑方向是正北方向，而初始缉私艇的位置在x轴正向。如果放宽这个假定，也就是当这个夹角是任意角度时，如何建立方程进行求解。 2，如果有多个走私艇在一个位置上进行交易，而缉私艇向该方向追赶。这些走私艇向不同方向四散逃走，问如何安排追赶路线？ [1] 追击路线问题 返 回 范例：[2] 地中海鲨鱼问题 意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究，他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中，发现鲨鱼的比例有明显增加（见下表）。 但为何鲨鱼的比例大幅增加呢？生物学家Ancona无法解释这个现象，于是求助于著名的意大利数学家V.Volterra，希望建立一个食饵—捕食系统的数学模型，定量地回答这个问题. 因为战争使捕鱼量下降，食用鱼增加，显然鲨鱼也随之增加。 战争为什么使鲨鱼数量增加？是什么原因？ 想 [2] 地中海鲨鱼问题 ①x1(t), x2(t)分别是食饵、捕食者在t时刻的数量； ②r1, r2是食饵、捕食者的固有增长率； λ1是捕食者掠取食饵的能力， λ2是食饵对捕食者的供养能力； 1、符号说明： 食饵—捕食系统的数学模型 ① 捕食者的存在使食饵的增长率降低，假设降低的程度与捕食者数量成正比，即 2、基本假设： ②　食饵对捕食者的数量x2起到增长的作用， 其程度与食饵数量x1成正比，即： 综合以上①和②，得到如下模型： 模型(一)：不考虑人工捕获的情况 食饵—捕食系统的数学模型 该 模型反映了在没有人工捕获的自然环境中食饵与捕食者之间的制约关系，没有考虑食饵和捕食者自身的阻滞作用，是Volterra提出的最简单的模型. 给定一组具体数据，用matlab软件求解。 食饵： r1= 1, λ1= 0.1, x10=