cost function计算:
function J = computeCost(X, y, theta)
%COMPUTECOST Compute cost for linear regression
% J = COMPUTECOST(X, y, theta) computes the cost of using theta as the
% parameter for linear regression to fit the data points in X and y
% Initialize some useful values
m = length(y); % number of training examples
% You need to return the following variables correctly
J = 0;
% ====================== YOUR CODE HERE ======================
% Instructions: Compute the cost of a particular choice of theta
% You should set J to the cost.
<span style="color:#cc0000;"><strong>h=X*theta;
sError=(h-y).^2;%注意点之前有空格 表示每个元素取平方 如果不加点则表示矩阵的平方
J=sum(sError)/(2*m);</strong></span>
% =========================================================================
end
gradience descent计算
function [theta, J_history] = gradientDescent(X, y, theta, alpha, num_iters)
%GRADIENTDESCENT Performs gradient descent to learn theta
% theta = GRADIENTDESENT(X, y, theta, alpha, num_iters) updates theta by
% taking num_iters gradient steps with learning rate alpha
% Initialize some useful values
m = length(y); % number of training examples
J_history = zeros(num_iters, 1);
for iter = 1:num_iters
% ====================== YOUR CODE HERE ======================
% Instructions: Perform a single gradient step on the parameter vector
% theta.
%
% Hint: While debugging, it can be useful to print out the values
% of the cost function (computeCost) and gradient here.
%
<span style="color:#ff0000;"> p=theta(1)-alpha/m*(sum((X*theta-y).*X(:,1)));
q=theta(2)-alpha/m*(sum((X*theta-y).*X(:,2)));
theta(1)=p; %要保证两个参数同时改变!!
theta(2)=q;</span>
% ============================================================
% Save the cost J in every iteration
J_history(iter) = computeCost(X, y, theta);
end
end