记录一下第一个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)如图: