eg10_1

code1:

clear

clc

a=[0.1 1;0.2 1.1;0.3 1.2;0.4 1.3;0.5 1.1;0.6 0.9;0.7 0.7;0.8 1.0;0.9 1.2;1 1.3];

plot(a(:,1),a(:,2),'*');%fig(b)

axis([0.1 1 0.5 1.5]);

pause

p=polyfit(a(:,1),a(:,2),4);

fprintf('order of polynomial is 4');

p

t=[0.1:0.001:0.95];

nev=polyval(p,t);

plot(t,nev,a(:,1),a(:,2),'*');%fig(c)

axis([0 1.1 0.6 1.4]);

grid

title('order of polynomial is 4');

xlabel('input data');

ylabel('output data');

pause

fig(a)

fig(b)

fig(c)