3.7 Исходные тексты сценариев

// вычисляет значения полинома с коэффициентами a на массиве x

function[pol]=polyval(x, a);

pol=0;

for i=1:1:length(a);

pol=pol+a(i)*x.^(i-1);

end;

endfunction;

//вычисляет массив разностей полинома c набором коэффициентов a на //массиве данных z

function[zr]=R(a, z);

pol=0;

for i=1:1:length(a);

pol=pol+a(i)*z(1).^(i-1);

end;

zr=z(2)-pol;

endfunction;

//Вычисляет набор коэффициентов полинома МНК s по массивам данных x, y

// err - ошибка аппроксимации

function[s, err]=MNK(x, y, a);

z=[x;y];

exec("R.sce");

[s,err]=datafit(R,z,a);

endfunction;

//вычисляет СКО полинома с набором коэффициентов a на массивах данных x и y

function[err1]=SKO(x, y, a);

err=0;

exec("polyval.sce");

fori=1:1:length(x);

err=err+(y-polyval(x,a)).^2;

end;

err1=err*err';

err1=sqrt(err1/length(x));

endfunction

//Исходный текст сценария PrimerMNK.sce

//По данным x и y получает полиномы 1, 2, 3, 4, 5 степеней по МНК

// вычисляет СКО этих полиномов и в графическом окне строит исходные //данные и графики полиномов

clf

z=[x;y];

plot(x,y,'red*');

a1=[0;0];

exec("MNK.sce");

[s1,err1]=MNK(x,y,a1);

exec("SKO.sce");

sko1=SKO(x,y,s1);

exec("polyval.sce");

a1=s1'

sko1

xx=x(1):(x(length(x))-x(1))/100:x(length(x));

y1=polyval(xx,s1);

plot(xx,y1,'blue');

sleep(1000);

a2=[0;0;0];

[s2,err2]=MNK(x,y,a2);

sko2=SKO(x,y,s2);

a2=s2'

sko2

y2=polyval(xx,s2);

plot(xx,y2,'green');

sleep(1000);

a3=[0;0;0;0];

[s3,err3]=MNK(x,y,a3);

sko3=SKO(x,y,s3);

a3=s3'

sko3

y3=polyval(xx,s3);

plot(xx,y3,'cyan');

sleep(1000);

a4=[0;0;0;0;0];

[s4,err4]=MNK(x,y,a4);

sko4=SKO(x,y,s4);

a4=s4'

sko4

y4=polyval(xx,s4);

plot(xx,y4,'magenta');

sleep(1000);

a5=[0;0;0;0;0;0];

[s5,err5]=MNK(x,y,a5);

sko5=SKO(x,y,s5);

a5=s5'

sko5

y5=polyval(xx,s5);

plot(xx,y5,'grey');

sleep(1000);