一复杂函数的优化设计——matlab程序

博大广阔

function aaaaa=aaa()
clc
clear all
  
F0=10000000000;
D=rand(1,20)*20+70;
d=rand(1,20)*2+4;
n=rand(1,20)*18+30;
phi=rand(1,20)*(120.3931-107.8064)+107.8064;
for i=1:1:length(d)
    for j=1:1:length(D)
        for k=1:1:length(n)
            for L=1:1:length(phi)                                 
                    XX=[d(i),D(j),n(k),phi(L)];
                  F=F_f(XX);                    
                    if F<=F0
                        F0=F;
                        sa=XX;
                    end        
                              
            end
        end
    end
end

XX_0=sa;
F_0=res(XX_0)


LB=[4;70;30;101.5];  LU=[6;90;48;126.7];
  option= optimset('Algorithm','active-set');%比较五个算法
% must be 'active-set', 'trust-region-reflective', 'interior-point', 'interior-point-convex', 'levenberg-marquardt',
% 'trust-region-dogleg', 'lm-line-search', or 'sqp'.
X0=sa;
[XX_1,value]=fmincon(@F_f,X0,[],[],[],[],LB,LU,@contion,option)
F_1=res(XX_1)



%XX_2=[5.03033,80.8013,47.5583,113.5];
% F_2=res(XX_2)  
end



