Error message: Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 2.078506e-18.

Question

Ashi Khajotia 2023-5-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1973589-error-message-warning-matrix-is-close-to-singular-or-badly-scaled-results-may-be-inaccurate-rcon

回答： Abhijeet 2023-5-30

Hello, I am trying to run this code, the values for D_g are going infinity at some points. Is there anything I can do?

th = 0:10:90;
lam = 3:0.01:7;
lp = 0.168; lc = 8.932;
ep = 1 - ((lam.^2.*lc)./(lp.^2*(lc - 1i.*lam)));
er = real(ep);
eim =  imag(ep);
ng = sqrt(0.5*(sqrt((er.^2) + (eim.^2)) + er));
kg = sqrt(0.5*(sqrt((er.^2) + (eim.^2)) - er));
nkg = ng + 1i*kg;
dg = 0.01;
for z = 1:length(th)
    for j = 1:numel(lam)
        phig(z,j) = 2*pi.*(dg./lam(j)).*sqrt((ng(j) - 1i*kg(j) ).^2 - sind(th(z).^2)); 
        kgx(z,j) = (2.*pi./lam(j)).*sqrt((ng(j) - 1i*kg(j) ).^2 - sind(th(z).^2));
        
        Dg = [kgx(z,j)/nkg(j) -kgx(z,j)/nkg(j); nkg(j) -nkg(j)];
        D_g = inv(Dg)
    end
end

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

请先登录，再进行评论。

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

Answer 1

Animesh 2023-5-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1973589-error-message-warning-matrix-is-close-to-singular-or-badly-scaled-results-may-be-inaccurate-rcon#answer_1246749

Hello! The code is likely giving infinite values for D_g because the Dg matrix is singular for certain values of ng and kg. In other words, the determinant of Dg is zero for these values, meaning that the inverse of Dg does not exist and D_g cannot be computed.

To avoid this issue, you can check whether the determinant of Dg is zero before attempting to compute the inverse. If the determinant is zero, you can skip the calculation and set D_g to NaN or some other appropriate value to indicate that it is undefined.

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

请先登录，再进行评论。

Answer 2

Abhijeet 2023-5-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1973589-error-message-warning-matrix-is-close-to-singular-or-badly-scaled-results-may-be-inaccurate-rcon#answer_1246869

在 MATLAB Online 中打开

Hi,

The values of D_g is going to infinity at some points because of the fact that the determinant of the matrix `Dg` is equalt to or close to zero.

You can fix this by adding a conditional check before calculating `D_g`, to check if the determinant of the Dg matrix is zero/close to zero. If it is zero, assign `D_g` to NaN or some other suitable value that indicates/fixes the error.

Here is the modified code for your reference:

th = 0:10:90;
lam = 3:0.01:7;
lp = 0.168; lc = 8.932;
ep = 1 - ((lam.^2.*lc)./(lp.^2*(lc - 1i.*lam)));
er = real(ep);
eim =  imag(ep);
ng = sqrt(0.5*(sqrt((er.^2) + (eim.^2)) + er));
kg = sqrt(0.5*(sqrt((er.^2) + (eim.^2)) - er));
nkg = ng + 1i*kg;
dg = 0.01;
for z = 1:length(th)
    for j = 1:numel(lam)
        phig(z,j) = 2*pi.*(dg./lam(j)).*sqrt((ng(j) - 1i*kg(j) ).^2 - sind(th(z).^2)); 
        kgx(z,j) = (2.*pi./lam(j)).*sqrt((ng(j) - 1i*kg(j) ).^2 - sind(th(z).^2));
        
        Dg = [kgx(z,j)/nkg(j) -kgx(z,j)/nkg(j); nkg(j) -nkg(j)];
        if abs(det(Dg)) < 1e-6  % check for singularity
           % handle or skip the calculation for this value of j
           continue;
        else
            D_g = inv(Dg);
        end
    end
end

Alternatively, You can try adding a small constant to the diagonal of the matrix `Dg`. You can achieve this by adding a small value `epsilon` to the diagonal matrix.

This can be achieved by modifying the code like:

Dg = [kgx(z,j)/nkg(j) -kgx(z,j)/nkg(j); nkg(j) -nkg(j)];
epsilon = 1e-6;  % small constant value
Dg = Dg + epsilon*eye(size(Dg));  % add small value to diagonal
D_g = inv(Dg);

Thanks !!