[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;
未定义与 'double' 类型的输入参数相对应的函数 'UnBalance_'。
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
1条回答 默认 最新
悬赏问题
- ¥15 素材场景中光线烘焙后灯光失效
- ¥15 请教一下各位,为什么我这个没有实现模拟点击
- ¥15 执行 virtuoso 命令后,界面没有,cadence 启动不起来
- ¥50 comfyui下连接animatediff节点生成视频质量非常差的原因
- ¥20 有关区间dp的问题求解
- ¥15 多电路系统共用电源的串扰问题
- ¥15 slam rangenet++配置
- ¥15 有没有研究水声通信方面的帮我改俩matlab代码
- ¥15 ubuntu子系统密码忘记
- ¥15 保护模式-系统加载-段寄存器