13-优化设计模型的几何描述-2.m

MATLAb 机械源码文件好东西啊

% 二维约束优化几何描述(例1组合)
% 按等间隔矢量产生二维网格矩阵
xx1=linspace(-1,3,40);
xx2=linspace(0,4,40);
[x1,x2]=meshgrid(xx1,xx2);
% 数学模型
f=x1.^2+x2.^2-4*x1+4;          % 目标函数f
g1=x1-x2+2;                    % 约束函数g1
g2=-x1.^2+x2-1;                % 约束函数g2
g3=x1;                         % 约束函数g3
% 设计空间
figure(1);
surfc(x1,x2,f);
title('\bf 目标函数和约束函数曲面');
xlabel('设计变量 \bf X1');
ylabel('设计变量 \bf X2');
zlabel('目标函数值和约束函数值');
hold on;                       % 保持图形
surfc(x1,x2,g1);
surfc(x1,x2,g2);
surfc(x1,x2,g3);
% 设计平面
figure(2);
h=contour(x1,x2,f);
clabel(h);
axis equal                     % 两坐标轴的定标因子相等
title('\bf 设计平面');
xlabel('设计变量 \bf X1');
ylabel('设计变量 \bf X2');
hold on;
h=contour(x1,x2,g1);
clabel(h);
h=contour(x1,x2,g2);
clabel(h);
h=contour(x1,x2,g3);
clabel(h);