电力系统潮流计算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 %求修正后的ｗ

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程序。相关推荐

edge robert matlab,哪位熟悉matlab的大神路过瞄一眼哈
cxfx(believe truth believe me) UID240430 帖子100 精华积分1755 蛋蛋币1755 枚威望0 BT积分0 阅读权限60 性别男在线时间125 小时注册 ...
php输出后以密文显示,加密解密,_求大神帮看一下，如何根据下面的程序，把明文生成密文谢谢，加密解密 - phpStudy...
求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 def parserLicense(self,encCode): data = base ...
用java编写圆锥_求java大神帮忙求大神帮助！Java
导航:网站首页 > 求java大神帮忙求大神帮助!Java 求java大神帮忙求大神帮助!Java 相关问题: 匿名网友: 普通类 public class CircularA { //求圆 ...
交流潮流matlab程序,大神们，求个电力系统潮流计算的matlab程序。
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 n=input('Please input n\n'); %n表示系统的节点数 d0=input('Please input d0\n'); %d0表示系 ...
matlab ktrlink,大神们，怎么设置滑动滑动条然后出来的图形也跟着变化？
本帖最后由我叫W儿于 2016-5-25 20:56 编辑大神们!我是初学者.在GUI设计中,我把传递函数里的变量以滑动条的形式表达了,但是,运行出图后,再改变滑动条中滑块位置,图不会跟着改变. ...
不同cpu matlab,求助大神，为何不同机器运行MATLAB结果不同
求助:不同机器运行MATLAB结果不同我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...
求解，某M1水卡数据计算分析/大神们求指导！
标题-某M1水卡数据计算分析洗澡水真的是好贵啊! 用Proxmark3解了一下水卡的数据,发现只有6扇区和7扇区数据有变化. 看了些资料发现 6扇区0块第1 2位表示金额.现在卡里有2位数的钱没 ...
抖音黑客帝国代码雨——程序大神如何攻占抖音？简直不可思议！
在普通人的脑补中,程序员可能是这样的不知道为什么,这张图总是能莫名触发技哥的笑点. 亦或是这样的每天无休止的与代码和Bug打着交道-- 等等,这好像不对啊.当初看黑客帝国的时候可不是这么心酸的,看 ...
evcdf matlab,求助大神
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...

电力系统潮流计算matlab程序,大神们，求个电力系统潮流计算的matlab程序。

电力系统潮流计算matlab程序,大神们，求个电力系统潮流计算的matlab程序。相关推荐

最新文章

热门文章