標簽[differentialequations.jl]

[英]Exclamation mark in Julia "!" after function in scripts

我正在嘗試將我的 ODE 的多個解保存在一個數組中。現在這就是我得到的：如您所見，如果使用特定的 ode 求解器，我只想將 append 添加到此數組。但是，“追加”。語句忽略此語句並運行循環的每次迭代，我嘗試預分配數組 sol: 使用這樣的語句： sols[i] = solution但這里 ...

julia微分方程解如何優化

[英]How to optimize the julia differential equation solution

我正在運行以下 Julia 代碼來求解一組微分方程，但它需要太多時間，有時我的機器會被壓碎。有什么解決方案可以減少運行時間和 memory 使用率嗎？我不太熟悉 Julia。因此非常感謝您的幫助。我只能用 RK4 方法解決它。 ...

Julia 中的解耦 ODE

[英]Decoupled ODEs in Julia

Lorenz 吸引子是 Julia 文檔中的原型示例，用於使用DifferentialEquations.jl求解 ODE 系統。假設我想針對一千種不同的初始條件求解該系統，同時保持吸引子參數固定。什么是解決和存儲這 1000 個解耦 ODE 數據的好（即有效）方法？ ...

向 DifferentialEquations.jl 提供稀疏模式時出現 BoundsError

[英]BoundsError when providing sparsity pattern to DifferentialEquations.jl

我正在嘗試使用 DifferentialEquations.jl 求解一組剛性常微分方程。當我只提供帶或不帶 Jacobian 的 ODE 函數時，它工作正常，但是當我嘗試使用jac_prototype關鍵字提供稀疏模式時，它總是導致 BoundsError。在為稀疏模式提供示例 Jacobi ...

Distributions.jl 中的“期望”命令最近是否不起作用？

[英]Does the command "expectation" in Distributions.jl not work recently?

我嘗試了以下命令當然，第二個命令是不正確的，並且仍然缺少許多參數，但它仍然適用於 nWME（不是工作最小示例）。我收到的錯誤消息是“錯誤：UndefVarError：未定義期望”。我是否缺少其他一些軟件包來運行“期望”命令？請幫忙。謝謝你。 ...

使用 DiffEqFlux 的 Sciml 沖突

[英]Sciml conflict using DiffEqFlux

我在使用 sciml 訓練 DiffEqFlux 教程時遇到錯誤。去年夏天，我使用了完全相同的代碼並進行了一些修改，沒有任何問題，但似乎找不到解決方案。我正在使用 Julia v.1.7.1 錯誤是：我不確定它是否是我正在使用的特定版本的 Julia。這直接來自https://doc ...

Julia BVP 求解器中邊界條件代碼的令人費解的結果

[英]Puzzling result from boundary condition code in Julia BVP solver

我正在嘗試按照此處找到的示例，使用 BoundaryValueDiffEq 包來解決 Julia 中的邊界值問題。在邊界條件函數中，示例需要一個 for 循環來單獨更新每個索引，à la 我想使用以下代碼，它應該更有效：盡管兩個版本的代碼的殘差結果值相同，但求解器在第一種情況下給出了正確 ...

在 Julia 微分方程中用向量求解 IVP

[英]Solving IVPs with vectors in Julia DifferentialEquations

我對朱莉婭很陌生，所以提前道歉。我正在嘗試同時解決 Julia 中的兩個微分方程並且遇到了麻煩。我在 Python 中有一個 MWE，即：這可以工作並且表現出我想要的方式。我嘗試在 Julia 中執行以下操作（同樣的問題），這給我帶來了麻煩：無論我如何編輯它，我都會在 Julia ...

DifferentialEquations.jl 多線程

[英]DifferentialEquations.jl multithreading

下面的代碼正在求解具有不同變量的 ODE，而PA s 對應於不同的Δ s，這是我在每次循環后唯一感興趣的數據。由於循環是獨立的，我假設我可以對循環進行多線程處理，所以我在for循環前面添加了Threads.@threads ：但它會導致julia終端崩潰（使用vscode），似乎與GR有關。 ...

Julia 微分方程抑制檢測到的不穩定性的警告

[英]Julia Differential Equations suppress warning of detected instabilities

我有一個程序，它使用 Julia 的微分方程 package 來模擬粒子的路徑。模擬允許粒子撞擊設備 - 為了防止繼續模擬此類粒子，我使用了求解器的unstable_check檢查（特別是EulerHeun求解器）。但是，這會導致如下警告：當我模擬數千個粒子時，這可能會很煩人（而且很慢）。我 ...

