查看: 3448|回复: 5

求解方程组，万分感谢！！！

字体大小: 正常放大

15 主题	3 听众	391 积分

升级 30.33%

TA的每日心情

	开心 2016-4-13 15:19

签到天数: 120 天

[LV.7]常住居民III

群组: 学术交流A

电梯直达

1^#

发表于 2011-5-9 23:36 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

syms R r x;
f1=sym('12*r*(R+r)-(R^2+r^2+R*r)*(x^2*r^2+4)');
f2=sym('x^2*(R^2+R+1)-6*(R+1)');
f3=('((R+r)*(R^2+R+1)-R*(R+1)*(R^2+r^2+R*r))^2*(-2*R^2+R+1)/3/(R+1)-((1-R)*(R+r)*(R^2+r^2+R*r)-r*(1-R^2)*(R^2+r^2+R*r))^2');
[R,r,x]=solve(f1,f2,f3)

为什么matlab一直算不出来，永远busy？求解，谢谢了！！！

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

小企鹅

7 主题	3 听众	364 积分

升级 21.33%

TA的每日心情

	郁闷 2011-10-16 19:41

签到天数: 29 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 武汉大学

群组: 华中科技大学

群组: 数与学-建和模

2^#

发表于 2011-5-10 08:05 |只看该作者 |招呼Ta 关注Ta

>> syms R r x
>> f1=sym('12*r*(R+r)-(R^2+r^2+R*r)*(x^2*r^2+4)');
>> f2=sym('x^2*(R^2+R+1)-6*(R+1)');
>> f3=('((R+r)*(R^2+R+1)-R*(R+1)*(R^2+r^2+R*r))^2*(-2*R^2+R+1)/3/(R+1)-((1-R)*(R+r)*(R^2+r^2+R*r)-r*(1-R^2)*(R^2+r^2+R*r))^2');
>> [R,r,x]=solve(f1,f2,f3)
Warning: Warning, solutions may have been lost

R =

-.14304308155148888972226061285614

r =

-1.1565759480515754524712804644619

x =

-2.4207603103462084953299164205756
呵呵我的算出来了，你的定义变量 syms 语句后面不用分号