function [c,ceq] = mycon(x)
%MYCON 建立不等式约束的等式和不等式的多项式
% [c,ceq] = mycon(x)返回的是关于变量矩阵x的非线性等式和不等式约束
%% 需要引入的电网已知量,包括负荷有功Pd、负荷无功Qd、导纳矩阵Y、线路的有功限制Pl,并分解Y得电导G和电纳B矩阵,这些量设为全局变量
global Pd Qd G B ;
%% 将矩阵x分解成四个不同量的矩阵
global n k; %% 全局变量n、k
Pg=x(1:k);Qg=x(k+1:2*k);V=x(2*k+1:2*k+n);a=x(2*k+n+1:2*(k+n));
%% 扩展Pg、Qg向量成完整的PG、QG向量
global C;
PG=zeros(n,1);QG=zeros(n,1);
for j=1:k
i=C.gen(j,1);
PG(i)=Pg(j);
QG(i)=Qg(j);
end
%% 建立非线性等式约束ceq(x)
ceq=zeros(2*n,1); %% 赋予向量ceq运算空间
ceqp=ceq(1:n);ceqq=ceq(n+1:2*n); %% 将向量分解成有功ceqp和无功ceqq两个向量
for i=1:n
sp=0; sq=0;
for j=1:n
sp0=V(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)));
sp=sp+sp0; %% 计算有功多项式中的求和量
sq0=V(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)));
sq=sq+sq0; %% 计算无功多项式中的求和量
end
ceqp(i)=PG(i)/100-Pd(i)/100-V(i)*sp; %% 计算有功多项式
ceqq(i)=QG(i)/100-Qd(i)/100-V(i)*sq; %% 计算无功多项式
end
ceq=[ceqp;ceqq]; %% 将有功向量ceqp和无功向量ceqq合成为ceq
%% 建立非线性不等式约束c(x)
c=zeros(n,1);
end
opf.zip_opf 分布_变压器_最优控制_潮流约束_负荷参数
版权申诉
5星 · 超过95%的资源 81 浏览量
2022-07-15
04:21:30
上传
评论 1
收藏 4KB ZIP 举报
刘良运
- 粉丝: 70
- 资源: 1万+
最新资源
- EXCEL从图片链接获取图片
- Untitled-1 (1).dwt
- ChatGPT是一款由OpenAI研发的人工智能聊天机器人程序,以下是对其特点、发展历程和功能的详细归纳: 特点 技术驱动:C
- 基于STM32的毕业设计项目可以涵盖多个领域和应用,以下是一个典型的基于STM32的毕业设计项目框架,并结合参考文章中的相关数字
- 在植物大战僵尸杂交版中,代码资源或秘籍指令可以提供各种游戏内的效果 以下是一些常用的代码资源和秘籍指令,以及它们的描述和用途:
- C语言题库-信.doc
- 对于端午节代码资源,你可以考虑以下几个方向: ### 1. 端午节主题的小游戏 你可以创建一个端午节主题的小游戏,比如"捞粽
- 如果你是在寻找编程相关的节日主题代码资源,我可以为你提供一些常见的做法和示例 以下是一些可能会有帮助的方法: ### 1. 制
- NX二次开发uc6496 函数介绍
- YOLO(You Only Look Once)是一个流行的实时目标检测算法,它通过单个神经网络模型同时预测图像中多个目标的位置
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