matlab解全微分方程,通过MATLAB求二阶全微分方程解析解

1.对于二阶全微分方程a，不同的a,b,c取值会求出不同的解析

解，解析解又是由齐次解和特解组成。其中，齐次解由特征方程决定，而特解的决定因素则比较复杂。

2.对于二阶全微分方程的分析，我们大致分为三种情况：

b^2-4ac>0(两个不同的实根)

b^2-4ac=0(两个相同的重根)

b^2-4ac<0(两个不同的复数根)

对三种情况进行MATLAB编程，分析齐次解和特解后，再改变W的值，观察解析解的变化

3.b^2-4ac>0的情况

STEP1：求解析解

s1=dsolve('D2y+3*Dy+2*y=0','y(0)=2,Dy(0)=0','t');

s2=dsolve('D2y+3*Dy+2*y=sin(t)','y(0)=2,Dy(0)=0','t');

s3=dsolve('D2y+3*Dy+2*y=sin(2*t)','y(0)=2,Dy(0)=0','t');

s4=dsolve('D2y+3*Dy+2*y=sin(5*t)','y(0)=2,Dy(0)=0','t');

s5=dsolve('D2y+3*Dy+2*y=sin(13*t)','y(0)=2,Dy(0)=0','t');

s6=dsolve('D2y+3*Dy+2*y=sin(25*t)','y(0)=2,Dy(0)=0','t');

STEP2：绘制图形

(1)求w=1情况下的通解和齐次解

t=1:0.1:10;

s1=4*exp(-t)-2*exp(-2*t) %general solution

s2=-3/10*cos(t)+1/10*sin(t)-11/5*exp(-2*t)+9/2*exp(-t) %special solution

subplot(2,1,1);

plot(t,s2);

xlabel('t')

ylabel('y(t)')

title('general solution ')

subplot(2,1,2);

plot(t,s1);

xlabel('t')

ylabel('y(t)')

title('special solution')