查看: 354|回复: 9|关注: 0

[已解决] 多元方程求解

[复制链接]

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
多元方程组,这些方程均等于0,求x,怎么求解?方程如下

function y=k(x)
y1=8/3 - (16*x(1)*x(6))/(4*x(1)^2 - 4*x(6)^2);
y2=26 - (24*x(1)*x(6) + 16*x(1)*x(7) + 16*x(2)*x(6))/(4*x(1)^2 + 8*x(2)*x(1) - 8*x(6)^2 - 8*x(7)*x(6));
y3=- (57*x(1)*x(6) + 24*x(1)*x(7) + 24*x(2)*x(6) + 16*x(1)*x(8) + 16*x(2)*x(7) + 16*x(3)*x(6))/(9*x(1)^2 + 8*x(1)*x(2) + 8*x(3)*x(1) + 4*x(2)^2 - 20*x(6)^2 - 16*x(6)*x(7) - 8*x(8)*x(6) - 4*x(7)^2) - 273/20;
y4=- (36*x(1)*x(6) + 57*x(1)*x(7) + 57*x(2)*x(6) + 24*x(1)*x(8) + 24*x(2)*x(7) + 24*x(3)*x(6) + 16*x(1)*x(9) + 16*x(2)*x(8) + 16*x(3)*x(7) + 16*x(4)*x(6))/(18*x(1)*x(2) + 8*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) - 40*x(6)*x(7) - 16*x(6)*x(8) - 8*x(6)*x(9) - 8*x(7)*x(8) + 4*x(1)^2 + 4*x(2)^2 - 16*x(6)^2 - 8*x(7)^2) - 593/84;
y5=- (36*x(1)*x(6) + 36*x(1)*x(7) + 36*x(2)*x(6) + 57*x(1)*x(8) + 57*x(2)*x(7) + 57*x(3)*x(6) + 24*x(1)*x(9) + 24*x(2)*x(8) + 24*x(3)*x(7) + 24*x(4)*x(6) + 16*x(1)*x(10) + 16*x(2)*x(9) + 16*x(3)*x(8) + 16*x(4)*x(7) + 16*x(5)*x(6))/(8*x(1)*x(2) + 18*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) - 32*x(6)*x(7) - 40*x(6)*x(8) - 16*x(6)*x(9) - 16*x(7)*x(8) - 8*x(6)*x(10) - 8*x(7)*x(9) + 4*x(1)^2 + 9*x(2)^2 + 4*x(3)^2 - 16*x(6)^2 - 20*x(7)^2 - 4*x(8)^2) - 2211/464;
y6=- (36*x(1)*x(7) + 36*x(2)*x(6) + 36*x(1)*x(8) + 36*x(2)*x(7) + 36*x(3)*x(6) + 57*x(1)*x(9) + 57*x(2)*x(8) + 57*x(3)*x(7) + 57*x(4)*x(6) + 24*x(1)*x(10) + 24*x(2)*x(9) + 24*x(3)*x(8) + 24*x(4)*x(7) + 24*x(5)*x(6) + 16*x(2)*x(10) + 16*x(3)*x(9) + 16*x(4)*x(8) + 16*x(5)*x(7))/(8*x(1)*x(2) + 8*x(1)*x(3) + 18*x(1)*x(4) + 18*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) - 32*x(6)*x(7) - 32*x(6)*x(8) - 40*x(6)*x(9) - 40*x(7)*x(8) - 16*x(6)*x(10) - 16*x(7)*x(9) - 8*x(7)*x(10) - 8*x(8)*x(9) + 4*x(2)^2 + 4*x(3)^2 - 16*x(7)^2 - 8*x(8)^2) - 1465/287;
y7=- (36*x(1)*x(8) + 36*x(2)*x(7) + 36*x(3)*x(6) + 36*x(1)*x(9) + 36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 57*x(1)*x(10) + 57*x(2)*x(9) + 57*x(3)*x(8) + 57*x(4)*x(7) + 57*x(5)*x(6) + 24*x(2)*x(10) + 24*x(3)*x(9) + 24*x(4)*x(8) + 24*x(5)*x(7) + 16*x(3)*x(10) + 16*x(4)*x(9) + 16*x(5)*x(8))/(8*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) + 18*x(1)*x(5) + 18*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) + 8*x(3)*x(5) - 32*x(6)*x(8) - 32*x(6)*x(9) - 32*x(7)*x(8) - 40*x(6)*x(10) - 40*x(7)*x(9) - 16*x(7)*x(10) - 16*x(8)*x(9) - 8*x(8)*x(10) + 4*x(2)^2 + 9*x(3)^2 + 4*x(4)^2 - 16*x(7)^2 - 20*x(8)^2 - 4*x(9)^2) - 3037/636;
y8=- (36*x(1)*x(9) + 36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 36*x(1)*x(10) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 57*x(2)*x(10) + 57*x(3)*x(9) + 57*x(4)*x(8) + 57*x(5)*x(7) + 24*x(3)*x(10) + 24*x(4)*x(9) + 24*x(5)*x(8) + 16*x(4)*x(10) + 16*x(5)*x(9))/(8*x(1)*x(4) + 8*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) + 18*x(2)*x(5) + 18*x(3)*x(4) + 8*x(3)*x(5) + 8*x(4)*x(5) - 32*x(6)*x(9) - 32*x(7)*x(8) - 32*x(6)*x(10) - 32*x(7)*x(9) - 40*x(7)*x(10) - 40*x(8)*x(9) - 16*x(8)*x(10) - 8*x(9)*x(10) + 4*x(3)^2 + 4*x(4)^2 - 16*x(8)^2 - 8*x(9)^2) - 1887/248;
y9=- (36*x(1)*x(10) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 57*x(3)*x(10) + 57*x(4)*x(9) + 57*x(5)*x(8) + 24*x(4)*x(10) + 24*x(5)*x(9) + 16*x(5)*x(10))/(8*x(1)*x(5) + 8*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) + 18*x(3)*x(5) + 8*x(4)*x(5) - 32*x(6)*x(10) - 32*x(7)*x(9) - 32*x(7)*x(10) - 32*x(8)*x(9) - 40*x(8)*x(10) - 16*x(9)*x(10) + 4*x(3)^2 + 9*x(4)^2 + 4*x(5)^2 - 16*x(8)^2 - 20*x(9)^2 - 4*x(10)^2) - 20665/2618;
y10=- (36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 57*x(4)*x(10) + 57*x(5)*x(9) + 24*x(5)*x(10))/(8*x(2)*x(5) + 8*x(3)*x(4) + 8*x(3)*x(5) + 18*x(4)*x(5) - 32*x(7)*x(10) - 32*x(8)*x(9) - 32*x(8)*x(10) - 40*x(9)*x(10) + 4*x(4)^2 + 4*x(5)^2 - 16*x(9)^2 - 8*x(10)^2) - 3345/79;
y11=579/13 - (36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 36*x(4)*x(10) + 36*x(5)*x(9) + 57*x(5)*x(10))/(4*x(4)^2 + 8*x(4)*x(5) + 9*x(5)^2 + 8*x(3)*x(5) - 16*x(9)^2 - 32*x(9)*x(10) - 20*x(10)^2 - 32*x(8)*x(10));
y12=315/31 - (36*x(4)*x(10) + 36*x(5)*x(9) + 36*x(5)*x(10))/(4*x(5)^2 + 8*x(4)*x(5) - 16*x(10)^2 - 32*x(9)*x(10));
y13=15/7 - (36*x(5)*x(10))/(4*x(5)^2 - 16*x(10)^2);
y=[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13]';
end



