资源描述：
《电力系统分析潮流计算的计算机算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、潮流计算的计算机算法实验报告姓名：学号：班级：一、实验目的掌握潮流计算的计算机算法。熟悉MATLAB，并掌握MATLAB程序的基本调试方法。二、实验准备根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法)的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。三、实验要求每人一组，在实验课时内，用MATLAB调试和修改运行程序，用算例计算输出潮流结果。四、实验程序clear;%清空内存n=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=inp
2、ut('请输入平衡节点号:isb=');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、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;fori=1:n%对PQ节点的处理ifi~=isb&B2(i,6)==2h=h+1;
4、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))*(real(
5、Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend%三：对PV节点的处理，注意这时不可再将h初始化为0fori=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
6、)));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;%四：创建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,
7、用于存储有功功率、无功功率和电压幅值的不平衡量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))-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(
