华南理工大学数学实验实验三.docVIP

实验三 微分方程 地 点： 4 号楼 4104 房； 实验台号： 66 实验日期与时间： 2017 年 04 月 15 日 评 分： 预 习 检 查 纪 录 ： 实验教师： 刘小兰 电子文档存放位置： 电子文档文件名： 卓越班-66-陈艺荣实验三 批改意见： 1 实验目的 ? 了解求微分方程解析解的方法 ? 了解求微分方程数值解的方法 ? 学习单自由度阻尼系统，观察阻尼系数对系统的影响 ? 模拟弹簧振动，学习实时动画编程的原理 ? 了解 dsolve，ode45 指令的使用方法 ? 了解 Simulink 仿真的设计思想 2 问题 1 2.1 问题描述 用 dsolve 函数求解下列微分方程 ?y??(x) ? y?(x) ? 2y(x) （2） ? ? ? ? y(0) 1, y (0) 0 ? 2.2 实验原理 1、微分方程的解析解 1 解析解指在一定条件下能够以数学表达式直接表达出来的的解。 2、dsolve 命令 用法：dsolve(‘equation’,’condition’,’v’) 说明： （1） equation 是方程式，condition 是条件，v 是自变量（缺省为 t） （2）若不带条件，则解中带积分常数 （3）如果没有显式解，则系统尝试给出隐式解 （4）如果无隐式解，则返回空符号。 格式： （1）y’表示为 Dy，y’’表示为 D2y，依次类推 （2）有多个方程或多个条件时，写多个相应的参数即可。 2.3 算法与编程 2.3.1 编程描述 在 dsolve(‘equation’,’condition’,’v’)函数中，分别编辑输入微分方程、 约束条件、自变量，然后运行 M 文件。 2.3.2 实现代码 %shiyan3_1.m 用 dsolve 函数求解微分方程 %题目说明 用 dsolve 函数求解微分方程 % y 输出的函数因变量 % x 输出的函数自变量 % Chen Yirong 修改于 2017-04-15 % 代码编辑 matlab 版本：MATLAB R2014a clear,clc; %刷新工作区 %--------函数主体--------% 2 y=dsolve(D2y=Dy+2*y,y(0)=1,Dy(0)=0,x) %D2y=Dy+2*y 为方程式，y(0)=1,Dy(0)=0 为约束条件 %--------函数主体--------% 2.4 实验结果 在 matlab 中运行 shiyan3_1，得到图 1 所示的结果 图 1 问题 1 的求解结果 因此，微分方程 y??(x) ? y?(x) ? 2y(x) ? ? ? ? ? y(0) 1, y (0) 0 ? 的特解为： y ? 2e?x ? e2x 3 3 3 问题 2 3.1 问题描述 我缉私雷达发现，距离 d 处有一走私船正以匀速 a 沿直线行驶，缉私舰立即 以最大速度（匀速 v）追赶。若用雷达进行跟踪，保持船的瞬时速度方向始终指 向走私船，则缉私舰的运动轨迹是怎么的？是否能够追上走私船？如果能追上， 需要多长时间？ y M0 M(x, y) d ? S0 S x 图 2 缉私船与走私船的运动模型 3.2 实验原理 在生产和科研中所处理的微分方程往往很复杂且大多得不出一般解。而在实 际上对初值问题，一般是要求得到解在若干个点上满足规定精确度的近似值，或 者得到一个满足精确度要求的便于计算的表达式。 一阶微分方程的数值解 ode45 命令 用法：[t,Y]=ode45(odefun,tspan,y0) 说明： （1） odefun 是待求解一阶微分方程或方程组的句柄，对应一个 M 文件 （2）tspan 求解区间，y0 为初值 （3）返回值 t 为自变量的数据列 （4）返回值 Y 一般是矩阵，每列对应一个待解变量的数据列 （5）对方程组，待解变量，导数，初始值等，全部用数组表示 4 3.3 算法与编程 3.3.1 编程描述 对缉私船追上走私船的问题，为了便于分析，我们用数值法模拟，然后用 MATLAB 求解。 1、基本假设 假设 1 走私船正以匀速v 沿 x 轴正方向直线行驶 1 假设 2 缉私船的速度方向始终指向走私船并立即以最大速度（匀速v ）追赶 0 假设 3 不考虑风浪的影响，走私船和缉私船在这个过程中没有发生故障 2、符号说明 表 1 问题 2 符号说明 符号 符号说明 单位 M(x,y) 缉私船的位置函数 km S 走私船的初始位置，与坐标原点的距离 km 0 S(x,y) 走私船的位置函数 km d 初始时，缉私船与走私船的距离 km ? 缉私船运动方向与走私船运动方向的夹角 rad t 缉私船追上走私船用到的时间 h v 缉私船的最大运动速度值 km / h 0 v 走私船的运动速度值 km / h