kurs_ras4et3

clear all;
yi=[4.7950 7.3150 5.200];
S0sr=0.004;
m1=sqrt(3/2);
m2=sqrt(1/2);
b0=sum(yi)/3;
b1=(yi(2)-yi(1))/(2*m1);
b2=(2*yi(3)-yi(2)-yi(1))/(6*m2);
Sb=S0sr/3;
X1=[-m1 m1 0];
X2=[-m2 -m2 2*m2];
X0=[1 1 1];
F=X0.*b0+X1.*b1+X2.*b2;
C=[0.93 0.93 0.93];
fill3(X1,X2,F,C);
axis([-1.23 1.23 -1 1.5 5 30]);
grid on;
title('Поверхность отклика');
xlabel('X1');
ylabel('X2');
zlabel('F(X1,X2)');

S=Sb.*X0.^2+Sb.*X1.^2+Sb.*X2.^2;
Fn=F-2.*S.^0.5;
Fv=F+2.*S.^0.5;
figure;
c=[1 1 1];
fill3(X1,X2,F,C,X1,X2,Fn,c,X1,X2,Fv,c);
axis([-1.23 1.23 -1 1.5 5 30]);
grid on;
title('Поверхность отклика');
xlabel('X1');
ylabel('X2');
zlabel('F(X1,X2)');