matlab潮流计算求节点自导纳,大神们,求个电力系统潮流计算的matlab程序。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
n=input('Please input n\n'); %n表示系统的节点数
d0=input('Please input d0\n'); %d0表示系统的误差
m=input('Please inpit m\n'); %m表示系统的pq节点的个数
p0=ones(n-1,1);
q0=ones(m,1);
d=1;
dpv=ones(n-1,1);
dqv=ones(m,1);
dp=zeros(n-1,1);
dq=zeros(m,1);
for a=1:n
for b=1:n
G(a,b)=real(Y(a,b));
B(a,b)=imag(Y(a,b));
end
end %把Y矩阵的实部和虚部分别导出
for a=1:n-1
for b=1:n-1
B1(a,b)=B(a,b);
end
end %表示出B1
for a=1:m
for b=1:m
B2(a,b)=B(a,b);
end
end %表示出B2
while(d>d0)
for a=1:n-1
p0(a)=0;
for b=1:n
p0(a)=p0(a)+v(a)*v(b)*(G(a,b)*cos(w(a)-w(b))+B(a,b)*sin(w(a)-w(b)));
end
end %利用初值求得相应的功率
for a=1:m
q0(a)=0;
for b=1:n
q0(a)=q0(a)+v(a)*v(b)*(G(a,b)*sin(w(a)-w(b))-B(a,b)*cos(w(a)-w(b)));
end
end %利用初值求得相应的无功
for a=1:n-1
dp(a)=p(a)-p0(a);
end %求dp
for a=1:m
dq(a)=q(a)-q0(a);
end %求dq
for a=1:n-1
dpv(a)=dp(a)/v(a);
end %利用初值求出dp/v
for a=1:m
dqv(a)=dq(a)/v(a);
end %利用初值求出dq/v
dw=-1*B1\dpv;
for a=1:n-1
dw(a)=dw(a)/v(a);
end %求出dw的修正值
dv=-1*B2\dqv; %求出v的修正值
d=max(max(abs(dw)),max(abs(dv)));%误差衡量
for a=1:n-1
w(a)=w(a)+dw(a);
end %求修正后的w
for a=1:m
v(a)=v(a)+dv(a);
end %求修正后的v
end
这个是基于PQ分解法解潮流方程,我是按照自己的习惯编写的==。在运行本程序前,需要自行建立p(有功)、q(无功)、Y(节点导纳矩阵)、w(相角)、v(电压幅值)矩阵。
其中,节点编号请按照PQ节点、PV节点、平衡节点的次序编号。为了获得更好的收敛,可以将w统一设置为0,v统一设置为1。PQ解耦法适用于大电力系统,对于小型配电网等电力系统可能会存在较大的计算误差。
matlab潮流计算求节点自导纳,大神们,求个电力系统潮流计算的matlab程序。相关推荐
- 用c语言求元素移动的次数T(n),悬赏!求C语言大神帮忙,写一个数字推盘游戏的程序(N-puzzle) 具体的要求看图片...
满意答案 zcloj 2013.10.14 采纳率:56% 等级:12 已帮助:5505人 这是我在网上收藏在本地的数独游戏代码.看能否符合你的要求 #include int fun(int a ...
- matlab里exist是什么意思,通达信exist函数用法,exist()函数在matlab中怎么用?求教大神~...
Q1:exist()函数在matlab中怎么用?求教大神~ exist name 等价于 r=exist(name) ,在程序里面这样更加实用 0 不存在则返回值 1 name 可以是变量名,如果存在 ...
- matlab 五连珠,我是初学者,想求哪位大神帮帮我把这个小游戏程序简单写写注释,五连珠游戏...
我是初学者,想求哪位大神帮帮我把这个小游戏程序简单写写注释,五连珠游戏 五连珠游戏 游戏规则: 移动小球,五个或五个以上相同颜色的组成一行或一列就会消失 移动的时候只能移动到相连通的格子里去 func ...
- 用java编写圆锥_求java大神帮忙 求大神帮助!Java
导航:网站首页 > 求java大神帮忙 求大神帮助!Java 求java大神帮忙 求大神帮助!Java 相关问题: 匿名网友: 普通类 public class CircularA { //求圆 ...
- 饭店菜单选择c语言,求C语言大神,帮忙做一下这个题,要源代码设计实现一个餐厅点餐系统? 爱问知识人...
求C语言大神,帮忙做一下这个题,要源代码 设计实现一个餐厅点餐系统,实现电子点餐.要求实现二级菜单: 一级菜单: 请选择: 1.点餐 2.查看已点菜品 3.添菜 求C语言大神,帮忙做一下这个题,要源代 ...
- php输出后以密文显示,加密解密,_求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢,加密解密 - phpStudy...
求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 def parserLicense(self,encCode): data = base ...
- 求解,某M1水卡数据计算分析/大神们求指导!
标题-某M1水卡数据计算分析 洗澡水真的是好贵啊! 用Proxmark3解了一下水卡的数据,发现只有6扇区和7扇区数据有变化. 看了些资料发现 6扇区0块 第1 2位表示金额.现在卡里有2位数的钱 没 ...
- matlab计算npv,哪位大神看看这个程序呀,参数下面给了
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 K=0; N=input('请输入节点数:N='); NPQ =input('请输入节点数:NPQ='); NPV =input('请输入节点数:NPV= ...
- python加权求和_Python计算加权平均分的问题?大神求解答。。。。。
def getSubjectStats(subject, weights): ''' e.g. subject:[['Bart','geology'],[10,5,85]]--grades即是[10, ...
最新文章
- 技术的顶点,是普惠人类还是超越人类
- 201521123029《Java程序设计》第1周学习总结
- windows10风格 springboot activiti 整合项目框架源码 shiro 安全框架 druid 数据库连接池...
- 做一名真正的软件工程师
- 【驱动】使用结构体 file_operations封装驱动设备的操作 | 结构体初始化
- mysql+firewall_mysql - ERROR 1123(HY000):无法初始化函数'mysql_firewall'; 插件初始化功能失败 - 堆栈内存溢出...
- .NET Core微服务之基于EasyNetQ使用RabbitMQ消息队列
- tomcat7 java_java开发环境配置(windows下JDK7+tomcat7)
- 第六次作业阶段一 学生成绩录入系统设计与实现
- Android 通知栏Notification的整合 全面学习 (一个DEMO让你完全了解它)
- unity NGUI下载 支持unity高版本
- 《政府采购货物和服务招标投标管理办法》指明的方向
- 【转载】MIXLY二次开发(傻瓜版)
- CatBoost快速入门
- thinkphp 检测上传的图片中是否含有木马脚本
- the sdk seems invalid pycharm
- linux安装程序出现了问题,linux安装extundelete以及对遇到问题的解决
- ubuntu18.04安装roboware studio
- C语言ctime函数
- linux系统怎么两个盘对拷,Linux分区对拷等小技巧
热门文章
- 中小企业的管理,难题很可能出在这三个方面
- 美国付费电视行业2015年财报大集合
- Hello MySQL(三)——SQL语言
- OpenOfficeException:not obtain service: com.sun.star.ucb.FileContentProvider 解决
- Java使用POI写(创建)Excel
- FreeModbus V1.6 主机
- Ubuntu 22.04 LTS 中文桌面/服务器正式版发布 (内含 ISO 镜像下载)
- 服务器的信息失败怎么回事,发送小程序信息 本机正常 服务器就发送失败 这是怎么回事?...
- 计算机选购配置项目活动,BIM电脑配置专辑之如何选择一台项目专业级别笔记本电脑...
- 如何正确高效准确的使用搜索引擎?