微分方程的基本概念
微分方程:一般的,凡表示未知函数、未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程。
微分方程的阶:微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的阶。
微分方程的解:使得微分方程成立的函数称为微分方程的解。不含任意常数的解称为微分方程的特解。若微分方程的解中所含的相互独立的任意常数的个数与微分方程的阶数相等,称这个解为微分方程的通解。
微分方程的通解:如果微分方程的解中含有任意常数,且任意常数的个数与微分方程的阶数相同,这样的解叫做微分方程的通解。
微分方程的解析解
一般只有比较简单的微分方程才能求出解析解,后面做建模一般都是求数值解。
这里我们讲的基本都是直接用MATLAB程序算的
MATLAB求解命令为
dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’)
记号:在表达微分方程时,用字母D表示求微分,D2、D3等表示求高阶微分.任何D后所跟的字母为因变量,自变量可以指定或由系统规则选定为确省.
例一:求 d u d t = 1 + u 2 \frac{du}{dt}=1+u^2dtdu=1+u2通解
dsolve('Du=1+u^2','t')
tan(C1 + t)%通解
1i%虚数+i、-i数值解
-1i
例二:
{ d 2 y d x 2 + 4 d y d x + 29 y = 0 y ( 0 ) = 0 y ′ ( 0 ) = 15 \begin{cases} \frac{d^2y}{dx^2}+4\frac{dy}{dx}+29y=0 \\ y(0)=0\\ y'(0)=15 \end{cases}⎩⎪⎨⎪⎧dx2d2y+4dxdy+29y=0y(0)=0y′(0)=15
y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')
y =
3*sin(5*x)*exp(-2*x)
[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z', 't');
x=simplify(x) % 将x化简
y=simplify(y)
z=simplify(z)
x =
C1*exp(-t) + C3*exp(2*t)
y =
exp(-2*t)*(C2 + C1*