用隐式有限差分法实现一维扩散方程的类无穷边界条件

文章/答案/技术大牛

发布

1回答

matlab、numerical-methods

我使用下面的MATLAB代码，用隐式有限差分法实现了杆上的一维扩散方程。; t1imp = L\R; 我想把右边界条件改为无限边界条件，这样10m处边界的存在对10m区域内的热流和温度流

浏览 8提问于2017-07-24得票数 2

2回答

Fick第一扩散定律在模拟细胞膜水运动中的应用

python

我想用fick定律做上面的模拟，用有限差分法。我的问题是整个问题的开始。物理特性： 边界条件： C = sc.zeros(5) C[0Fick<em

浏览 4提问于2012-01-14得票数 1

回答已采纳

3回答

用Theano求解反应扩散系统

python、scipy、theano

我是提亚诺的新手，我试图实现一个反应扩散系统的数值积分器--这个版本的：目前，我的表达方式是：import theano.tensor asT.dscalar('e')a1 = T.dscalar('a1') dvdt = e*(u - a1*v - a0) 所以我还没有实现

浏览 7提问于2016-02-21得票数 8

回答已采纳

1回答

带黑匣子的伴随方法

openmdao

如果将多个求解器作为外部代码组件(完整的黑匣子)附加到模型上，是否有一种实现伴随方法的方法？因为它也不太可能实现一个复杂的步骤方法。如果要使用基于梯度的优化器，唯一的方法是有限差分？

浏览 0提问于2018-04-12得票数 1

回答已采纳

3回答

二维热传导方程的数值快速算法

python、python-3.x、numpy、heat

我正在寻找一种用python求解2D热方程的方法。我已经实现了有限差分法，但它是慢动作(进行100,000次模拟需要30分钟)。在我的代码的第一种形式中，我使用了有限差分的2D方法，我的栅格是5000x250 (x，y)。现在我想降低计算速度，我的想法是找到其中U是热函数。对于实现，我在2D情况下使

浏览 1提问于2015-05-31得票数 4

2回答

IIR滤波器是否使用卷积运算？

signal-processing

根据这个卷积运算的定义，我是否正确地假设它是而不是，通常用于IIR滤波器，因为IIR滤波器可以在右边引用y[]吗？

浏览 11提问于2013-02-07得票数 0

1回答

Drake中的自动微分

drake

我有一个关于Drake中自动区分(AD)的问题。如果有人用python为系统编写DAE，Drake如何计算用于优化的导数？它在后端是否有AD包，或者它使用有限差分。另外，它如何实现四元数表示的积分和微分。非常感谢您提前抽出时间。

浏览 3提问于2021-05-20得票数 1

2回答

关于数值方法的概念混淆

numerical-methods

我有一个关于数值方法的概念性问题。有限元法、连续有限元法、间断有限元法、连续galerkin法和间断galerkin法有什么区别？他们中有些人也是一样的吗？提前感谢

浏览 5提问于2016-02-04得票数 1

回答已采纳

1回答

Python中的显式(前向)和隐式(向后) Euler方法

python、jupyter、ode

我试图实现显式和隐式Euler方法，以近似于以下ODE的解决方案：欧拉方法用有限差分近似导数：前向方法基于先前的解x(t)显式计算x(t+dt)。x(t+dt) = x(t) + f(x，t)dt 向后方法是隐式<e

浏览 11提问于2022-01-20得票数 2

回答已采纳

1回答

有限域上椭圆曲线上随机点的生成

elliptic-curves、implementation

为了应用ECC算法，我编写了一个椭圆曲线的实现程序。然而，在大多数情况下，Alice需要选择给定曲线上的点P。选择这样一点的一般程序是什么？给出一个小的例子，例如y^2 = x^3 + x + 1 on F_{25}，有算法在曲线上生成随机点吗？在我的实现点上，这一领域的代表是多项式，如果这是相关的。

浏览 0提问于2019-04-07得票数 6

回答已采纳

1回答

解Laplace方程

code-golf、math、calculus

偏微分方程(偏微分方程)。拉普拉斯或扩散方程在物理学中经常出现，如热方程、变形、流体力学等。现实生活是3D的，但我们想说“你好，世界！”不唱"99瓶啤酒“这项任务是在1D中完成的。你可以把它解释为用某种力量绑在两端墙壁上的橡胶长袍。有限差分将-u_i''

浏览 0提问于2016-11-06得票数 15

回答已采纳

7回答

GPGPU上的金融应用

cuda、finance、gpgpu

我想知道使用GPGPU可以实现什么样的金融应用程序。我知道使用CUDA在GPGPU上使用蒙特卡洛模拟进行期权定价/股票价格估计。有人能列举出在金融领域使用GPGPU的各种可能性吗？

浏览 4提问于2010-05-11得票数 16

1回答

有效地表示FiPy中的“零阶”饱和汇，并且没有负浓度。

python、numerical-methods、pde、fipy

我试图用FiPy来求解一个简单的二维粒子扩散场景的瞬态行为，这个场景中有一个小的圆形水槽(粒子浓度= n)。我希望的接收器是可饱和的，也就是说，应该有一个最大的可能损失率内的接收器。用FiPy?表示这一点的最有效的方法是什么？"]) 在这种情况下，n(t)的求解效果很好，而且我可以在不牺牲太多精度的情况下采取大的<em

浏览 2提问于2022-08-22得票数 0

7回答

从字符串求解方程到C++

c++、c、string、linear-algebra、solver

我想知道是否有人有关于如何做一些听起来简单但看起来不像编程的事情的信息或经验。其想法是:给出一个包含等式的字符串，例如："2*x = 10“(这很简单，但它可能会变得非常复杂，如sqrt(54)*35=x^2；以此类推……)程序将返回x=5，并可能给出他是如何到达那里的日志。对于信息，有这个网站()，它用PHP做同样的事情，但不是开源的。谢谢你的帮助！

浏览 0提问于2010-05-17得票数 3

回答已采纳

1回答

含动力学反应的隐式差分热方程的Matlab解

arrays、matlab、for-loop、matrix、finite-element-analysis

我试图用隐式有限差分法来模拟木材圆柱体内的热传导。对于圆柱形和球形形状，我使用的一般热方程是：其中p为形状因子，p=1为圆柱，p=2为球面。边界条件包括表面对流。有关该模型的更多细节，请参见下面Matlab代码中的评论。主要的m文件是：rhow = 650; % density of wood, kg

浏览 7提问于2013-09-13得票数 1

回答已采纳

2回答

C++中的多项式类

c++、object-oriented

这是我朋友的Java OOP课堂作业之一。他们的任务是用Java编写一个具有乘法功能的多项式类。我在C++中这样做是一种挑战，也是为了提高我的编程技能。我还添加了一些额外的功能。else} main()函数只是作为一个测试，所以它不是最优的。让我知道你的</em

