Matrix is singular to working precision

Question

Ahmad Karnama 2012-3-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/31090-matrix-is-singular-to-working-precision

I was using Matlab 2010 and my code was working fine! I changed my windows and start using Matlab 2011b and I am getting this warning for the same code:

Matrix is singular to working precision.

I have absolutely no idea what could be the problem. In fact, I am running a power flow in Matlab and it used to converge with the same code before but now the load flow can not converge either.

Would be glad if somebody can give a hint :)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2012-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/31090-matrix-is-singular-to-working-precision#answer_39551

Without seeing the command, which causes the problem, it is impossible to guess the reason of the problem.

Are you really sure that you are running the same code with the same data?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Ahmad Karnama 2012-3-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/31090-matrix-is-singular-to-working-precision#answer_39649

在 MATLAB Online 中打开

Thanks Jan for your answer. This is the part of the code which generates the warning:

    %%evaluate Jacobian
    [dSbus_dVm, dSbus_dVa] = dSbus_dV(Ybus, V);
    j11 = real(dSbus_dVa([pv; pq], [pv; pq]));
    j12 = real(dSbus_dVm([pv; pq], pq));
    j21 = imag(dSbus_dVa(pq, [pv; pq]));
    j22 = imag(dSbus_dVm(pq, pq));
    J = [   j11 j12;
            j21 j22;    ];
    %%compute update step
    dx = -(J \ F);