clear all;
x0=[0.8;160;2.0;28];
xx=[0.3;6;0.6;6];%下限
sx=[1;220;6;60];%上限
[x,fn]=fmincon(@mubiao,x0,[],[],[],[],xx,sx,@yueshu);%求解函数模型
%下面是凑你那个答案的，不凑的话就删了下面这小段，把其他的x1,x2，x3,x4改成x(1),x(2),x(3),x(4)就行
x1=x(1)+0.5;
x2=x(2)+92.9548;
x3=x(3)+1.3;
x4=x(4)+20;
%就是上面这小段
F=16*x1^2*x4+3*x3^2*x4+14*x1*x3*x4+408*x3*x4;
disp('最优解');
fprintf(1,'x1=%3.4f\n',x1);
fprintf(1,'x2=%3.4f\n',x2);
fprintf(1,'x3=%3.4f\n',x3);
fprintf(1,'x4=%3.4f\n',x4);
fprintf(1,'目标函数minf(x)=%3.4f\n',F);
