声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2102|回复: 1

[综合讨论] MATLAB中dsolve问题请教

[复制链接]
发表于 2010-5-18 15:49 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
截图00.jpg
上面的问题我在MATLAB中编程老是提示我错误,错误如下:
f=dsolve('Dx=-8/3*x+y*z','Dy=-10*y+10*z','Dz=-x*y+28*y-z','x(0)=0','y(0)=0','z(0)=1e-10','t')
??? Error using ==> dsolve
Error, (in dsolve/IC) The 'implicit' option is not available when giving Initial Conditions.

>>  [x,y,z]=dsolve('Dx=-8/3*x+y*z','Dy=-10*y+10*z','Dz=-x*y+28*y-z')
??? Error using ==> dsolve
Too many output arguments.

dsolve('Dx=-8/3*x+y*z','Dy=-10*y+10*z','Dz=-x*y+28*y-z')

ans =

[z(t)=0,x(t)=C1*exp(-8/3*t),y(t)=0,[y(t)=&where(_a,[diff(_b(_a),$(_a,2))*_b(_a)^2-1/3*(-3*diff(_b(_a),_a)^2*_b(_a)*_a-41*_a*diff(_b(_a),_a)*_b(_a)-3*_a^3*_b(_a)-30*_a^4-88*_a*_b(_a)+2160*_a^2+3*_b(_a)^2*diff(_b(_a),_a)+33*_b(_a)^2)/_a=0,_b(_a)=diff(y(t),t),_a=y(t),t=Int(1/_b(_a),_a)+C1,y(t)=_a]),z(t)=1/10*diff(y(t),t)+y(t),x(t)=1/10*(-diff(y(t),$(t,2))-11*diff(y(t),t)+270*y(t))/y(t)]]

提示我在初始化条件上有问题,却不知道问题出在哪里?
回复
分享到:

使用道具 举报

发表于 2010-5-18 16:55 | 显示全部楼层

回复 14楼 zxax1912 的帖子

用dsolve找不到解析解,用数值方法解答,试试ode楼上有的

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-16 15:27 , Processed in 0.059751 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表