回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

6

主题

1440

帖子

307

最佳答案
  • 关注者: 70
发表于 2020-6-26 16:50:47 | 显示全部楼层 |此回复为最佳答案
仅供参考
  1. clc;
  2. clear;
  3. syms x [1,10]
  4. y1=8/3 - (16*x(1)*x(6))/(4*x(1)^2 - 4*x(6)^2);
  5. y2=26 - (24*x(1)*x(6) + 16*x(1)*x(7) + 16*x(2)*x(6))/(4*x(1)^2 + 8*x(2)*x(1) - 8*x(6)^2 - 8*x(7)*x(6));
  6. y3=- (57*x(1)*x(6) + 24*x(1)*x(7) + 24*x(2)*x(6) + 16*x(1)*x(8) + 16*x(2)*x(7) + 16*x(3)*x(6))/(9*x(1)^2 + 8*x(1)*x(2) + 8*x(3)*x(1) + 4*x(2)^2 - 20*x(6)^2 - 16*x(6)*x(7) - 8*x(8)*x(6) - 4*x(7)^2) - 273/20;
  7. y4=- (36*x(1)*x(6) + 57*x(1)*x(7) + 57*x(2)*x(6) + 24*x(1)*x(8) + 24*x(2)*x(7) + 24*x(3)*x(6) + 16*x(1)*x(9) + 16*x(2)*x(8) + 16*x(3)*x(7) + 16*x(4)*x(6))/(18*x(1)*x(2) + 8*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) - 40*x(6)*x(7) - 16*x(6)*x(8) - 8*x(6)*x(9) - 8*x(7)*x(8) + 4*x(1)^2 + 4*x(2)^2 - 16*x(6)^2 - 8*x(7)^2) - 593/84;
  8. y5=- (36*x(1)*x(6) + 36*x(1)*x(7) + 36*x(2)*x(6) + 57*x(1)*x(8) + 57*x(2)*x(7) + 57*x(3)*x(6) + 24*x(1)*x(9) + 24*x(2)*x(8) + 24*x(3)*x(7) + 24*x(4)*x(6) + 16*x(1)*x(10) + 16*x(2)*x(9) + 16*x(3)*x(8) + 16*x(4)*x(7) + 16*x(5)*x(6))/(8*x(1)*x(2) + 18*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) - 32*x(6)*x(7) - 40*x(6)*x(8) - 16*x(6)*x(9) - 16*x(7)*x(8) - 8*x(6)*x(10) - 8*x(7)*x(9) + 4*x(1)^2 + 9*x(2)^2 + 4*x(3)^2 - 16*x(6)^2 - 20*x(7)^2 - 4*x(8)^2) - 2211/464;
  9. y6=- (36*x(1)*x(7) + 36*x(2)*x(6) + 36*x(1)*x(8) + 36*x(2)*x(7) + 36*x(3)*x(6) + 57*x(1)*x(9) + 57*x(2)*x(8) + 57*x(3)*x(7) + 57*x(4)*x(6) + 24*x(1)*x(10) + 24*x(2)*x(9) + 24*x(3)*x(8) + 24*x(4)*x(7) + 24*x(5)*x(6) + 16*x(2)*x(10) + 16*x(3)*x(9) + 16*x(4)*x(8) + 16*x(5)*x(7))/(8*x(1)*x(2) + 8*x(1)*x(3) + 18*x(1)*x(4) + 18*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) - 32*x(6)*x(7) - 32*x(6)*x(8) - 40*x(6)*x(9) - 40*x(7)*x(8) - 16*x(6)*x(10) - 16*x(7)*x(9) - 8*x(7)*x(10) - 8*x(8)*x(9) + 4*x(2)^2 + 4*x(3)^2 - 16*x(7)^2 - 8*x(8)^2) - 1465/287;
  10. y7=- (36*x(1)*x(8) + 36*x(2)*x(7) + 36*x(3)*x(6) + 36*x(1)*x(9) + 36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 57*x(1)*x(10) + 57*x(2)*x(9) + 57*x(3)*x(8) + 57*x(4)*x(7) + 57*x(5)*x(6) + 24*x(2)*x(10) + 24*x(3)*x(9) + 24*x(4)*x(8) + 24*x(5)*x(7) + 16*x(3)*x(10) + 16*x(4)*x(9) + 16*x(5)*x(8))/(8*x(1)*x(3) + 8*x(1)*x(4) + 8*x(2)*x(3) + 18*x(1)*x(5) + 18*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) + 8*x(3)*x(5) - 32*x(6)*x(8) - 32*x(6)*x(9) - 32*x(7)*x(8) - 40*x(6)*x(10) - 40*x(7)*x(9) - 16*x(7)*x(10) - 16*x(8)*x(9) - 8*x(8)*x(10) + 4*x(2)^2 + 9*x(3)^2 + 4*x(4)^2 - 16*x(7)^2 - 20*x(8)^2 - 4*x(9)^2) - 3037/636;
  11. y8=- (36*x(1)*x(9) + 36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 36*x(1)*x(10) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 57*x(2)*x(10) + 57*x(3)*x(9) + 57*x(4)*x(8) + 57*x(5)*x(7) + 24*x(3)*x(10) + 24*x(4)*x(9) + 24*x(5)*x(8) + 16*x(4)*x(10) + 16*x(5)*x(9))/(8*x(1)*x(4) + 8*x(2)*x(3) + 8*x(1)*x(5) + 8*x(2)*x(4) + 18*x(2)*x(5) + 18*x(3)*x(4) + 8*x(3)*x(5) + 8*x(4)*x(5) - 32*x(6)*x(9) - 32*x(7)*x(8) - 32*x(6)*x(10) - 32*x(7)*x(9) - 40*x(7)*x(10) - 40*x(8)*x(9) - 16*x(8)*x(10) - 8*x(9)*x(10) + 4*x(3)^2 + 4*x(4)^2 - 16*x(8)^2 - 8*x(9)^2) - 1887/248;
  12. y9=- (36*x(1)*x(10) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 57*x(3)*x(10) + 57*x(4)*x(9) + 57*x(5)*x(8) + 24*x(4)*x(10) + 24*x(5)*x(9) + 16*x(5)*x(10))/(8*x(1)*x(5) + 8*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) + 18*x(3)*x(5) + 8*x(4)*x(5) - 32*x(6)*x(10) - 32*x(7)*x(9) - 32*x(7)*x(10) - 32*x(8)*x(9) - 40*x(8)*x(10) - 16*x(9)*x(10) + 4*x(3)^2 + 9*x(4)^2 + 4*x(5)^2 - 16*x(8)^2 - 20*x(9)^2 - 4*x(10)^2) - 20665/2618;
  13. y10=- (36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 57*x(4)*x(10) + 57*x(5)*x(9) + 24*x(5)*x(10))/(8*x(2)*x(5) + 8*x(3)*x(4) + 8*x(3)*x(5) + 18*x(4)*x(5) - 32*x(7)*x(10) - 32*x(8)*x(9) - 32*x(8)*x(10) - 40*x(9)*x(10) + 4*x(4)^2 + 4*x(5)^2 - 16*x(9)^2 - 8*x(10)^2) - 3345/79;
  14. y11=579/13 - (36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 36*x(4)*x(10) + 36*x(5)*x(9) + 57*x(5)*x(10))/(4*x(4)^2 + 8*x(4)*x(5) + 9*x(5)^2 + 8*x(3)*x(5) - 16*x(9)^2 - 32*x(9)*x(10) - 20*x(10)^2 - 32*x(8)*x(10));
  15. y12=315/31 - (36*x(4)*x(10) + 36*x(5)*x(9) + 36*x(5)*x(10))/(4*x(5)^2 + 8*x(4)*x(5) - 16*x(10)^2 - 32*x(9)*x(10));
  16. y13=15/7 - (36*x(5)*x(10))/(4*x(5)^2 - 16*x(10)^2);
  17. y=[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13]';
  18. s=solve(y==0,[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10])