浏览 0提问于2020-03-22得票数 4

回答已采纳

4回答

为什么有限域在密码学中如此重要？

finite-field

我只是进入密码学，目前通过尝试实现一些密码算法来学习。我只是不明白为什么它们是相关的。我看到的一件事是，他们可以确保你的结果没有十进制，所以没有四舍五入的错误，但我强烈怀疑这就是它们如此重要的原因。如果有人能给我一个关于为什么需要他们的直觉，那就太好了。而且，它们是有限的，这是否会损害安全？

浏览 0提问于2021-06-19得票数 11

回答已采纳

10回答

C/C++中的导数？

c++、c、math、numerical、symbolic-math

我想做的一件事，就是取表达式的偏导数。所以如果是f(x,y) = x^2 + y^2，那么f相对于x的偏导数就是2x，相对于y的偏导数就是2y。我用有限差分方法写了一个小函数，但是我遇到了很多浮点精度的问题。例如，我最终使用1.99234而不是2。有没有支持符号区分的库？还有其他建议吗？

浏览 0提问于2009-01-07得票数 21

3回答

随机(Java 7)中的181783497276652981和8682522807148012是什么？

java、random

new Random()将使用当前的“种子统一器”，并使用System.nanoTime()进行XOR。一开始，这个评论给了我一个简单的线索。在线搜索这篇文章会产生。8682522807148012没有出现在报纸上，但是181783497276652981确实出现了--作为另一个数字的子字符串1181783497276652981，它是带有1前缀的181783497276652981本文认为1181783497276652981是一个对线性同余生成元产生好的“优点”的数字。这个数字只是错

浏览 3提问于2013-08-06得票数 114

回答已采纳

14回答

线性回归和逻辑回归的区别是什么？

machine-learning、data-mining、linear-regression

当我们必须预测 (或离散)结果的值时，我们使用。我相信我们也可以使用来预测给定输入值的结果的值。那么，这两种方法有什么不同呢？

浏览 4提问于2012-08-28得票数 269

回答已采纳

点击加载更多