用python如何求解非齐次线性方程组的一组正值特解

wby951

问下像这种方程怎么得到一组正解？
2*x_1 + 3*x_2 + 9*x_3 +...+ 13*x_n = 12345.56384
3*x_1 + 5*x_2 + 6*x_3 +...+ 9*x_n = 93834.798235
6*x_1 + 45*x_2 + 8*x_3 +...+ 0*x_n = 24125.73205
4*x_1 + 2*x_2 + 1*x_3 +...+ 0*dx_n = 124325.90537

一共四个方程，系数啥的是随机的，未知数的量大于4
怎么求得一组匹配方程并都是正数的解
目前个人是用矩阵，调用numpy.linalg.lstsq模块
可以解出来一组解，但是可能有负数。
不胜感激

mswsg

搞不懂数学问题

20032007

问题是怎么约束结果必须是正数吗？
我也遇到这样的问题。。。

foxwb

这是代数问题，和Python没多少关系。对这类欠定方程，我自己遇到的一种类似情况，是获得通解（对于你的这种情况，只有4个方程，应该更好解决，我的情况是，方程数都未定），最后定一组数量很少的系数。不过，我认为有更严谨的解法，我不会。

neodreamerus

这是个线性规划问题吧