复制代码

可得结果
s.x1

ans =

5
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-26 17:01:56 | 显示全部楼层
20141303 发表于 2020-6-26 16:50
仅供参考

可得结果

输入s.x1
出现未定义变量 "s" 或类 "s.x1"。
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-26 17:20:07 | 显示全部楼层
阁楼上的小火箭 发表于 2020-6-26 17:01
输入s.x1
出现未定义变量 "s" 或类 "s.x1"。

吧x(1)改为x1统一变量就好了
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

6

主题

1440

帖子

307

最佳答案
  • 关注者: 70
发表于 2020-6-26 17:20:50 | 显示全部楼层
在2019b版本运行正常啊
1.PNG
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

6

主题

1440

帖子

307

最佳答案
  • 关注者: 70
发表于 2020-6-26 17:21:50 | 显示全部楼层
阁楼上的小火箭 发表于 2020-6-26 17:20
吧x(1)改为x1统一变量就好了

那可能是版本问题,在我这可以运算,就没改
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

中级

1204 麦片

财富积分


5001500


0

主题

2767

帖子

258

最佳答案
  • 关注者: 170
发表于 2020-6-26 21:56:18 | 显示全部楼层
多解吧:
1:
x1        25.541834603311
x2        25.5418346033201
x3        25.5418346033247
x4        25.5418346033062
x5        25.5418346032994
x6        12.770917301656
x7        51.0836692066343
x8        12.7709173016592
x9        25.5418346033056
x10        5.1083669206596

