surfacearea.m

isight调用matlab接口进行优化的例子

% Function OgiveSurfaceArea(Length, BaseDiameter)
%
%This function will estimate the SurfaceArea of an Ogive nose with the specified length and
% base diameter
%
% Functions/Subroutines Required
%       OgiveRadius
%
%       Length           - Length of Ogive Nose
%       BaseDiameter     - Diameter of the base of the Ogive Nose
%
%       Total Surface area of the Ogive nose

% Length = 25.0
% BaseDiameter = 10.0

NumberOfSlices = 25

DeltaX = Length / NumberOfSlices

RMid = BaseDiameter / 2

FinenessRatio = Length / BaseDiameter
RBar = (FinenessRatio ^ 2) + 0.25
RadiusOfCurvature = RBar * BaseDiameter
XBar = DeltaX / Length

TermA = (FinenessRatio ^ 2) / (RBar ^ 2)
TermB = (1 - TermA * ((1 - XBar) ^ 2)) ^ 0.5
TermC = 1 - (2 * RBar * (1 - TermB))
Radius = RMid * TermC
if (Radius < 0) 
    OgiveRadius = 0
else
    OgiveRadius = RMid * TermC
end ;





SurfaceArea = 3.14 * OgiveRadius * sqrt((OgiveRadius ^ 2) + (DeltaX ^ 2))
XOld = DeltaX
RadiusOld = OgiveRadius

i = 1

while ( i <= (NumberOfSlices - 1) )
    X = XOld + DeltaX
    
     
    
    RMid = BaseDiameter / 2
    
    FinenessRatio = Length / BaseDiameter
    RBar = (FinenessRatio ^ 2) + 0.25
    RadiusOfCurvature = RBar * BaseDiameter
    XBar = DeltaX / Length
    
    TermA = (FinenessRatio ^ 2) / (RBar ^ 2)
    TermB = (1 - TermA * ((1 - XBar) ^ 2)) ^ 0.5
    TermC = 1 - (2 * RBar * (1 - TermB))
    Radius = RMid * TermC
    if (Radius < 0) 
        Radius = 0
    else
        Radius = RMid * TermC
   end ;
    
    
    
    SArea = 3.14 * (Radius + RadiusOld) * sqrt(((Radius - RadiusOld) ^ 2) + DeltaX ^ 2)
    SurfaceArea = SurfaceArea + SArea
    XOld = X
    RadiusOld = Radius
    
    i = i + 1
end 

OgiveSurfaceArea = SurfaceArea