matlab平面抛射方程,MATLAB 数学实验 第七章 微分方程与计算机模拟 PPT注记

第七章 微分方程与计算机模拟 PPT 注记

(2009-5-24)

只有三个内容，有一定难度和深度。三个内容是：

常微分方程初值问题求数据解和蝴蝶效应(洛仑兹模型)的动态仿真、追击曲线动态仿真、有阻力的抛射曲线实验中阻力系数(电影集结号中的火炮数据)

一、一阶常微分方程

MATLAB 求常微分方程初值问题

???=='0

0)(),(y x y y x f y 数值解方法是首先创建一个函数文件，用以描述一阶微分方程右端的二元函数，然后用ode23()求解对应的微分方程初值问题。由于命令行较少，可以在命令窗口实现

第一条命令

ode23('fun1',[1994,2020],12)

执行后，在图形窗口显示解曲线；

第二条命令

[T,N]=ode23('fun1',[1994,2020],12)执行后，在命令窗口显示数值解的数据

与练习题有关的两个例题 马尔萨斯模型：N dt

dN 15.0=，]3,0[]3,0[),(?∈y x

clear

[x,y]=meshgrid(0:0.2:3);

k=0.15*y;d=sqrt(1+k.^2);

px=1./d;py=k./d;

quiver(x,y,px,py)

axis image

蛇形曲线模型：y x

y 2112-+=' [x,y]=meshgrid(-5:.5:5,-1:.1:1);