РГР / metod_naimenshih_kvadratov_m_fayl_dlya_matlab

function [z]=MinKvad(x,y)

[m,n]=size(x);

c=size(n,4);

for i=1:n

c(i,1:4)=[1 x(i) x(i)^2 y(i)];

end

k=1;

r=ones(3,4);

while k<=3

if k==1

for i=1:4

r(k,i)=sum(c(1:n,i));

end

elseif k==2

d=c;

for i=1:n

d(i,1:4)=c(i,1:4)*x(i);

end

for i=1:4

r(k,i)=sum(d(1:n,i));

end

elseif k==3

for i=1:n

d(i,1:4)=c(i,1:4)*(x(i)^2);

end

for i=1:4

r(k,i)=sum(d(1:n,i));

end

end

k=k+1;

end

A=size(3,n-1);

for j=1:3

for i=1:4

if i~=4

A(j,i)=r(j,i);

else

B(j,1)=r(j,i);

end

end

end

A

B

C=A\B

Q=x;

h=1/(n-1);

for i=1:n

Q(1,i)=C(1,1)+C(2,1)*x(1,i)+C(3,1)*x(1,i)^2;

end

Q

Prognoz=C(1,1)+C(2,1)*(x(1,i)+1)+C(3,1)*(x(1,i)+1)^2;

sum=0;

for i=1:n

Q(1,i)=((Q(1,i)-y(1,i))^2);

sum=sum+Q(1,i);

end

'Serednokvadratu4ne vidxulenya'

s=sqrt(sum*h)

'Prognoz'

z=Prognoz;