clc
clear
orig=dlmread（‘fund4.txt‘）;
har=dlmread（‘harmonic4.txt‘）;
p_h=1;???????????????????%諧波源節點
a=-0.5+j*3^0.5/2;????????%基準項
n_node=orig（11）;
n_b=orig（12）;
n_t=orig（13）;
n_g=orig（14）;
n_l=orig（15）;
sb=100;vb=110;??????????????%基準標幺值
zb=vb*vb/sb;????????????????%阻抗基準值
%讀入數據
for?i=1:n_b????????%支路
????b（i）.i=orig（i+11）;
????b（i）.j=orig（i+12）;
????b（i）.r=orig（i+13）;
????b（i）.x=orig（i+14）;
????b（i）.b=orig（i+15）;
????b（i）.l=orig（i+16）;
end
for?i=1:n_t
????t（i）.i=orig（n_b+i+11）;
????t（i）.j=orig（n_b+i+12）;
????t（i）.r=orig（n_b+i+13）;
????t（i）.x=orig（n_b+i+14）;
????t（i）.k=orig（n_b+i+15）;
end
for?i=1:n_g???????????????%將平衡節點號改為最后一個（重要）！！！
????g（i）.i=orig（n_b+n_t+i+11）;
????g（i）.type=orig（n_b+n_t+i+12）;
????g（i）.p=orig（n_b+n_t+i+13）;
????g（i）.q=