
示例代码:
import numpy as np
import pandas as pd
# 参数设定
params = {
'm': 8, # 零配件数量
'n': 2, # 半成品数量
'p': [0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10], # 零配件次品率
'c': [2, 8, 12, 2, 8, 12, 8, 12], # 购买单价
'd': [1, 1, 2, 1, 1, 2, 1, 2], # 检测成本
'a': [8, 8], # 装配成本
'pf': 0.10, # 成品次品率
's': 200, # 市场售价
'l': 40, # 调换损失
'r': 6 # 拆解费用
}
# 决策变量
decisions = [0, 1]
# 利润计算函数
def calculate_profit(D, C, R, params):
cost_parts = sum(params['c'][i] + params['d'][i] * D[i] for i in range(params['m']))
cost_subassemblies = sum(params['a'][j] + params['d'][j + params['m']] * C[j] for j in range(params['n&#