octave - 如何绘制线性回归成本函数的等高线图？

Question

我正在关注 machine learning course在 Coursera 上，其中一个讲座提供了一个变量的线性回归成本函数的等值线图:

来源:https://www.coursera.org/learn/machine-learning/lecture/nwpe2/cost-function-intuition-ii

我认为从教育的角度来看，能够重现这张图表会很有用。我没有任何 Octave 音阶经验，因此我需要可以将其粘贴到 Octave 音阶命令窗口中的分步说明。

这里有人可以帮忙吗？

---

更新:

我得到了以下结果:

function cost = calc_cost (theta0, theta1)
   x = 1:10;
   y = x.*2;
   cost = arrayfun( @(t0, t1) ( 1/(length(x)) * sum( ((t0 + t1*x) - y).^2 )), theta0, theta1);
endfunction

[xx, yy] = meshgrid( -3000:50:3000, -3000:50:3000) ;
zz = calc_cost(xx, yy);
contour(xx, yy, zz )

Answer 1

如果您无法重写成本函数使其接受矩阵输入，您可以使用 arrayfun:

function cost = calc_cost (theta0, theta1)
   x = 1:10;
   y = x.*2;
   cost = ( 1/(length(x)) * sum( ((theta0 + theta1*x) - y).^2 ) );
 endfunction

[x,y] = meshgrid (linspace(-5000,5000,20), linspace(-500,500,20));
z = arrayfun (@calc_cost, x, y);
contour (x, y, z)
print out.png

给予