阻滞增长模型求解_浅谈用Matlab求解微分方程

这里总结一点matlab在微分方程中的应用，解微分方程有两种解，一种是解析解，一种是数值解，这两种分别对应不同的解法，下面就粗略的介绍一下两种解的解法。

解析解

利用dsolve函数进行求解

syms  x;
s = dsolve('eq1,eq2,...', ’cond1,cond2,...', 'v');
%eq：微分方程
%cond：条件
%v：独立变量
%形如：方程：y'= f(t,y),初值：y(t0) = y0

1.求解析解

求

的解析解

s = dsolve('D2y=3*y+2*x','x');
% D2y用以表示y的二阶导数，默认是以t为自变量的，所以最好指明自变量为x.
syms y(x);
s = dsolve([diff(y,x,2) == 3*y+2*x], [y(0) == 5]) 
% diff内依次是函数、自变量、微分阶数，方程用==表示相等而不是赋值

2.初值问题

求初值问题

s