资源描述：
《MATLAB 极坐标求解潮流计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、电力系统分析大作业——潮流计算班级：XXX姓名：XXX学号：XXX-8-程序清单：%直角坐标法求解潮流计算clear;close;clc;Branch=[120.1+0.4j0.01528j10310.3jinf1.11140.12+0.50j0.01920j10240.08+0.40j0.01413j10]%Branch矩阵：1、支路首端号；2、支路末端号；3、支路阻抗；4、支路对地导纳；5、支路的变化；6、支路首端处于K侧为1,1侧为0Y=zeros(4);%节点导纳矩阵fori=1:4ifBranch(i,6)==0%不含变压器的支路j=Branch(i,1);k=Branc
2、h(i,2);Y(j,k)=Y(j,k)-1/Branch(i,3);Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3)+Branch(i,4);Y(k,k)=Y(k,k)+1/Branch(i,3)+Branch(i,4);elsej=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-Branch(i,5)/Branch(i,3);-8-Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3);Y(k,k)=Y(k,k)+Branch(i,5)^2/Branch(i,3);endenddisp
3、('节点导纳矩阵');YG=real(Y);B=imag(Y);V=[1;1;1.1;1.05];%给定V的初始计算值e=real(V);f=imag(V);disp('节点电压的实部：')edisp('节点电压的虚部：')fdisp('节点注入有功功率：')Ps=[-0.3;-0.55;0.5;0]disp('节点注入无功功率：')Qs=[-0.18;-0.13;0;0]%由各节点电压向量（状态变量）可得各节点注入功率：P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G*f+B*e);del_W=[-0.30-P(1);-0.18-Q
4、(1);-0.55-P(2);-0.13-Q(2);0.5-P(3);1.1^2-e(3)^2-f(3)^2];n=0;%辅助循环计数变量while(any(del_W>1e-5))&(n<5)n=n+1;disp('迭代次数：');disp(n)%--------------------------------------------------------%修正方程式：%[-0.30-P(1)][del_e(1)]%[-0.18-Q(1)][del_f(1)]%[-0.55-P(2)][del_e(2)]%[-0.13-Q(2)]+J*[del_f(2)]=0%[0.5-P(3
5、)][del_e(3)]%[1.1^2-e(3)^2-f(3)^2][del_f(3)]%雅克比矩阵为：J=[P1_e1P1_f1P1_e2P1_f2P1_e3P1_e3;%Q1_e1Q1_f1Q1_e2Q1_f2Q1_e3Q1_f3%P2_e1P2_f1P2_e2P2_f2P2_e3P2_f3%Q2_e1Q2_f1Q2_e2Q2_f2Q2_e3Q2_f3%P3_e1P3_f1P3_e2P3_f2P3_e3P3_f3%V3_e1V3_f1V3_e2V3_f2V3_e3V3_f3]%---------------------------------------------------
6、-----%-------------------求雅可比矩阵的参数-----------------------%-8-GeBf=G*e-B*f;%辅助计算函数GfBe=G*f+B*e;%辅助计算函数fori=1:2forj=1:3ifi==jP_e(i,j)=-GeBf(i)-G(i,j)*e(i)-B(i,j)*f(i);P_f(i,j)=-GfBe(i)-G(i,j)*f(i)+B(i,j)*e(i);Q_e(i,j)=GfBe(i)+B(i,j)*e(i)-G(i,j)*f(i);Q_f(i,j)=-GeBf(i)+G(i,j)*e(i)+B(i,j)*f(i);else
7、P_e(i,j)=-G(i,j)*e(i)-B(i,j)*f(i);Q_f(i,j)=-P_e(i,j);P_f(i,j)=B(i,j)*e(i)-G(i,j)*f(i);Q_e(i,j)=P_f(i,j);endendendforj=1:2P_e(3,j)=-G(3,j)*e(3)-B(3,j)*f(3);P_f(3,j)=B(3,j)*e(3)-G(3,j)*f(3);V_e(3,j)=0;V_f(3,j)=0;endP_e(3,3)=-GeBf(3)-G(3,3
