Matlab微分方程的求解

求解常微分方程的通解

试解常微分方程：
$$x^2+y+\left(x-2y\right)y^{{}^{\prime }}=0$$
尝试用手算了一下，没算出来。。。

syms  y(x);
f=x^2+y+(x-2*y)*diff(y)
dsolve(f==0)

结果为：
$$\frac{x\pm \sqrt{\frac{4x^3}{3}+x^2+C_1}}{2}$$

求解常微分方程的初边值问题

$$y^{{}^{\prime \prime \prime }}-y^{{}^{\prime \prime }}=x,y^{{}^{\prime }}\left(1\right)=8,y^{{}^{\prime \prime }}\left(2\right)=4$$

syms  y(x);
dy=diff(y);d2y=diff(y,2)
f=diff(y,3)-diff(y,2)-x
y=dsolve(f==0,dy(1)==7,d2y(2)==4)
simplify(y)

simplify(y)可以简化结果
结果为：
$$7e^{x-2}-\left(\frac{7}{e}-\frac{17}{2}\right)x-\frac{x^2}{2}-\frac{x^3}{6}+\frac{1}{6}$$

求解常微分方程组：

$$\{\begin{array}{l}{f}^{{}^{\prime \prime }}+3g=\sin x\\ g^{{}^{\prime }}+f^{{}^{\prime }}=\cos x\end{array}\text{通解和初边值条件为}{f}^{{}^{\prime }}\left(2\right)=0,f\left(3\right)=3,g\left(5\right)=1\text{的解}$$

% 求解常微分方程组
syms f(x) g(x);
df=diff(f)
[f1,g1]=dsolve(diff(f,2)+3*g==sin(x),diff(g)+df==sin(x))
f1=simplify(f1),g1=simplify(g1)
% 带有边值条件
[f2,g2]=dsolve(diff(f,2)+3*g==sin(x),diff(g)+df==sin(x),df(2)==0 ,...
    f(3)==3,g(5)==1)
f2=simplify(f2),g2=simplify(g2)