使用 UPFC 计算电力系统网络潮流（Matlab代码实现）

news2024/9/29 7:20:54

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

本文为电力系统网络执行潮流计算，同时集成了统一潮流控制器 （UPFC） 的功能。UPFC是一种灵活的交流输电系统（FACTS）设备，可以调节输电线路参数，如阻抗和相位角，以控制功率流和电压曲线。

关键步骤：

  1. 初始化：该代码使用与发电机设置、负载需求、输电线路和 UPFC 相关的数据初始化电力系统网络。
  2. 潮流方程的制定：潮流方程是根据基尔霍夫电流定律（KCL）和基尔霍夫电压定律（KVL）建立的，考虑到UPFC的存在来控制特定区域的潮流。
  3. UPFC控制：UPFC调节功率流和电压曲线，以提高系统稳定性。它动态调整传输线的阻抗和相位角。
  4. 迭代求解：采用迭代方法利用UPFC模型求解非线性潮流方程。该代码使用雅可比矩阵迭代更新电压幅度和角度（VM 和 DL），直到实现收敛。
  5. 收敛检查：通过比较连续迭代的结果来验证收敛性。潮流计算将继续进行，直到满足指定的收敛标准。

📚2 运行结果


itn =

     1


itn =

     2


itn =

     3


itn =

     4

THE SOLUTION IS ------------------------

minsing =

    0.2198


Pmk =

   -0.4000


Qmk =

   -0.2500


VMcr =

    0.0489


DLcrD =

 -129.6673


VMvr =

    0.9708


DLvrD =

   -6.4230


Pbb =

  -9.6091e-08


Pcr =

   5.1958e-04


Pvr =

  -5.1968e-04


Qcr =

    0.0230


Qvr =

   -0.2837


VM =

    1.0500
    1.0450
    1.0100
    1.0500
    1.0100
    1.0500
    0.9993
    1.0171
    1.0000
    0.9939
    1.0087
    1.0330
    1.0054
    1.0144
    1.0064
    1.0085
    0.9932
    0.9894
    0.9825
    0.9846
    0.9820
    0.9830
    0.9894
    0.9754
    0.9844
    0.9661
    0.9990
    1.0010
    0.9786
    0.9669
    1.0019


DLD =

         0
   -1.7871
   -5.5002
   -4.1493
   -6.4175
   -6.3774
   -6.1804
   -3.5264
   -6.4199
   -8.3436
   -4.2953
   -7.8566
   -5.1138
   -8.8013
   -8.8296
   -8.2868
   -8.5811
   -9.4083
   -9.5479
   -9.3068
   -8.8566
   -8.8489
   -9.2266
   -9.3918
   -9.5838
  -10.0326
   -9.4077
   -5.6022
  -10.6993
  -11.6280
   -6.3644


P =

    0.9919
    0.5830
   -0.1000
    0.2000
   -0.4420
    0.2000
   -0.2280
   -0.0240
   -0.0000
   -0.0580
   -0.0760
   -0.1120
    0.0000
   -0.0620
   -0.0820
   -0.0350
   -0.0900
   -0.0320
   -0.0950
   -0.0220
   -0.1750
   -0.0000
   -0.0320
   -0.0870
    0.0000
   -0.0350
    0.0000
    0.0000
   -0.0240
   -0.1060
    0.0000


Q =

   -0.0351
    0.3490
    0.1583
    0.2564
    0.0058
    0.1297
   -0.1090
   -0.0120
    0.0000
   -0.0200
   -0.0160
   -0.0750
    0.0000
   -0.0160
   -0.0250
   -0.0180
   -0.0580
   -0.0090
   -0.0340
   -0.0070
   -0.1120
    0.0000
   -0.0160
   -0.0670
    0.0000
   -0.0230
    0.0000
    0.0000
   -0.0090
   -0.0190
   -0.0000


Ploss =

    0.0579


Qloss =

    0.2191


Lj =

    0.0237
    0.0168
    0.0615
    0.0909
    0.0198
    0.0548
    0.0182
    0.0799
    0.0864
    0.0784
    0.0941
    0.1073
    0.1148
    0.1102
    0.1066
    0.1058
    0.1045
    0.1180
    0.1131
    0.1333
    0.1004
    0.0237
    0.1328
    0.1547
    0.0594


t =

    0.1190

>>

部分代码：

JM3;
JM4=zeros(nb-g+1,nb-g+1);
for im4=1:nb-g
    for jm4=1:nb-g
        JM4(im4,jm4)=JM4(im4,jm4)+J4(im4,jm4);
    end
end
JM4;

