您需要創(chuàng)建一個函數(shù)來一次性處理輸入、積分和繪圖 ( sir_interactive_func),見下文:# For integration.import scipy.integrate # For arrays (Python does not have native arrays).import numpy as np# For graphing.import matplotlib.pyplot as plt # Prevents the pop-up graphs in a separate window.get_ipython().run_line_magic('matplotlib', 'inline')# Allows for an interactive widget bar.from ipywidgets import interactive S0 = 0.95I0 = 0.05R0 = 0.0def SIR_model(y, t, beta, gamma): S, I, R = y dS_dt = -beta*S*I dI_dt = beta*S*I - gamma*I dR_dt = gamma*I return([dS_dt, dI_dt, dR_dt,]) def sir_interactive_func(beta, gamma): # Graph from 0 to 100, include 10000 points. t = np.linspace(0, 100, 10000) solution = scipy.integrate.odeint(SIR_model, [S0, I0, R0], t, args=(beta, gamma)) solution = np.array(solution) plt.figure(figsize=[8, 5]) plt.plot(t, solution[:, 0], label="S(t)") plt.plot(t, solution[:, 1], label="I(t)") plt.plot(t, solution[:, 2], label="R(t)") plt.grid() plt.legend() plt.title("SIR Model") plt.xlabel("Time") plt.ylabel("Proportions of Populations") interactive_plot = interactive(sir_interactive_func, beta=(0.35,1,0.01), gamma=(0.1,1,0.01))interactive_plot
在python中使用子進程找不到文件錯誤
函數(shù)式編程
2023-02-22 16:04:32