matlab练习程序(常微分方程向量场)

十年热恋 提交于 2020-02-10 18:17:10

通过向量场能很直观的看到微分方程所有解的变化规律。

这里随便设了个方程:dx/dt = sin(t)*cos(x)+sin(t)。

由于方程本身就代表了x在t处的斜率,所以:

vt = cos(atan(f));

vx = sin(atan(f));

matlab代码如下:

clear all;
close all;
clc;

t = -5:0.3:5;
x = -5:0.3:5;

f = @(t,x) sin(t).*cos(x)+sin(t);

[t,x] = meshgrid(t,x);

vt = cos(atan(f(t,x)));
vx = sin(atan(f(t,x)));

quiver(t,x,vt,vx);

hold on
[t,x] = ode45(f,[-5 5],[3;-2;-4]);
plot(t,x,'r');

结果如下:

红线为微分方程的三个特解。

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!