吴恩达机器学习作业2

plotData

pos = find(y == 1)
neg = find(y == 0)
% 属于分类1的样本赋值为pos，属于分类2的样本赋值为neg

plot(X(pos,1),X(pos,2),'k+','LineWidth',2,'MarkerSize',7)
plot(X(neg,1),X(neg,2),'ko','MarkerFaceColor','y','MarkerSize',7)
% X(pos, 1), X(pos, 2)表示分类1样本的横纵坐标，'k+'表示线性为加号型、黑色，LineWidth', 2表示线宽为2，'MarkerSize', 7表示形状大小为7
% MarkerEdgeColor','y'表示圆框边缘颜色为蓝色

sigmoid

g = zeros(size(z))
for i = 1:size(z,1)
  for j = 1:size(z,2)
    g(i,j) = 1 / (1 + e ^ (-z(i,j)))
  endfor
endfor

costFunction

h = zeros(m,1);
hx = X * theta;
h = sigmoid(hx);

for i = 1:m
  J = J + (-y(i)) * log(h(i)) - (1.-y(i)) * log(1.-h(i))
endfor
J = J / m;

sum = zeros(size(X,2),1);
for i  = 1:m
  sum = sum + (h(i) - y(i)) .* X(i,:)'
endfor
grad = sum / m  

predict

h = sigmoid(X * theta)
for i = 1:m
  if h(i) >= 0.5 
    P = 1
  else
    p = 0
  endif
endfor

costFunctionReg

h = sigmoid(X * theta)
for i = 1:m
  J = J + (-y(i)) * log(h(i)) - (1-y(i)) * log(1-h(i))
endfor
J = J / m
reg = 0
for j = 2 : size(X,2)
  reg = reg + theta(j) * theta(j)
endfor
reg = lambda * reg / (2 * m )
J = J + reg

sum = zeros(size(X,2),1)
for i = 1:m
  sum = sum + (h(i) - y(i)) .* X(i,:)'
endfor
sum = sum / m
grad(1) = sum(1)
for j = 2:size(X,2)
  grad(j) = sum(j) + lambda * theta(j) / m 
endfor