MATLAB与高等数学--常微分方程(ODE)的求解

常微分方程(ODE)的求解

常微分方程的求解使用dsolve命令

>> clear;
>> syms y(t) a;
>> q = diff(y,t,1)==a*y;
>> s = solve(q)
s =
 
C1*exp(a*t)

我们还可以指定常数值，带入：

>> C1 = 2;a = 4;
>> f  = subs(s)
 
f =
 
2*exp(4*t)

初始条件作为参数放在后面，注意要使用赋值号。：

>> syms f(t) ;
>> q = diff(f,t,1)==f*t/(t-5);
>> s = dsolve(q,f(0)==2)
 
s =
 
-(2*exp(t + 5*log(- 5 + t)))/3125

二阶或者更高阶的处理方法类似：

• 使用syms声明变量:
  由于求解的是微分方程，所以f跟自变量是函数关系，但具体的表达式不清楚，所以我们可以直接声明：
  syms f(t)，后面也可以跟上其它未知数。
• 根据阶数使用diff命令：
  固定用法:

diff(f,t,n)

f代表因变量，t代表自变量，也就是对谁求导，n是方程的阶数。

• 然后写出等式：
  固定格式：

q = 方程左边==方程右边

• 然后使用dsolve命令，并将结果返回给一个变量

当初始条件中含有一阶导数的条件时，我们先要把一阶导数表示出来，然后将条件带入：

>>