基于MATLAB开发的电力系统潮流计算程序，采用极坐标计算方式

function [pnum,linenum,pnod,benchmarkP,Kmax,center,et,type,line,branch,linetran,PGm,Qlimit,ppram,lflimit,qnum]=read11(name) fid=fopen(name); pnum=fscanf(fid,'%g',1); linenum=fscanf(fid,'%g',1); pnod=fscanf(fid,'%g',1); benchmarkP=fscanf(fid,'%g',1); Kmax=fscanf(fid,'%g',1); center=fscanf(fid,'%g',1); et=fscanf(fid,'%g',1); type=fscanf(fid,'%g',1); i=1; while fscanf(fid,'%g',[1,1])>0 line(i,:)=fscanf(fid,'%g',[1,5]); i=i+1; end i=1; while 1 branch(i,1)=fscanf(fid,'%g',[1,1]); if branch(i,1)~=0 branch(i,2)=fscanf(fid,'%g',[1,1]); i=i+1; else branch(i,:)=[]; break; end end i=1; while fscanf(fid,'%g',[1,1])>0 linetran(i,:)=fscanf(fid,'%g',[1,7]); i=i+1; end i=1; while 1 PGm(i,1)=fscanf(fid,'%g',[1,1]); if PGm(i,1)~=0 PGm(i,2:5)=fscanf(fid,'%g',[1,4]); i=i+1; else PGm(i,:)=[]; break; end end i=1; while 1 Qlimit(i,1)=fscanf(fid,'%g',[1,1]); if Qlimit(i,1)~=0 Qlimit(i,2:4)=fscanf(fid,'%g',[1,3]); i=i+1; else Qlimit(i,:)=[]; break; end end i=1; while 1 ppram(i,1)=fscanf(fid,'%g',[1,1]); if ppram(i,1)~=0 ppram(i,2:6)=fscanf(fid,'%g',[1,5]); i=i+1; else ppram(i,:)=[]; break; end end i=1; while 1 if fscanf(fid,'%g',[1,1])>0 lflimit(i,:)=fscanf(fid,'%g',[1,4]); i=i+1; else break; end end PGm(:,2:5)=PGm(:,2:5)./benchmarkP; [qnum,mm]=size(Qlimit);