3

x=[1; 5; 9; 17; 22; 27; 30; 37; 44; 48];
y=[21.74747774;32.10382622;43.24388629;118.1762691;180.6262061;163.4765075;223.4358954;261.5011088;298.3508735;331.5182586];
n=10
x1=sort(x);

% function y1=ypory(x,x1,y)
%Функция ypory упорядочивает массив значений функции y по порядку
% возрастания элементов массива аргументов x1. Входные параметры:
% x и y – исходные массивы аргументов x и функции y, x1 – вариационный ряд аргументов
% Выходной параметр y1 содержит значения функции y, соответствующие значениям
% вариационного ряда x1
[k, k1]=size(x);
if k1>k
for i=1:k
y1(i)=0;
end;
return;
end;
for i=1:k
for j=1:k
if abs(x1(i)-x(j))<10^-5
y1(i)=y(j);
end;
end;
end;
y1=y1';
a=0
for i=1:n
a=a+y1(i);
end;
a=a/n
%Вычислена предварительная оценка параметра а.
xmean=mean(x1);
ymean=mean(y1);
b=0;
c=0;
for i=1:n
b=b+(x1(i)-xmean)*y1(i);
c=c+(x1(i)-xmean)^2;
end;
b=b/c
% Найдена оценка параметра b.
a=a-b*xmean
% Исправлена оценка параметра а.
for i=1:n
yr(i)=a+b*x1(i);
end;
plot(x1,yr,'-*r',x1,y1,'-+g'),grid
bet=0.95;
alf=1-bet;
t=tinv(1-alf/2,n-2)
D1=0;
for i=1:n
D1=D1+(y1(i)-yr(i))^2;
end;
D1=D1/(n-2);
aleft=a-t*sqrt(D1/(n-2));
aright= a+t*sqrt(D1/(n-2));
sum=0;
for i=1:n
sum=sum+(x1(i)-xmean)^2;
end;
epsb=sqrt(n/(n-2)*D1/sum);
bleft = b-t*epsb;
bright=b+t*epsb;
for i=1:n
eps(i)=y1(i)-yr(i);
yleft(i)=yr(i)-t*sqrt(D1/(n-2)*(1+n*(x1(i)-xmean)^2/sum));
yright(i)=yr(i)+t*sqrt(D1/(n-2)*(1+n*(x1(i)-xmean)^2/sum));
end;
plot(x1, eps,'-*r'),grid;
plot(x1, yr,'-*r',x1,yleft,'-+g',x1,yright,'-+b'),grid