Matlab第二课——常微分方程初步

ODE--常微分方程

ODE是常微分方程的英文缩写，即ordinary diffrential equation，如果在微分方程中，自变量的个数只有一个，这就是ODE方程，例如形如F(x,y,y',y")=0的方程就是一个二阶ODE方程；
PDE方程指偏微分方程，即：partial differential equation，指的是自变量的个数为两个或两个以上的微分方程，如y"(t)+y'(x)=0(这里的导数指的是y对t及x的偏导数)。

‘等差’数列

linspace(start, end, difference)

向量（矩阵）的运算

加减乘除之前一定要有点运算符号啊！！！否则会报错的！

分号的用法

命令行里面输入命令不立即执行，用分号

画图函数

plot（x,y,'style'）

画图

参数为自变量，因变量（不确定，后面说好像是xyz轴的变量），样式（）

hold on

保持图像

legend(‘m1 desciption’,‘m2  desciption’,...)

参数是单引号加图例文字，貌似要按照画图的顺序输入

不过这样的话如果我不想要某个画好的函数的图例怎么办呢？？？凉拌！

用Matlab ODE45解ODE方程

[t,y] = ode45(@odefun,tspan,y0)

t是自变量

y是算出来的因变量

ode45是计算模式

@后面是方程名称

tspan是自变量的计算取值区间，是一个1*2的矩阵

y0是因变量的初始值

定义区间

直接[0, 4*pi]; '''π在matlab里面用pi替换

--------------------

计算一个点的取值

首先得到解 sol = ode45(@reactor,tRange,C0)

然后用deval函数  C16 = deval(sol,1.6);

画出图像

tspan = [0 5];

t = linspace(0,5,100)

C = deval (sol,t)

plot(t,C)

定义微分方程

在脚本中声明

function dxdt = myODEfun(t,x)

enter formula here

end 

一定要写end！！！

用Matlab解微分方程（组）

解微分方程组，要将方程组改写为以上形式

来源：CSDN

作者：雨亮风芳

链接：https://blog.csdn.net/qq_41145832/article/details/79168080