ex2_31.m

good code for matlab by mili , i than you

%	Chapter 2: Optimization with Matlab
%  Dr. P.Venkataraman
%  Example 2.3 Sec.2.3
%
%	graphical solution using matlab (two design variables)
%  Optimal design of a Flag Pole for high winds
%  Ref. 2.4
%----------------------------------------------------
% global statement is used to share same information between 
% various m-files
global ELAS SIGALL TAUALL GAM FS GRAV
global RHO CD FLAGW SPEED LP L DELT
%-----------------------------------------------------
% Initialize values
ELAS = 200e+09;   		% Pa
SIGALL = 250E+06; 		% Pa
TAUALL = 145e+06;			%Pa
GAM	= 7860;				% kg/m3
FS = 2.5;					% factor of safety
GRAV = 9.81;				% gravitational acceleration
%-------------------
RHO = 1.225;				% kg/m3
CD = 1.0	;					% drag coefficient
FLAGW = 5000;				% N
SPEED = 156.46	;			% m/s
%-------------------
LP = 8;						% m
L = 10;
DELT = 0.05;				% m
%------------------
g1val = SIGALL/FS	% right hand side values
g2val = TAUALL/FS
g3val = DELT
g4val = 0.001
%-----------------------------------------------------------
x1=0.02:0.01:1;	% the semi-colon at the end prevents the echo 
x2=0.025:0.01:1;	% these are also the side constraints
% x1 and x2 are vectors filled with numbers 
% note a way to avoid x1 = x2
[X1 X2] = meshgrid(x1,x2);
% generates matrices X1 and X2 correspondin
% vectors x1 and x2


f1 = obj_ex3(X1,X2);
% the objecive function is evaluated over the entire mesh

% Constraints are evaluated 
ineq1 = ineq1_ex3(X1,X2);

ineq2 = ineq2_ex3(X1,X2);

ineq3 = ineq3_ex3(X1,X2);

ineq4 = ineq4_ex3(X1,X2);

%contour(x1,x2,f1,10);

[C1,h1] = contour(x1,x2,f1,[0,10000,50000,100000, ...
           150000, 200000, 250000, 300000],'g-');
%   [0 0.1 0.3 0.6 0.8 1.0 1.2 1.5 1.8 2.0 2.2 2.4 2.6 2.8 3.0]);
clabel(C1,h1);
set(gca,'xtick',[0 0.2 0.4 0.6 0.8 1.0])
set(gca,'ytick',[0 0.2 0.4 0.6 0.8 1.0])
xlabel('outside diameter','FontName','times','FontSize',12);
ylabel('inside diameter','FontName','times','FontSize',12)
grid
hold on


[C2,h2] = contour(x1,x2,ineq1,[g1val,g1val],'r-');
set(h2,'LineWidth',1)
k2 = gtext('g1');
set(k2,'FontName','Times','FontWeight','bold','FontSize',14,'Color','red')
grid


[C3,h3] = contour(x1,x2,ineq2,[g2val,g2val],'m-');
%clabel(C3,h3);
set(h3,'LineWidth',1)
k3 = gtext('g2');
set(k3,'FontName','Times','FontWeight','bold','FontSize',14,'Color','m')
grid


[C4,h4] = contour(x1,x2,ineq3,[g3val,g3val],'b-');
%clabel(C4,h4);
set(h4,'LineWidth',1)
k4 = gtext('g3');
set(k4,'FontName','Times','FontWeight','bold','FontSize',14,'Color','b')
grid

[C5,h5] = contour(x1,x2,ineq4,[g4val,g4val],'k-');
%clabel(C5,h5);
set(h5,'LineWidth',1)
k5 = gtext('g4');
set(k5,'FontName','Times','FontWeight','bold','FontSize',14,'Color','k')
grid

% the equality and inequality constraints are not written 
% with 0 on the right hand side. If you do write them that way
% you would have to include [0,0] in the contour commands

%xlabel('outside diameter','FontName','times','FontSize',12);	% label for x-axes
%ylabel('inside diameter','FontName','times','FontSize',12);
title('Flag Pole problem', ...
   'FontName','times','FontSize',10)
%grid

hold off