MATLAB牛拉法,MATLAB潮流程序(IEEE14直角坐标牛拉法).doc

MATLAB潮流程序(IEEE14直角坐标牛拉法)

MATLAB潮流程序(IEEE14 直角坐标牛拉法)

clear

baseMVA=100; %功率基值

%%读Data1中数据

load Data1.txt

Bus=Data1(:,1); %节点号

Vtype=Data1(:,5); %节点类型

Pload=Data1(:,8); %负载有功

Qload=Data1(:,9); %负载无功

Pgen=Data1(:,10); %发电机发出有功

Qgen=Data1(:,11); %发电机发出无功

Vset=Data1(:,13); %电压设定点

Qsh=Data1(:,17); %并联电容电纳标幺值

%%读Data2中数据

load Data2.txt

II=Data2(:,1);

JJ=Data2(:,2); %支路端点号

Ltype=Data2(:,5); %线路类型

R=Data2(:,6); %两点间电阻

X=Data2(:,7); %两点间电抗

B=Data2(:,8)/2; %线路对地电纳

K=Data2(:,14); %变压器非标准变压比

%%求导纳矩阵Y

y1=zeros(14);

y2=zeros(14);

y3=zeros(14);

lin=length(II); %支路数

for x=1:lin

switch Ltype(x)

case 1

y1(II(x),JJ(x))=1/(R(x)+i*X(x));

y1(JJ(x),II(x))=y1(II(x),JJ(x));

y3(II(x),JJ(x))=i*B(x);

y3(JJ(x),II(x))=i*B(x);

case 2

y1(II(x),JJ(x))=1/((R(x)+i*X(x))*K(x));

y1(JJ(x),II(x))=y1(II(x),JJ(x));

y2(II(x),JJ(x))=(1-K(x))/((R(x)+i*X(x))*K(x)^2);

y2(JJ(x),II(x))=(K(x)-1)/((R(x)+i*X(x))*K(x));

end

clear x

Y=zeros(14);

for x=1:14

Y(x,x)=sum(y1(x,:))+sum(y2(x,:))+sum(y3(x,:))+i*Qsh(x);

end

clear x;

Y=Y-y1;

G=real(Y);

B=imag(Y);

%%设电压初值

U=Vset;

e=real(U);

f=imag(U);

Ps=zeros(1,14);

Qs=zeros(1,14);

D=ones(26,1);

for x=1:14

Ps(x)=(Pgen(x)-Pload(x))/baseMVA;

Qs(x)=(Qgen(x)-Qload(x))/baseMVA;

end

clear x;

N=0;

Jacbi=zeros(26);

while max(abs(D))>0.000001

for x=2:14 %节点功率及电压不平衡量

switch Vtype(x)

case 1 %PQ节点

D(2*x-3)=Ps(x)-e(x)*(G(x,:)*e-B(x,:)*f)-f(x)*(G(x,:)*f+B(x,:)*e);

D(2*x-2)=Qs(x)-f(x)*(G(x,:)*e-B(x,:)*f)+e(x)*(G(x,:)*f+B(

MATLAB牛拉法,MATLAB潮流程序(IEEE14直角坐标牛拉法).doc相关推荐

matlab程序里nbus=x的意思,MATLAB潮流程序(IEEE14 直角坐标牛拉法)
MATLAB潮流程序(IEEE14 直角坐标牛拉法)clear baseMVA=100; %功率基值 %%读Data1中数据 load Data1.txt Bus=Data1(:,1); %节点号 V ...
小程序中scroll-view的下拉刷新和小程序页面的下拉刷新开启方法
scroll-view的下拉刷新 video.wxml中 <scroll-view class="videoScroll" scroll-y="true" ...
弦截法c语言程序,The Secant Method(正割法、弦截法) 附C语言代码
弦截法是一种求方程根的基该方法,在计算机编程中经常使用. 他的思路是这种:任取两个数x1.x2,求得相应的函数值f(x1).f(x2).假设两函数值同号,则又一次取数.直到这两个函数值异号为止. 连接 ...
MATLAB牛拉法计算潮流,matlab潮流计算
<matlab潮流计算>由会员分享,可在线阅读,更多相关<matlab潮流计算(14页珍藏版)>请在装配图网上搜索. 1.附录1使用牛顿拉夫逊法进行潮流计算的Matlab程序代 ...
用matlab计算潮流牛拉法,matlab潮流计算
<matlab潮流计算>由会员分享,可在线阅读,更多相关<matlab潮流计算(14页珍藏版)>请在人人文库网上搜索. 1.附录 1使用牛顿拉夫逊法进行潮流计算的Matlab ...
【极坐标下牛顿—拉夫逊潮流计算(matlab版+python版)】
程序名称## 极坐标下牛顿-拉夫逊潮流计算(matlab版+python版) 程序功能(对象) 适用于任意大小的纯交流电网,支持节点和支路的增删: 适用于接入多个风电.光伏等分布式电源: 子函数包含: ...
matlab潮流程序,IEEE33节点matlab潮流程序.doc
IEEE33节点matlab潮流程序 Bus [1,0,0 ; 2,100, 60; 3,90,40; 4,120,80; 5,60,30; 6,60,20; 7, 200, 100 ; 8,200, ...
matlab灰关联度,[转载]灰色关联和加权关联分析法matlab程序
%by allen %灰色关联分析法和灰色加权关联分析法matlab程序,其区别主要在于求关联度是使用平均值法还是加权法 %平均值法为一般普通求法,加权值法可以根据侧重点不同进行分析 %matlab6 ...
powerworld电力系统仿真，潮流计算，短路计算，电力系统分析。潮流计算对比，牛拉法，PQ分解法对比
powerworld电力系统仿真,潮流计算,短路计算,电力系统分析.潮流计算对比,牛拉法,PQ分解法对比编号:1710662437866344电气女博士

MATLAB牛拉法,MATLAB潮流程序(IEEE14直角坐标牛拉法).doc

MATLAB牛拉法,MATLAB潮流程序(IEEE14直角坐标牛拉法).doc相关推荐

最新文章

热门文章