示例代码:
import numpy as np
import pandas as pd
# 参数设定
params = {
'p1': 0.10, 'p2': 0.10, 'c1': 4, 'c2': 2, 'd1': 2, 'd2': 3,
'pf': 0.10, 'a': 6, 'df': 3, 's': 56, 'l': 6, 'r': 5
}
# 决策变量
decisions = [0, 1]
# 利润计算函数
def calculate_profit(D1, D2, C, R, params):
cost_parts = (params['c1'] + params['d1'] * D1) + (params['c2'] + params['d2'] * D2)
cost_assembly = params['a'] + params['df'] * C
revenue = params['s'] * (1 - params['pf']) * (1 - C * params['df'])
cost_rework = params['r'] * R * params['pf'] * C
loss_replacement = params['l'] * params[