所以我有一個(gè)問題，此時(shí)我有點(diǎn)迷失。因此，任何意見都將不勝感激，因?yàn)槲椰F(xiàn)在真的很掙扎！我有一個(gè)模型，想使用我得到的一些實(shí)驗(yàn)數(shù)據(jù)來檢查/優(yōu)化。一般來說，我的模型需要兩個(gè)輸入（我們稱之為：時(shí)間和溫度）并有 8 個(gè)變量 (x0-x7)。該模型生成兩個(gè)輸出（out1 和 out2）。我的每組實(shí)驗(yàn)數(shù)據(jù)都為我提供了 4 組可用于優(yōu)化的信息：2 個(gè)輸入（時(shí)間和溫度）和 2 個(gè)實(shí)驗(yàn)結(jié)果（結(jié)果 1 和結(jié)果 2）。最終我想最小化 result1 和 out1 以及 result2 和 out2 之間的差異。因此，基本上用幾組數(shù)據(jù)最小化兩個(gè)殘差，這些數(shù)據(jù)受 8 個(gè)共同參數(shù) (x0-x7) 的影響。我對(duì)參數(shù) x0-x7 有一些限制，這可以有所幫助，但除此之外沒有真正的限制。到目前為止，我已經(jīng)嘗試使用scipy.minimize我的實(shí)驗(yàn)結(jié)果數(shù)據(jù)集進(jìn)行迭代，如下所示（非常示意性）：import numpy as npfrom scipy.optimize import minimize    Experiment=[['Set 1','Set 2',             'Set 3','Set 4'],                   [Out 1-1,Out 1-2,                    Out 1-3,Out 1-4],                   [Out 2-1,Out 2-2,                    Out 2-3,Out 2-4],            ]global curr_casecurr_case=0 #just for debugging in the first place    def objective_fcn(x):            SetFitParameters(x) #x0-x7            #---------probably totally dumb: iteration-----------    global curr_case    #number of experimental set        curr_case=curr_case+1    if curr_case==len(Experiment):        curr_case=0    #----------------------------------------------------            getTemp(curr_case) # function that gets time and temperature from experimental data as two arrays - time and temperature            RefVariables(x) #sets some global variabales needed for ModelCal using x0-x7            ModelCal(time,Temperature)  #gives Out1 and Out2            f1 = abs(Out1[Upper_index-1]-Experiment[1][curr_case]) #compares Out1 with result1 (from experimental data)    f2 = abs(Out2[Upper_index-1]-Experiment[2][curr_case]) #compares Out2 with result2 (from experimental data)            # some weighting factors for the future - maybe?    A=1    B=1           return A*f1+B*f2       bounds_x1=(1450,1700) #upper and lower bonds of x0bounds_x2=(0.1,1)bounds_x3=(1450,1700)bounds_x4=(0.1,7)bounds_x5=(1450,1700)bounds_x6=(0.1,7)bounds_x7=(1450,1700)bounds_x8=(0.1,7)