短路计算matlab程序.docxVIP

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Power System Analysis Project: Matlab Program %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear all close all clc a=exp(j*2*pi/3); A=[1 1 1 ; 1 a^2 a ; 1 a a^2]; Sref=900e6; Uref=300e3; Iref=Sref/(sqrt(3)*Uref); % data from the different documents: Positive; Negative; Zero; % run power flow: opt = mpoption(ENFORCE_Q_LIMS, 1); results=runpf(case9gs,opt); Vp=[results.bus(:,8).*(cos(pi/180*results.bus(:,9))+j*sin(pi/180*result s.bus(:,9)))]; clc; % Menu: Bucle=1; while (Bucle==1) % Menu: disp( ) disp(%%%%%% Program to calcula