24小时热门版块排行榜    

查看: 1557  |  回复: 0

youhave

铁虫 (正式写手)

[求助] 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 youhave 的主题更新
信息提示
请填处理意见