1、Y（m,m）=Y（m,m）+1/（B（a,3）+j*B（a,4）;Y（ n,n ）=Y（ n,n ）+1/（B（a,3）+j*B（a,4）/（B（a,5）A2）;Y（m,n）=-1/（B（a,3）+j*B（a,4）/B（a,5）;Y（n,m）=Y（m,n）;elseif B（a,5）Y（m,m）=Y（m,m）+1/（B（a,3）+j*B（a,4）-j*B（a,5）;Y（n,n）=Y（n,n）+1/（B（a,3）+j*B（a,4）-j*B（a,5）;Y（m,n）=-1/（B（a,3）+j*B（a,4）;end%=Y（2,2）=Y（2,2）+j*0.01;Z=inv（Y）;If=1/Z（3,3）;
2、 %接金属短路， Zf=0for i=1:5V（1,i）=1-Z（i,3）/Z（3,3）;%= %计算各支路电流for a=1:0 I（m,n）=（B（a,5）*V（1,m）-V（1,n）/（B（a,3）+j*B（a,4）;I（m,n）=（V（1,m）-V（1,n）/（B（a,3）+j*B（a,4）;%=ans.txtw fprintf（fid,Y_matrixn%= %输出导纳阵for j=1:k=Y（i,j）;re=real（k）;fprintf（fid,%8.4f,re）;im=imag（k）;if im=0+%8.4fjtn%=%= %输出阻抗阵，导纳阵求逆Z_matrixn5 k=Z
3、（i,j）; re=real（k）; im=imag（k）;0 fprintf（fid,=0 fprintf（fid,end fprintf（fid,%=%= %打印输出短路电流 fprintf（fid,If= re=real（If）; im=imag（If）; end fprintf（fid,nVn%=%= %打印输出节点电压5 fprintf（fid,Note%d:V%d=,i,i）; k=V（1,i）;%输出打印支路电流fprin tf（fid,lbrn4k=l（i,j）；if（re=0|im=0）Bra nch%d-%d:l%d%d=,i,j,i,j）;%8.4f,re）;%8.4fj
4、n+%8.4fjn % 关闭文件 附：1-20. 00. 105L 05230. 0240.0650. 016q0. 0180. 050. 0130.030. 080. 02_q0- 00. 1840. S02、程序输出结果 节点导纳阵：0. DDfla -9. 5238j1 00004$.07D1jQ. ODOtH-n oooojo.oooc*o.aocoj3.C00C4O.OOOOj0. ODCO? D7O3j9. 10B5-33. D90lj-4. 9O9D+13.5288jy mgs+L0.958ja.tooc+D.MOOj0. OOCQ4-0. 1000-4. 99B0+13.53B
5、8J11.3729-31.T.3F3 舛.tcoc+O.OOOftj0. QOCLH-n. DOriQj7. 109C+in.EKRCj-a. 3?3SF17 7C54j10.4835-34,528flju.tmc-FS.6612Jn. Qdooj1.00004n.tnnijjn. OQODfn nrnnjn.n. conn -51.节点阻抗阵:l_jh.it rixD. 0040-S.282SjCL 0042-8. BC70j-0. J028-0. 033d-s.asrej-0.0036-fl.1954：d 0042-S.8070JQL 0044-氏 3473j-0. CIQ29 -S.26
6、71J-0, 0030-工西紂-0.0038-9.鸽 BL;-0. C02B-6.S25j-0, 0023-9. 271j0.0067 -S,2436j-o. wio-p. 362（ j-0,0010-9.6430j-0. 0034-E,827Sj-a oo3s-9r !90j-0. a010 -0.252OJ0. OOGD2431j0.00-0. 0036-0,1064j-a 0033-0. B651j-a loio -o+ei30jO.OJ61-S.62S6j0.0054-9+ S4呂巧短路电流：If= 0.0001+ 0.1082j节点电压：四、实验总结这是我的第一次上机实验，感觉稍微有
7、点难，主要还是在工具软件一 C语言或者MATLAB的运用上，但是我相信，以后学习中，我会努力掌握的，这是我把理论 应用于实际中的必要桥梁！一、 实验目的根据所给的电力系统， 编制牛顿法潮流计算程序， 通过计算机进行调试， 最后完成一个切实可行的电力系统计算应用程序。 通过自己设计电力系统计算程序加深对电力系统分析的 理解，同时加强计算机实际应用能力。二、 实验内容电力系统潮流计算的计算程序设计及编制和调试。三、 程序框图四、实验过程1、实验程序及说明 clear%打开输入数据%读 8 个数%得4X 4的0阵node4.txt,8）; C=fscanf（fid,3,（A（1,1）-1）;C=CB
8、（2,5）=1/B（2,5）;%=%与前一实验同法求导纳阵 for a=1:A（2,1） m=B（a,1）; n=B（a,2）; if B（a,5）0 Y（m,m）=Y（m,m）+1/（B（a,3）+j*B（a,4）; Y（n,m）=Y（m,n）;else if B（a,5）Q（2*i-1,1）=C（i,3）;Q（2*i,1）=0;elseQ（2*i-1,1）=1;Q（2*A（3,1）-1,1）=A（4,1）;Q（2*A（3,1）,1）=0;end fid=fopen（answer.txt= 节点电压 V=n 迭代计数 tt V1=e1+jf1tttt V2=e2+jf2ttttV3=e3+j
9、f3ttn%= 求 W 阵= W=zeros（2*（A（1,1）-1）,50）;for x=1:50 % 设置迭代次数为 50 次A（7,1）k=C（i,1）;p=0;q=0;m=0;n=0;A（1,1）g=real（Y（k,j）;b=imag（Y（k,j）;e=Q（2*j-1,1）;f=Q（2*j,1）;p=p+g*e-b*f;q=q+g*f+b*e;m=m+g*e-b*f;n=n+g*f+b*e;W（2*k,x）=C（k,2）-Q（2*k-1,1）*p-Q（2*k）*q;W（2*k-1,x）=C（k,3）-Q（2*k）*m+Q（2*k-1,1）*n;End %=PQ 节点%=PV 节点fo
10、r l=1:A（8,1）k=C（l+A（7,1）,1）;A（1,1） g=real（Y（k,j）; b=imag（Y（k,j）; e=Q（2*j-1,1）; f=Q（2*j,1）; p=p+g*e-b*f; q=q+g*f+b*e;W（2*k-1,x）=（C（k,3）A2-（Q（2*k-1,1）A2+Q（2*k）A2）;End%=%比较是否符合条件Max=0;2*（A（1,1）-1）Max=max（abs（W（i,x）,Max）;if Max%8.4f+%8.4fjtt ,Q（2*k-1,1）,Q（2*k,1）;%8.4f%8.4fjttend %=%= 平衡点功率 =k=A（3,1）;v=0
11、;j=sqrt（-1）;for b=1:m=conj（Y（k,b）;p=Q（2*b-1,1）+Q（2*b,1）*j;n=conj（p）;v=v+m*n;end Wp=（Q（2*k-1,1）+j*Q（2*k,1）*v;%= 节点电压 V= %fid=fopen（ 节点电压 Vn node%d:%8.4f+%8.4fjn,k,Q（2*k-1,1）,Q（2*k,1）; else%8.4f%8.4fjn end平衡点功率 P+jQ=%8.4f+%8.4fj,real（Wp）,imag（Wp）;附注：n ode4.txt %输入数据4 4 4 1. 05 0.00001 50 2 12 0, 10.4-
12、0. 015281 00. 31. 14 0. L20. 5-0. 01S204 0. 080. 4-CL 01413-0. 3-0. 18-0, 55Q. 13CL 5L 12、程序输出结果迭代计数V2=e2+jV3=e3+jf30. 093& -0. 0088J0.9763 -0. 1078j1. 1000+ D. 12575-O.OOSfij0. 9590 -0. 1084j1.0924+ 0. 1289j0. 9846 -0. 0080 j0. 95S7 -0. 1084j1.09244- 0. lSgiOj节点电压node!: 0.9846 -0.0083 node2: 0.9587 -0. 1084jnode3： 1.0924+ 0, 12P0;i平衡点功率P+jQ= 0. 36 79+ 0.五、实验总结本次实验是我把应用理论知识的重要实践，经过实验，我有两点感想，首先，作为一名工科学生，应该能熟练运用 C语言和MATLAB等工具，其次，理论如果不用于实践，就永远不知道理论是用来干什么的，学到头一直是满脑子的浆糊，所以， 以我的切身经历建议，把这门实验放在跟课程平行的时间上进行， 这样不仅有利于实验开展，也有利于学生更加深刻地学习！