有沒有更快的方法來求解具有多維數組的微分方程

[英]Is there any faster method to solve differential equations with multidimensional array

我想要解決一個復雜的網絡系統，該系統涉及通過多維數組起作用的高階交互項。我已經編寫了相應的代碼，但是獲得結果需要花費太多時間。有什么方法可以更快地求解微分方程嗎？代碼如下： #= 我發現的問題是在我的代碼中計算廣義訂單參數需要太多時間。我是 Julia 的新手。誰能向我建議一種更快的方法 ...

如何在 Julia 中評估和獲取神經網絡的導數

[英]How can I evaluate and take the derivative of a neural net in Julia

我已經用神經網絡求解了一個微分方程。我在下面留下代碼示例。我希望能夠計算該神經網絡相對於其輸入“x”的一階導數，並針對任何“x”評估該導數。 1-請注意，我計算der = discretize.derivative 。這是神經網絡對“x”的導數嗎？使用這個表達式，如果我輸入[first(d ...

訓練神經微分方程時如何保證收斂？

[英]How to guarantee convergence when training a neural differential equation?

我目前正在完成 Julia 語言的 SciML 教程研討會練習（ https://tutorials.sciml.ai/html/exercises/01-workshop_exercises.ZFC35FDC70D5FC67D23E8C ）具體來說，我堅持練習 6 第 3 部分，其中涉及訓練神 ...

FFTW.jl 用於 2D 陣列：擴散只發生在 1D

[英]FFTW.jl for 2D array: Diffusion only happening in 1D

根據我的閱讀，當 A 是 2D 數組時，使用 FFTW.jl / AbstractFFTs.jl 的 fft(A) 應該在 2D 中執行 fft，而不是按列執行。知道為什么當（我認為）我將縮放的二階空間導數添加到 u(t,x) 時，我只看到列方向的擴散，就好像使用顯式求解器來計算時間一樣？謝謝 ...

解決 Julia 中的 DDE 的問題

[英]Issue with solving a DDE in Julia

我正在嘗試使用 Julia 的DifferentialEquations.jl package 來解決 DDE 系統。我能夠使用下面的代碼解決兩個對象的問題當我試圖使我的代碼（ MultiDelayedBuckle.jl ）可擴展以解決任意數量的對象時，我遇到了一個我無法理解的錯誤。下面你可以 ...

Julia 中的二階 ODE 給出錯誤結果

[英]Second order ODE in Julia giving wrong results

我正在嘗試使用 julia 提供的 DifferentialEquations.jl，在我嘗試在二階 ODE 上使用它之前它工作正常。例如考慮二階 ODE x''(t) = x'(t) + 2* x(t)，帶初始條件 x'(0) = 0, x(0) = 1 它的解析解為：x(t) = 2/3 e ...

.= Julia 中的運算符

[英].= operator in Julia

在代碼中取自： https://tutorials.sciml.ai/html/models/01-classical_physics.html ，如下所示：我不明白 function 諧波振盪器中 .= 運算符的用法。使用 = 給了我錯誤的答案。所以，我想知道.=與=有何不同？它不是向量化 ...

Julia：何時在微分方程系統中鍵入真或假？

[英]Julia: When to type true or false in DifferentialEquations systems?

我從DiffEqFinancial.jl的源代碼中汲取靈感來正確設置我自己的系統。我false true目的：那么{true}和{false}在這種情況下做了什么？我想知道我是否也應該在我自己的用戶代碼中使用它。我隱約意識到“嚴格類型類型，松散類型功能”的格言，所以我想知道為什么要在這里打 ...

具有離散強迫的非齊次 ODE function 在 Julia 中使用微分方程 package

[英]Nonhomogeneous ODE with discrete forcing function using DifferentialEquations package in Julia

在Julia中，我想用DifferentialEquations.jl package來解決 \ddot{u} + f(u,\dot{u},p) = g(t) 其中 g(t) 是在時間 t 的等距離瞬間給出的值向量。這種情況與https://diffeq.sciml.ai/stable/tuto ...

是否可以使用回調來訪問 Julia 的微分方程集合問題中的單個軌跡？

[英]Is it possible to use callbacks to access a single trajectory in Julia's DifferentialEquations Ensemble Problems?

我是 Julia 的新手，並嘗試使用 Julia package 微分方程同時求解同一組耦合 ODE 的多個條件。我的系統是一個實驗的 model，在其中一種情況下，我在過程的中途增加了一個因變量的數量。我希望能夠調整這個單一軌跡的條件，但是到目前為止我只能一次調整所有軌跡。是否可以使用回調 ...