使用Juila語言求解Unit Commitment 問題

1. 安裝Julia

image.png
  • 安裝完成后運行界面如下


    image.png

2. 配置IDE

  • VScode Julia + VSCode - 知乎 (zhihu.com)
    • 安裝插件:Julia 和 code runner
    • 在擴展列表中,右擊 Julia,選擇擴展設置(Extension Settings)
    • 將 Julia: Executable Path (julia.executablePath)的值改為 julia.exe 的絕對路徑。
      • 如果路徑使用反斜線則需要轉義(\ -> \),建議使用正斜線(/)
      • 如果之前添加了環境變量,可以只寫 julia // 或不寫

] 進入pkg模式,可以下載所需的包
? 進入help模式
;進入shell模式,可以執行shell命令

3. 安裝所需包

  • 打開Julia, 進入pkg模式(])
  • 輸入以下命令
add Cbc
add JuMP
add UnitCommitment

4. 求解Unit Commitment 問題,示例如下

新建sample.jl文件,運行即可

using Cbc
using JuMP
using UnitCommitment

import UnitCommitment:
    Formulation,
    KnuOstWat2018,
    MorLatRam2013,
    ShiftFactorsFormulation

# Read benchmark instance
instance = UnitCommitment.read_benchmark(
    "matpower/case118/2017-02-01",
)

# Construct model (using state-of-the-art defaults)
model = UnitCommitment.build_model(
    instance = instance,
    optimizer = Cbc.Optimizer,
)

# Construct model (using customized formulation)
model = UnitCommitment.build_model(
    instance = instance,
    optimizer = Cbc.Optimizer,
    formulation = Formulation(
        pwl_costs = KnuOstWat2018.PwlCosts(),
        ramping = MorLatRam2013.Ramping(),
        startup_costs = MorLatRam2013.StartupCosts(),
        transmission = ShiftFactorsFormulation(
            isf_cutoff = 0.005,
            lodf_cutoff = 0.001,
        ),
    ),
)

# Modify the model (e.g. add custom constraints)
@constraint(
    model,
    model[:is_on]["g3", 1] + model[:is_on]["g4", 1] <= 1,
)

# Solve model
UnitCommitment.optimize!(model)

# Extract solution
solution = UnitCommitment.solution(model)
UnitCommitment.write("output.json", solution)
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容