MATLAB: Resolve inf result

Question

Я пытаюсь вычислить следующее. Проблема в том, что я получаю inf для третьего элемента, который вычисляется. Как я могу разрешить это, чтобы больше не было inf-value?

num_labels = 2; 
htheta = [0.2 0.5; 0.8 0.1; 0.6 1; 0.7 0.4] 
y =  [1; 0; 0; 1;] 
m = 4; 


wes = ones(m, num_labels); 
J = 0; 


for k = 1:num_labels 
     for r = 1:m 
     if (htheta(r, k) < 0.5) && y(r, 1) ~= 0 
      wes(r, k) = 2; 
     else if (htheta(r, k) > 0.5) && y(r, 1) ~= 1 
       wes(r, k) = 4; 
      else wes(r, k) = 1; 
      end 
     end 
    end 
    yk = y == k; 
    hthetak = htheta(:, k); 
    Jk1 = (-yk .* log(hthetak) - (1 - yk) .* log(1 - hthetak)); 
    Jk2 = (-yk .* log(hthetak) - (1 - yk) .* log(1 - hthetak)) .* wes(:, k); 
    Jk3 = 1/m * sum((-yk .* log(hthetak) - (1 - yk) .* log(1 - hthetak)) .* wes(:, k)); 
    J = J + Jk3; 
end 

Jk1 
Jk2 
Jk3 
J 

Answer 1

htheta равно 1 в одном из своих пар индексов (3,2).

Jk3 рассчитывается с использованием значения log(1 - hthetak), что не определено для hthetak == 1. В то же время, yk имеет нулевое значение, так как k=2 и y никогда не равен 2.

Это делает продукт -Inf, и вы вычитая это, следовательно, вы получите Inf.

Эта страница MathWorks о Infinity and NaN может быть полезна.