kurs_ras4et

clear all;
y=[5.024 3.562 3.240 2.632 1.863 1.731 1.610
   5.112 3.500 2.992 2.348 2.157 1.796 1.779
   4.669 3.804 2.616 2.139 2.173 1.615 1.533
   4.939 3.936 2.796 2.989 1.710 1.958 1.626
   5.168 3.514 3.164 2.449 2.249 1.979 1.560];
x=[0.0 0.5 1.0 1.5 2.0 2.5 3.0];
yi=sum(y)/5;
plot(x,y,'.',x,yi,'*');
grid on;
S0i=var(y);
G0=max(S0i)/sum(S0i);
S0=sum(S0i)/7;
S0sr=S0/6;
dx=x(1,2)-x(1,1);
z=(x-(x(1,1)-dx))/dx;
f0(1:7)=1;
b0=sum(yi)/7;
S1=sum((yi-b0.*f0).^2)/6;
Fp1=S1/S0sr;
f1=z-4;
b1=sum(yi.*f1)/sum(f1.^2);
S2=sum((yi-b0.*f0-b1.*f1).^2)/5;
Fp2=S2/S0sr;
figure;
plot(z,yi,'.');
hold on;
grid on;
F0=(f0.*b0)';
plot(z,F0,'g');
F1=(f0.*b0+f1.*b1);
plot(z,F1,'r');
f2=f1.*f1-f0.*(7^2-1)/(4*(4-1));
b2=sum(yi.*f2)/sum(f2.^2);
F2=(f0.*b0+f1.*b1+f2.*b2);
plot(z,F2);
S3=sum((yi-b0.*f0-b1.*f1-b2.*f2).^2)/4;
Fp3=S3/S0sr;
f3=f2.*f1-f1.*(7^2-2^2)/(4*(4*2^2-1))*2^2;
b3=sum(yi.*f3)/sum(f3.^2);
F3=(f0.*b0+f1.*b1+f2.*b2+f3.*b3);
plot(z,F3,'m');
S4=sum(yi-b0.*f0-b1.*f1-b2.*f2-b3.*f3).^2/3;
Fp4=S4/S0sr;
Sb0=S0sr/sum(f0.^2);
Sb1=S0sr/sum(f1.^2);
Sb2=S0sr/sum(f2.^2);
Sb3=S0sr/sum(f3.^2);
S=Sb0*f0.^2+Sb1*f1.^2+Sb2*f2.^2+Sb3*f3.^2;
Fv=F3+2.*S.^0.5;
Fn=F3-2.*S.^0.5;
figure;
plot(x,Fv,':',x,Fn,':',x,F3);
grid on;
St0=b0/Sb0;
St1=b1/Sb1;
St2=b2/Sb2;
St3=b3/Sb3;