ode

    1熱度

    1回答

    我想用sympy來求解一個複雜的微分方程組。我使用sympy快速計算時間導數,之後我有兩個包含衍生物的衍生物方程。這些方程不是線性的,也不符合sympy識別的方程式,所以它會拋出未實現的錯誤。有沒有更簡單的方法來解決這些方程(甚至是數值),並獲得它們各自的運動定律(隨時間變化的值)?這可能是非常低效的,所以如果有人知道更高效的流程,我就會全神貫注,我主要開始使用sympy,因爲計算衍生工具相當快(

    -1熱度

    1回答

    由於之前不清楚,我對此表示歉意。我現在更瞭解這個功能,但可以在幾個方面使用一些幫助。 我想找回轉換(X)與音量(V)的關係,或者反過來也可以。在我看來,傳統的「時代」術語是我想用0-1的X序列取代的,X是轉換,記得如此界定0和1.0 下面,rw是反應速率,並且是任何給定時刻的分壓的函數,它們被描述爲Pw,Px,Py和Pz,它們本身是初始條件(P.w0,v.0)的函數,轉換也是X. 預先感謝您 rm

    0熱度

    1回答

    我想運行求解&將我的ode解決方案繪製爲以下三個不同的時間,每次只改變一個初始條件。我想解決q.0 = 16,q.0 = 20,q.0 = 18 (這將是很好,如果這都繪製在矩陣..!在此先感謝) rm(list = ls()) par(mfrow = c(1,1)) D = 6 ; A = pi* D^2/4 ; # ODE solver library(deSolve) p

    0熱度

    1回答

    我想在我的ode45函數中使用一個變量值,該函數在另一個函數中進行計算。爲此,我使用global來訪問所有函數中的變量,但不以這種方式工作。我的主要功能代碼如下: function window_sine %%%%%%% Here is the global variable global vgth; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    2熱度

    2回答

    我目前有一個與時間有關的常數系統。例如。 def fun(u, t, a, b, c): x = u[0] y = u[1] z = u[2] dx_dt = a * x + y * z dy_dt = b * (y-z) dz_dt = -x*y+c*y-z return [dx_dt, dy_dt, dz_dt] 常數是

    1熱度

    1回答

    我試圖編譯庫, http://www.tbi.univie.ac.at/~raim/odeSolver/ 這一個特別,我已經得到了這個必要的庫。 但是,無論何時我嘗試使用「make」進行編譯,我都會收到以下錯誤消息。 Making all in src make[1]: Entering directory `/home/user/SOSLib/buildSOSlib/src' /usr/bi

    0熱度

    1回答

    我想找到具有2個自由度的簡單動力系統的響應,由ODE定義:m ddx(t) + k ddx(t) = f(t)。該系統被激勵力f1 = f01 * sin(omega*t)在點1: import numpy as np m1 = 10 m2 = 8 k1 = 1e6 k2 = 3e6 k3 = 7e6 m = np.diag([m1,m2]) k = np.array([[k1

    -5熱度

    1回答

    The following system of ODEs can be written in the format dx/dt = Ax where A is a matrix. What is A? dx1/dt = 4(x2-x1) dx2/dt = -2*x1 - x2 - 9*x3 dx3/dt = x2 + 5*x3 輸入以下格式的答案: [0 0 0 0 0 0 0 0

    2熱度

    1回答

    我試着用scipy.integrate.ode解決一個1階的ode方程。 DH/DX = - S 0 *(1 - (HN /小時)^ 3)/(1 - (HC /小時)^ 3) 的初始條件爲x = 0和h = 10。 S0 = 0.001當x < 15000, S0 = 0.0005當x> = 15000。 HN =(F * q^2/8 * G * S0)^(1/3) HC =( (1/3) f,q

    0熱度

    1回答

    ode的主要功能如下所示。 function dxdt = state(t,x,vgth,vgval1,vgval2) vgval=vgval1+vgval2; p=1; k=10^0.7; window1=1-((2*x)-1).^(2*p); dxdt=k*(vgval-vgth+1.2)*window1; end 腳本如下。 step=0.01; t = 0:step:10