IEEE 14 节点系统

function mloss=yue（tt） %chaoliu
mloss=0;
b1=tt（1）;%tt（1）-tt（5）发电机端电压
u2=tt（2）;
u3=tt（3）;
u4=tt（4）;
u5=tt（5）;
t1=tt（6）;%tt（7）-tt（9）变压器抽头
t2=tt（7）;
t3=tt（8）;
q9=tt（9）;%无功补偿容量
qc9=0.1*round（q9）;
ta1=1+0.0125*round（t1）;%变压器变比
ta2=1+0.0125*round（t2）;%变压器变比
ta3=1+0.0125*round（t3）;%变压器变比
%n=input（‘请输入节点数:n=‘）;
%n1=input（‘输入支路数:n1=‘）;
%isb=input（‘输入平衡母线接线号:isb=‘）;
%pr=input（‘输入误差精度:pr=‘）;
%B1=input（‘输入由支路参数形成的矩阵:B1=‘）;
%B2=input（‘输入由节点参数形成的矩阵:B2=‘）;
%X=input（‘输入由节点号及其对地阻抗形成的矩阵:X=‘）;
n=14;
n1=20;  
isb=1;
pr=0.00001;
B1=[1 2 0.01938+0.05917i 0.0264 1 0;1 5 0.05403+0.22304i 0.01700 1 0;2 3 0.04699+0.19797i 0.02190 1 0;2 4 0.05811+0.17632i 0.01870 1 0;2 5 0.05695+0.17388 0.01730 1 0;3 4 0.06701+0.17103i 0.00640 1 0;4 5 0.01335+0.04211i 0 1 0;5 6 0.25202i 0 ta3 0;4 7 0.20912i