%% 三门水电站三相短路程序
%% 电气一班 刘永康 0805010114
%% 短路数据读入部分(读入数据为自定应格式)
clear
clc
format long
fid=fopen('C:\Users\john\Desktop\三门水电站程序\三门水电站设计短路部分\短路读入数据.txt','r');
tline=fgetl(fid);
node1=[];node2=[];node3=[];R=[];X1=[]; X2=[]; %初始化,或者后面加个变量a,用node1(a),就不用初始定义了
while tline(1)~='L'
tline=fgetl(fid);
end
while tline(1)=='L'
node1=[node1 eval(tline(5))]; %平衡结点为1,PV结点为2,PQ结点为3
node2=[node2 eval(tline(13))];
R=[R eval(tline(19))];
X1=[X1 eval(tline(25:30))];
tline=fgetl(fid);
end
while tline(1)~='H'
tline=fgetl(fid);
end
while tline(1)=='H'
node3=[node3 eval(tline(5))];
X2=[X2 eval(tline(