未定义与 'double' 类型的输入参数相对应的函数 'UnBalance_'。

[n,baseVolume,kmax,eps,BalanceNode,lineI,lineJ,lineR,lineX,lineB,groundI,groundB,transI,transJ,transR,transX,transk,operation,P,Q,PVNode,U,angU] = readDate(filename); %打开读取文件
tic
[Y] = formY(n,lineI,lineJ,lineR,lineX,lineB,groundI,groundB,transI,transJ,transR,transX,transk); %生成节点导纳矩阵
format long
for k=0:kmax
[deltpq,P,Q,U,angU] = UnBalance_(BalanceNode,Y,P,Q,U,angU); %计算不平衡量
pre=max(abs(deltpq)); %计算最大不平衡量，修改精度
if pre <=eps

disp('潮流收敛');
break;
else
[J]= Jac_(BalanceNode,n,angU,U,P,Q); %形成雅可比矩阵
y = Jac(-deltpq); %计算电压幅值和相角的修正量
U=U+y(n+1:2*n);
angU=angU+y(1:n); %重新计算电压幅值和电压相角

if k>=Kmax
disp('潮流计算不收敛'); %超过最大迭代次数则输出潮流不收敛
break;
end
end
end
toc;

楼上的回答看不懂.....