|
|
[求助]
matlab中调用dsolve函数求解微分方程出现多解的问题
求解方程如下:

调用dsolve函数代码如下:
dsolve('D2Ti-DTi*q11/(m*cp)+l*k*q11*Ti/(m*cp)^2=l*k*(q11*T0-q1)/(m*cp)^2','Ti(0)=T0','y')
求得的解:
ans =
C4*exp((y*(q11 - (q11^2 - 4*k*l*q11)^(1/2)))/(2*cp*m)) - exp((y*(q11 + (q11^2 - 4*k*l*q11)^(1/2)))/(2*cp*m))*(C4 - T0 + (2*(k*l*q1 - T0*k*l*q11))/((q11 + (q11^2 - 4*k*l*q11)^(1/2))*(q11^2 - 4*k*l*q11)^(1/2)) - (2*(k*l*q1 - T0*k*l*q11))/((q11 - (q11^2 - 4*k*l*q11)^(1/2))*(q11^2 - 4*k*l*q11)^(1/2))) + (2*exp((y*(q11 + (q11^2 - 4*k*l*q11)^(1/2)))/(2*cp*m))*(k*l*q1 - T0*k*l*q11))/(exp((y*(q11^2 - 4*k*l*q11)^(1/2))/(2*cp*m) + (q11*y)/(2*cp*m))*(q11 + (q11^2 - 4*k*l*q11)^(1/2))*(q11^2 - 4*k*l*q11)^(1/2)) - (2*exp((y*(q11^2 - 4*k*l*q11)^(1/2))/(2*cp*m) - (q11*y)/(2*cp*m))*exp((y*(q11 - (q11^2 - 4*k*l*q11)^(1/2)))/(2*cp*m))*(k*l*q1 - T0*k*l*q11))/((q11 - (q11^2 - 4*k*l*q11)^(1/2))*(q11^2 - 4*k*l*q11)^(1/2))
但是观察微分方程的形式不止一个解,而matlab只返回了一个,刚学matlab,对函数不太懂,请教通晓matlab的虫友!
惠览预谢!
[ Last edited by youhave on 2013-5-4 at 20:24 ] |
|