有誰知道為什么我的python代碼中出現這個錯誤以及如何糾正它？

當我運行下面的代碼時，出現此錯誤：

yr = solve_ivp(sec_lotes,t_sim,V_inicial[0],V_inicial[1],V_inicial[2],V_inicial[3],method='RK45')
TypeError: solve_ivp() got multiple values for argument 'method' 

這是代碼：

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
#When I run this code I get this error
#yr=
#solve_ivp(sec_lotes,t_sim,V_inicial[0],V_inicial[1],V_inicial[2],V_inicial[3],method='RK45')
#TypeError: solve_ivp() got multiple values for argument 'method' 
###Parámetros y constantes del modelo
kg = (0.00626 * 60) / ((1000) / (60 ** 2))  # coeficiente de transferencia de masa en 
kg/m2*s*kpa
pws = 12.352 * ((1000) / (60 ** 2))  # a 50°C presión de vapor de saturación del agua en kpa
pwinf = 2.33 * ((1000) / (60 ** 2))  # presión de vapor del agua en kpa
h = 0.0718 * 60  # coeficiente de transferencia de calor convectivo en kw/m2*°C
pa = 1.092  # densidad del aire en kg/m3
ps = 530  # densidad del sólido en kg/m3
Cpha = 1.006 * ((1000) / (60 ** 2))  # calor del aire húmedo en kj/kg*°C
k = 0.0005641 * 60  # conductividad térmica de la papa en kw/m*°C
Lc = 0.05  # espesor de la papa en m
delta_Hw = 2382 * ((1000) / (60 ** 2))  # calor de vaporización en kj/kg
V = 0.00000392  # volumen real en m3
V1o = 0.000549  # volumen de lecho inicial en m3
avo = 200  # área de la partícula/ volumen inicial (rodaja de papa) en m-1
c1 = 0.0267  # constante adimensional
c3 = -1.656  # constante adimensional
q1 = 0.0107  # constante adimensional
q2 = 1.287  # constante adimensional
q3 = -1.513  # constante adimensional
A = 1.612  # constante adimensional
B = 3.114  # constante adimensional
C = 2.175  # constante adimensional
D = -3.194  # constante adimensional
E = 3.661  # constante adimensional
F = -1.661  # constante adimensional
G = 0.339  # constante adimensional
H = 1.246  # constante adimensional
J = -1.385  # constante adimensional
K = 0.792  # constante adimensional
Bio = 0.291
Gs = 0.0183 * 60  # flujo de aire kg/s
Lo = 0.035  # longitud del lecho inicial en m
S = 0.07625  # área de sección transversal
pso = 500  # densidad del sólido inicial
eo = 0.9  # porisidad inicial
# nodos de ***delta_A***
M = 100
# FRONTERAS : 1s=inicial     2s=final
# Condiciones de frontera de y
Y0 = 0.01578
# Condiciones de frontera de Tg
Tg0 = 50

# Condición inicial de x
x0 = 1
# Condición inicial de Ts
Ts0 = 25
# *******************************************************************************
#*********************************Ecuaciones adicionales*********************************

def sec_lotes(t,Y,x,Ts,Tg):
# *********************************Ecuaciones derivadas*********************************
    for i in range(0, M):
        plt.clf()
        delta_A = i / M

        if i == 0:
            Y = Y0
            Tg= Tg0
            x=x0
            Ts=Ts0
        else:
            ##Ecuaciones adicionales
            aw = np.exp(-c1 * x[i] ** (c3))   np.exp(q1 * np.exp(-q2 * x[i]) * x[i] ** (q3) * 
            np.log(pws))
            fi = 0.00532 * ((x[i] / x0) ** -1.079)
            nw = (kg * (aw * pws * (Ts) - pwinf)) / ((1   (fi / x0) * Bio))
            delta_Hs = (1   q1 * np.exp(-q2 * x[i]) * x[i] ** (q3)) * delta_Hw
            Cphs = A   B * (x[i] / (1   x[i]))
            av = (C   D * (x[i] / x0)   E * (x[i] / x0) ** (2)   F * (x[i] / x0) ** (3)) * avo
            e = 1 - V / ((G   H * (x[i] / x0)   J * (x[i] / x0) ** (2)   K * (x[i] / x0) ** (3)) * V1o)
        ###############################
            dYdt = ((nw * aw * (1 - e)) / (pa * e)) - (1 / (S * Lo)) * (Gs / (pa * e)) * ((ps * (1 - e)) / (pso * (1 - eo))) * ((Y[i] - Y[i - 1]) / delta_A)
            dxdt = -((nw * av) / ps)

            dTsdt= ((av / (ps * (1   x[i]) * Cphs)) * (h * (Tg - Ts) - nw * delta_Hs))

            dTgdt = -(((h * av * (1 - e)) / (pa * e * Cpha)) * (Tg - Ts)) - (1 / S * Lo) * (Gs / (pa * e)) * (ps * (1 - e) / pso * (1 - e)) * ((Tg[i] - Tg[i - 1]) / delta_A)

            sol = [dYdt, dxdt, dTsdt, dTgdt]

            return sol


Y0 = 0.01578  # Humedad del sólido
x0 = 1.0  # Humedad del aire
Ts0 = 20.0  # Temperatura inicial del sólido en °C
Tg0 = 50.0  # Temperatura inicial del aire en °C
V_inicial = [Y0, x0, Ts0, Tg0]
# Intervalo de tiempo para la simualción (min)
t_sim = np.linspace(0, 1000, 5)  # 10000,10
# Integración del modelo

yr = solve_ivp(sec_lotes,t_sim,V_inicial[0],V_inicial[1],V_inicial[2],V_inicial[3],method='RK45')

uj5u.com熱心網友回復：

V_inicial[0],V_inicial[1],V_inicial[2],V_inicial[3] 應該是一個陣列，我從來沒有使用過這個函式，但我一直在尋找你可以在這里閱讀的檔案：

https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.solve_ivp.html

該函式將 V_inicial[0] 標識為 y0 引數，將 V_inicial[1] 標識為方法引數，當您撰寫 method='RK45' 時，您正在為方法引數定義另一個值。嘗試使用陣列，并告訴我 :))。

yr = solve_ivp(sec_lotes,t_sim,{V_inicial[0],V_inicial[1],V_inicial[2],V_inicial[3]},method='RK45')

標籤：Python 蟒蛇-3.x scipy

上一篇：猜謎游戲不會停止

下一篇：我們如何在資料框中的列中的字串周圍添加引號