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程序。相关推荐
- 用java编写圆锥_求java大神帮忙 求大神帮助!Java
导航:网站首页 > 求java大神帮忙 求大神帮助!Java 求java大神帮忙 求大神帮助!Java 相关问题: 匿名网友: 普通类 public class CircularA { //求圆 ...
- php输出后以密文显示,加密解密,_求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢,加密解密 - phpStudy...
求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 def parserLicense(self,encCode): data = base ...
- 求解,某M1水卡数据计算分析/大神们求指导!
标题-某M1水卡数据计算分析 洗澡水真的是好贵啊! 用Proxmark3解了一下水卡的数据,发现只有6扇区和7扇区数据有变化. 看了些资料发现 6扇区0块 第1 2位表示金额.现在卡里有2位数的钱 没 ...
- matlab lsa,omlsa cohen大神的OM-LSA仿真代码
MATLAB soft 联合开发网 - pudn.com
omlsa 所属分类:其他 开发工具:matlab 文件大小:4KB 下载次数:22 上传日期:2019-01-16 10:01:04 上 传 者:wraII 说明: cohen大神的OM-LSA仿 ...
- matlab lognpdf用法,求助大神
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...
- matlab在电力系统潮流计算程序,大神们,求个电力系统潮流计算的matlab程序。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 n=input('Please input n\n'); %n表示系统的节点数 d0=input('Please input d0\n'); %d0表示系 ...
- 采用morison方程基于matlab计算大直径波浪力,用Morison方程计算分析悬浮隧道所受波浪力初探...
用Morison方程计算分析悬浮隧道所受波浪力初探 用Morison方程计算分析悬浮隧道所受波浪力初探 第16卷第3期石家庄铁道学院第卷第期后豕饮迥r子阮子很 2003年9月JOURNALOFSHIJ ...
- matlab 代码开头,诸位大神请帮我讲解一下每个程序开头这些行的意思
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function varargout = Image_Main(varargin) % IMAGE_MAIN MATLAB code for Image_ ...
- 计算机演示题打不开,大神为你演示win7系统计算机上右键管理打不开的还原技巧...
当我们经常使用的电脑工作的时候,时间久了难免会遇到win7系统计算机上右键管理打不开的问题,如今就有用户反映在使用电脑的时候遇到win7系统计算机上右键管理打不开的情况不知怎么解决,针对这个问题小编就 ...
最新文章
- Ubuntu root 密码 sudo passwd
- JavaScript中变量提升是什么?如何实现?
- 关于网页的一些小知识点
- 电脑pin码忘了登录不进系统_想要大屏打王者,家里电视用不了?投屏到电脑试试...
- 华科计算机考研2022年分数线,2022年华中科技大学软件工程考研分数线、参考书、上岸前辈初复试经验...
- ST-Link刷成J-Link
- leetcode59题:滑动窗口的最大值(用双端队列实现)
- 一般纳税人税额计算_一般纳税人企业缴税计算
- burpsuite基本设置与重复提交比较(一)
- 使用PPT编写一个序号抽奖
- 2018年微信数据报告
- java8中Function函数
- 以图搜图在线网站汇总(共7个)
- Ubuntu 16.10 解决 DNS 解析失败
- 一些有用的英语学习资料
- 【Spring】Spring IOC学习笔记
- 欧几里得算法扩展欧几里得
- 【重磅整理】提前看287篇ICLR-2021 深度强化学习领域论文得分汇总列表
- css 动画 呼吸 呼吸灯 效果
- 2017计算机考证结果