谐波潮流程序

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=