实验7 微分方程建模与计算机动态仿真.pptVIP

武汉理工大学数学系 数学实验与数学软件 实验7 微分方程建模与计算机动态仿真 问题的某些特性随时间或空间而变化，要求 分析其变化规律，预测其未来动态及其控制手 段等，通常需要建立对象的动态或稳态模型， 如：运动轨迹模型、反应扩散等。 本节主要介绍常微分方程建模方法、及确 定性系统的时间步长仿真方法。 前言 计算机仿真是利用计算机对一个实际系统的结构和行为进行动态演示，以评价或预测该系统的行为效果。它是解决较复杂的实际问题的一条有效途径。 计算机仿真在航空、机电、冶金等工业部门及社会经济、交通运输、生态系统等方面有着广泛的应用，已成为分析、研究和设计各种系统的重要手段。 实验7 概率统计与计算机仿真 计算机仿真及其matlab实现 确定性系统的计算机仿真 随机系统的计算机仿真 仿真类型 实验7 概率统计与计算机仿真 从系统特性上，可以分为： 从仿真方法上，可以分为 时间步长法 事件表法 Monte Carlo 方法 例 四人追逐问题 试验1 四人追逐问题（确定性系统的时间步长法仿真） 如图所示，正方形ABCD的四个顶点各有一个人， 在t = 0时刻，四人同时出发以匀速按顺时针方向走向 下一个人。 在行进过程中，他们每个人都始终保持对准各自的目标。 试用计算机模拟出每个人的行进轨迹。 A B C D 实验7 概率统计与计算机仿真 AB位置坐标的连线确定了 A的追逐方向? 假设在某一时刻t，A在(x1,y1)，B在(x2,y2), A(x1,y1) B(x2,y2) 前进方向 时刻，A的新位置 (x1’,y1’)： 求解分析： 依次类推可得出B、C、D的新时刻的位置。 实验7 概率统计与计算机仿真 若t时刻，A的坐标(x1,y1)， 则 时刻，A的坐标(x1’,y1’)： 在时间间隔无穷小时，两边求极限，将此时差分--微分，则有： 因涉及到四人，本题适宜建立极坐标下的微分方程，不做详细介绍。 仿真过程： ⑴初始化：仿真总时间t，采样时间间隔dt，行进速度v， 各点的起始位置； ⑵循环次数：k为循环变量，n = t/dt为循环次数； ⑶位置纪录：生成两个2维数组XX[k, i] 和YY[k，i] 存储第i个人的位置横坐标为XX(:,i)，纵坐标为YY(:,i)， (4)在每个新的时刻，记录下各点的位置，并添加到数组中 的最后一个元素。 解释和运行程序：zhuiji.m 实验7 概率统计与计算机仿真 思考： 在两人追逐问题中，思考一些逃跑策略或追逐策略，提 出一个实际问题，并建立数学模型． 若四个人速度不同，程序又应作如何变化？ 试验2 导弹追踪问题（确定性系统的时间步长法仿真） 设位于坐标原点的导弹发现正北方向(1,0)处海面上有敌艇一艘 以v0米/秒的速度向正东方向行驶。该基地立即发射导弹跟 踪追击敌艇，导弹速率为5v0米/秒，自动导航系统使导弹在 任一时刻都能对准敌艇。 （1）设导弹在t时刻的位置为 ，试推导 的理论表达式及x,y的初值 （2）分别给出在t时刻及在t+?t时刻，敌舰和导弹各自的位置； （3）现要求以时间间隔?t为0.01秒，利用计算机进行仿真得到 导弹和敌舰的运行轨迹，直到导弹击中敌艇为止。 一. 微分方程建模 由导弹的速度为v2，则： 由导弹的方向对准敌舰，有： 联立(1),(2)，令v0=1,可得导弹运动轨迹的参数方程为： （2）分别给出在t时刻及在t+?t时刻，敌舰和导弹各自的位置； 设t时刻，敌舰在(x２,y２)处，t时刻，导弹的位置为(x１,y１) 此时它的运动方向为：(x1,y1)? (x2,y2)， t+?t时，敌舰新位置为(1, y2+v0*?t ) ， 此时导弹的位置为： 其中 ， 随着时间的增加，直到两者的坐标相等。 问题的理论求解（通过联立两个方程）： 得其解析解为：