昆明理工大学计算机辅助分析课程设计及matlab自动实现潮流调节与计算

1、重要声明：1、本文只是提供一定参考，有错勿怪。2.、 本文采用的潮流调节及计算只针对本实验， 其中有很多固化的数据与参数，以及特殊的公式、请根据要求的不同进行一定的调整 （潮流计算的程序适用任何情况、 不同的只是系统参数） 。3、本文潮流调节为自动调节生成，调节变量有发电机2、 3 电压，三台变压器分接头，没有进行无功补偿，请根据具体情况加入或删减变量。4、自动调节潮流程序中，控制变量为各点电压首次迭代值、最终各点电压。最终选择标准为网损最低的一组数据。5、本实验报告为草版，有错误，但懒得改。6、鉴于个人所有文档免费得原则，本文为免费文档，请尽量下载。7、最后！都免费了，要求不要太高。要求及原

2、始数据电力系统潮流计算是电力系统分析中的一种最基本的计算， 是对复杂电力系统正常和故障条件下稳态运行状态的计算。 潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷. 各点电压是否满足要求， 功率的分布和分配是否合理以及功率损耗等。 对现有电力系统的运行和扩建， 对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。 潮流计算结果可用如电力系统稳态研究， 安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。 产生短路故障的主要原因是电力设备绝

3、缘损坏。 短路故障分为三相短路、 两相短路、 单相接地短路及两相接地短路。 其中三相短路时三相电流仍然对称， 其余三类短路统成为不对称短路。 短路故障大多数发生在架空输电线路。 电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。 短路计算可以为设备的选择提供原始数据。一、课程设计的目的与任务:该课程设计为电气工程及其自动化专业的主要专业课设计之一， 是学习电力系统分析基础和电力系统计算机辅助分析课程后的一个重要的实践性教学环节。 其目的在于巩固和加深对电力系统潮流和短路电流计算基本原理的理解， 学习和掌握应用计算机进行电力系统设计和计算的方法，培养学生独立分析和解决问题的能力。二、

4、课程设计的基本要求:掌握电力系统等值模型和参数计算， 以及潮流和短路计算的基本原理， 学会应用计算机计算系统潮流分布和短路电流的方法。三、课程设计选题原则 :该课程设计是根据电力系统分析基础和电力系统计算机辅助分析课程内容， 结合 实际工程和科研的电力系统网络进行系统的潮流和短路电流计算。四、课程设计内容:1、课程设计地点：机房和教室2、设计内容（ 1）电力系统潮流计算（2）电力系统短路电流计算具体内容包括：1）系统等值模型和参数计算。2）利用现有的潮流和短路电流计算程序进行计算，或编制电力系统潮流和短路计算程序。3）计算结果分析。4）编写说明书。3、时间安排1、系统等值模型和参数计算（2 天

5、）2、读懂现有的算法程序或视需要编制算法程序并进行计算（4 天）3、计算结果整理和分析（1 天）4、编写说明书（2 天）5、成绩考核（1 天）五、课程设计主要参考资料：1、李光琦.电力系统暂态分析.北京：水利电力出版社.第二版2、陈琦.电力系统稳态分析.北京：水利电力出版社.第二版3、孟祥萍.基于matlab的电力系统设计与计算.长春：吉林人民出版社4、孟祥萍.电力系统分析.北京：高等教育出版社5、电力系统潮流和短路计算课程设计指导书.自编六、课程设计成绩考核：根据学生的平时表现，说明书，以及设计完成后的口试成绩等进行综合评定。课程结束时，要求学生有课程设计说明书。课程设计成绩以平时成绩（占5

6、0%） 、设计说明书（占50%）进行综合评定。原始数据与思考题目标电网单线图busbus-6bt-is-jbus-1系统参数表1.线路参数表线路编号线路型号线路长度（km）线路电阻 q/km线路正序电抗q/kn）线路容纳之半s/km4-5lgj-240/301040.0470.41.78 x 10 64-6lgj-120/701100.0741.47x 10 65-7lgj-120/25980.0791.60 x 10 66-9lgj-95/551630.0921.80x 10 67-8lgj-240/30880.0471.78x 10 68-9lgj-240/301220.0471.78x