JM1(ku-1,ku-1)=JM1(ku-1,ku-1)+PkDLk;
JM1(ku-1,nb)=JM1(ku-1,nb)+PkDLCR;
JM1(ku-1,nb+1)=JM1(ku-1,nb+1)+PkDLVR;
JM1(mu-1,mu-1)=JM1(mu-1,mu-1)+PmDLm;
JM1(mu-1,nb)=JM1(mu-1,nb)+PmDLCR;
JM1(nb,ku-1)=JM1(nb,ku-1)+PmkDLk;
JM1(nb,mu-1)=JM1(nb,mu-1)+PmkDLm;
JM1(nb,nb)=JM1(nb,nb)+PmkDLCR;
JM1(nb+1,ku-1)=JM1(nb+1,ku-1)+PBBDLk;
JM1(nb+1,mu-1)=JM1(nb+1,mu-1)+PBBDLm;
JM1(nb+1,nb)=JM1(nb+1,nb)+PBBDLCR;
JM1(nb+1,nb+1)=JM1(nb+1,nb+1)+PBBDLVR;

JM2(ku-1,ku-g)=JM2(ku-1,ku-g)+PkVk;
JM2(ku-1,nb-g+1)=JM2(ku-1,nb-g+1)+PkVCR;
JM2(mu-1,mu-g)=JM2(mu-1,mu-g)+PmVm;
JM2(mu-1,nb-g+1)=JM2(mu-1,nb-g+1)+PmVCR;
JM2(nb,ku-g)=JM2(nb,ku-g)+PmkVk;
JM2(nb,mu-g)=JM2(nb,mu-g)+PmkVm;
JM2(nb,nb-g+1)=JM2(nb,nb-g+1)+PmkVCR;
JM2(nb+1,ku-g)=JM2(nb+1,ku-g)+PBBVk;
JM2(nb+1,mu-g)=JM2(nb+1,mu-g)+PBBVm;
JM2(nb+1,nb-g+1)=JM2(nb+1,nb-g+1)+PBBVCR;

JM3(ku-g,ku-1)=JM3(ku-g,ku-1)+QkDLk;
JM3(ku-g,nb)=JM3(ku-g,nb)+QkDLCR;
JM3(ku-g,nb+1)=JM3(ku-g,nb+1)+QkDLVR;
JM3(mu-g,mu-1)=JM3(mu-g,mu-1)+QmDLm;
JM3(mu-g,nb)=JM3(mu-g,nb)+QmDLCR;
JM3(nb-g+1,ku-1)=JM3(nb-g+1,ku-1)+QmkDLk;
JM3(nb-g+1,mu-1)=JM3(nb-g+1,mu-1)+QmkDLm;
JM3(nb-g+1,nb)=JM3(nb-g+1,nb)+QmkDLCR;

JM4(ku-g,ku-g)=JM4(ku-g,ku-g)+QkVk;
JM4(ku-g,nb-g+1)=JM4(ku-g,nb-g+1)+QkVCR;
JM4(mu-g,mu-g)=JM4(mu-g,mu-g)+QmVm;
JM4(mu-g,nb-g+1)=JM4(mu-g,nb-g+1)+QmVCR;
JM4(nb-g+1,ku-g)=JM4(nb-g+1,ku-g)+QmkVk;
JM4(nb-g+1,mu-g)=JM4(nb-g+1,mu-g)+QmkVm;
JM4(nb-g+1,nb-g+1)=JM4(nb-g+1,nb-g+1)+QmkVCR;

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]陈淮金,白中.含UPFC电力系统的潮流计算研究[J].电力系统自动化, 1996, 20(3):5.DOI:CNKI:SUN:DLXT.0.1996-03-005.

[2]郎兵.含UPFC电力系统的潮流计算研究[J].华北电力技术, 1998(7):5.DOI:10.3969/j.issn.1003-9171.1998.07.004.

[3]武历忠,徐诚.含有UPFC的电力系统潮流计算[C]//2015年云南电力技术论坛论文集(下册).2015.DOI:CNKI:SUN:YNDJ.0.2016-04-040.

🌈4 Matlab代码实现

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/984520.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

Ubuntu yolov5 环境配置

Ubuntu yolov5 环境配置

查看Ubuntu版本 $ cat /proc/version Linux version 5.4.0-150-generic (builddbos03-amd64-012) (gcc version 7.5.0 (Ubuntu 7.5.0-3ubuntu1~18.04)) #167~18.04.1-Ubuntu SMP Wed May 24 00:51:42 UTC 2023虚拟机磁盘扩容 因为在环境搭建过程中遇到了磁盘空间不足的问题&a…
阅读更多...
索尼 toio™ 应用创意开发征文｜探索创新的玩乐世界——索尼 toio™

