机器学习第二次编程作业（week3）

这次作业完成加修改花了我五个小时作业（人还是太傻了）
这是要完成的作业

这是我完成的情况（我第一次才三十分，共上传了五次左右吧，做错不要紧，知道自己错哪了才重要）

plotData.m文件

function plotData(X, y)
%PLOTDATA Plots the data points X and y into a new figure 
%   PLOTDATA(x,y) plots the data points with + for the positive examples
%   and o for the negative examples. X is assumed to be a Mx2 matrix.

% Create New Figure
figure; hold on;

% ====================== YOUR CODE HERE ======================
% Instructions: Plot the positive and negative examples on a
%               2D plot, using the option 'k+' for the positive
%               examples and 'ko' for the negative examples.
%以特征1,2分别为坐标轴，在图上标出各点，y=1和y=0分别用‘+’和‘o’标出

pos = find(y == 1);#返回一个列向量，向量的元素是在y向量中y=1元素的位置
neg = find(y == 0);
plot(X(pos,1),X(pos,2),'k+')#X(pos,1)返回一个列向量，元素是以pos的元素为行，然后取第一列，
plot(X(neg,1),X(neg,2),'ko')






% =========================================================================



hold off;

end

sigmoid.m

function g = sigmoid(z)
%SIGMOID Compute sigmoid function
%   g = SIGMOID(z) computes the sigmoid of z.

% You need to return the following variables correctly 
g = zeros(size(z));

% ====================== YOUR CODE HERE ======================
% Instructions: Compute the sigmoid of each value of z (z can be a matrix,
%               vector or scalar).
g = 1./(1+exp(-z));



% =============================================================

end

costFunction.m

function [J, grad] = costFunction(theta, X, y)
%COSTFUNCTION Compute cost and gradient for logistic regression
%   J = COSTFUNCTION(theta, X, y) computes the cost of using theta as the
%   parameter for logistic regression and the gradient of the cost
%   w.r.t. to the parameters.

% Initialize some useful values
m = length(y); % number of training examples

% You need to return the following variables correctly 
J = 0;
grad = zeros(size(theta));

% ====================== YOUR CODE HERE ======================
% Instructions: Compute the cost of a particular choice of theta.
%               You should set J to the cost.
%               Compute the partial derivatives and set grad to the partial
%               derivatives of the cost w.r.t. each parameter in theta
%
% Note: grad should have the same dimensions as theta
%
h = sigmoid(X * theta)

J = 1/m * sum(-y.*log(h) - (1-y).*log(1-h)) #想不出来的话，用笔在纸上把矩阵写出来

for i = 1:size(theta,1)
  sum1=0
  for j = 1:m
       sum1 = sum1 + (h(j)-y(j))*X(j,i)
  end
  grad(i) = 1/m*sum1
  
end



% =============================================================

end

注意公式

predict.m

function p = predict(theta, X)
%PREDICT Predict whether the label is 0 or 1 using learned logistic 
%regression parameters theta
%   p = PREDICT(theta, X) computes the predictions for X using a 
%   threshold at 0.5 (i.e., if sigmoid(theta'*x) >= 0.5, predict 1)

m = size(X, 1); % Number of training examples

% You need to return the following variables correctly
p = zeros(m, 1);

% ====================== YOUR CODE HERE ======================
% Instructions: Complete the following code to make predictions using
%               your learned logistic regression parameters. 
%               You should set p to a vector of 0's and 1's
%
h = sigmoid(X *theta)
for i = 1:m
    if h(i)>=0.5
      p(i) = 1
    else
      p(i)= 0
    end
end

costFunctionReg.m

function [J, grad] = costFunctionReg(theta, X, y, lambda)
%COSTFUNCTIONREG Compute cost and gradient for logistic regression with regularization
%   J = COSTFUNCTIONREG(theta, X, y, lambda) computes the cost of using
%   theta as the parameter for regularized logistic regression and the
%   gradient of the cost w.r.t. to the parameters. 

% Initialize some useful values
m = length(y); % number of training examples

% You need to return the following variables correctly 
J = 0;
grad = zeros(size(theta));

% ====================== YOUR CODE HERE ======================
% Instructions: Compute the cost of a particular choice of theta.
%               You should set J to the cost.
%               Compute the partial derivatives and set grad to the partial
%               derivatives of the cost w.r.t. each parameter in theta

h = sigmoid(X * theta)

J = 1/m * sum(-y.*log(h) - (1-y).*log(1-h))+lambda/(2*m)*sum(theta(2:end).^2)#这里的theta是从第二行开始，第一行是theta0，一般来说不包括进去

for i = 1:length(theta)
  sum1= 0
  for j = 1:m
       sum1 = sum1 + (h(j)-y(j))*X(j,i);
  end
  if i == 1 #注意这里要分情况讨论
    grad(i) = 1/m * sum1 ;
  else
    grad(i) = 1/m * sum1 + (lambda/m) * theta(i);
  end
  
end




% =============================================================

end

用到的公式

（多说一句，if语句那里我一开始是写等于号的，找了好久都没发现自己的错误，后面执行submit的时候，会执行所有文件，我看到了命令行窗口上的warning才发现了自己的错误）

如果以上有什么错误的地方或者你有什么建议或想法，欢迎找我讨论