VUP 脚本编程(一):内嵌脚本

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
x_Iij_square=sdpvar(32,N,'full');%电流平方 x_ui_square=sdpvar(33,N,'full');%电压平方 x_pij=sdpvar(32,N,'full');%线路有功功率 x_qij=sdpvar(32,N,'full');%线路无功功率 x_Pdg=sdpvar(3,N,'full');%dg功率 u_Pdg=binvar(3,N,'full');%dg出力状态 x_Pil=sdpvar(2,N,'full');%切负荷功率 u_Pil=binvar(2,N,'full');%切负荷状态 x0_ug=binvar(3,N,'full');%机组启动 x0_vg=binvar(3,N,'full');%机组停止 x_R=sdpvar(1,N,'full');%备用负荷 %% 定义约束条件 Constraints=[]; P_dg=sdpvar(33,24,'full');%将dg出力情况放入对应的节点 for i=1:33 if i==7 Constraints=[Constraints,P_dg(i,:)==x_Pdg(1,:)]; elseif i==24 Constraints=[Constraints,P_dg(i,:)==x_Pdg(2,:)]; elseif i==25 Constraints=[Constraints,P_dg(i,:)==x_Pdg(3,:)]; else Constraints=[Constraints,P_dg(i,:)==0]; end end P_il=sdpvar(33,24,'full');%将IL合同约束放入对应的节点 for i=1:33 if i==8 Constraints=[Constraints,P_il(i,:)==x_Pil(1,:)]; elseif i==25 Constraints=[Constraints,P_il(i,:)==x_Pil(2,:)]; else Constraints=[Constraints,P_il(i,:)==0]; end end Constraints=[Constraints,x_ui_square(1,:)==12.66^2]; %平衡节点每小时电压平方 Constraints=[Constraints,x_Iij_square>=0]; % Constraints=[Constraints,x_pij(1,:)>=0]; %% 支路欧姆定律 for r=1:32 Constraints=[Constraints,x_ui_square(Branch(r,2),:)-... x_ui_square(Branch(r,3),:)+(r_ij(r)^2+x_ij(r)^2)*x_Iij_square(r,:)-... 2*(r_ij(r)*x_pij(r,:)+x_ij(r)*x_qij(r,:))==0]; end %dg爬坡约束 for t=1:23 for i=1:3 Constraints=[Constraints,x_Pdg(i,t+1)-x_Pdg(i,t)<=vup]; Constraints=[Constraints,x_Pdg(i,t)-x_Pdg(i,t+1)<=vdn]; end end %IL合同约束 Constraints=[Constraints,0<=x_Pil<=Pilmax.*u_Pil]; %dg最小启停时间约束 for i=1:3 Constraints=[Constraints,consequtiveON([onoffhistory(i,:) u_Pdg(i,:)],Tmup)]; Constraints=[Constraints,consequtiveON(1-[onoffhistory(i,:) u_Pdg(i,:)],Tmd)]; end for t=2:N%公式7这段代码什么意思
06-11

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值