NF=1;
for i1=1:NF
clear;
n1=9;
n2=9;
n0=9;
nl1=15;
nl2=15;
nl0=15;
Lf=input('请输入短路类型(单相接地短路输入‘1’,两相短路输入‘2’,两相短路接地输入‘3’):Lf=');
f=7;
If=0;
zf=0;
zg=0;
B1= [0 1 0.0553i 0 1 0 ;
0 2 0.1228i 0 1 0;
0 3 0.1785i 0 1 0 ;
1 4 0.0476i 0 1.050 0;
0 5 0.6514-0.2413i 0 1 0;
0 6 0.9174-0.2752i 0 1 0;
2 7 0.0474i 0 1.025 0;
0 8 1.0492-0.4590i 0 1 0 ;
3 9 0.0484i 0 1.025 0;
4 5 0.0854i 0 1 1;
4 6 0.0907i 0 1 1;
5 7 0.1248i 0 1 1;
6 9 0.1255i 0 1 1;
7 8 0.0696i 0 1 1;
8 9 0.0922i 0 1 1];
B2=[0 1 0.0553i 0 1 0 ;
0 2 0.1228i 0 1 0;
0 3 0.1785i 0 1 0 ;
1 4 0.0476i 0 1.050 0;
0 5 0.6514-0.2413i 0 1 0;
0 6 0.9174-0.2752i 0 1 0;
2 7 0.0474i 0 1.025 0;
0 8 1.0492-0.4590i 0 1 0 ;
3 9 0.0484i 0 1.025 0;
4 5 0.0854i 0 1 1;
4 6 0.0907i 0 1 1;
5 7 0.1248i 0 1 1;
6 9 0.1255i 0 1 1;
7 8 0.0696i 0 1 1;