kurs_ras4et2

clear all;
y=[4.75 4.84 4.81 4.78
7.21 7.72 7.15 7.18
5.16 5.25 5.18 5.21
10.74 10.79 10.86 10.82];
yi=sum(y,2)/5;
S0i=var(y,0,2);
G0=max(S0i)/sum(S0i);
S0=sum(S0i)/4;
S0sr=S0/5;
b0=sum(yi)/4;
b1=(-yi(1)+yi(2)-yi(3)+yi(4))/4;
b2=(-yi(1)-yi(2)+yi(3)+yi(4))/4;
b12=(yi(1)-yi(2)-yi(3)+yi(4))/4;
Sb0=S0sr/4;
X1=-1:0.1:1;
X0=ones(length(X1));
X2=-1:0.1:1;
X12=[1;-1;-1;1];
[X1,X2]=meshgrid(X1,X2);
F=b0.*X0+b1.*X1+b2.*X2+b12.*X1.*X2;
S=Sb0*4;
Fn=F-2.*S.^0.5;
Fv=F+2.*S.^0.5;
figure;
meshc(X1,X2,F);
axis([-1.1 1.1 -1.1 1.1 5 30]);
colormap(copper);
title('Поверхность отклика');
xlabel('X1');
ylabel('X2');
zlabel('F(X1,X2)');
figure;
mesh(X1,X2,F);
colormap(winter);
title('Поверхность отклика при p=0.95');
xlabel('X1');
ylabel('X2');
zlabel('F(X1,X2)');
hold on;
mesh(X1,X2,Fn);
mesh(X1,X2,Fv);
axis([-1.1 1.1 -1.1 1.1 5 30]);