潮流计算的计算机算法实验报告姓名： 学号： 班级： 一、实验目的掌握潮流计算的计算机算法。熟悉 MATLAB，并掌握 MATLAB 程序的基本调试方法.二、实验准备根据课程内容，熟悉 MATLAB 软件的使用方法,自行学习 MATLAB 程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法（或 PQ 分解法） 的计算程序，用相应的算例在 MATLAB 上进行计算、调试和验证。三、实验要求每人一组，在实验课时内，用 MATLAB 调试和修改运行程序，用算例计算输出潮流结果。四、实验程序clear; ％清空内存n=input（’请输入节点数：n=’)；n1=input(’请输入支路数:n1=’)；isb=input（’请输入平衡节点号：isb=’）；pr=input('请输入误差精度：pr=’);B1=input(’请输入支路参数:B1=')；B2=input（'请输入节点参数：B2=’）;X=input（’节点号和对地参数：X=')；Y=zeros（n）; Times=1； ％一：创建节点导纳矩阵for i=1：n1 if B1（i,6）==0 ％不含变压器的支路 p=B1(i，1)； q=B1（i,2）; Y（p，q)=Y（p,q)-1/B1(i,3); Y(q，p）=Y(p，q); Y(p,p）=Y(p，p)+1/B1（i,3)+0。5＊B1（i，4）； Y(q，q)=Y（q,q）+1/B1（i，3)+0.5＊B1（i,4）； else ％含有变压器的支路 p=B1（i,1)； q=B1(i，2); Y(p,q）=Y（p，q）—1/(B1（i，3)＊B1(i,5)）; Y(q，p）=Y(p,q); Y(p,p)=Y（p，p）+1/B1（i，3）； Y(q,q）=Y（q,q）+1/(B1(i，5）^2＊B1(i,3）); endendY； %将 OrgS、DetaS 初始化 OrgS=zeros（2*n-2，1）; DetaS=zeros（2*n—2,1); %二：创建 OrgS，用于存储初始功率参数h=0;j=0;for i=1：n ％对 PQ 节点的处理 if i~=isb&B2（i，6)==2 h=h+1; for j=1:n OrgS(2＊h—1,1)=OrgS(2*h-1，1)+real（B2（i,3)）＊（real(Y（i,j))*real(B2（j,3）)-imag（Y(i,j))*I mag(B2(j，3)）)+imag(B2（i，3))*(real（Y（i,j)）*imag（B2（j,3）)+imag(Y（i,j）)＊real（B2（j,3)））; OrgS(2*h,1）=OrgS(2*h,1）+imag（B2(i，3))＊(real(Y（i,j)）*real(B2(j，3）)—imag（Y(i,j)）＊imag（B2(j,3))）—real(B2（i,3))*(real（Y(i，j））＊imag(B2(j，3))+imag（Y（i,j）)＊real(B2(j,3)）)； end endend％三：对 PV 节点的处理，注意这时不可再将 h 初始化为 0for i=1：n if i~=isb＆B2（i，6)==3 h=h+1; for j=1：n OrgS（2*h-1,1)=OrgS(2*h-1,1)+real(B2（i，3)）*(real（Y(i,j））＊real(B2(j，3）)-imag（Y...
