CS229作业一第二题

之前没学过Matlab，所以第一题的代码直接看的答案，看完之后熟悉了Matlab的写作习惯，第二题顺利自己搞定。

题目地址：http://cs229.stanford.edu/materials.html

Matlab代码：

x=load('q2x.dat');
y=load('q2y.dat');
x=[ones(size(x,1),1) x];
theta=pinv(x'*x)*x'*y;
figure; hold on;
m=size(x,1);
plot(x(:,2),y,'b.');
x1=min(x(:,2)):.01:max(x(:,2));
y1=theta(1)+theta(2)*x1;
plot(x1,y1);
tau=[ 0.1 0.3 0.8 2 10];
colors=['r' 'g' 'm' 'y' 'k'];
for i=1:5;
    for j=1:size(x1,2);
        w=zeros(m,m);
        for l=1:m;
            w(l,l)=exp(-(x1(j)-x(l,2))*(x1(j)-x(l,2))/(2*tau(i)*tau(i)))/2;
        end;
        theta=pinv(x'*w*x)*x'*w*y;
        y1(j)=theta(1)+theta(2)*x1(j);
    end;
    plot(x1,y1,colors(i));
end;
legend('trainingdata','linear','tau=0.1','tau=0.3','tau=0.8','tau=2',...
'tau=10');