索尼 toio™ 应用创意开发征文｜探索创新的玩乐世界——索尼 toio™

导语： 在技术的不断进步和发展中，玩具也逐渐融入了智能化的潮流。索尼 toio™作为一款前沿的智能玩具，给孩子和成人带来了全新的游戏体验。本文将介绍索尼 toio™的特点、功能和应用场景，让读者了解这个令人兴奋的创新产品。 1. 了…
阅读更多...
抢先一步，新华三携手HPE正在重塑AI存储新格局

抢先一步，新华三携手HPE正在重塑AI存储新格局

当前智能IT的创新 正呈现出蓬勃发展之势 【全球存储观察 ｜ 热点关注】 当前，数字经济在中国获得了长足进步，并驱动着算力快速进化。为此，智能IT的创新也呈现出了前所未有的蓬勃之势。 从智能存储中枢发布到现在，新华…
阅读更多...
LAMP搭建wordpress并使用reids加速网页

LAMP搭建wordpress并使用reids加速网页

L linux A apache hhtpd M mysql/maridb P PHP1、 安装php rpm -ivh http://rpms.famillecollet.com/enterprise/remi-release-7.rpm yum install -y --enablereporemi --enablereporemi-php72 php php-opcache php-devel php-mysqlnd php-gd php-redis2、 安装mysql5.7 2.1、…
阅读更多...
考完试家长怎么查看孩子成绩和等级？

考完试家长怎么查看孩子成绩和等级？

考试结束了，孩子们的成绩和等级也出来了，对于家长来说，如何快速方便地查看孩子的成绩和等级呢？今天，我要向大家介绍一个非常实用的工具——易查分，让家长们便捷高效了解孩子的学习成果。 好消息&#xff01…
阅读更多...
如何封装自动化测试框架？（超详细~）

如何封装自动化测试框架？（超详细~）

封装自动化测试框架，测试人员不用关注框架的底层实现，根据指定的规则进行测试用例的创建、执行即可，这样就降低了自动化测试门槛，能解放出更多的人力去做更深入的测试工作。 本篇文章就来介绍下，如何封装自动化测试框…
阅读更多...
【已解决】ognl.PropertyAccessor

【已解决】ognl.PropertyAccessor

在Spring boot2.x用TemplateEngine处理数据得时候&#xff0c;出现以下错误&#xff1a; 定位到代码行&#xff1a; 解决办法&#xff1a;修改thymeleaf的依赖&#xff1a; <!-- thymeleaf --><dependency><groupId>org.thymeleaf</groupId><…
阅读更多...
低功耗窗帘电机解决方案成功应用并通过 Matter 1.1 认证

低功耗窗帘电机解决方案成功应用并通过 Matter 1.1 认证

Nordic Semiconductor官方宣布与HooRii Tech&#xff08;和众科技&#xff09;携手合作&#xff0c;基于 Nordic nRF52840 芯片平台打造的 HRN71模组&#xff0c;成功赋能低功耗窗帘电机品牌发布Matter产品。低功耗窗帘电机获得 Matter 1.1 认证意味着它具有与其他 Matter 认证…
阅读更多...
C. Ntarsis‘ Set

C. Ntarsis‘ Set

Problem - C - Codeforces 思路&#xff1a;这个题求一个满足条件的最小的&#xff0c;我们可以想到二分可以求满足条件的最小值&#xff0c;我们考虑二分答案&#xff0c;当当前的枚举的为mid时&#xff0c;我们考虑它会怎样变化&#xff0c;首先一开始mid的排名就是mid&#…
阅读更多...
Python，Bytetrack 源码解读，参数，源码解释，逐句分析代码，目标追踪

Python，Bytetrack 源码解读，参数，源码解释，逐句分析代码，目标追踪

文章目录 1、得到索引2、高得分框参与匹配&#xff0c;可能会留下有匹配不了的框3、低得分框参与匹配4、处理 unconfirmed 匹配5、创建新的【STrack对象】6、扔掉太久没匹配到框的【STrack对象】7、输出追踪框 1、得到索引 self.args.track_thresh是轨迹阈值。轨迹的得分是iou…
阅读更多...
RabbitMQ: topic 结构

RabbitMQ: topic 结构

生产者 package com.qf.mq2302.topic;import com.qf.mq2302.utils.MQUtils; import com.rabbitmq.client.Channel; import com.rabbitmq.client.Connection;public class Pubisher {public static final String EXCHANGE_NAME"mypubilisher";public static void ma…
阅读更多...
c++通过tensorRT调用模型进行推理

c++通过tensorRT调用模型进行推理

模型来源&#xff1a; 算法工程师训练得到的onnx模型 c对模型的转换&#xff1a; 拿到onnx模型后&#xff0c;通过tensorRT将onnx模型转换为对应的engine模型&#xff0c;注意&#xff1a;训练用的tensorRT版本和c调用的tensorRT版本必须一致。 如何转换&#xff1a; 算法工…
阅读更多...
Json“牵手”亚马逊商品详情数据方法，亚马逊商品详情API接口，亚马逊API申请指南

Json“牵手”亚马逊商品详情数据方法，亚马逊商品详情API接口，亚马逊API申请指南

亚马逊平台是美国最大的一家网络电子商务公司&#xff0c;亚马逊公司是1995年成立&#xff0c;刚开始只做网上书籍售卖业务&#xff0c;后来扩展到了其他产品。现在已经是全世界商品品种最多的网上零售商和第二互联网公司&#xff0c;亚马逊是北美洲、欧洲等地区的主流购物平台…
阅读更多...
为什么5G 要分离 CU 和DU？（4G分离RRU 和BBU）

为什么5G 要分离 CU 和DU？（4G分离RRU 和BBU）

在 Blog 一文中&#xff0c;5G--BBU RRU 如何演化到 CU DU&#xff1f;_5g rru_qq_38480311的博客-CSDN博客 解释了4G的RRU BBU 以及 5G CU DU AAU&#xff0c;主要是讲了它们分别是什么。但是没有讲清楚 为什么&#xff0c;此篇主要回答why。 4G 为什么分离基站为 RRU 和 BBU…
阅读更多...
什么是原生IP？原生IP与住宅IP有何区别？

什么是原生IP？原生IP与住宅IP有何区别？

相信许多做跨境的都会接触到IP代理&#xff0c;比如电商平台、社媒平台、收款平台等等&#xff0c;都会检测IP。那也会经常听到一些词汇&#xff1a;原生IP、住宅IP&#xff0c;这两者之间有什么区别呢&#xff1f;什么业务需要用到呢&#xff1f;接下来带大家具体了解一下。 什…
阅读更多...
React Antd可编辑单元格，非官网写法，不使用可编辑行和form验证

React Antd可编辑单元格，非官网写法，不使用可编辑行和form验证

antd3以上的写法乍一看还挺复杂&#xff0c;自己写了个精简版 没用EditableRowCell的结构&#xff0c;也不使用Context、高阶组件等&#xff0c;不使用form验证 最终效果&#xff1a; class EditableCell extends React.Component {state {editing: false};toggleEdit () &…
阅读更多...
SFUD固件移植

SFUD固件移植

SFUD作用 SFUD 是一款开源的串行 SPI Flash 通用驱动库。由于现有市面的串行 Flash 种类居多&#xff0c;各个 Flash 的规格及命令存在差异&#xff0c; SFUD 就是为了解决这些 Flash 的差异现状而设计&#xff0c;让我们的产品能够支持不同品牌及规格的 Flash&#xff0c;提高…
阅读更多...
Android 修改代码后不生效问题的终极方案

Android 修改代码后不生效问题的终极方案

前言&#xff1a; 最近遇到几个项目&#xff0c;都出现了代码修改后&#xff0c;直接点studio上的run&#xff0c;跑起来后代码没生效&#xff0c;如果重新clean rebuild可以生效&#xff0c;但是这太浪费时间了。网上找了各种方案&#xff0c;前面几个项目&#xff0c;有的是可…
阅读更多...
手写Spring：第19章-JDBC功能整合

手写Spring：第19章-JDBC功能整合

文章目录 一、目标&#xff1a;JDBC功能整合二、设计&#xff1a;JDBC功能整合三、实现&#xff1a;JDBC功能整合3.1 工程结构3.2 整合JDBC功能核心类图3.3 数据源操作3.3.1 数据源操作抽象类3.3.2 JDBC 工具类 3.4 数据库执行3.4.1 语句处理器接口3.4.2 结果处理器接口3.4.3 行…
阅读更多...
嵌入式Linux驱动开发（LCD屏幕专题）（四）

嵌入式Linux驱动开发（LCD屏幕专题）（四）

单Buffer的缺点与改进方法 1. 单Buffer的缺点 如果APP速度很慢&#xff0c;可以看到它在LCD上缓慢绘制图案 即使APP速度很高&#xff0c;LCD控制器不断从Framebuffer中读取数据来显示&#xff0c;而APP不断把数据写入Framebuffer 假设APP想把LCD显示为整屏幕的蓝色、红色 很…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023