Scipy.optimize.linprog文件:值错误inpu无效

2024-10-01 09:24:15 发布

您现在位置:Python中文网/ 问答频道 /正文

我试图解决一个小问题,只是为了在游戏中增加一些单位的产量,其中Alpha是一个变化系数(它设置了变量之间的差异):

import numpy as np
import scipy.optimize as opti
alpha = 0.05
C = np.array([-1,-1,-1,-1,-15,-3,-3,-4,0,0,0,0,0,0])
B = np.array([1600,0,0,0,0,0,0,0,0,0,0,0,0,0])
MatriceC = np.array([\
np.array([14-((1-alpha)*8),7-((1-alpha)*8),7-((1-alpha)*25),18-((1- 
alpha)*12),30-((1-alpha)*30),40-((1-alpha)*40),18-((1-alpha)*1),76-((1- 
alpha)*16),-1,0,0,0,0,0]),\
np.array([14-((1+alpha)*8),7-((1+alpha)*8),7-((1+alpha)*25),18- 
((1+alpha)*12),30-((1+alpha)*30),40-((1+alpha)*40),18-((1+alpha)*1),76- 
((1+alpha)*16),0,-1,0,0,0,0])*(-1),\
np.array([14-((1-alpha)*30),7-((1-alpha)*2),7-((1-alpha)*13),18-((1- 
alpha)*7),30-((1-alpha)*30),40-((1-alpha)*40),18-((1-alpha)*24),76-((1- 
alpha)*56),0,0,-1,0,0,0]),\
np.array([14-((1+alpha)*30),7-((1+alpha)*2),7-((1+alpha)*13),18- 
((1+alpha)*7),30-((1+alpha)*30),40-((1+alpha)*40),18-((1+alpha)*24),76- 
((1+alpha)*56),0,0,0,-1,0,0])*(-1),\
np.array([8-((1-alpha)*30),8-((1-alpha)*2),25-((1-alpha)*13),12-((1- 
alpha)*7),30-((1-alpha)*30),40-((1-alpha)*40),1-((1-alpha)*24),16-((1- 
alpha)*56),0,0,0,0,-1,0]),\
np.array([8-((1+alpha)*30),8-((1+alpha)*2),25-((1+alpha)*13),12- 
((1+alpha)*7),30-((1+alpha)*30),40-((1+alpha)*40),1-((1+alpha)*24),16- 
((1+alpha)*56),0,0,0,0,0,-1])*(-1)])
#print(help(opti.linprog))
print(np.shape(MatriceC))
print(np.shape(B))
opti.linprog(C,A_eq=MatriceC,b_eq=B) #This causes the error...

我得到的结果是:

(6, 14)
(14,)
ValueError: Invalid input for linprog with method = 'simplex'.  The number 
of rows in A_eq must be equal to the number of values in b_eq

考虑到我得到的矩阵的形状。我不明白我做错了什么。你知道吗

附言:

我试过添加

MatriceC = MatriceC.T

就在linprog调用之前,它仍然输出相同的错误。它确实把(6,14)的形状改成了(14,6)(这是合乎逻辑的)


Tags: oftheinimportalphanumberasnp
1条回答
网友
1楼 · 发布于 2024-10-01 09:24:15

在把你的MatriceC传递给linprog之前,先用MatriceC.T转置它

^{}根据他们的文档:

Minimize: c^T * x

Subject to: A_ub * x <= b_ub A_eq * x == b_eq

为了满足上述方程,矩阵的维数必须一致。阅读Matrix Multiplication。你知道吗

相关问题 更多 >