我正在努力使发电资产的每小时利润最大化。 据我所知,我可能会使用quadprog :: solve.QP。
我做了大部分所需的数据准备,已经给了我一个96 x 5的数据框。
这些列包含以下信息:
前两列已完成,这导致我的二次优化。目标函数如下:
max Profit[i] = Volume[i] * (Price[i] - Cost[i])
主要问题是发电成本是生产量的函数(我已预先确定,而且还取决于各种静态值)。
除此之外,某一季度的生产量不得与先前的产量相差不多,假设超过20兆瓦。 生产量不得超过最大产量,不得低于最低产量。
我尝试按如下方式实现优化问题:
体积= x1
价格= x2
费用= x3
Profit = x1 * (x2 - x3) --> max
Profit = x1*x2 - x1*x3 --> max
带
x3 = f(x1)
受
的约束x1(t) >= x1(t-1) - 20
x1(t) <= x1(t-1) + 20
x1 <= max(x1)
x1 >= min(x1)
从quadprog手册中我读到了我需要使用
solve.QP(Dmat, dvec, Amat, bvec)
但老实说,我不知道如何填写两个矩阵和两个向量。
有人可以帮忙吗?
我希望所提供的信息足够。
干杯, 蒂尔曼