《电力系统中的潮流计算程序》由会员分享，可在线阅读，更多相关《电力系统中的潮流计算程序(7页珍藏版)》请在人人文库网上搜索。

1、B2=input(请输入各节点参数形成的矩阵：B2=);%本程序的功能是用牛顿-拉夫逊法进行潮流计算n=input(请输入节点数:n=);nl=input(请输入支路数：nl=);isb=input(请输入平衡母线节电号:isb=);pr=input(请输入误差精度:pr=);B1=input(请输入由支路参数形成的矩阵：B仁);变压器侧为1，否则为0B2=input(请输入各节点参数形成的矩阵:B2=);X=input(请输入由节点号及其对地阻抗形成的矩阵:X=);X=i nput(请输入由节点号及其对地阻抗形成的矩阵:X=);Y=zeros( n);U=zeros(1, n);cta=ze。

2、ros(1, n);V=zeros(1, n);O=zeros(1, n);S1=zeros (n l);for i=1: nif X(i,2)=0;p=X(i,1);Y(p,p)=X(i,2);endendfor i=1: nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5F2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,。

3、4)./2;end %求导纳矩阵G=real(Y);B=imag(Y);for i=1: ncta(i)=a ngle(B2(i,3);U(i)=abs(B2(i,3);%V(i)=B2(i,4);endfor i=1: nS(i)=B2(i,1)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=0;IT2=1;while IT2=0IT2=0;t1=1;t2=1;for i=1: nif i=isbC(i)=O;D(i)=O; for j1=1: nC(i)=C(i)+U(i)*U(j1)*(G(i,j1)*cos(cta(。

4、i)-cta(j1)+B(i,j1)*si n(cta(i)-cta(j1);D(i)=D(i)+U(i)*U(j1)*(G(i,j1)*si n(cta(i)-cta(j1)-B(i,j1)*cos(cta(i)-cta(j1);endDP(t1)=P(i)-C(i);t1=t1+1;if B2(i,6)=2DQ(t2)=Q(i)-D(i);t2=t2+1;endendendt1=t1-1;t2=t2-1;DPQ=DP;DQ; % 求 DP,DQfor i=1:t1+t2if abs(DPQ(i)prIT2=IT2+1;endendH=zeros(t1,t1);N=zeros(t1,t2);。

5、K=zeros(t2,t1);L=zeros(t2,t2);for i=1:t1for j1= 1:t1if j1=isb&j1=iH(i,j1)=0-U(i)*U(j1)*(G(i,j1)*si n(cta(i)-cta(j1)-B(i,j1)*cos(cta(i)-cta(j1);elseif j1=isb&j1=iH(i,j1)=U(i)A2*B(i,j1)+D(i);endendendfor i=1:t1for j1= 1:t2if j1=isb&j1=iN(i,j1)=0-U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1)+B(i,j1)*si n(cta(。

6、i)-cta(j1);elseif j1=isb&j1=iN(i,j1)=0-U(i)A2*G(i,j1)-C(i);endendendfor i=1:t2for j1= 1:t1if j1=isb&j1=iK(i,j1)=U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1)+B(i,j1)*si n(cta(i)-cta(j1); elseif j1=isb&j1=iK(i,j1)=U(i)A2*G(i,j1)-C(i);endendendfor i=1:t2for j1= 1:t2if j1=isb&j1=iL(i,j1)=0-U(i)*U(j1)*(G(i,j1)。

7、*si n(cta(i)-cta(j1)-B(i,j1)*cos(cta(i)-cta(j1);elseif j1=isb&j1=iL(i,j1)=U(i)A2*B(i,j1)-D(i);endendendJ=H,N;K,L;%求雅可比矩阵modify=-JDPQ;Dcta=modify(1:t1,:);t3=U(:,1:t2);DU=diag(t3,0)*modify(t1+1:t1+12,:);t4=1;for i=1:t1if B2(i,6 )=1cta(1,i)=cta(1,i)+Dcta(t4,1);t4=t4+1;endendt5=1;for i=1:t2if B2(i,6)=2U。

8、(1,i)=U(1,i)+DU(t5,1);t5=t5+1;endendICT1=ICT1+1;end %修正原值for i=1: nUU(i)=U(i)*cos(cta(i)+1i*U(i)*si n(cta(i);endfor p=1: nc(P)=0；for q=1: nc(p)=c(p)+conj(丫(p,q)*conj(UU(q); ends(p)=UU(p)*c(p);enddisp();disp(各节点电压U为(节点从小到大排列):);disp(UU);disp();disp(各节点电压相角为(节点从小到大排列):);disp(180*a ngle(UU)/pi);disp();。

9、disp(按公式计算全部线路功率，结果如下:);for i=1: nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(p,q)=UU(p)*(conj(UU(p)*conj(B1(i,4)./2)+(conj(UU(p)*B1(i,5)-conj(UU(q) *conj(1./(B1(i,3)*B1(i,5);% 各条支路首端功率 Sif=p,q,Si(p,q);disp(f);endfor i=1: nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1)。

10、;endSj(q,p)=UU(q)*(conj(UU(q)*conj(B1(i,4)./2)+(conj(UU(q)./B1(i,5)-conj(UU(p) *conj(1./(B1(i,3)*B1(i,5);% 各条支路末端功率 Sjf=q,p,Sj(q,p);disp(f);enddisp();disp(各条支路的功率损耗DS为(顺序同您输入B1时一样):);for i=1: nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);% 各条支路功率损耗 DS disp(DS(i)。

11、;endSp=0;for i=1: n Sp=Sp+UU(isb)*conj(Y(isb,i)*conj(UU(i);enddisp(平衡节点的功率:);disp(Sp);土豪，但是我们缺少贵族。富而不骄，莫若富而好礼。”如今我们不缺富不贵只能是土豪，你可以一夜暴富，但是贵气却需要三代以上的培养。孔子说高贵是大庇天下寒士俱欢颜的豪气与悲悯之怀，高贵是位卑未敢忘忧国的壮志与担当之志高贵是先天下之忧而忧的责任之心。精神的财富和高贵的内心最能养成性格的高贵，以贵为美，在不知不觉中营造出和气的氛围；以贵为高，在潜移默化中提升我们的素质。以贵为尊，在创造了大量物质财富的同时，精神 也提升一个境界。一个心灵高贵的人举手投足间都会透露出优雅的品质，一个道德高贵的社会大街小巷都会留露出和谐的温馨，一个气节高贵的民族一定是让人尊崇膜拜的民族。别让富而不贵成为永久的 痛。分享一段网上流传着改变内心的风水的方法，让我们的内心高贵起来：喜欢付出，福报就越来越多；喜欢感恩，顺利就越来越多；喜欢助人，贵人就越来越多；喜欢知足，快乐就越来越多；喜欢逃避，失败就越来越多；喜欢分享，朋友就越来越多。喜欢生气，疾病就越来越多；喜欢施财，富贵就越来越多；喜欢享福，痛苦就越来越多；喜欢学习，智慧就越来越多。