我已经为 cost函数编写了代码,但给出的答案不正确。
我已经阅读了很多次代码,但是我找不到错误。
这是我的代码:-
function J = computeCost(X, y, theta)
m = length(y); % number of training examples
s = 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
我已经检查了矩阵乘法,并且代码正在正确计算。
答案 0 :(得分:2)
似乎您误解了操作员评估顺序。实际上
1/2*m ~= 1/(2*m)
考虑到这一点,看来您正在计算平均值。除了重新设计轮子之外,通常最好使用内置函数来完成这项工作,从而实现更清晰(且更不易出错)的实现:
function J = computeCost(X, y, theta)
h = X * theta;
sq = (h - y).^2;
J = 1/2 * mean(sq);
end
computeCost( [1,2;1,3;1,4;1,5], [7;6;5;4], [0.1;0.2] )
% ans = 11.9450