问题标签 [differential-equations]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
2 回答
2323 浏览

wolfram-mathematica - 在 Mathematica 的两个独立绘图命令中使用微分方程的解

我在尝试在两个单独的绘图命令中使用NDSolve的答案时遇到了问题。为了说明这个问题,我将使用一个简单的微分方程和一个绘图命令。如果我写这样的东西:

它求解方程并计算 x[4] 没有问题,但情节变空,我不知道为什么。

在我的实际问题中,我的方程是一个包含多个函数的相当复杂的系统,而不是x[4],我绘制了已求解函数的参数图。我最终打算将所有这些都包含在Manipulate语句中,因此我不希望NDSolve语句出现多次(花费太长时间)并且我不能只提前计算它(因为它有很多参数)。


编辑:我想澄清和扩展我的问题:我真正想做的是通过以下方式将我的绘图语句包含在Manipulate语句中:

因为只有Manipulate语句给参数a赋值,所以我不能事先计算NDSolve的答案。另外,由于我的实际方程系统非常复杂且非线性,我不能使用符号函数DSolve

抱歉,如果之前不清楚。

0 投票
4 回答
8718 浏览

c# - 使用 C# 求解偏微分方程

我正在做一个项目(C# 和 .NET Framework),它需要我解决一些偏微分方程。是否有任何基于 .NET Framework 的特定库可供我查看并使我的工作更简单?

我使用过 MATLAb,在那里求解偏微分方程非常简单。我怎么解决这个问题?

0 投票
1 回答
627 浏览

maple - 在微分方程组中分配

当我在数值上求解一个由两个微分方程组成的系统时:

关于“t”,则“t=4”的解(例如)为 sol(4) 形式:

现在,如何将 n1(t) 和 n2(t) 的值用于另一个方程中的所有“t”,即“p”,其中涉及 n1(t) 或 n2(t)(例如:{p =a+n1(t)*n2(t)+f(t)},其中定义了“a”和“f(t)”),并为“t”的间隔绘制“p”?

0 投票
1 回答
7201 浏览

matlab - 如何用 MATLAB 求解微分方程

我有个问题。那是 :

我想用 MATLAB 解决这个问题,但我不能。

你能帮助我吗 ?

0 投票
3 回答
820 浏览

python - 在Python中并行化微分方程的分辨率

我正在使用odeint函数求解一个常微分方程组。是否有可能(如果是的话如何)轻松地并行化这类问题?

0 投票
3 回答
4718 浏览

math - 帮助辛积分器

我正在尝试开发物理模拟,并且我想实现四阶辛积分方法。问题是我一定弄错了数学,因为我的模拟在使用辛积分器时根本不起作用(与四阶龙格-库塔积分器相比,它对模拟工作得相当好)。我一直在谷歌上搜索这个,我能找到的都是关于这个主题的科学文章。我试图调整文章中使用的方法,但没有运气。我想知道是否有人有使用辛积分器的模拟源代码,最好是模拟引力场,但任何辛积分器都可以。源是什么语言并不重要,但我会欣赏使用 C 风格语法的语言。谢谢!

0 投票
1 回答
550 浏览

matlab - 使用 matlab 求解 ODE

我有以下 ODE:

其中 k15 和 k16 是常数。

关于如何解决它的任何想法?

谢谢!阿米特

0 投票
2 回答
1265 浏览

math - 抛物线偏微分方程

我的应用程序有一些抛物线偏微分方程......它们是相互关联的,并使用用户通过桌面应用程序的 UI 输入的一些变量。
您能否指导我了解哪种软件或库或特定语言最适合上述目的?

0 投票
2 回答
5219 浏览

animation - Mathematica 2D 热方程动画

我正在绘制二维的温度梯度图,遇到了很多麻烦。我目前的方法是定义一个插值函数,然后尝试多次绘制它,然后为该图表设置动画。这是我到目前为止所拥有的:

所以前两个变量只是控制变化率。我正在求解的方程是基本的二维热方程,其中 dT/dt=a(d^2T/dx^2+d^2T/dy^2)。初始条件将所有内容设置为 0,然后将边缘定义为热变化的来源。现在它从 t=0 扫描到 t=6 的 9x9 块。

第二部分尝试为功能工作设置动画。

不幸的是,这不起作用,我疯了想找出原因。我最初认为它与插值函数有关,但现在我对动画代码的工作也不太有信心。有人有想法么?

0 投票
1 回答
16289 浏览

java - Java中的微分方程

我正在尝试在java中创建一个简单的SIR-epidemics模型模拟程序。

基本上,SIR 由三个微分方程组定义:
S'(t) = - l(t) * S(t)
I'(t) = l(t) * S(t) - g(t) * I(t)
R'(t) = g(t) * I(t)

S - 易感人群,I - 感染者,R - 康复者。

l(t) = [c * x * I(t)] / N(T)

c - 接触人数,x - 传染性(与病人接触后生病的概率),N(t) - 总人口(恒定)。

如何在 Java 中求解这样的微分方程?我不认为我知道任何有用的方法来做到这一点,所以我的实现产生了垃圾。

}

我将不胜感激任何帮助,非常感谢提前!