ADMM算法代码的一些问题

qq_42887487

于 2024-06-18 19:23:34 发布

阅读量259

点赞数 1

文章标签：算法

本文链接：https://blog.csdn.net/qq_42887487/article/details/139779124

版权

% 负荷平衡约束
for j = 1:t
    Constraints1 = [Constraints1;sum(p_g1(:,j))-p_ex12(:,j)+p_ex13(:,j)+p_b1(:,j) == load1(:,j)+p_s1(:,j)];
end

背景介绍下：即有三个微电网，微网可以从主电网购电，或者电力共享，达到负荷平衡

ops = sdpsettings('verbose',2,'debug',1,'solver','gurobi');

% 微网1目标函数
obj1 = sum(sum(a1'*p_g1.^2))+sum(sum(b1'*p_g1))+sum(c1)-lambda(1,:)*p_ex12'+lambda(2,:)*p_ex13'+pri_b*p_b1'-pri_s*p_s1'+sum(0.5*omega(1,:).*(norm(p_ex12-np_ex12,2).^2))+sum(0.5*omega(2,:).*(norm(p_ex13-np_ex13,2).^2));
ans1 = optimizer(Constraints1,obj1,ops,{lambda,np_ex12,np_ex13,omega},{p_ex12,p_ex13,p_g1,p_b1,p_s1});

为什么lambda是与p_ex12，p_ex13相乘，而不是乘以负荷平衡约束

p_ex12是微电网12交互功率