Octave：梯度下降 线性回归编程

记录一下第一个Octave的代码。

传入的参数分别为自变量x，因变量（真实值）y，参数(随便给一个初值就行，但要注意向量的大小要和x，y匹配)

function J=function3(x, y, T)  %代价函数

J=1/(2*size(T,1))*sum((x*T-y).^2)
j=0;
m=0;
t1=0;
while 1,
	if(t1>=1000000),
		break;
	end

	for i=1:1:size(T,1),
		T(i)=T(i)-0.01*(1/size(T,1))*sum(x(:,i)'*(x*T-y));  %'
	end

	m=1/(2*size(T,1))*sum((x*T-y).^2);
	J=[J;j];                            %用来记录每次下降后的代价函数值
	if(abs(m-j)<=0.000000000000001),
		break;
	end

	j=m;

	t1=t1+1;
end

收敛效果(代价函数J)如图：