function F=F_f(XX)
%等价的性能函数
d=XX(1);D=XX(2);n=XX(3);phi=XX(4);
F =((4.131*d^8*n^2)/D^2 - (14944.0*d^4*n)/D - (3.743*d^8*n^2*phi)/D^2 + (0.8479*d^8*n^2*phi^2)/D^2 + (6769.0*d^4*n*phi)/D + 1.834e7)^(1/2) + ((6739.0*d^4*n)/D + (3.297*d^8*n^2)/D^2 - (3.706*d^8*n^2*phi)/D^2 + (1.041*d^8*n^2*phi^2)/D^2 - (3787.0*d^4*n*phi)/D + 5.761e6)^(1/2) + ((5.017*d^8*n^2)/D^2 - (24333.0*d^4*n)/D - (4.568*d^8*n^2*phi)/D^2 + (1.04*d^8*n^2*phi^2)/D^2 + (11088.0*d^4*n*phi)/D + 3.357e7)^(1/2) + ((5.264*d^8*n^2)/D^2 - (26388.0*d^4*n)/D - (4.816*d^8*n^2*phi)/D^2 + (1.102*d^8*n^2*phi^2)/D^2 + (12077.0*d^4*n*phi)/D + 3.744e7)^(1/2) + ((4.419*d^8*n^2)/D^2 - (18566.0*d^4*n)/D - (3.997*d^8*n^2*phi)/D^2 + (0.9037*d^8*n^2*phi^2)/D^2 + (8396.0*d^4*n*phi)/D + 2.347e7)^(1/2) + ((11422.0*d^4*n)/D + (4.373*d^8*n^2)/D^2 - (4.852*d^8*n^2*phi)/D^2 + (1.346*d^8*n^2*phi^2)/D^2 - (6337.0*d^4*n*phi)/D + 9.647e6)^(1/2) + ((6570.0*d^4*n)/D + (2.585*d^8*n^2)/D^2 - (2.889*d^8*n^2*phi)/D^2 + (0.8072*d^8*n^2*phi^2)/D^2 - (3671.0*d^4*n*phi)/D + 9.107e6)^(1/2) + ((5.533*d^8*n^2)/D^2 - (28444.0*d^4*n)/D - (5.095*d^8*n^2*phi)/D^2 + (1.173*d^8*n^2*phi^2)/D^2 + (13100.0*d^4*n*phi)/D + 4.144e7)^(1/2) + ((7042.0*d^4*n)/D + (2.628*d^8*n^2)/D^2 - (2.928*d^8*n^2*phi)/D^2 + (0.8158*d^8*n^2*phi^2)/D^2 - (3923.0*d^4*n*phi)/D + 1.007e7)^(1/2) + ((13222.0*d^4*n)/D + (4.74*d^8*n^2)/D^2 - (5.233*d^8*n^2*phi)/D^2 + (1.444*d^8*n^2*phi^2)/D^2 - (7297.0*d^4*n*phi)/D + 1.146e7)^(1/2) + ((10244.0*d^4*n)/D + (2.978*d^8*n^2)/D^2 - (3.25*d^8*n^2*phi)/D^2 + (0.8868*d^8*n^2*phi^2)/D^2 - (5591.0*d^4*n*phi)/D + 1.72e7)^(1/2) + ((5.822*d^8*n^2)/D^2 - (30399.0*d^4*n)/D - (5.403*d^8*n^2*phi)/D^2 + (1.253*d^8*n^2*phi^2)/D^2 + (14100.0*d^4*n*phi)/D + 4.53e7)^(1/2) + ((4.015*d^8*n^2)/D^2 - (13144.0*d^4*n)/D - (3.648*d^8*n^2*phi)/D^2 + (0.8286*d^8*n^2*phi^2)/D^2 + (5968.0*d^4*n*phi)/D + 1.638e7)^(1/2) + ((5149.0*d^4*n)/D + (2.664*d^8*n^2)/D^2 - (3.007*d^8*n^2*phi)/D^2 + (0.8485*d^8*n^2*phi^2)/D^2 - (2906.0*d^4*n*phi)/D + 5.484e6)^(1/2) + ((5237.0*d^4*n)/D + (2.742*d^8*n^2)/D^2 - (3.095*d^8*n^2*phi)/D^2 + (0.8735*d^8*n^2*phi^2)/D^2 - (2956.0*d^4*n*phi)/D + 5.322e6)^(1/2) + ((10955.0*d^4*n)/D + (3.064*d^8*n^2)/D^2 - (3.325*d^8*n^2*phi)/D^2 + (0.9022*d^8*n^2*phi^2)/D^2 - (5944.0*d^4*n*phi)/D + 1.9e7)^(1/2) + ((5144.0*d^4*n)/D + (2.605*d^8*n^2)/D^2 - (2.939*d^8*n^2*phi)/D^2 + (0.8291*d^8*n^2*phi^2)/D^2 - (2902.0*d^4*n*phi)/D + 5.73e6)^(1/2) + ((5419.0*d^4*n)/D + (2.843*d^8*n^2)/D^2 - (3.207*d^8*n^2*phi)/D^2 + (0.9046*d^8*n^2*phi^2)/D^2 - (3057.0*d^4*n*phi)/D + 5.25e6)^(1/2) + ((5214.0*d^4*n)/D + (2.564*d^8*n^2)/D^2 - (2.891*d^8*n^2*phi)/D^2 + (0.8149*d^8*n^2*phi^2)/D^2 - (2939.0*d^4*n*phi)/D + 6.058e6)^(1/2) + ((6.127*d^8*n^2)/D^2 - (32000.0*d^4*n)/D - (5.738*d^8*n^2*phi)/D^2 + (1.343*d^8*n^2*phi^2)/D^2 + (14988.0*d^4*n*phi)/D + 4.846e7)^(1/2) + ((5714.0*d^4*n)/D + (2.967*d^8*n^2)/D^2 - (3.344*d^8*n^2*phi)/D^2 + (0.9425*d^8*n^2*phi^2)/D^2 - (3221.0*d^4*n*phi)/D + 5.283e6)^(1/2) + ((7.28*d^8*n^2)/D^2 - (26888.0*d^4*n)/D - (7.149*d^8*n^2*phi)/D^2 + (1.755*d^8*n^2*phi^2)/D^2 + (13200.0*d^4*n*phi)/D + 3.522e7)^(1/2) + ((5352.0*d^4*n)/D + (2.539*d^8*n^2)/D^2 - (2.86*d^8*n^2*phi)/D^2 + (0.8054*d^8*n^2*phi^2)/D^2 - (3014.0*d^4*n*phi)/D + 6.47e6)^(1/2) + ((11622.0*d^4*n)/D + (3.15*d^8*n^2)/D^2 - (3.398*d^8*n^2*phi)/D^2 + (0.9163*d^8*n^2*phi^2)/D^2 - (6266.0*d^4*n*phi)/D + 2.082e7)^(1/2) + ((11122.0*d^4*n)/D + (7.173*d^8*n^2)/D^2 - (7.524*d^8*n^2*phi)/D^2 + (1.973*d^8*n^2*phi^2)/D^2 - (5829.0*d^4*n*phi)/D + 8.627e6)^(1/2) + ((6144.0*d^4*n)/D + (3.117*d^8*n^2)/D^2 - (3.51*d^8*n^2*phi)/D^2 + (0.9879*d^8*n^2*phi^2)/D^2 - (3459.0*d^4*n*phi)/D + 5.442e6)^(1/2) + ((3505.0*d^4*n)/D + (3.611*d^8*n^2)/D^2 - (3.485*d^8*n^2*phi)/D^2 + (0.8407*d^8*n^2*phi^2)/D^2 - (1691.0*d^4*n*phi)/D + 1.831e7)^(1/2) + ((5557.0*d^4*n)/D + (2.53*d^8*n^2)/D^2 - (2.846*d^8*n^2*phi)/D^2 + (0.8002*d^8*n^2*phi^2)/D^2 - (3125.0*d^4*n*phi)/D + 6.97e6)^(1/2) + ((4.264*d^8*n^2)/D^2 - (16744.0*d^4*n)/D - (3.858*d^8*n^2*phi)/D^2 + (0.8727*d^8*n^2*phi^2)/D^2 + (7575.0*d^4*n*phi)/D + 2.072e7)^(1/2) + ((9855.0*d^4*n)/D + (4.046*d^8*n^2)/D^2 - (4.508*d^8*n^2*phi)/D^2 + (1.256*d^8*n^2*phi^2)/D^2 - (5490.0*d^4*n*phi)/D + 8.176e6)^(1/2) + ((7.035*d^8*n^2)/D^2 - (30844.0*d^4*n)/D - (6.818*d^8*n^2*phi)/D^2 + (1.652*d^8*n^2*phi^2)/D^2 + (14944.0*d^4*n*phi)/D + 4.391e7)^(1/2) + ((12200.0*d^4*n)/D + (3.233*d^8*n^2)/D^2 - (3.465*d^8*n^2*phi)/D^2 + (0.9282*d^8*n^2*phi^2)/D^2 - (6536.0*d^4*n*phi)/D + 2.258e7)^(1/2) + ((5228.0*d^4*n)/D + (3.608*d^8*n^2)/D^2 - (3.514*d^8*n^2*phi)/D^2 + (0.8555*d^8*n^2*phi^2)/D^2 - (2546.0*d^4*n*phi)/D + 1.998e7)^(1/2) + ((6.747*d^8*n^2)/D^2 - (32722.0*d^4*n)/D - (6.458*d^8*n^2*phi)/D^2 + (1.546*d^8*n^2*phi^2)/D^2 + (15666.0*d^4*n*phi)/D + 4.884e7)^(1/2) + ((5828.0*d^4*n)/D + (2.535*d^8*n^2)/D^2 - (2.847*d^8*n^2*phi)/D^2 + (0.799*d^8*n^2*phi^2)/D^2 - (3272.0*d^4*n*phi)/D + 7.569e6)^(1/2) + ((6851.0*d^4*n)/D + (3.606*d^8*n^2)/D^2 - (3.545*d^8*n^2*phi)/D^2 + (0.8713*d^8*n^2*phi^2)/D^2 - (3368.0*d^4*n*phi)/D + 2.167e7)^(1/2) + ((12655.0*d^4*n)/D + (3.311*d^8*n^2)/D^2 - (3.523*d^8*n^2*phi)/D^2 + (0.9374*d^8*n^2*phi^2)/D^2 - (6730.0*d^4*n*phi)/D + 2.417e7)^(1/2) + ((8835.0*d^4*n)/D + (2.817*d^8*n^2)/D^2 - (3.104*d^8*n^2*phi)/D^2 + (0.8551*d^8*n^2*phi^2)/D^2 - (4868.0*d^4*n*phi)/D + 1.39e7)^(1/2) + ((8346.0*d^4*n)/D + (3.603*d^8*n^2)/D^2 - (3.576*d^8*n^2*phi)/D^2 + (0.8872*d^8*n^2*phi^2)/D^2 - (4141.0*d^4*n*phi)/D + 2.329e7)^(1/2) + ((3.772*d^8*n^2)/D^2 - (7640.0*d^4*n)/D - (3.474*d^8*n^2*phi)/D^2 + (0.8001*d^8*n^2*phi^2)/D^2 + (3519.0*d^4*n*phi)/D + 1.327e7)^(1/2) + ((7532.0*d^4*n)/D + (3.51*d^8*n^2)/D^2 - (3.936*d^8*n^2*phi)/D^2 + (1.104*d^8*n^2*phi^2)/D^2 - (4223.0*d^4*n*phi)/D + 6.284e6)^(1/2) + ((3.719*d^8*n^2)/D^2 - (5771.0*d^4*n)/D - (3.448*d^8*n^2*phi)/D^2 + (0.799*d^8*n^2*phi^2)/D^2 + (2675.0*d^4*n*phi)/D + 1.318e7)^(1/2) + ((3.68*d^8*n^2)/D^2 - (3889.0*d^4*n)/D - (3.435*d^8*n^2*phi)/D^2 + (0.8015*d^8*n^2*phi^2)/D^2 + (1815.0*d^4*n*phi)/D + 1.354e7)^(1/2) + ((12933.0*d^4*n)/D + (3.382*d^8*n^2)/D^2 - (3.572*d^8*n^2*phi)/D^2 + (0.9431*d^8*n^2*phi^2)/D^2 - (6827.0*d^4*n*phi)/D + 2.55e7)^(1/2) + ((6165.0*d^4*n)/D + (2.554*d^8*n^2)/D^2 - (2.861*d^8*n^2*phi)/D^2 + (0.8014*d^8*n^2*phi^2)/D^2 - (3454.0*d^4*n*phi)/D + 8.277e6)^(1/2) + ((16822.0*d^4*n)/D + (5.572*d^8*n^2)/D^2 - (6.074*d^8*n^2*phi)/D^2 + (1.655*d^8*n^2*phi^2)/D^2 - (9165.0*d^4*n*phi)/D + 1.516e7)^(1/2) + ((9678.0*d^4*n)/D + (3.596*d^8*n^2)/D^2 - (3.603*d^8*n^2*phi)/D^2 + (0.9026*d^8*n^2*phi^2)/D^2 - (4848.0*d^4*n*phi)/D + 2.473e7)^(1/2) + ((7.461*d^8*n^2)/D^2 - (20733.0*d^4*n)/D - (7.427*d^8*n^2*phi)/D^2 + (1.848*d^8*n^2*phi^2)/D^2 + (10322.0*d^4*n*phi)/D + 2.428e7)^(1/2) + ((13000.0*d^4*n)/D + (3.443*d^8*n^2)/D^2 - (3.608*d^8*n^2*phi)/D^2 + (0.9451*d^8*n^2*phi^2)/D^2 - (6809.0*d^4*n*phi)/D + 2.647e7)^(1/2) + ((3.651*d^8*n^2)/D^2 - (2006.0*d^4*n)/D - (3.433*d^8*n^2*phi)/D^2 + (0.8073*d^8*n^2*phi^2)/D^2 + (943.6*d^4*n*phi)/D + 1.429e7)^(1/2) + (1.616e-65*abs(5.885e67*D - 1.31e64*d^4*n + 7.353e63*d^4*n*phi))/abs(D) + ((3.837*d^8*n^2)/D^2 - (9492.0*d^4*n)/D - (3.516*d^8*n^2*phi)/D^2 + (0.8052*d^8*n^2*phi^2)/D^2 + (4348.0*d^4*n*phi)/D + 1.383e7)^(1/2) + ((9530.0*d^4*n)/D + (2.895*d^8*n^2)/D^2 - (3.176*d^8*n^2*phi)/D^2 + (0.8709*d^8*n^2*phi^2)/D^2 - (5227.0*d^4*n*phi)/D + 1.548e7)^(1/2) + ((10822.0*d^4*n)/D + (3.584*d^8*n^2)/D^2 - (3.625*d^8*n^2*phi)/D^2 + (0.9166*d^8*n^2*phi^2)/D^2 - (5470.0*d^4*n*phi)/D + 2.589e7)^(1/2) + (2.964e-64*abs(3.6e66*D - 1.911e62*d^4*n + 1.078e62*d^4*n*phi))/abs(D) + ((4358.0*d^4*n)/D + (7.412*d^8*n^2)/D^2 - (7.68*d^8*n^2*phi)/D^2 + (1.989*d^8*n^2*phi^2)/D^2 - (2258.0*d^4*n*phi)/D + 6.172e6)^(1/2) + ((12833.0*d^4*n)/D + (3.494*d^8*n^2)/D^2 - (3.631*d^8*n^2*phi)/D^2 + (0.9432*d^8*n^2*phi^2)/D^2 - (6665.0*d^4*n*phi)/D + 2.701e7)^(1/2) + ((11733.0*d^4*n)/D + (3.564*d^8*n^2)/D^2 - (3.638*d^8*n^2*phi)/D^2 + (0.9285*d^8*n^2*phi^2)/D^2 - (5989.0*d^4*n*phi)/D + 2.67e7)^(1/2) + ((12411.0*d^4*n)/D + (3.535*d^8*n^2)/D^2 - (3.641*d^8*n^2*phi)/D^2 + (0.9376*d^8*n^2*phi^2)/D^2 - (6391.0*d^4*n*phi)/D + 2.708e7)^(1/2) + ((17888.0*d^4*n)/D + (6.014*d^8*n^2)/D^2 - (6.503*d^8*n^2*phi)/D^2 + (1.758*d^8*n^2*phi^2)/D^2 - (9669.0*d^4*n*phi)/D + 1.597e7)^(1/2) + ((3.918*d^8*n^2)/D^2 - (11322.0*d^4*n)/D - (3.573*d^8*n^2*phi)/D^2 + (0.8146*d^8*n^2*phi^2)/D^2 + (5163.0*d^4*n*phi)/D + 1.487e7)^(1/2) + ((4.595*d^8*n^2)/D^2 - (20433.0*d^4*n)/D - (4.16*d^8*n^2*phi)/D^2 + (0.9414*d^8*n^2*phi^2)/D^2 + (9245.0*d^4*n*phi)/D + 2.655e7)^(1/2) + ((17699.0*d^4*n)/D + (6.447*d^8*n^2)/D^2 - (6.909*d^8*n^2*phi)/D^2 + (1.851*d^8*n^2*phi^2)/D^2 - (9477.0*d^4*n*phi)/D + 1.512e7)^(1/2) + ((8561.0*d^4*n)/D + (3.759*d^8*n^2)/D^2 - (4.203*d^8*n^2*phi)/D^2 + (1.175*d^8*n^2*phi^2)/D^2 - (4786.0*d^4*n*phi)/D + 7.069e6)^(1/2) + ((7.553*d^8*n^2)/D^2 - (12799.0*d^4*n)/D - (7.622*d^8*n^2*phi)/D^2 + (1.923*d^8*n^2*phi^2)/D^2 + (6453.0*d^4*n*phi)/D + 1.407e7)^(1/2) + ((7580.0*d^4*n)/D + (2.682*d^8*n^2)/D^2 - (2.978*d^8*n^2*phi)/D^2 + (0.827*d^8*n^2*phi^2)/D^2 - (4209.0*d^4*n*phi)/D + 1.119e7)^(1/2) + ((4.795*d^8*n^2)/D^2 - (22344.0*d^4*n)/D - (4.35*d^8*n^2*phi)/D^2 + (0.9865*d^8*n^2*phi^2)/D^2 + (10144.0*d^4*n*phi)/D + 2.992e7)^(1/2) + ((7.54*d^8*n^2)/D^2 - (4003.0*d^4*n)/D - (7.712*d^8*n^2*phi)/D^2 + (1.972*d^8*n^2*phi^2)/D^2 + (2047.0*d^4*n*phi)/D + 7.599e6)^(1/2) + ((1709.0*d^4*n)/D + (3.618*d^8*n^2)/D^2 - (3.46*d^8*n^2*phi)/D^2 + (0.8274*d^8*n^2*phi^2)/D^2 - (817.3*d^4*n*phi)/D + 1.675e7)^(1/2) + ((15111.0*d^4*n)/D + (5.143*d^8*n^2)/D^2 - (5.644*d^8*n^2*phi)/D^2 + (1.549*d^8*n^2*phi^2)/D^2 - (8292.0*d^4*n*phi)/D + 1.343e7)^(1/2) + ((8180.0*d^4*n)/D + (2.745*d^8*n^2)/D^2 - (3.037*d^8*n^2*phi)/D^2 + (0.8403*d^8*n^2*phi^2)/D^2 - (4526.0*d^4*n*phi)/D + 1.246e7)^(1/2) + ((6.439*d^8*n^2)/D^2 - (32944.0*d^4*n)/D - (6.093*d^8*n^2*phi)/D^2 + (1.442*d^8*n^2*phi^2)/D^2 + (15588.0*d^4*n*phi)/D + 5.003e7)^(1/2) + (5.063e-62*abs(5.506e64*D + 1.962e61*d^4*n - 1.061e61*d^4*n*phi))/abs(D) + (4.92e-65*abs(5.951e67*D + 2.327e64*d^4*n - 1.247e64*d^4*n*phi))/abs(D) + (1.043e-63*abs(1.917e66*D + 5.873e62*d^4*n - 3.242e62*d^4*n*phi))/abs(D) + (2.195e-66*abs(2.288e68*D + 4.325e65*d^4*n - 2.284e65*d^4*n*phi))/abs(D) + (1.888e-61*abs(1.387e64*D - 9.999e60*d^4*n + 4.977e60*d^4*n*phi))/abs(D) + (3.992e-63*abs(2.323e65*D + 2.871e62*d^4*n - 1.478e62*d^4*n*phi))/abs(D) + (1.627e-61*abs(6.779e63*D + 8.254e60*d^4*n - 4.02e60*d^4*n*phi))/abs(D) + (4.493e-66*abs(2.719e68*D + 3.601e64*d^4*n - 2.027e64*d^4*n*phi))/abs(D) + (8.82e-65*abs(1.657e67*D + 3.926e63*d^4*n - 2.195e63*d^4*n*phi))/abs(D) + ((15588.0*d^4*n)/D + (6.843*d^8*n^2)/D^2 - (7.259*d^8*n^2*phi)/D^2 + (1.925*d^8*n^2*phi^2)/D^2 - (8262.0*d^4*n*phi)/D + 1.237e7)^(1/2) + (1.798e-64*abs(3.126e67*D - 9.98e63*d^4*n + 4.776e63*d^4*n*phi))/abs(D) + (1.119e-61*abs(8.594e63*D - 3.718e59*d^4*n + 1.681e59*d^4*n*phi))/abs(D) + (1.684e-62*abs(5.607e64*D + 7.784e61*d^4*n - 3.723e61*d^4*n*phi))/abs(D) + (2.55e-65*abs(6.374e67*D + 6.295e64*d^4*n - 3.261e64*d^4*n*phi))/abs(D) + (4.51e-63*abs(2.295e65*D + 2.959e62*d^4*n - 1.428e62*d^4*n*phi))/abs(D) + (6.076e-64*abs(1.885e66*D + 2.202e63*d^4*n - 1.082e63*d^4*n*phi))/abs(D) + (1.01e-62*abs(1.401e65*D - 2.529e61*d^4*n + 1.144e61*d^4*n*phi))/abs(D) + (5.558e-63*abs(1.034e66*D - 3.018e62*d^4*n + 1.428e62*d^4*n*phi))/abs(D) + (1.034e-63*abs(4.521e65*D + 1.022e63*d^4*n - 4.735e62*d^4*n*phi))/abs(D) + (1.619e-65*abs(5.615e67*D - 1.561e64*d^4*n + 8.745e63*d^4*n*phi))/abs(D) + (9.657e-64*abs(4.291e65*D + 3.639e62*d^4*n - 1.649e62*d^4*n*phi))/abs(D) + (3.174e-66*abs(1.31e68*D + 5.458e65*d^4*n - 2.791e65*d^4*n*phi))/abs(D) + (7.598e-65*abs(1.787e67*D + 3.654e63*d^4*n - 2.049e63*d^4*n*phi))/abs(D) + (3.899e-64*abs(1.832e66*D + 3.122e63*d^4*n - 1.468e63*d^4*n*phi))/abs(D) + (2.35e-64*abs(1.796e67*D - 4.857e63*d^4*n + 2.236e63*d^4*n*phi))/abs(D) + (2.95e-64*abs(9.082e66*D - 2.379e63*d^4*n + 1.083e63*d^4*n*phi))/abs(D) + (4.879e-64*abs(2.106e66*D + 2.464e63*d^4*n - 1.258e63*d^4*n*phi))/abs(D) + (2.749e-64*abs(4.165e66*D + 4.695e63*d^4*n - 2.352e63*d^4*n*phi))/abs(D) + (8.005e-63*abs(1.471e65*D + 1.399e61*d^4*n - 7.885e60*d^4*n*phi))/abs(D) + (6.813e-65*abs(1.702e67*D + 1.938e64*d^4*n - 9.614e63*d^4*n*phi))/abs(D) + (8.105e-65*abs(1.973e67*D + 5.231e63*d^4*n - 2.914e63*d^4*n*phi))/abs(D) + (6.491e-65*abs(1.566e67*D - 2.055e63*d^4*n + 1.157e63*d^4*n*phi))/abs(D) + (3.011e-64*abs(9.467e66*D + 4.325e63*d^4*n - 2.294e63*d^4*n*phi))/abs(D) + (4.08e-66*abs(6.402e67*D + 1.679e65*d^4*n - 7.654e64*d^4*n*phi))/abs(D) + (7.031e-65*abs(9.384e66*D + 1.448e64*d^4*n - 7.585e63*d^4*n*phi))/abs(D) + (1.023e-61*abs(8.166e63*D + 1.243e61*d^4*n - 5.895e60*d^4*n*phi))/abs(D) + (1.38e-62*abs(2.902e65*D - 1.379e62*d^4*n + 6.77e61*d^4*n*phi))/abs(D) + (6.923e-65*abs(7.826e64*D + 1.071e64*d^4*n - 5.775e63*d^4*n*phi))/abs(D) + (1.944e-63*abs(2.54e66*D - 6.913e62*d^4*n + 3.207e62*d^4*n*phi))/abs(D) + ((3.63*d^8*n^2)/D^2 - (134.7*d^4*n)/D - (3.443*d^8*n^2*phi)/D^2 + (0.8161*d^8*n^2*phi^2)/D^2 + (63.79*d^4*n*phi)/D + 1.538e7)^(1/2) + (1.703e-65*abs(1.727e67*D - 3.589e64*d^4*n + 1.958e64*d^4*n*phi))/abs(D) + (1.632e-65*abs(6.39e67*D - 5.826e63*d^4*n + 3.284e63*d^4*n*phi))/abs(D) + (6.87e-65*abs(2.186e66*D - 9.823e63*d^4*n + 5.33e63*d^4*n*phi))/abs(D) + (2.045e-63*abs(2.463e66*D - 9.149e62*d^4*n + 4.434e62*d^4*n*phi))/abs(D) + (1.036e-63*abs(5.682e65*D + 1.106e63*d^4*n - 5.161e62*d^4*n*phi))/abs(D) + (1.651e-63*abs(6.924e65*D + 4.005e61*d^4*n - 2.259e61*d^4*n*phi))/abs(D) + (8.23e-66*abs(3.607e67*D + 6.395e64*d^4*n - 2.905e64*d^4*n*phi))/abs(D) + (1.784e-62*abs(3.538e64*D + 9.056e60*d^4*n - 4.097e60*d^4*n*phi))/abs(D) + (2.515e-64*abs(2.18e67*D - 6.062e63*d^4*n + 2.839e63*d^4*n*phi))/abs(D) + (3.796e-65*abs(9.063e67*D - 2.438e64*d^4*n + 1.115e64*d^4*n*phi))/abs(D) + (1.107e-64*abs(1.156e67*D + 1.957e63*d^4*n - 1.099e63*d^4*n*phi))/abs(D) + (5.015e-65*abs(2.129e67*D - 3.648e64*d^4*n + 1.841e64*d^4*n*phi))/abs(D) + (5.823e-64*abs(6.238e65*D + 1.632e63*d^4*n - 7.518e62*d^4*n*phi))/abs(D) + (4.852e-65*abs(2.033e67*D - 3.546e63*d^4*n + 1.994e63*d^4*n*phi))/abs(D) + (2.81e-64*abs(2.859e66*D + 3.858e63*d^4*n - 2.004e63*d^4*n*phi))/abs(D) + (5.08e-65*abs(4.805e67*D + 2.871e64*d^4*n - 1.506e64*d^4*n*phi))/abs(D) + (4.5e-65*abs(5.024e67*D + 1.614e64*d^4*n - 8.854e63*d^4*n*phi))/abs(D) + (2.007e-64*abs(2.686e66*D - 2.458e63*d^4*n + 1.354e63*d^4*n*phi))/abs(D) + (7.593e-65*abs(5.583e66*D - 7.25e63*d^4*n + 3.976e63*d^4*n*phi))/abs(D) + (1.056e-65*abs(1.056e68*D + 2.204e63*d^4*n - 1.244e63*d^4*n*phi))/abs(D) + (7.184e-65*abs(2.775e67*D - 6.64e63*d^4*n + 3.012e63*d^4*n*phi))/abs(D) + (7.028e-63*abs(9.091e64*D - 6.251e61*d^4*n + 3.458e61*d^4*n*phi))/abs(D) + (9.735e-65*abs(2.97e66*D + 8.489e63*d^4*n - 3.888e63*d^4*n*phi))/abs(D) + (1.773e-61*abs(6.208e63*D + 7.055e60*d^4*n - 3.568e60*d^4*n*phi))/abs(D) + (6.533e-65*abs(1.22e67*D - 5.215e63*d^4*n + 2.906e63*d^4*n*phi))/abs(D) + (2.22e-64*abs(3.264e66*D - 1.75e63*d^4*n + 9.719e62*d^4*n*phi))/abs(D) + (5.205e-66*abs(2.095e68*D - 3.341e63*d^4*n + 1.886e63*d^4*n*phi))/abs(D) + (3.249e-65*abs(2.641e67*D - 9.097e63*d^4*n + 5.083e63*d^4*n*phi))/abs(D) + (2.791e-64*abs(9.097e66*D + 3.057e63*d^4*n - 1.666e63*d^4*n*phi))/abs(D) + (2.801e-64*abs(1.202e66*D + 3.14e63*d^4*n - 1.671e63*d^4*n*phi))/abs(D) + (3.357e-64*abs(5.296e66*D + 1.526e63*d^4*n - 8.467e62*d^4*n*phi))/abs(D) + (2.613e-64*abs(6.478e65*D + 3.099e63*d^4*n - 1.66e63*d^4*n*phi))/abs(D);
end

博大广阔

程序目的：性能函数由五个因变量。求性能函数的最小值下，对应的五个因变量的值。约束函数是非线性的：function [c ceq]=contion(XX)
%约束条件
d=XX(1);D=XX(2);n=XX(3);phi=XX(4);
c(1)=126.7-phi-3.466*D*(D-d)/(d*(4*D-d));
c(2)=phi-101.5-3.466*D*(D-d)/(d*(4*D-d));
ceq(1)=0;
ceq(2)=0;
end
求解思路：先进行随机最有搜索法，获得较优的初始值。而后用函数fmincon .求解最有值。