odeint

0熱度

1回答

我想使用的代碼 tmax,dt=5000,1.68 t=linspace(0,tmax,num=round(tmax/dt)+1) def d_displ(x,t,a,b,c,d,e): #x'=y, y'=a*x+(b/(x**2.0+c))*x+d*y+e y=x[0] dy=x[1] dx=[dy,a*x+(b/(x**2.0+c))*x+d*y

0熱度

1回答

犰狳的cx_mat和Boost的odeint編譯錯誤

我打算解決幾個矩陣微分方程，形式爲d/dt (X) = F(X)，其中X是一個大的複數矩陣，F表示它的一些函數。我試圖用Boost的odeint和state_type作爲Armadillo的cx_mat。但它會爲受控步進器類型生成編譯錯誤。我的示例代碼如下 #include <armadillo> #include <iostream> #include <boost/numeric/odei

1熱度

1回答

從odeint scipy使用的函數提取值python

我有以下腳本來使用odeint來計算dRho。 P_r = 10e5 rho_r = 900 L = 750 H = 10 W = 150 A = H * W V = A * L fi = 0.17 k = 1.2e-13 c = 12.8e-9 mu = 2e-3 N = 50 dV = V/N dx = L/N P_in = P_r rho_in = rho

0熱度

1回答

無法使odeINT接受由循環生成的方程

我想使用Python的odeINT來集成由循環生成的多組方程。方程式都是耦合的，所以必須通過一次調用odeINT同時進行集成。問題是初始條件（「y0」）必須是列表或矩陣列表（而不僅僅是列表）。 odeINT給出了這個錯誤：「初始條件y0必須是一維的」。我想知道如何解決這個問題。這是一個代碼示例;非常感謝任何想法。 class network: def __init__(self):

-3熱度

1回答

我是否需要安裝Boost來構建odeint？

我想通過使用odeint包來解決ODE問題。 odeint網站提供他們的軟件包下載。我下載了包含在我的項目中的zip文件，但它不起作用。據我所知，如果我下載整個boost軟件包，它會自動包含odeint。我已經成功地做到了這一點。但我的問題是「我可以直接從odeint網站下載odeint文件，並從那裏建立它，而無需安裝整個boost庫？」在我們的情況下，我們沒有足夠的時間讓學生在課堂上下載和

1熱度

2回答

如何使在matplotlib

兩個滑塊我想使兩個滑塊在matplotlib手動更改在我的捕食模型N和P的值： import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint def lotka(x,t,params): N, P = x alpha, beta, gamma, delta =

0熱度

1回答

錯誤：索引2超出軸1的大小1的範圍

我使用Python（odeint）解決激光速率方程（常微分方程第一階）的以下小程序存在問題。當我運行該程序，錯誤總是顯示出來： index 2 is out of bounds for axis 0 with size 1 的公式是正確的;我不知道這個錯誤是什麼意思。任何人都可以解釋這是什麼意思，以及如何解決它？ def Rate(y,t): D = y[0] P =

1熱度

2回答

描述化學方程式的Scipy ODE錯誤

我有一大組ODE代表生物系統中的化學通量。分子正在反應，隔離和循環。我試圖讓這個功能以這樣一種方式發揮作用，讓我知道在一系列條件下可以生產出多少某些產品。我使用這些包 import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt import math import pylab

0熱度

1回答

boost odeint：自定義類和向量空間代數的受控步進器

我使用boost odeint來演化微分方程。時間和值由一個包裝類在一個普通數字（做高斯錯誤傳播）給出。 using Time = Number<double>; using Value = Number<double>; 的狀態和它的衍生物是由自定義類給定的與聲明 class State : boost::addable<State>, boost::multipliable<State

0熱度

1回答

Python-使用兩個方程的odeint

我試圖繪製一段時間（1800秒）的溫度。我有兩個初始值問題，其中func2依賴於func1和兩個初始值。這裏有兩個功能： rate = dA/dt = -(3.083e8*np.exp(-56000/(8.314*Temp))*A*0.033) dT/dt = (-0.45*-98000*rate+5.7431*(273.15-Temp))/(2018.94) 其中A是物質的濃度和溫度是