《《电力系统分析潮流计算的计算机算法》》由会员分享，可在线阅读，更多相关《《电力系统分析潮流计算的计算机算法》（9页珍藏版）》请在金锄头文库上搜索。
1、.、三、四、潮流计算的计算机算法实验报告姓名：学号：班级：实验目的掌握潮流计算的计算机算法。熟悉MATLAB，并掌握MATLAB程序的基本调试方法。实验准备根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法（或PQ分解法）的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。实验要求每人一组，在实验课时内，用MATLAB调试和修改运行程序，用算例计算输出潮流结果。实验程序clear;%活空内存n=input（请输入节点数:n=）;n1=input（请输入支路数:n1=）;isb=input（请输入平衡节点号：isb=
2、）;pr=input（请输入误差精度:pr=）;B1=input（请输入支路参数:B1=）;B2=input（请输入节点参数:B2=）;X=input（节点号和对地参数:X=）;Y=zeros(n);Times=1;%一：创建节点导纳矩阵fori=1:n1ifB1(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)
3、;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)A2*B1(i,3);endendY;%将OrgS、DetaS初始化OrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%二：创建OrgS,用于存储初始功率参数h=0;j=0;fori=1:n%对PQ节点的处理ifi=isb&B2(i,6)=2h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(r
4、eal(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*Imag(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);endendend%三：对PV节点的处理，注意这时不可再将h初始化为0
5、fori=1:nifi=isb&B2(i,6)=3h=h+1;forj=1:nOrgS(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)*imag(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)*(rea
6、l(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);endendendOrgS;%四：创建PVU用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;fori=1:nifB2(i,6)=3t=t+1;PVU(t,1)=B2(i,3);endendPVU;%五：创建DetaS用于存储有功功率、无功功率和电压幅值的不平衡量h=0;fori=1:n%对PQ节点的处理ifi=isb&B2(i,6)=2h=h+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2
7、)-OrgS(2*h,1);endendt=0;fori=1:n%六：对PV节点的处理，注意这时不可再将h初始化为0ifi=isb&B2(i,6)=3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1)A2+imag(PVU(t,1)A2-real(B2(i,3)A2-imag(B2(i,3)A2;endendDetaS;%七：创建I,用于存储节点电流参数i=zeros(n-1,1);h=0;fori=1:nifi=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(
8、2*h,1)*sqrt(-1)/conj(B2(i,3);endendI;%八：创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;fori=1:n%对PQ节点的处理ifB2(i,6)=2h=h+1;forj=1:nifj=isbk=k+1;ifi=j%对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+re
9、al(I(h,1);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1);else%非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*
10、h,2*k)=Jacbi(2*h-1,2*k-1);endifk=(n-1)%将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendk=0;fori=1:n%对PV节点的处理ifB2(i,6)=3h=h+1;forj=1:nifj=isbk=k+1;ifi=j%对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i
11、,3)+real(I(h,1);Jacbi(2*h,2*k-1)=2*imag(B2(i,3);Jacbi(2*h,2*k)=2*real(B2(i,3);else%非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endifk=(n-1)%将用于内循环的指针置于初始值以确保雅可比矩阵换
12、行k=0;endendendendendJacbi;%九：求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU;%修正节点电压j=0;fori=1:n%对PQ节点处理ifB2(i,6)=2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfori=1:n%对PV节点的处理ifB2(i,6)=3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2;%十：开
13、始循环*whileabs(max(DetaU)prOrgS=zeros(2*n-2,1);%!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0；fori=1:nifi=isb&B2(i,6)=2h=h+1;forj=1:nOrgS(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)*imag(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
14、(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);endendendfori=1:nifi=isb&B2(i,6)=3h=h+1;forj=1:nOrgS(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)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j
15、)*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);endendendOrgS;%创建DetaSh=0;fori=1:nifi=isb&B2(i,6)=2h=h+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1);endendt=