2:
x1        48.0689085822221
x2        48.0689085822258
x3        48.0689085822241
x4        48.0689085822202
x5        48.0689085822193
x6        24.0344542911117
x7        96.1378171644467
x8        24.0344542911116
x9        48.0689085822202
x10        9.61378171644389
x11        62.1853201650083
x12        23.7014923011884
x13        61.7574368703609

3:
x1        45.4000618726738
x2        45.4000618724733
x3        45.4000618719761
x4        45.4000618728415
x5        45.4000618734146
x6        22.7000309363709
x7        90.800123744825
x8        22.700030936159
x9        45.4000618730238
x10        9.08001237467787

4:
x1: 44.1343730434392
x2: 44.1343730446726
x3: 44.1343730448411
x4: 44.1343730421187
x5: 44.1343730400391
x6: 22.0671865218445
x7: 88.2687460882428
x8: 22.0671865220386
x9: 44.1343730413197
x10: 8.82687460806033
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-28 13:20:35 | 显示全部楼层
20141303 发表于 2020-6-26 17:20
在2019b版本运行正常啊

前辈,我的实际模型数据都不是整数,然后这样处理后发现一直处于运行状态,下面是数据,麻烦您帮我看下

clc;
clear;
syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10
y1= - (1434141.9807292360491942333490217*x1*x6)/(350130.94404664958219363331291123*x1^2 - 367039.72503822869462774846694618*x6^2) - 2125253851959049/36028797018963968;
y2=1749361779637073/4503599627370496 - (8004.9372223610097635664182183229*x1*x6 + 1434141.9807292360491942333490217*x1*x7 + 1434141.9807292360491942333490217*x2*x6)/(1974.7022761497755957556966456574*x1^2 + 700261.88809329916438726662582246*x2*x1 - 2027.836419934238523100305885565*x6^2 - 734079.45007645738925549693389236*x7*x6);
y3=- (23.11993999196850561454972711548*x1*x6 + 8004.9372223610097635664182183229*x1*x7 + 8004.9372223610097635664182183229*x2*x6 + 1434141.9807292360491942333490217*x1*x8 + 1434141.9807292360491942333490217*x2*x7 + 1434141.9807292360491942333490217*x3*x6)/(5.7300759489573745164830483604488*x1^2 + 3949.4045522995511915113932913148*x1*x2 + 700261.88809329916438726662582246*x3*x1 + 350130.94404664958219363331291123*x2^2 - 5.8300590373004836233593479131794*x6^2 - 4055.67283986847704620061177113*x6*x7 - 734079.45007645738925549693389236*x8*x6 - 367039.72503822869462774846694618*x7^2) - 4692846113499865/144115188075855872;
y4=832770614242175/4503599627370496 - (0.033349678458066127068533966807495*x1*x6 + 23.11993999196850561454972711548*x1*x7 + 23.11993999196850561454972711548*x2*x6 + 8004.9372223610097635664182183229*x1*x8 + 8004.9372223610097635664182183229*x2*x7 + 8004.9372223610097635664182183229*x3*x6 + 1434141.9807292360491942333490217*x1*x9 + 1434141.9807292360491942333490217*x2*x8 + 1434141.9807292360491942333490217*x3*x7 + 1434141.9807292360491942333490217*x4*x6)/(11.460151897914749032966096720898*x1*x2 + 3949.4045522995511915113932913148*x1*x3 + 700261.88809329916438726662582246*x1*x4 + 700261.88809329916438726662582246*x2*x3 - 11.660118074600967246718695826359*x6*x7 - 4055.67283986847704620061177113*x6*x8 - 734079.45007645738925549693389236*x6*x9 - 734079.45007645738925549693389236*x7*x8 + 0.0083069946274930885583100599963803*x1^2 + 1974.7022761497755957556966456574*x2^2 - 0.0083678982736249750473689346108586*x6^2 - 2027.836419934238523100305885565*x7^2);
y5=- (0.000024891988409476623903992286914176*x1*x6 + 0.033349678458066127068533966807495*x1*x7 + 0.033349678458066127068533966807495*x2*x6 + 23.11993999196850561454972711548*x1*x8 + 23.11993999196850561454972711548*x2*x7 + 23.11993999196850561454972711548*x3*x6 + 8004.9372223610097635664182183229*x1*x9 + 8004.9372223610097635664182183229*x2*x8 + 8004.9372223610097635664182183229*x3*x7 + 8004.9372223610097635664182183229*x4*x6 + 1434141.9807292360491942333490217*x1*x10 + 1434141.9807292360491942333490217*x2*x9 + 1434141.9807292360491942333490217*x3*x8 + 1434141.9807292360491942333490217*x4*x7 + 1434141.9807292360491942333490217*x5*x6)/(0.016613989254986177116620119992761*x1*x2 + 11.460151897914749032966096720898*x1*x3 + 3949.4045522995511915113932913148*x1*x4 + 3949.4045522995511915113932913148*x2*x3 + 700261.88809329916438726662582246*x1*x5 + 700261.88809329916438726662582246*x2*x4 - 0.016735796547249950094737869221717*x6*x7 - 11.660118074600967246718695826359*x6*x8 - 4055.67283986847704620061177113*x6*x9 - 4055.67283986847704620061177113*x7*x8 - 734079.45007645738925549693389236*x6*x10 - 734079.45007645738925549693389236*x7*x9 + 0.0000061960526636098115098586266230578*x1^2 + 5.7300759489573745164830483604488*x2^2 + 350130.94404664958219363331291123*x3^2 - 0.00000625*x6^2 - 5.8300590373004836233593479131794*x7^2 - 367039.72503822869462774846694618*x8^2) - 3345647445825393/144115188075855872;
y6=4240998265180857/36028797018963968 - (0.000024891988409476623903992286914176*x1*x7 + 0.000024891988409476623903992286914176*x2*x6 + 0.033349678458066127068533966807495*x1*x8 + 0.033349678458066127068533966807495*x2*x7 + 0.033349678458066127068533966807495*x3*x6 + 23.11993999196850561454972711548*x1*x9 + 23.11993999196850561454972711548*x2*x8 + 23.11993999196850561454972711548*x3*x7 + 23.11993999196850561454972711548*x4*x6 + 8004.9372223610097635664182183229*x1*x10 + 8004.9372223610097635664182183229*x2*x9 + 8004.9372223610097635664182183229*x3*x8 + 8004.9372223610097635664182183229*x4*x7 + 8004.9372223610097635664182183229*x5*x6 + 1434141.9807292360491942333490217*x2*x10 + 1434141.9807292360491942333490217*x3*x9 + 1434141.9807292360491942333490217*x4*x8 + 1434141.9807292360491942333490217*x5*x7)/(0.000012392105327219623019717253246116*x1*x2 + 0.016613989254986177116620119992761*x1*x3 + 11.460151897914749032966096720898*x1*x4 + 11.460151897914749032966096720898*x2*x3 + 3949.4045522995511915113932913148*x1*x5 + 3949.4045522995511915113932913148*x2*x4 + 700261.88809329916438726662582246*x2*x5 + 700261.88809329916438726662582246*x3*x4 - 0.0000125*x6*x7 - 0.016735796547249950094737869221717*x6*x8 - 11.660118074600967246718695826359*x6*x9 - 11.660118074600967246718695826359*x7*x8 - 4055.67283986847704620061177113*x6*x10 - 4055.67283986847704620061177113*x7*x9 - 734079.45007645738925549693389236*x7*x10 - 734079.45007645738925549693389236*x8*x9 + 0.0083069946274930885583100599963803*x2^2 + 1974.7022761497755957556966456574*x3^2 - 0.0083678982736249750473689346108586*x7^2 - 2027.836419934238523100305885565*x8^2);
y7=- (0.000024891988409476623903992286914176*x1*x8 + 0.000024891988409476623903992286914176*x2*x7 + 0.000024891988409476623903992286914176*x3*x6 + 0.033349678458066127068533966807495*x1*x9 + 0.033349678458066127068533966807495*x2*x8 + 0.033349678458066127068533966807495*x3*x7 + 0.033349678458066127068533966807495*x4*x6 + 23.11993999196850561454972711548*x1*x10 + 23.11993999196850561454972711548*x2*x9 + 23.11993999196850561454972711548*x3*x8 + 23.11993999196850561454972711548*x4*x7 + 23.11993999196850561454972711548*x5*x6 + 8004.9372223610097635664182183229*x2*x10 + 8004.9372223610097635664182183229*x3*x9 + 8004.9372223610097635664182183229*x4*x8 + 8004.9372223610097635664182183229*x5*x7 + 1434141.9807292360491942333490217*x3*x10 + 1434141.9807292360491942333490217*x4*x9 + 1434141.9807292360491942333490217*x5*x8)/(0.000012392105327219623019717253246116*x1*x3 + 0.016613989254986177116620119992761*x1*x4 + 0.016613989254986177116620119992761*x2*x3 + 11.460151897914749032966096720898*x1*x5 + 11.460151897914749032966096720898*x2*x4 + 3949.4045522995511915113932913148*x2*x5 + 3949.4045522995511915113932913148*x3*x4 + 700261.88809329916438726662582246*x3*x5 - 0.0000125*x6*x8 - 0.016735796547249950094737869221717*x6*x9 - 0.016735796547249950094737869221717*x7*x8 - 11.660118074600967246718695826359*x6*x10 - 11.660118074600967246718695826359*x7*x9 - 4055.67283986847704620061177113*x7*x10 - 4055.67283986847704620061177113*x8*x9 - 734079.45007645738925549693389236*x8*x10 + 0.0000061960526636098115098586266230578*x2^2 + 5.7300759489573745164830483604488*x3^2 + 350130.94404664958219363331291123*x4^2 - 0.00000625*x7^2 - 5.8300590373004836233593479131794*x8^2 - 367039.72503822869462774846694618*x9^2) - 8068061030834815/576460752303423488;
y8=2524576098304283/72057594037927936 - (0.000024891988409476623903992286914176*x1*x9 + 0.000024891988409476623903992286914176*x2*x8 + 0.000024891988409476623903992286914176*x3*x7 + 0.000024891988409476623903992286914176*x4*x6 + 0.033349678458066127068533966807495*x1*x10 + 0.033349678458066127068533966807495*x2*x9 + 0.033349678458066127068533966807495*x3*x8 + 0.033349678458066127068533966807495*x4*x7 + 0.033349678458066127068533966807495*x5*x6 + 23.11993999196850561454972711548*x2*x10 + 23.11993999196850561454972711548*x3*x9 + 23.11993999196850561454972711548*x4*x8 + 23.11993999196850561454972711548*x5*x7 + 8004.9372223610097635664182183229*x3*x10 + 8004.9372223610097635664182183229*x4*x9 + 8004.9372223610097635664182183229*x5*x8 + 1434141.9807292360491942333490217*x4*x10 + 1434141.9807292360491942333490217*x5*x9)/(0.000012392105327219623019717253246116*x1*x4 + 0.000012392105327219623019717253246116*x2*x3 + 0.016613989254986177116620119992761*x1*x5 + 0.016613989254986177116620119992761*x2*x4 + 11.460151897914749032966096720898*x2*x5 + 11.460151897914749032966096720898*x3*x4 + 3949.4045522995511915113932913148*x3*x5 + 700261.88809329916438726662582246*x4*x5 - 0.0000125*x6*x9 - 0.0000125*x7*x8 - 0.016735796547249950094737869221717*x6*x10 - 0.016735796547249950094737869221717*x7*x9 - 11.660118074600967246718695826359*x7*x10 - 11.660118074600967246718695826359*x8*x9 - 4055.67283986847704620061177113*x8*x10 - 734079.45007645738925549693389236*x9*x10 + 0.0083069946274930885583100599963803*x3^2 + 1974.7022761497755957556966456574*x4^2 - 0.0083678982736249750473689346108586*x8^2 - 2027.836419934238523100305885565*x9^2);
y9=- (0.000024891988409476623903992286914176*x1*x10 + 0.000024891988409476623903992286914176*x2*x9 + 0.000024891988409476623903992286914176*x3*x8 + 0.000024891988409476623903992286914176*x4*x7 + 0.000024891988409476623903992286914176*x5*x6 + 0.033349678458066127068533966807495*x2*x10 + 0.033349678458066127068533966807495*x3*x9 + 0.033349678458066127068533966807495*x4*x8 + 0.033349678458066127068533966807495*x5*x7 + 23.11993999196850561454972711548*x3*x10 + 23.11993999196850561454972711548*x4*x9 + 23.11993999196850561454972711548*x5*x8 + 8004.9372223610097635664182183229*x4*x10 + 8004.9372223610097635664182183229*x5*x9 + 1434141.9807292360491942333490217*x5*x10)/(0.000012392105327219623019717253246116*x1*x5 + 0.000012392105327219623019717253246116*x2*x4 + 0.016613989254986177116620119992761*x2*x5 + 0.016613989254986177116620119992761*x3*x4 + 11.460151897914749032966096720898*x3*x5 + 3949.4045522995511915113932913148*x4*x5 - 0.0000125*x6*x10 - 0.0000125*x7*x9 - 0.016735796547249950094737869221717*x7*x10 - 0.016735796547249950094737869221717*x8*x9 - 11.660118074600967246718695826359*x8*x10 - 4055.67283986847704620061177113*x9*x10 + 0.0000061960526636098115098586266230578*x3^2 + 5.7300759489573745164830483604488*x4^2 + 350130.94404664958219363331291123*x5^2 - 0.00000625*x8^2 - 5.8300590373004836233593479131794*x9^2 - 367039.72503822869462774846694618*x10^2) - 3879526205558635/144115188075855872;
y10=4986590317127023/36028797018963968 - (0.000024891988409476623903992286914176*x2*x10 + 0.000024891988409476623903992286914176*x3*x9 + 0.000024891988409476623903992286914176*x4*x8 + 0.000024891988409476623903992286914176*x5*x7 + 0.033349678458066127068533966807495*x3*x10 + 0.033349678458066127068533966807495*x4*x9 + 0.033349678458066127068533966807495*x5*x8 + 23.11993999196850561454972711548*x4*x10 + 23.11993999196850561454972711548*x5*x9 + 8004.9372223610097635664182183229*x5*x10)/(0.000012392105327219623019717253246116*x2*x5 + 0.000012392105327219623019717253246116*x3*x4 + 0.016613989254986177116620119992761*x3*x5 + 11.460151897914749032966096720898*x4*x5 - 0.0000125*x7*x10 - 0.0000125*x8*x9 - 0.016735796547249950094737869221717*x8*x10 - 11.660118074600967246718695826359*x9*x10 + 0.0083069946274930885583100599963803*x4^2 + 1974.7022761497755957556966456574*x5^2 - 0.0083678982736249750473689346108586*x9^2 - 2027.836419934238523100305885565*x10^2);
y11=8494363061721879/72057594037927936 - (0.000024891988409476623903992286914176*x3*x10 + 0.000024891988409476623903992286914176*x4*x9 + 0.000024891988409476623903992286914176*x5*x8 + 0.033349678458066127068533966807495*x4*x10 + 0.033349678458066127068533966807495*x5*x9 + 23.11993999196850561454972711548*x5*x10)/(0.0000061960526636098115098586266230578*x4^2 + 0.016613989254986177116620119992761*x4*x5 + 5.7300759489573745164830483604488*x5^2 + 0.000012392105327219623019717253246116*x3*x5 - 0.00000625*x9^2 - 0.016735796547249950094737869221717*x9*x10 - 5.8300590373004836233593479131794*x10^2 - 0.0000125*x8*x10);
y12=- (0.000024891988409476623903992286914176*x4*x10 + 0.000024891988409476623903992286914176*x5*x9 + 0.033349678458066127068533966807495*x5*x10)/(0.0083069946274930885583100599963803*x5^2 + 0.000012392105327219623019717253246116*x4*x5 - 0.0083678982736249750473689346108586*x10^2 - 0.0000125*x9*x10) - 4515855507282783/4503599627370496;
y13=- (0.000024891988409476623903992286914176*x5*x10)/(0.0000061960526636098115098586266230578*x5^2 - 0.00000625*x10^2) - 1478743747115975/17592186044416;

y=[y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13]';
s=solve(y==0,[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10])
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-28 13:30:16 | 显示全部楼层
20141303 发表于 2020-6-26 17:20
在2019b版本运行正常啊

原问题可见链接
多元方程求解,传递函数
https://www.ilovematlab.cn/thread-600111-1-1.html
(出处: MATLAB中文论坛)
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-28 13:33:07 | 显示全部楼层
shihe 发表于 2020-6-26 21:56
多解吧:
1:
x1        25.541834603311

对,其实上面的解约分后基本一样,不过您是怎么求出的? 还有您能帮我看下这个问题的实际情况吗,

多元方程求解,传递函数
https://www.ilovematlab.cn/thread-600111-1-1.html
(出处: MATLAB中文论坛)


这个问题的方程组可以看我给楼下的回复,都不是整数,求解一直处于运行状态,怎么回事?
回复此楼 已获打赏: 0 积分

举报

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

本版积分规则

关闭

站长推荐上一条 /4 下一条

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