基于前推回代法的电力系统潮流计算

16-6-12下午7:34MATLAB命令行窗口第1页，共3页function[z,s,v]=sback(z,s,v)%%%{函数功能：逆着功率传输方向计算功率分布输入：节点阻抗矩阵z，负荷矩阵初值s，电压初值v输出：节点阻抗矩阵z，负荷矩阵修正值s，电压初值v备注：s上三角为末端功率，下三角为首端功率%}%%n=length(z);%计算节点数forp=1:ni=n+1-p;forq=1:nj=n+1-q;if(ij&&z(i,j)~=0)s(i,j)=sum(s(:,j))-s(i,j);%计算ij支路末端功率s(j,i)=s(i,j)+...((real(s(i,j))).^2+(imag(s(i,j))).^2)./...((v(j)).^2).*z(i,j);%计算ij支路首段功率endendend%%function[z,s,v]=vforward(z,s,v)%%%{函数功能：顺着功率传输方向计算节点电压输入：节点阻抗矩阵z，负荷矩阵s，电压初值v输出：节点阻抗矩阵z，负荷矩阵s，电压修正值v备注：s上三角为末端功率，下三角为首端功率%}%%n=length(z);%计算节点数fori=1:nforj=1:nif(ij&&z(i,j)~=0)v(j)=sqrt((v(i)-((real(s(j,i)).*real(z(i,j))+...imag(s(j,i)).*imag(z(i,j)))./v(i))).^2+...((real(s(j,i)).*imag(z(i,j))-...imag(s(j,i)).*real(z(i,j)))./...v(i)).^2);%由始节点电压计算终节点电压功能函数sback.m：功能函数vforward.m基于前推回代法的电力系统潮流计算16-6-12下午7:34MATLAB命令行窗口第2页，共3页endendenda=z;%返回计算后的zb=s;%返回计算后的sc=v;%返回计算后的v%%%%%{程序功能：基于前推回代法的电力系统潮流计算%}%%%clear,clc;%清除工作区，清除命令窗口loaddata_experiment4.mat%加载原始数据z，s，vv=[110,110,110,110,110];%输入电压初值z1=z;%备份原始数据s1=s;%s上三角为末端功率，下三角为首端功率v1=v;n=3;%输入迭代次数nfori=1:n-1[z1,s1,v1]=sback(z1,s1,v1);[z1,s1,v1]=vforward(z1,s1,v1);end%完成第1至n-1次迭代[z2,s2,v2]=sback(z1,s1,v1);[z2,s2,v2]=vforward(z2,s2,v2);%完成第n次迭代delta_v=v2-v1;%计算v的误差s2%输出s2v2%输出v2delta_v%输出v的误差%%主函数main.m：16-6-12下午7:34MATLAB命令行窗口第3页，共3页s2=2.0000+2.0000i14.0294+14.0294i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i14.0621+14.0621i2.0000+2.0000i4.0000+4.0000i8.0054+8.0027i0.0000+0.0000i0.0000+0.0000i4.0027+4.0053i4.0000+4.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i8.0267+8.0241i0.0000+0.0000i4.0000+4.0000i4.0000+4.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i4.0054+4.0027i4.0000+4.0000iv2=110.0000109.7443109.6349109.4518109.3421delta_v=1.0e-09*0-0.2493-0.2690-0.5509-0.5940运行主函数输出结果：