coffeetu 发表于 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.

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

                              

lamborbikini 发表于 2012-3-18 13:06:41

用mathematica试试?

tukaiwu 发表于 2012-3-18 14:27:27

可以用matlab里的solve函数,具体这个函数的用法,你在command window里输下“doc solve”就知道了~~

coffeetu 发表于 2012-3-22 18:44:38

lamborbikini 发表于 2012-3-18 13:06 static/image/common/back.gif
用mathematica试试?

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

coffeetu 发表于 2012-3-22 18:45:55

tukaiwu 发表于 2012-3-18 14:27 static/image/common/back.gif
可以用matlab里的solve函数,具体这个函数的用法,你在command window里输下“doc solve”就知道了~~ ...

我用的就是solve函数求解的,但是就是一直处于Busy状态,给不出结果,也不报错。。。
页: [1]
查看完整版本: 求助:如何解六个方程六个未知数的方程组???