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系统计算机上右键管理打不开的情况不知怎么解决,针对这个问题小编就 ...
最新文章
- 2014 Super Training #8 C An Easy Game --DP
- 洛谷P1194 买礼物
- 软工实践第三次作业-原型设计
- Winform中使用NPOI实现导出Excel并文件另存为
- W-D-S-UART编程
- linux 中文文件名不能下载不了,linux 64位系统 mod_encoding解决中文文件名不能下载问题...
- 老齐python-基础7(文件操作、迭代)
- leetcode python3 简单题168. Excel Sheet Column Title
- HDU 4389——X mod f(x)(数位DP)
- git使用puttygen生成公钥私钥
- 互联网常见的大数据分析指标
- python怎么编辑图片中的文字_Python在图片中添加文字的两种方法
- Mac操作系统-软件安装
- 阿里云技术大咖分享新内容新交互时代下的新技术、新机会
- 计算机考研就业率,这8个考研专业,就业率一直居高不下,报了就是赚了!
- ios中在app应用内刷新小组件数据
- 条码打印软件如何制作渐变色二维码
- 微信8.0,初心改不改?
- 老大威武,afterDelay,
- 【观察】明略数据吴明辉:AI不是一个技术,而是一个时代
热门文章
- Cpp / 通用引用、引用折叠与完美转发问题
- TCP/IP / 状态转换图
- 启明云端分享|ESP32-C3-MINI-1和ESP32-C3-MINI-1U的区别是什么?
- int类型存小数 mysql_MySQL面试题-数据类型
- 0xbc指令 st75256_DDOS终极加速列车算法
- mysql query结果集_如何解决PHP使用mysql_query查询超大结果集超内存问题
- proj4经纬度bl转换xy_多种坐标系之间的转换之Proj.NET_转载
- STM32L0开发——ADC多通道采集,IDE和IAR开发注意事项
- nginx与php-fpm通讯方式
- Python——你应该知道这些