需要用dsolve解一个形如A*Dy*Dy+B*Dy+C*y+E=0的方程,但是符号变量名似乎对通解的形式有影响 >> syms Tw Oe Zr Cs t y >> dsolve('Tw*Dy*Dy+Oe*Dy+Zr*y+Cs=0','t') ans = -((Oe + Oe*lambertw(0, exp(-(Oe + Zr*(C1 + t))/Oe)/Oe))^2 - Oe^2 + 4*Cs*Tw)/(4*Tw*Zr) -((Oe + Oe*lambertw(0, -exp(-(Oe + Zr*(C2 + t))/Oe)/Oe))^2 - Oe^2 + 4*Cs*Tw)/(4*Tw*Zr) -Cs/Zr >> syms A B C E >> dsolve('A*Dy*Dy+B*Dy+C*y+E=0','t') ans = -(4*A*exp(1) - B^2 + (B + B*lambertw(0, exp(-(B + C*(C1 + t))/B)/B))^2)/(4*A*C) -(4*A*exp(1) - B^2 + (B + B*lambertw(0, -exp(-(B + C*(C2 + t))/B)/B))^2)/(4*A*C) -exp(1)/C 两次利用dsolve解的常微分方程形式上是一样的,只是符号名称有变化 但是两次运行程序得到的解似乎是不一样的,比如上面没有exp(1)下面有,那么是否有一种是错的? 如何修改才能让dsolve得到正确的解呢? 希望各位大佬赐教呀 |
1 条回复