Selasa, 17 April 2012

regresi polinom

%regresi polinom
clear; clc;

%input data
fin=fopen('data.dat');
data=textscan(fin,'%f %f','r');
fclose(fin);

x=data{1};
Y=data{2};

G1=[ones(length(x),1) x];
G2=[ones(length(x),1) x x.^2];
G3=[ones(length(x),1) x x.^2 x.^3];
G4=[ones(length(x),1) x x.^2 x.^3 x.^4];

next page