machine-learning ex1_2

%Load Data

data = load('ex1data2.txt');

X = data(:, 1:2);
y = data(:, 3);
m = length(y);

[X mu sigma] = featureNormalize(X);

function [X_norm, mu, sigma] = featureNormalize(X)
X_norm = X;
mu = zeros(1, size(X, 2));
sigma = zeros(1, size(X, 2));
%disp(X);
mu =mean(X);
%fprintf('%f %f\n',mu(1),mu(2));
sigma=std(X);
%fprintf('%f %f\n',sigma(1),sigma(2))
n=size(X,2);

for iter=1:n
    X_norm(:,iter)=(X_norm(:,iter)-mu(iter))/sigma(iter);
end

X = [ones(m, 1) X];

alpha = 0.01;
num_iters = 400;

theta = zeros(3, 1);
[theta, J_history] = gradientDescentMulti(X, y, theta, alpha, num_iters);

function [theta, J_history] = gradientDescentMulti(X, y, theta, alpha, num_iters)
m = length(y); % number of training examples
J_history = zeros(num_iters, 1);
for iter = 1:num_iters
    delta=((theta' * X' - y')* X)';
    theta=theta - alpha / m * delta;
    J_history(iter) = computeCostMulti(X, y, theta);
    fprintf('%f \n',J-history(iter));

end

end

figure;
plot(1:numel(J_history), J_history, '-b', 'LineWidth', 2);
xlabel('Number of iterations');
ylabel('Cost J');
fprintf('Theta computed from gradient descent: \n');
fprintf(' %f \n', theta);
fprintf('\n');

price = 0;

fprintf(['Predicted price of a 1650 sq-ft, 3 br house ' ...
         '(using gradient descent):\n $%f\n'], price);


fprintf('Program paused. Press enter to continue.\n');
pause;

fprintf('Solving with normal equations...\n');

%% Load Data
data = csvread('ex1data2.txt');
X = data(:, 1:2);
y = data(:, 3);
m = length(y);
% Add intercept term to X
X = [ones(m, 1) X];
% Add intercept term to X
X = [ones(m, 1) X];

fprintf('Theta computed from the normal equations: \n');
fprintf(' %f \n', theta);
fprintf('\n');

price = 0;

fprintf(['Predicted price of a 1650 sq-ft, 3 br house ' ...
         '(using normal equations):\n $%f\n'], price);