在Python中优化函数的多个输出变量

Question

我目前正在研究一种算法来确定风力涡轮机支撑结构的成本。我写的算法需要优化初始输入支撑结构的重量，这样应力水平不会超过，而是接近所用材料性能的失效标准。另一个要求是结构的固有频率需要限定在2个值之间。为了优化结构，可以改变4个变量。在

我可以使用斯皮比。优化使用几个设计参数优化结构重量的库，但同时考虑了支撑结构的固有频率和最大应力值？在

我正在优化的函数如下所示：

def func(self, x):
    self.properties.D_mp = x[0]                    # Set a new diameter for the monopile
    self.properties.Dtrat_tower = x[1]             # Set a new thickness ratio for the tower
    self.properties.Dtrat_tp = x[2]                # Set a new thickness ratio for the transition piece  
    self.properties.Dtrat_mud = x[3]               # Set a new thickness ratio for the mudline region of the monopile

    self.UpdateAll()                               # Update the support structure based on the changes in variables above

    eig = self.GetEigenFrequency()                 # Get the natural frequency
    maxUtil = self.GetMaximumUtilisationFactor()   # Get the maximum utilisation ratio on the structure (more than 1 means stress is higher than maximum allowed)

    # Natural frequency of 0.25 and utilisation ratio of 1 are ideal
    # Create some penalty...
    penalty = (100000 * abs((eig - 0.25)))
    penalty +=  (100000 * abs(maxUtil - 1)) 

    return self.GetTotalMass() + penalty

提前谢谢！在

Answer 1

您可以使用的leatsq函数辛辣。优化. 在

在我的例子中，它是用两个变量拟合指数函数：

def func_exp(p, x, z):                                     
# exponential function with multiple parameters  
    a, b, c, d, t, t2 = p[0], p[1], p[2], p[3], p[4], p[5]
    return a*np.exp(b + pow(x,c)*t + pow(z,d)*t2)

但是要使用leatsq函数，你需要创建一个error函数，这个函数需要优化。在

^{pr2}$

要使用它：

p0=[1e4,-1e-3,1,1,-1e-2, -1e-6]                           
# First parameters 
pfit_fin, pcov, infodict, errmsg, success = leastsq(err, p0, args=(X,Y,Z), full_output=1, epsfcn=0.000001)

因此，返回最佳参数，以生成结果：

Y_2= func_exp(pfit_fin, X,Y)        

我希望这对你有帮助

克里斯。在