matlab曲面分析回归方程

matlab曲面分析回归方程 Y1=2.167-0.018x1+0.01x2-0.009x3+0.012x4+0.085x1^2+0.074x3^2 Y2=39.872-0.198x1+0.085x2-0.264x3-0.105x4+0.617x2^2 y121670018?x1001?x20009?x30012?x40085?x120074?x32=.-.+.-.+.+.+. y2398720198?x10085?x20264?x30105?x40617?x22=.-.+.-.-.+. clear all; %响面应应应 x1=50:0.01:70; x2=30:0.01:50; [X1,X2]=meshgrid(x1,x2);X3=50;X4=70; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X1,X2,Y1) shading flat rotate3d clear all; %响面应应应 x2=30:0.01:50; x3=40:0.01:60; [X2,X3]=meshgrid(x2,x3);X1=60;X4=70; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X2,X3,Y1) shading flat rotate3d clear all; %响面应应应 x1=50:0.01:70; x3=40:0.01:60; [X1,X3]=meshgrid(x1,x3);X2=40;X4=70; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X1,X3,Y1) shading flat rotate3d clear all; %响面应应应 x1=50:0.01:70; x4=60:0.01:80; [X1,X4]=meshgrid(x1,x4); X2=40;X3=50; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X1,X4,Y1) shading flat rotate3d clear all; %响面应应应 x2=30:0.01:50; x4=60:0.01:80; [X2,X4]=meshgrid(x2,x4);X1=60;X3=50; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X2,X4,Y1) shading flat rotate3d clear all; %响面应应应 x3=40:0.01:60; x4=60:0.01:80; [X3,X4]=meshgrid(x3,x4);X1=60;X2=40; Y1=2.167-0.018.*X1+0.01.*X2-0.009.*X3+0.012.*X4+0.085.*X1.^2+0.074.*X3.^2; surfc(X3,X4,Y1) shading flat rotate3d clear all; %响面应应应 x1=50:0.01:70; x2=30:0.01:50; [X1,X2]=meshgrid(x1,x2);X3=50;X4=70; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X1,X2,Y2) shading flat rotate3d clear all; %响面应应应 x1=50:0.01:70; x3=40:0.01:60; [X1,X3]=meshgrid(x1,x3);X2=40;X4=70; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X1,X3,Y2) shading flat rotate3d clear all; %响面应应应 x1=50:0.01:70; x4=60:0.01:80; [X1,X4]=meshgrid(x1,x4);X2=40;X3=50; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X1,X4,Y2) shading flat rotate3d clear all; %响面应应应 x2=30:0.01:50; x3=40:0.01:60; [X2,X3]=meshgrid(x2,x3);X1=60;X4=70; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X2,X3,Y2) shading flat rotate3d clear all; %响面应应应 x2=30:0.01:50; x4=60:0.01:80; [X2,X4]=meshgrid(x2,x4);X1=60;X3=50; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X2,X4,Y2) shading flat rotate3d clear all; %响面应应应 x3=40:0.01:60; x4=60:0.01:80; [X3,X4]=meshgrid(x3,x4);X1=60;X2=40; Y2=39.872-0.198*X1+0.085*X2-0.264*X3-0.105*X4+0.617*X2.^2; surfc(X3,X4,Y2) shading flat rotate3d x = 70 30 60 80 fval = -64.4480 exitflag = 1 f1 = 684.5270 f2 = 559.6220 function f=ff12(x) f=4.647-0.54*x(1)+0.225*x(2)-0.766*x(3)-0.351*x(4)-0.043*x(1)^2+0.309*x(2)^2-0.037*x(3)^2 lb=[-1.546;-1.546;-1.546;-1.546];ub=[1.546;1.546;1.546;1.546]; x0=[-2;-2;-2;-2]; [x,fval,exitflag]=fmincon(@ff12,x0,[],[],[],[],lb,ub)f1=2.04-0.106*x(1)+0.062*x(2)+0.072*x(4)-0.053*x(3)+0.085*x(1)^2+0.074*x(3)^2 f2=11.333-1.586*x(3)-1.186*x(1)-0.63*x(4)+0.512*x(2)+0.617*x(2)^2 x = 1.5460 -0.3641 1.5460 1.5460 fval = 1.8531exitflag = 1 f1 = 2.2630f2 = 5.9689