ex2_4.m

good code for matlab by mili , i than you

%	Chapter 2: Optimization with Matlab
%  Dr. P.Venkataraman
%  Example 2.4 Sec.2.3
%
%	graphical solution using matlab (two design variables)
%  Minimum fin volume for efficient heat transfer
%  material Aluminum
%----------------------------------------------------
% global statement is used to share same information between 
% various m-files
global N H K W AREA 
%-----------------------------------------------------
% Initialize values  
N = 20						% number of fins
W = 0.1						% width of fins
H = 50.0	% convection coefficient W/m*m
K = 177.0	% thermal conductivity W/m-K
AREA = 0.015 % available fin foot print area 
%-----------------------------------------------------------
% right hand limits for the functions
h1val = AREA;
g1val = 0.95;
g2val = 0.94;
%---------------------------------------------------------
x1=0.001:0.0001:0.005;	 
x2=0.01:0.001:0.03;% x1 and x2 are vectors filled with numbers 
[X1 X2] = meshgrid(x1,x2);
% generates matrices X1 and X2 corresponding to
% vectors x1 and x2


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

% Constraints are evaluated 
eq1 = eq1_ex4(X1,X2);

[ineq1, ineq2]  = ineq_ex4(X1,X2);


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

[C1,h1] = contour(x1,x2,f1,[0.00001, 0.00002, 0.00004, ...
      0.00006, 0.00008, 0.0001],'g-');
clabel(C1,h1);
set(gca,'xtick',[0.001 0.0015 0.002 0.0025 0.003 0.0035 ...
      0.004 0.0045 0.005])
set(gca,'ytick',[0.010 0.015 0.02 0.025 0.03 0.125]);
xlabel('fin length','FontName','times','FontSize',12);
ylabel('fin width','FontName','times','FontSize',12)
grid
%hold on

figure
contour(x1,x2,eq1,[h1val,h1val],'r-');
set(gca,'xtick',[0.001 0.0015 0.002 0.0025 0.003 0.0035 ...
      0.004 0.0045 0.005])
set(gca,'ytick',[0.010 0.015 0.02 0.025 0.03 0.125]);

xlabel('fin length','FontName','times','FontSize',12);
ylabel('fin width','FontName','times','FontSize',12)
grid


figure
contour(x1,x2,ineq1,[0.968,0.968],'r-');
hold on
contour(x1,x2,ineq1,[0.969,0.969],'b-');
set(gca,'xtick',[0.001 0.0015 0.002 0.0025 0.003 0.0035 ...
      0.004 0.0045 0.005])

set(gca,'ytick',[0.010 0.015 0.02 0.025 0.03 0.125]);
xlabel('fin length','FontName','times','FontSize',12);
ylabel('fin width','FontName','times','FontSize',12)
hold off
grid

figure
contour(x1,x2,ineq2,[g2val,g2val],'r-');
hold on
contour(x1,x2,ineq2,[1.01*g2val,1.01*g2val],'b-');
set(gca,'xtick',[0.001 0.0015 0.002 0.0025 0.003 0.0035 ...
      0.004 0.0045 0.005])
set(gca,'ytick',[0.010 0.015 0.02 0.025 0.03 0.125]);
xlabel('fin length','FontName','times','FontSize',12);
ylabel('fin width','FontName','times','FontSize',12)
hold off
grid