mathematica求解方程组怎么只返回正的解
mathematica吧
全部回复
仅看楼主
level 2
611715xs 楼主
Solve[{u1*v1 - x1*y1 == 1252, u1*v2 - x2*y1 == 1163,
u1*v3 - x3*y1 == -288, u1*v4 - x4*y1 == 801, u1*v5 - x5*y1 == 973,
u1*v6 - x6*y1 == 1851, u1*v7 - x7*y1 == 952, u1*v8 - x8*y1 == -255,
u1*v9 - x9*y1 == 2817, u1*v10 - x10*y1 == 645,
u1*v11 - x11*y1 == 2192, u1*v12 - x12*y1 == -386,
u1*v13 - x13*y1 == 2474, u1*v14 - x14*y1 == 2100,
u1*v15 - x15*y1 == 1430, u1*v16 - x16*y1 == -539,
u1*v17 - x17*y1 == -79, u2*v1 - x1*y2 == 805, u2*v2 - x2*y2 == 967,
u2*v3 - x3*y2 == -695, u2*v4 - x4*y2 == 506, u2*v5 - x5*y2 == 546,
u2*v6 - x6*y2 == 1618, u2*v7 - x7*y2 == 838, u2*v8 - x8*y2 == -733,
u2*v9 - x9*y2 == 2916, u2*v10 - x10*y2 == 150,
u2*v11 - x11*y2 == 2371, u2*v12 - x12*y2 == -969,
u2*v13 - x13*y2 == 2448, u2*v14 - x14*y2 == 2224,
u2*v15 - x15*y2 == 1463, u2*v16 - x16*y2 == -1016,
u2*v17 - x17*y2 == -281, u3*v1 - x1*y3 == -2228,
u3*v2 - x2*y3 == -1132, u3*v3 - x3*y3 == -1668,
u3*v4 - x4*y3 == -1464, u3*v5 - x5*y3 == -2072,
u3*v6 - x6*y3 == -1464, u3*v7 - x7*y3 == -728,
u3*v8 - x8*y3 == -1980, u3*v9 - x9*y3 == -288,
u3*v10 - x10*y3 == -2280, u3*v11 - x11*y3 == 212,
u3*v12 - x12*y3 == -2396, u3*v13 - x13*y3 == -736,
u3*v14 - x14*y3 == 0, u3*v15 - x15*y3 == -220,
u3*v16 - x16*y3 == -1904, u3*v17 - x17*y3 == -844,
u4*v1 - x1*y4 == 325, u4*v2 - x2*y4 == 689, u4*v3 - x3*y4 == -975,
u4*v4 - x4*y4 == 192, u4*v5 - x5*y4 == 112, u4*v6 - x6*y4 == 1236,
u4*v7 - x7*y4 == 646, u4*v8 - x8*y4 == -1071, u4*v9 - x9*y4 == 2682,
u4*v10 - x10*y4 == -300, u4*v11 - x11*y4 == 2267,
u4*v12 - x12*y4 == -1373, u4*v13 - x13*y4 == 2156,
u4*v14 - x14*y4 == 2088, u4*v15 - x15*y4 == 1331,
u4*v16 - x16*y4 == -1322, u4*v17 - x17*y4 == -427,
u5*v1 - x1*y5 == 1129, u5*v2 - x2*y5 == 959, u5*v3 - x3*y5 == -51,
u5*v4 - x4*y5 == 726, u5*v5 - x5*y5 == 910, u5*v6 - x6*y5 == 1494,
u5*v7 - x7*y5 == 766, u5*v8 - x8*y5 == 3, u5*v9 - x9*y5 == 2088,
u5*v10 - x10*y5 == 690, u5*v11 - x11*y5 == 1583,
u5*v12 - x12*y5 == -53, u5*v13 - x13*y5 == 1880,
u5*v14 - x14*y5 == 1536, u5*v15 - x15*y5 == 1067,
u5*v16 - x16*y5 == -212, u5*v17 - x17*y5 == 23,
u6*v1 - x1*y6 == -3501, u6*v2 - x2*y6 == -1839,
u6*v3 - x3*y6 == -2481, u6*v4 - x4*y6 == -2298,
u6*v5 - x5*y6 == -3234, u6*v6 - x6*y6 == -2418,
u6*v7 - x7*y6 == -1206, u6*v8 - x8*y6 == -2955,
u6*v9 - x9*y6 == -756, u6*v10 - x10*y6 == -3510,
u6*v11 - x11*y6 == 69, u6*v12 - x12*y6 == -3567,
u6*v13 - x13*y6 == -1392, u6*v14 - x14*y6 == -240,
u6*v15 - x15*y6 == -495, u6*v16 - x16*y6 == -2808,
u6*v17 - x17*y6 == -1263, u7*v1 - x1*y7 == -1349,
u7*v2 - x2*y7 == -574, u7*v3 - x3*y7 == -1269,
u7*v4 - x4*y7 == -891, u7*v5 - x5*y7 == -1295,
u7*v6 - x6*y7 == -669, u7*v7 - x7*y7 == -326,
u7*v8 - x8*y7 == -1488, u7*v9 - x9*y7 == 387,
u7*v10 - x10*y7 == -1515, u7*v11 - x11*y7 == 617,
u7*v12 - x12*y7 == -1817, u7*v13 - x13*y7 == -10,
u7*v14 - x14*y7 == 444, u7*v15 - x15*y7 == 143,
u7*v16 - x16*y7 == -1493, u7*v17 - x17*y7 == -628,
u8*v1 - x1*y8 == 1214, u8*v2 - x2*y8 == 950, u8*v3 - x3*y8 == 134,
u8*v4 - x4*y8 == 784, u8*v5 - x5*y8 == 1008, u8*v6 - x6*y8 == 1448,
u8*v7 - x7*y8 == 740, u8*v8 - x8*y8 == 214, u8*v9 - x9*y8 == 1836,
u8*v10 - x10*y8 == 840, u8*v11 - x11*y8 == 1346,
u8*v12 - x12*y8 == 210, u8*v13 - x13*y8 == 1704,
u8*v14 - x14*y8 == 1328, u8*v15 - x15*y8 == 946,
u8*v16 - x16*y8 == 20, u8*v17 - x17*y8 == 110,
u9*v1 - x1*y9 == 2651, u9*v2 - x2*y9 == 1929, u9*v3 - x3*y9 == 631,
u9*v4 - x4*y9 == 1718, u9*v5 - x5*y9 == 2254, u9*v6 - x6*y9 == 2878,
u9*v7 - x7*y9 == 1466, u9*v8 - x8*y9 == 845, u9*v9 - x9*y9 == 3276,
u9*v10 - x10*y9 == 2010, u9*v11 - x11*y9 == 2301,
u9*v12 - x12*y9 == 937, u9*v13 - x13*y9 == 3152,
u9*v14 - x14*y9 == 2320, u9*v15 - x15*y9 == 1705,
u9*v16 - x16*y9 == 488, u9*v17 - x17*y9 == 393,
u10*v1 - x1*y10 == 1769, u10*v2 - x2*y10 == 1516,
u10*v3 - x3*y10 == -111, u10*v4 - x4*y10 == 1137,
u10*v5 - x5*y10 == 1421, u10*v6 - x6*y10 == 2367,
u10*v7 - x7*y10 == 1214, u10*v8 - x8*y10 == -30,
u10*v9 - x9*y10 == 3339, u10*v10 - x10*y10 == 1065,
u10*v11 - x11*y10 == 2539, u10*v12 - x12*y10 == -127,
u10*v13 - x13*y10 == 2998, u10*v14 - x14*y10 == 2460,
u10*v15 - x15*y10 == 1705, u10*v16 - x16*y10 == -373,
u10*v17 - x17*y10 == 22, u11*v1 - x1*y11 == -2659,
u11*v2 - x2*y11 == -1356, u11*v3 - x3*y11 == -1979,
u11*v4 - x4*y11 == -1747, u11*v5 - x5*y11 == -2471,
u11*v6 - x6*y11 == -1757, u11*v7 - x7*y11 == -874,
u11*v8 - x8*y11 == -2350, u11*v9 - x9*y11 == -369,
u11*v10 - x10*y11 == -2715, u11*v11 - x11*y11 == 231,
u11*v12 - x12*y11 == -2843, u11*v13 - x13*y11 == -898,
u11*v14 - x14*y11 == -20, u11*v15 - x15*y11 == -275,
u11*v16 - x16*y11 == -2257, u11*v17 - x17*y11 == -1002,
u12*v1 - x1*y12 == 1805, u12*v2 - x2*y12 == 1453,
u12*v3 - x3*y12 == 105, u12*v4 - x4*y12 == 1164,
u12*v5 - x5*y12 == 1484, u12*v6 - x6*y12 == 2232,
u12*v7 - x7*y12 == 1142, u12*v8 - x8*y12 == 213,
u12*v9 - x9*y12 == 2934, u12*v10 - x10*y12 == 1200,
u12*v11 - x11*y12 == 2179, u12*v12 - x12*y12 == 179,
u12*v13 - x13*y12 == 2692, u12*v14 - x14*y12 == 2136,
u12*v15 - x15*y12 == 1507, u12*v16 - x16*y12 == -94,
u12*v17 - x17*y12 == 121}, {u1, u2, u3, u4, u5, u6, u7, u8, u9, u10,
u11, u12, v1, v2, v3, v4, v5, v6, v7, v8, v9, v10, v11, v12, v13,
v14, v15, v16, v17, y1, y2, y3, y4, y5, y6, y7, y8, y9, y10, y11,
y12, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14,
x15, x16, x17}]
方程是由uv-xy=c构造的 12*17,里面的数,是我自己选择的一组正的u v x y构造出来的,所以肯定会有一组正解
2020年07月21日 12点07分 1
吧务
level 10
限制求解域为PositiveReals
或者解出所有的以后自己用Select挑
2020年07月21日 14点07分 2
请问是在x17后面加上PositiveReals就行了吗
2020年07月22日 04点07分
吧务
level 15
……我觉得你应该先看这帖:https://tieba.baidu.com/p/4850440026
2020年08月01日 03点08分 3
1