7、10 6说明：线路零序电抗为正序电抗 3倍。表2.变压器参数表线路编号变压器型号变压器变比（kv）短路电压白分数（）2-7sspl-220000242±3x 2.5%/20103-9sspl-120000242±3x 2.5%/1561-4sspl-240000242±3x2.5%/17.511说明：变压器零序电抗与正序电抗相等，且均为a/y0接法。表3.发电机参数表1发电机额定功率mw额定电压kv额定功率因数120016.50.852180180.85310013.80.85表4.发电机参数表2发电机母线名xdq-, xdqtd0（s）xqq一'xqq一

8、'tq0 (s)tj (s)110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.80331.31250.18138.591.25780.25000.6006.02表5.负荷数据表节点号有功负荷（mw）无功负荷（mva）51255069030810035思考题1、牛顿拉夫逊法与pq分解法有哪些联系？有哪些区别？二者的计算性能如何？2、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情 况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分 的？哪些可调？哪些不

9、可调？3、发电机节点的注入无功为负值说明了什么？4、线路有功潮流最有可能的流向是？线路无功潮流最有可能的流向是？5、负荷功率因数对系统潮流有什么影响？6、潮流控制的主要手段有哪些？7、如何降低系统网损？8、如何校验主设备的短路容量？9、利用节点阻抗矩阵计算短路电流和利用导纳矩阵计算短路电流的算法有何异同？10、减小短路电流的措施有哪些?设计主要内容：1、完成调试潮流计算与短路计算程序调试；计算程序以 电力系统计算机辅助分析 所用教材中的例程为基础， 要求掌握算法原理，建议进行适当改写以利于本次设计应用。2、根据设计资料准备计算数据；整理原始数据，使其符合计算程序的输入条件。要求比较采用变压器n

10、型模型前后两种情况的差异。3、为指定电网安排合理的稳态运行方式；合理设置目标电网中的控制变量， 以 1 号机为平衡机， 将潮流分布调整为可行的运行方式（要求所有节点的电压不得低于 1.0 p.u. ，也不能高于1.05p.u. ）绘制潮流分布图。 分析该运行方式的特点， 提出进一步降低网损的改进的措施， 并利用计算结果简单证明之。具体内容参考设计资料所附思考题。4、列表分析节点7 发生各种短路故障时的短路电流，以及运行方式对故障电流的影响；5、完成设计说明书。要求结构清晰，即从电网模型建立，到初态潮流分布计算，到问题分析， 再到调整后的方案计算。 故障电流计算也类似。 不必逐一回答所有思考题。

