问题标签 [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 投票
3 回答
3825 浏览

cuda - 用 CUDA 求解二维扩散(热)方程

我正在尝试解决一些标准问题来学习 CUDA。例如,我正在使用以下代码求解二维扩散方程。但是我的结果与标准结果不同,我无法弄清楚。

有人可以告诉我这个过程是否有任何问题吗?任何帮助将不胜感激。

0 投票
1 回答
1370 浏览

matlab - 连续时间马尔可夫过程

求解 CT 马尔可夫过程的方法有哪些?

我知道对于诸如生死或某些队列之类的已知过程,可以分析解决问题吗?但是,如果无法解析可解,如何解决?

看来应该使用数值方法。但是,我还没有找到任何支持它的工具。MATLAB有ode方法,但是除了设置过程的(常微分方程)ode外,还要添加一个matlab不支持的线性方程?

如何解决这些问题?
有什么工具或其他方法吗?

0 投票
1 回答
1821 浏览

javascript - 为不同的网络屏幕尺寸加载不同的图像

我有这段代码,每次用户加载 index.html 时,它都会从数组中自动加载不同的图片。这是jquery代码:

我是新来的,可以将图像适应不同的屏幕分辨率。所以我认为如果有人打开我的网站,例如使用 2560/1440 像素的 imac,图像将使用此代码正确调整,但我想它会完全像素化。所以我认为,我必须创建一个更大的图像文件,以便那些计算机加载更大的文件以适应分辨率。我想避免其他屏幕正常的用户出于速度原因加载大文件。我可以在这段代码中添加什么来使更大的屏幕加载更大的文件,这样它就不会像素化?!?!

PD 如果您还知道不同屏幕尺寸组的最佳图像分辨率,那将非常有帮助!

谢谢!

0 投票
2 回答
777 浏览

numerical-methods - 使用 odeint 函数定义

很菜鸟的问题,所以请多多包涵。

我正在按照此处给出的示例-> http://www.codeproject.com/Articles/268589/odeint-v2-Solving-ordinary-differential-equations

特别是,我正在研究这个函数:

就我而言,R 采用一系列值(具有 100 个双精度值的向量)。

odeint 被称为:

我想为每个 R 值执行此操作。我该如何完成此操作?我对 C++/OOP 的了解有限,但我愿意学习。

谢谢你。

0 投票
1 回答
2128 浏览

numerical-methods - 数值求解两个耦合非线性二阶微分方程

我在拉格朗日力学中遇到了以下微分方程组。您能否建议一种数值方法,以及有关如何解决它的相关链接和参考。另外,在 Matlab 或 Mathematica 上是否有更短的实现?

mx (y 点)^2 + mgcosy - Mg - (M=m)(x 双点) =0

gsiny + 2(x 点)(y 点 + x(y 双点)=0

其中 (x 点) 或 (y 点)= dx/dt 或 dy/dt,双点表示 wrt 时间的双导数。

0 投票
2 回答
1646 浏览

wolfram-mathematica - 带参数的 NDSolve

我是 Mathematica 的新手,可能我正在寻找的东西很简单,我尝试了一些方法,但我无法成功。

我有一个 DE,我想要一个由 DE 参数参数化的解决方案。更准确地说,
让我们将 DE 定义为: x'=f[x,t,a]其中“a”是参数,我希望得到x作为时间和参数的函数的解,我将函数 g 命名如下。 x=g[t,a]

当 DE 是基本的时,它很容易使用DSolve和获取它,但是我找不到使用NDSolve. 我知道这NDSolve仅适用于数值争论,但我正在寻找一种方法来将结果与不同的参数值进行比较。
用不同的参数值显示/比较结果的简单方法也可以。

0 投票
4 回答
2417 浏览

haskell - Haskell - 优化微分方程求解器

我正在学习 Haskell,并尝试尽可能快地用 C 语言编写代码。对于这个练习,我正在为一个简单的一维物理系统编写欧拉积分器。

  • C 代码使用 GCC 4.5.4 和-O3. 它在1.166秒内运行。
  • Haskell 代码使用 GHC 7.4.1 和-O3. 它在21.3秒内运行。
  • 如果我用 编译 Haskell -O3 -fllvm,它会在4.022秒内运行。

那么,我是否缺少优化 Haskell 代码的内容?

PS.:我使用了以下论点:1e-8 5.

C代码:

哈斯克尔代码:

0 投票
0 回答
363 浏览

wolfram-mathematica - NDSolve 不返回任何内容

任何人都可以看看这个.nb,因为我很确定我的 7 个未知数和 3 个微分未知数的 10 个方程组是有效的?NDSolve 什么都不返回,它甚至不会给我一条错误消息。

链接到笔记本

0 投票
1 回答
1279 浏览

java - 如何将函数作为 Java 中 JTextField 的输入?

我正在尝试编写一个 java 程序,它采用两个方程 f(x) 和 g(x) 作为 JTextFields 中用户的输入。但是,我无法弄清楚我应该如何获取字符串输入并将其转换为可行的方程式。有人可以帮我吗?我实在想不通。

如果有人想知道,我正在尝试从用户输入中绘制欧拉法,但要做到这一点,我需要输入一个初始点和一个 h 值,这很容易,然后我需要这两个方程 f(x) 和g(x) 来自形成系统 F(x) 的用户输入。

谢谢。

0 投票
2 回答
4570 浏览

matlab - 如何在 Matlab 中求解这些耦合微分方程?

如图所示,我有三个偏微分方程 (PDE) 和一个变量的解析解。使用这些方程我想解决 \phi(x,y,t), p(x,y,t), C_{a}(x,y,t) 和 C_{b}(x,y,t)即在空间和时间方面。

我知道 Matlab 中有一个函数pdepe( )可以解决一维抛物线椭圆 PDE 的初始边界值问题。我想知道如何使用这个函数或 Matlab 中的其他函数来解决下面描述的二维和耦合问题。

问题:

以下两个方程分别代表两个物种 a 和 b 的 PDE: 在此处输入图像描述

其中 D_{h} 和 q 给出为:

在此处输入图像描述

在此处输入图像描述

这里,R_{a}=R_{b}=R,其中 R 为:

在此处输入图像描述

最后,最后一个方程为:

在此处输入图像描述

初始和边界条件:

在此处输入图像描述

总域大小为 10 cm x 5 cm,y 形子域的宽度为 0.5 cm。该子域的初始 phi 为 0.50,而在周围矩阵中则为 0.26。1 Pa 和 0 Pa 的常数 p 分别保持在边界 (1) 和 (2) 处,对应于大约 10^-3 mm^-1 的梯度。边界 (3) 和 (4) 上的 p 由边界 (1) 和 (2) 之间的线性梯度确定。C_{a} = 2 mol m^-3 和 C_{b} = 0.2302 mol m^-3 的常数 C 保持在边界 (3),而边界 (4) 的浓度设置为 C_{a} = 1 mol m^-3 和 C_{b} = 0.4603 mol m^-3。边界 (1) 处的浓度由边界 (3) 和 (4) 之间的恒定梯度确定,同时设置平流通量边界条件 $$(\frac{\partial C}{\partial x} = 0)$$在 (2) 处的出口处。