通过MATLAB求二阶全微分方程解析解.docVIP

对于二阶全微分方程a ，不同的a,b,c取值会求出不同的解析解，解析解又是由齐次解和特解组成。其中，齐次解由特征方程决定，而特解的决定因素则比较复杂。 对于二阶全微分方程的分析，我们大致分为三种情况： b^2-4ac0(两个不同的实根) b^2-4ac=0（两个相同的重根） b^2-4ac0（两个不同的复数根） 对三种情况进行MATLAB编程，分析齐次解和特解后，再改变W的值，观察解析解的变化 b^2-4ac0的情况 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) %gen