我编写了一个成本函数的代码,但它给出的答案是错误的。
我已经多次阅读了代码,但找不到错误所在。
这是我的代码:-
function J = computeCost(X, y, theta)m = length(y); % number of training exampless = 0;h = 0;sq = 0;J = 0;for i = 1:m h = theta' * X(i, :)'; sq = (h - y(i))^2; s = s + sq;end J = (1/2*m) * s;end
示例:-
computeCost( [1 2; 1 3; 1 4; 1 5], [7;6;5;4], [0.1;0.2] )ans = 11.9450
这里答案应该是11.9450,但我的代码给出的结果是:-
ans = 191.12
我已经检查了矩阵乘法,代码计算是正确的。
回答:
看起来你误解了运算符的计算顺序。实际上
1/2*m ~= 1/(2*m)
考虑到这一点,似乎你在计算平均值。通常情况下,与其重新发明轮子,不如使用内置函数来完成工作,这会得到一个更加清晰(且错误更少)的实现:
function J = computeCost(X, y, theta) h = X * theta; sq = (h - y).^2; J = 1/2 * mean(sq);endcomputeCost( [1,2;1,3;1,4;1,5], [7;6;5;4], [0.1;0.2] )% ans = 11.9450