11、参数计算、潮流调节、与潮流计算matlab 源程序% 参数计算的clc;clear;l=104,110,98,163,88,122;r=0.047,0.074,0.079,0.092,0.047,0.047;b=1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006;disp(线路电抗x=');x1=100/220a2*0.4*l; % 线路各段电抗(均为标幺值，取 ub=220kv,sb=100mva , 下同)disp(x1);for i=1:6r(i)=100/220a2*l(i)*r(i);%加分号避免出现过程项，

12、下同end%线路各段电阻disp('线路电阻 r='); disp(r); %输出 rfor i=1:6z(i)=r(i)+j*x1(i);z0(i)=r(i)+3*j*x1(i);enddisp(线路正序、负序阻抗z1=z2='); disp(z);for i=1:6bl(i)=220a2/100*b(i)*l(i);b(i)=0+j*bl(i);end%线路各容纳之半disp(线路容纳之半b/2='); disp(b);%输出bdisp('线路零序阻抗z0:');disp(z0);u0=11,10,6;sn=240,220,120;un=24

13、2,242,242;for i=1:3;xt(i)=j*u0(i)*242a2/sn(i)/220a2;xt0(i)=xt(i);end;disp('变压器正、负序电抗值xt： ');disp(xt);disp('变压器零序电抗值 xt0： ');disp(xt0);x=0.0608,0.1198,0.1813;%发电机暂态电抗有名值ugn=16.5,18,13.8;for i=1:3k(i)=242/ugn(i);endfor i=1:3xdg(i)=j*100/220a2*x(i)*k(i)*k(i);end%发电机暂态电抗disp('发电机暂态电抗

14、xdg='); disp(xdg);%输出xdgpl=125,90,100;ql=50,30,35;for i=1:3zl(i)=100/(pl(i)a2+ql(i)a2)*(pl(i)+j*ql(i);end% 负荷节点阻抗disp('负荷节点阻抗zl='); disp(zl); %输出zl%支数参数形成的矩阵b1、支数参数形成的矩阵b2、节点号及其对地阻抗形成 的矩阵 x 的形成b1=1,4,xt(1),0,1,0;2,7,xt(2),0,1,0;3,9,xt(3),0,1,0;4,5,r(1)+x1(1)*j,b(1),1,0;4,6, r(2)+ x1(2)*j

15、,b(2),1,0;5,7,r(3)+x1(3)*j,b(3),1,0;6,9,r(4)+x1(4)*j,b(4),1,0;7,8,r(5)+x1(5)*j, b(5),1,0;8,9,r(6)+x1(6)*j,b(6),1,0;disp('支数参数形成的矩阵：b1='); disp(b1);%输出支数参数形成的矩阵b1p=200,180,100;w=(p/100)+(p/0.85)*sqrt(1-0.85*0.85)/100)*j;b2=w(1),0,1.05,1,0,1;w(2),0,1,1,0,3;w(3),0,1,1,0,3;0,0,1,0,0,2;0,1.25+0.5

16、i,1,0,0,2;0, 0.9+0.3i,1,0,0,2;0,0,1,0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;disp('支数参数形成的矩阵：b2='); disp(b2);%输出支数参数形成的矩阵b2% 潮流自动调节的与自动计算的n=9;nl=9;isb=1;pr=0.0001;mm=0;x=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;disp('*disp('只考虑各点电压水平，得到的满足要求的数据如下：)；for aa=0:5b2(2,4)=1+0.01*aa;for ab=0:5b2(3,4

17、)=1+0.01*ab;for ba=0:5b1(1,5)=0.95+0.025*ba;for bb=0:5b1(2,5)=0.95+0.025*bb;for bc=0:5b1(3,5)=0.95+0.025*bc;y=zeros(n);e=zeros(1,n);f=zeros(1,n);v=zeros(1,n);o=zeros(1,n);s1=zeros(nl); for i=1:nif x(i,2) =0;p=x(i,1);y(p,p)=1./x(i,2);endendfor i=1:nlif b1(i,6)=0p=b1(i,1); q=b1(i,2);else p=b1(i,2); q=

18、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,5)a2)+b1(i,4)./2;y(p,p)=y(p,p)+1./b1(i,3)+b1(i,4)./2;endg=real(y);b=imag(y);for i=1:ne(i)=real(b2(i,3);f(i)=imag(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

19、);ict1=0;it2=1;n0=2*n;n=n0+1;a=0;while it2 =0it2=0;a=a+1;for i=1:nif i=isbc(i)=0;d(i)=0;for j1=1:nc(i)=c(i)+g(i,j1)*e(j1)-b(i,j1)*f(j1);d(i)=d(i)+g(i,j1)*f(j1)+b(i,j1)*e(j1);endp1=c(i)*e(i)+f(i)*d(i);q1=f(i)*c(i)-d(i)*e(i);v2=e(i)a2+f(i)a2;if b2(i,6)=3dp=p(i)-p1;dq=q(i)-q1;for j1=1:nif j1=isb &

20、j1=ix1=-g(i,j1)*e(i)-b(i,j1)*f(i);x2=b(i,j1)*e(i)-g(i,j1)*f(i);x3=x2;x4=-x1;p=2*i-1;q=2*j1-1;j(p,q)=x3;j(p,n)=dq;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x4;j(m,q)=x2;elseif j1=i & j1=isbx1=-c(i)-g(i,i)*e(i)-b(i,i)*f(i);x2=-d(i)+b(i,i)*e(i)-g(i,i)*f(i);x3=d(i)+b(i,i)*e(i)-g(i,i)*f(i);x4=-c(i)+g(i,

21、i)*e(i)+b(i,i)*f(i);p=2*i-1;q=2*j1-1;j(p,q)=x3;j(p,n)=dq;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x4;j(m,q)=x2;endendelsedp=p(i)-p1;dv=v(i)a2-v2;for j1=1:nif j1=isb&j1=ix1=-g(i,j1)*e(i)-b(i,j1)*f(i);x2=b(i,j1)*e(i)-g(i,j1)*f(i);x5=0;x6=0;p=2*i-1;q=2*j1-1;j(p,q)=x5;j(p,n)=dv;m=p+1;j(m,q)=x1;j(m,n)=

22、dp;q=q+1;j(p,q)=x6;j(m,q)=x2;elseif j1=i&j1=isbx1=-c(i)-g(i,i)*e(i)-b(i,i)*f(i);x2=-d(i)+b(i,i)*e(i)-g(i,i)*f(i);x5=-2*e(i);x6=-2*f(i);p=2*i-1;q=2*j1-1;j(p,q)=x5;j(p,n)=dv;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x6;j(m,q)=x2;endendendendendfor k=3:n0k1=k+1;n1=n;for k2=k1:n1j(k,k2)=j(k,k2)./j(k,k)

23、;endj(k,k)=1;if k=3k4=k-1;for k3=3:k4for k2=k1:n1j(k3,k2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endif k=n0,break;endfor k3=k1:n0for k2=k1:n1j(k3,k2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endelsefor k3=k1:n0for k2=k1:n1j(k3,k2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endendendfor k=3:2:n0-1l=(k+1)./2;e

24、(l)=e(l)-j(k,n);k1=k+1;f(l)=f(l)-j(k1,n);endfor k=3:n0det=abs(j(k,n);if det>=prit2=it2+1;endendict2(a)=it2;ict1=ict1+1;for k=1:ndy(k尸sqrt(e(k)a2+f(k)a2);endfor i=1:ndy(ict1,i)=dy(i);endfor vb=1:nif dy(1,vb)<=1.05&dy(1,vb)>=1.00;zy(vb)=1;else zy(vb)=0;endendend %用高斯消去法解"w=-j*v"

25、for k=1:nv(k)=sqrt(e(k)a2+f(k)a2);o(k)=atan(f(k)./e(k)*180./pi;ende=e+f*j;for va=1:nif v(va)<=1.05&v(va)>=1.000;yj(va)=1;else yj(va)=0;endendif yj(1)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)=1if zy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)=1;dispc发电机2出口电压：');disp(b2(2,4

26、);dispc发电机3出口电压：');disp(b2(3,4);disp('变压器1接头电压：');disp(b1(1,5);disp('变压器2接头电压：');disp(b1(2,5);disp('变压器3接头电压：');disp(b1(3,5);disp('各节点的电压大小v为(节点号从小到大排列)：)；disp(v);for p=1:nc(p)=0;for q=1:nc(p)=c(p)+conj(y(p,q)*conj(e(q);ends(p)=e(p)*c(p);endfor i=1:nlif b1(i,6)=0p=b1(

27、i,1);q=b1(i,2);else p=b1(i,2);q=b1(i,1);endsi(p,q)=e(p)*(conj(e(p)*conj(b1(i,4)./2)+(conj(e(p)*b1(i,5)-conj(e(q)*conj(1./( b1(i,3)*b1(i,5);endfor i=1:nlif b1(i,6)=0p=b1(i,1);q=b1(i,2);else p=b1(i,2);q=b1(i,1);endsj(q,p)=e(q)*(conj(e(q)*conj(b1(i,4)./2)+(conj(e(q)./b1(i,5)-conj(e(p)*conj(1./( b1(i,3)

28、*b1(i,5);endzds=0;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);zds=ds(i)+zds;endmm=mm+1;zdss(mm尸sqrt(real(zds)a2+imag(zds)a2);disp('*');for i=1:nzb224(mm)=b2(2,4);zb234(mm)=b2(3,4);zb115(mm)=b1(1,5);zb125(mm)=b1(2,5);zb135(mm)=b1(3,5);end end e

29、nd end endend end enddisp('*');disp('满足各点电压要求的数据组共有：);disp(mm);lol=min(zdss);for dota=1:mmif zdss(dota)=lol mm=dota;end end');disp('当加入线路损耗时，得到其中线路损耗最小的一组组号及关键数据为:disp(mm);disp。发电机2出口电压：');disp(zb224(mm);disp。发电机3出口电压：');disp(zb234(mm);disp('变压器1接头电压：');disp(zb115

30、(mm);disp('变压器2接头电压：');disp(zb125(mm);disp('变压器3接头电压：');disp(zb115(mm);b2(2,4)=zb224(mm);b2(3,4)=zb234(mm);b1(1,5)=zb115(mm);b1(2,5)=zb125(mm);b1(3,5)=zb135(mm);disp('用此组数据算的的潮流结果为：');disp('节点数：');disp(n);disp('节点支路数：');disp(nl);disp('平衡节点：');disp(isb)

31、;disp(精度：)；disp(pr);disp('b1=');disp(b1);disp('b2=');disp(b2);disp('x=');disp(x);y=zeros(n);e=zeros(1,n);f=zeros(1,n);v=zeros(1,n);o=zeros(1,n);s1=zeros(nl);disp('节点数：)；disp(n);disp('节点支路数：');disp(nl);disp('平衡节点：');disp(isb);disp(精度：)；disp(pr);mm=0;x=1,0;2

32、,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;for i=1:nif x(i,2) =0;p=x(i,1);y(p,p)=1./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,5)a2)+b1(i,4)./2;y(p,p)=y(p,p)+1./b1(i,3)+b1(i,4)./2;endg=real(

33、y);b=imag(y);for i=1:ne(i)=real(b2(i,3);f(i)=imag(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;n0=2*n;n=n0+1;a=0;while it2 -=0it2=0;a=a+1;for i=1:nif i=isbc(i)=0;d(i)=0;for j1=1:nc(i)=c(i)+g(i,j1)*e(j1)-b(i,j1)*f(j1);d(i)=d(i)+g(i,j1

34、)*f(j1)+b(i,j1)*e(j1);endp1=c(i)*e(i)+f(i)*d(i);q1=f(i)*c(i)-d(i)*e(i);v2=e(i)a2+f(i)a2;if b2(i,6)=3dp=p(i)-p1;dq=q(i)-q1;for j1=1:nif j1=isb & j1=ix1=-g(i,j1)*e(i)-b(i,j1)*f(i);x2=b(i,j1)*e(i)-g(i,j1)*f(i);x3=x2;x4=-x1;p=2*i-1;q=2*j1-1;j(p,q)=x3;j(p,n)=dq;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=

35、x4;j(m,q)=x2; elseif j1=i & j1=isbx1=-c(i)-g(i,i)*e(i)-b(i,i)*f(i);x2=-d(i)+b(i,i)*e(i)-g(i,i)*f(i);x3=d(i)+b(i,i)*e(i)-g(i,i)*f(i);x4=-c(i)+g(i,i)*e(i)+b(i,i)*f(i);p=2*i-1;q=2*j1-1;j(p,q)=x3;j(p,n)=dq;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x4;j(m,q)=x2; end end elsedp=p(i)-p1;dv=v(i)a2-v2;for j

36、1=1:nif j1=isb&j1=ix1=-g(i,j1)*e(i)-b(i,j1)*f(i);x2=b(i,j1)*e(i)-g(i,j1)*f(i);x5=0;x6=0;p=2*i-1;q=2*j1-1;j(p,q)=x5;j(p,n)=dv;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x6;j(m,q)=x2; elseif j1=i&j1=isbx1=-c(i)-g(i,i)*e(i)-b(i,i)*f(i);x2=-d(i)+b(i,i)*e(i)-g(i,i)*f(i);x5=-2*e(i);x6=-2*f(i);p=2*i-1;

37、q=2*j1-1;j(p,q)=x5;j(p,n)=dv;m=p+1;j(m,q)=x1;j(m,n)=dp;q=q+1;j(p,q)=x6;j(m,q)=x2; endendendendendfor k=3:n0k1=k+1;n1=n;for k2=k1:n1j(k,k2)=j(k,k2)./j(k,k);endj(k,k)=1;if k=3k4=k-1;for k3=3:k4for k2=k1:n1j(k3,k2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endif k=n0,break;endfor k3=k1:n0for k2=k1:n1j(k3,k

38、2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endelsefor k3=k1:n0for k2=k1:n1j(k3,k2)=j(k3,k2)-j(k3,k)*j(k,k2);endj(k3,k)=0;endendendfor k=3:2:n0-1l=(k+1)./2;e(l)=e(l)-j(k,n);k1=k+1;f(l)=f(l)-j(k1,n);endfor k=3:n0det=abs(j(k,n);if det>=prit2=it2+1;endendict2(a)=it2;ict1=ict1+1;for k=1:ndy(k尸sqrt(e(k)a

39、2+f(k)a2);endfor i=1:ndy(ict1,i)=dy(i);endend %用高斯消去法解"w=-j*v"disp(迭代次数);disp(ict1);disp。没有达到精度要求的个数');disp(ict2);for k=1:nv(k)=sqrt(e(k)a2+f(k)a2);o(k)=atan(f(k)./e(k)*180./pi;ende=e+f*j;disp('各节点的实际电压标幺值 e为(节点号从小到大排列)：);disp(e);disp('各节点的电压大小v为(节点号从小到大排列)：');disp(v);disp(

40、'各节点的电压角。为(节点号从小到大排列)：)；disp(o);for p=1:nc(p)=0;for q=1:nc(p)=c(p)+conj(y(p,q)*conj(e(q);ends(p)=e(p)*c(p);enddisp('各节点白功率s为(节点号从小到大排列)：');disp(s);disp('各条支路的首段功率si为(顺序同您输入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);endsi(p,q)=e(p)*(conj(e(p)*conj(b1(i

41、,4)./2)+(conj(e(p)*b1(i,5)-conj(e(q)*conj(1./( b1(i,3)*b1(i,5);disp(si(p,q);enddisp('各条支路的末段功率sj为(顺序同您输入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);endsj(q,p)=e(q)*(conj(e(q)*conj(b1(i,4)./2)+(conj(e(q)./b1(i,5)-conj(e(p)*conj(1./( b1(i,3)*b1(i,5);disp(sj(q,p);end

42、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);disp(ds(i);endfor i=1:ict1cs(i)=i;enddisp('以下是各节点每次迭代后的电压值(如图所示) ');plot(cs,dy),xlabel('迭代次数'),ylabel('电压'),title('电压迭代次数曲线');最后是课程设计报告课程设

43、计课程名称：电力系统分析设计题目：基于matlab计算程序的电力系统运行分析学院：电力工程学院专业：电气工程自动化年 级：2009学生姓名易军（你可以称我易大神，可惜 dota、lol都是菜鸟顶级水平）指导教师：宋琪、谢静、单节杉日 期： 2013-1-13教务处制电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障 条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求， 功率的分布和分配是否 合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对

44、电力 系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究, 安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短 路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两 相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大 多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。目标电网接线图g1、电网模型的建立设计中，选取sb =100mva， u b

45、=220kv ，将所有支路的参数都折算到 220kv 电压等级侧，本文利用 matlab 计算参数其源程序为：clc;clear;l=104,110,98,163,88,122;r=0.047,0.074,0.079,0.092,0.047,0.047;b=1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006;disp('线路电抗x=');x1=100/220a2*0.4*l;%线路各段电抗(均为标幺值，取 ub=220kv,sb=100mv a)disp(x1);for i=1:6r(i)=100/220a2*

46、l(i)*r(i);enddisp('线路电阻r=');disp(r);for i=1:6z(i)=r(i)+j*x1(i);z0(i)=r(i)+3*j*x1(i);enddisp('线路正序、负序阻抗 z1=z2=');disp(z);for i=1:6bl(i)=220a2/100*b(i)*l(i);b(i)=0+j*bl(i);enddisp('线路容纳之半 b/2='); disp(b);disp('线路零序阻抗z0:');disp(z0);u0=11,10,6;sn=240,220,120;un=242,242,24

47、2;for i=1:3;xt(i)=j*u0(i)*242a2/sn(i)/220a2;xt0(i)=xt(i);end;disp('变压器正、负序电抗值 xt： ');disp(xt);disp('变压器零序电抗值xt0： ');disp(xt0);x=0.0608,0.1198,0.1813;%加分号避免出现过程项，下同% 线路各段电阻%输出 r%线路各容纳之半%输出 b%发电机暂态电抗有名值ugn=16.5,18,13.8;for i=1:3k(i)=242/ugn(i);endfor i=1:3xdg(i月*100/220a2*x(i)*k(i)*k(i

48、);end%发电机暂态电抗disp('发电机暂态电抗xdg='); disp(xdg);pl=125,90,100;ql=50,30,35;for i=1:3zl(i)=100/(pl(i)a2+ql(if2)*(pl(i)+j*ql(i);end%负荷节点阻抗disp(' 负荷节点阻抗zl='); disp(zl);0.0101 + 0.0860i0.0168 + 0.0909i0.0160 + 0.0810i0.0310 + 0.1347i0.0085 + 0.0727i计算结果为：线路电抗x=0.08600.09090.08100.13470.07270.

49、1008线路电阻r=0.01010.01680.01600.03100.00850.0118线路正序、负序阻抗z1=z2=columns 1 through 5column 60.0118 + 0.1008i线路容纳之半b/2=columns 1 through 50 + 0.0896i0 + 0.0783i0 + 0.0759i0 + 0.1420i0 + 0.0758icolumn 60 + 0.1051i线路零序阻抗 z0:columns 1 through 50.0101 + 0.2579i0.0168 + 0.2727i0.0160 + 0.2430i0.0310 + 0.4041i

50、0.0085 + 0.2182icolumn 60.0118 + 0.3025i变压器正、负序电抗值xt：0 + 0.0555i0.0550i0.0605i变压器零序电抗值 xt0 ：0 + 0.0555i0.0550i0.0605i发电机暂态电抗 xdg=0 + 0.0270i负荷节点阻抗 zl=0.0447i0.1152i0.8909 + 0.3118i0.6897 + 0.2759i1.0000 + 0.3333i支数参数形成的矩阵： b1=columns 1 through 51.00004.00000+0.0555i01.00002.00007.00000+0.0550i01.000

51、03.00009.00000+0.0605i01.00004.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.0810i0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000column 6000000000支数参数形成的矩

52、阵： b2=columns 1 through 52.0000 + 1.2395i01.05001.000001.8000+1.1155i01.00001.000001.0000+0.6197i01.00001.00000001.00000001.2500 +0.5000i1.00000000.9000 +0.3000i1.000000001.00000001.0000 +0.3500i1.000000001.000000column 61.00003.00003.00002.00002.00002.00002.00002.00002.0000二.系统等值电路图的绘制根据以上计算结果，得到系统等值电路图如下:ja