MATLAB潮流计算仿真设计

附录 A MATLAB 程序%本程序的功能是用牛顿——拉夫逊法进行潮流计算% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号clear;n=10;%input('请输入节点数：n=');nl=11;%input('请输入支路数：nl=');isb=1;%input('请输入平衡母线节点号：isb=');pr=0.00001;%input('请输入误差精度：pr=');B1=[1 2 1.755e-2+4.155e-2i 0.26i 1 0; 1 4 3.159e-2+7.479e-2i 0.1215i 1 0; 1 6 3.159e-2+7.479e-2i 0.1215i 1 0; 2 3 3.68e-3+0.11135i 0 0.909 1; 4 5 3.68e-3+0.11135i 0 0.909 1; 4 6 2.808e-2+6.648e-2i 0.108i 1 0; 6 7 3.0865e-3+0.0833i 0 0.909 1; 6 8 3.159e-2+7.479e-2i 0.1215i 1 0; 6 10 2.457e-2+5.817e-2i 0.0945i 1 0; 8 9 3.0865e-3+0.0833i 0 0.909 1; 8 10 2.808e-2+6.648e-2i 0.108i 1 0];%input('请输入由支路参数形成的矩阵： B1=');B2=[0 0 1.05 1.05 0 1; 0 0 1 0 0 2; 0 0.6+0.3718i 1 0 0 2; 0 0 1 0 0 2; 0 0.4+0.247i 1 0 0 2; 0 0 1 0 0 2; 0 0.35+0.2169i 1 0 0 2; 0 0 1 0 0 2; 0 0.5+0.3099i 1 0 0 2; 0.8 0 1.05 1.05 0 3];%input('请输入各节点参数形成的矩阵： B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%－－－－－－－修改部分－－－－－－－－－－－－ym=0;SB=100;UB=220;%ym=input('您输入的参数是标么值？（若不是则输入一个不为零的数值）');if ym~=0 %SB=input('请输入功率基准值:SB='); %UB=input('请输入电压基准值:UB='); YB=SB./UB./UB; BB1=B1; BB2=B2; for i=1:nl B1(i,3)=B1(i,3)*YB; B1(i,4)=B1(i,4)./YB; end disp('B1矩阵B1='); disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB; B2(i,2)=B2(i,2)./SB; B2(i,3)=B2(i,3)./UB; B2(i,4)=B2(i,4)./UB; B2(i,5)=B2(i,5)./SB; end disp('B2矩阵B2='); disp(B2)end% % %---------------------------------------------------for i=1:nl %支路数 if B1(i,6)==0 %左节点处于低压侧 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B...