找回密码
 注册
Simdroid-非首页
查看: 73|回复: 4

求助:如何解六个方程六个未知数的方程组???

[复制链接]
发表于 2012-3-16 10:53:28 | 显示全部楼层 |阅读模式 来自 上海
我的方程组如下:
k1=(R2*(1+(v2*x1+v2*k1^2*(k2^2-1)*x2-2*y1+((k1^2-1)-v2*(k1^2-1)))/(E2*(k1^2-1))+a2*t))/(R1*(1+((1-v1)*y1-v1*x1)/E1+a1*t));
k2=(R3*(1+(-2*y2-v3*(k2^2-1)*x2)/(E3*(k2^2-1))+a3*t))/(R2*(1+(v2*x1+v2*k1^2*(k2^2-1)*x2-2*y1+((k1^2-1)-v2*(k1^2-1)))/(E2*(k1^2-1))+a2*t));
(x1-2*v1*y1)/E1+a1*t_1-(-x1-k1^2*(k1^2-1)*x2+2*v2*y1-2*v2*k1^2*y2)/(E2*(k1^2-1))+a2*t=0;
(-x1-k1^2*(k1^2-1)*x2+2*v1*y1-2*v1*k1^2*y2)/(E2*(k1^2-1))+a2*t=((k2^2-1)*x2+2*v3*y2)/(E3*(k2^2-1))+a3*t;
((1-v1)*y1-v1*x1)/E1+a1*t=((-(k1^2+1)+v2*(k1^2-1))*y1+2*k1^2*y2+v2*x1+v2*k1^2*(k2^2-1)*x2)/(E2*(k1^2-1))+a2*t;
(v2*x1+v2*k1^2*(k2^2-1)*x2-2*y1+((k1^2-1)-v2*(k1^2-1)))/(E2*(k1^2-1))+a2*t=((-(k2^2+1)-v3*(k2^2-1))*y2-v3*(k2^2-1)*x2)/(E3*(k2^2-1))+a3*t;
其中,E1 E2 E3 v1 v2 v3 a1 a2 a3 R1 R2都是已知参数
x1 x2 y1 y2 k1 k2是未知数
我尝试用solve的方法去解,程序没有报错,但是一直在算,算一个小时也没有出结果。
我也尝试用fsolve的方式去解,但是总是报错:fsolve stopped because the problem appears to be locally singular.

求高手赐教,有无办法可以解得次方程组。

                              

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2012-3-18 13:06:41 | 显示全部楼层 来自 美国
Simdroid开发平台
用mathematica试试?
回复 不支持

使用道具 举报

发表于 2012-3-18 14:27:27 | 显示全部楼层 来自 湖南长沙
可以用matlab里的solve函数,具体这个函数的用法,你在command window里输下“doc solve”就知道了~~
回复 不支持

使用道具 举报

 楼主| 发表于 2012-3-22 18:44:38 | 显示全部楼层 来自 上海
lamborbikini 发表于 2012-3-18 13:06
用mathematica试试?

谢谢你的提醒。
我用mathematica算了下,结果也是一直在running...不出结果,郁闷呀:'(
回复 不支持

使用道具 举报

 楼主| 发表于 2012-3-22 18:45:55 | 显示全部楼层 来自 上海
tukaiwu 发表于 2012-3-18 14:27
可以用matlab里的solve函数,具体这个函数的用法,你在command window里输下“doc solve”就知道了~~ ...

我用的就是solve函数求解的,但是就是一直处于Busy状态,给不出结果,也不报错。。。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-7-8 08:28 , Processed in 0.031290 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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