我正在尝试用有限差分法在python中编写2D Schrödinger方程的曲线图

文章/答案/技术大牛

发布

1回答

python、physics

我已经开始编写一个程序，用有限差分法计算二维Schrödingers方程的解。我想要使用等高线图或一些其他图形显示，通过接受用户输入的尺寸和网格点数量，以图形方式显示解决方案。我将hbar^2/2m设置为1，并将电位(V)设置为0，从而简化了Schrödinger方程，从而得出： -(dψ^2/dx^2 + dψ^2/

浏览 16提问于2020-04-27得票数 0

1回答

如何有效地对稀疏矩阵进行lambdify？

python、scipy、sparse-matrix、sympy、lambdify

我正在使用渐近来构建一个稀疏(N^2 x N^2)矩阵，并尝试将其转换为稀疏scipy矩阵。它被用于求解二维网格上依赖于k向量k= (kx，ky)的Schrödinger方程的有限差分方法。SciPy的稠密矩阵我目前的方法是使用lambdify从渐近矩阵A创建函数B，该矩阵为某个k向量提供有限差分矩阵，我将其转换为稀疏矩阵M。

浏览 52提问于2020-06-23得票数 0

1回答

欧拉-拉格朗日方程数值解的算法建议

python、numerical-methods、pde

在我的课程作业中，我被告知要用Python编写一个算法来求解具有Dirichlet边界条件的欧拉-拉格朗日方程。非常感谢!

浏览 0提问于2018-11-17得票数 0

4回答

用Python进行简单的物理模拟？

python、modeling、simulation

我想知道类似的，具体的模拟，作为模拟浇水一块田地。我知道小Simpy，Numpy和Pygame。我想要得到关于它们的例子。

浏览 3提问于2009-02-02得票数 8

回答已采纳

1回答

如何在pyhthon中初始化2D不连续的初始数据

python、numerical-methods、pde、piecewise

我刚刚用Python编写了单位平方热方程的有限差分求解器，为了检查方程的一个重要性质，我需要给出一个不连续的函数作为初始数据。我的域在[0,1] x [0,1]中，我希望有一个如下所示的函数，但在2D中，即变量(x，y)的分段函数。

浏览 23提问于2019-11-23得票数 2

3回答

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

python、python-3.x、numpy、heat

我正在寻找一种用python求解2D热方程的方法。我已经实现了有限差分法，但它是慢动作(进行100,000次模拟需要30分钟)。在我的代码的第一种形式中，我使用了有限差分的2D方法，我的栅格是5000x2

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

4回答

用Python进行波浪模拟

python、numpy、simulation、physics、wave

我想要模拟三维空间中一些物体上具有吸收和反射的传播波。我想用python来实现。我应该使用numpy吗？有没有一些我应该使用的特殊的库？提前谢谢。

浏览 3提问于2011-02-10得票数 10

回答已采纳

4回答

初学者指南自己的CFD代码-二维欧拉方程

math、fluid-dynamics

你知道一个好的，特别容易的指南来编写自己的计算流体动力学解算器，用于二维欧拉方程吗？我只是想了解像Fluent这样的商业软件在做什么。当它足够简单时，我想向一些朋友展示如何做到这一点并编写代码。不幸的是，我无法找到如何将此转换为数字应用程序。安德烈亚斯

浏览 0提问于2011-10-27得票数 0

2回答

波动方程二维数值解的作图问题

matlab、numerical-methods、pde

我正在尝试用MATLAB (R2020a)制作二维波动方程的动画。到目前为止，我相信有限差分法是正确的。但是，当我尝试使用matlab中的"surf“命令绘制这些值时，它不起作用。到目前为止，这是我的脚本： clearclose alltMax = 30; xM

浏览 22提问于2021-02-04得票数 2

回答已采纳

1回答

python反问题的拉普拉斯平滑

python、optimization、matrix-multiplication

我目前正在研究一个用python编写的反问题，在这个问题中，观测到的GPS数据被用来推断地球以下的位移(例如地震)。地球以下的位移可以在网格中离散，描述这个问题的方程可以写成： G(n*n) * S(n*1) = U(n*1) [1] 其中G是描述问题的模型矩阵，S是地球向量以下的位移，U是表面观测值向量。N是在反演中使用的</e

浏览 22提问于2019-05-19得票数 0

回答已采纳

1回答

PARFOR在有限差分格式中的应用

matlab、parallel-processing、pde、parfor

我正在尝试用有限差分法求解一个非线性扩散方程。它本质上是两个嵌套的for循环，内部的是网格上的循环，外部的是时间推进方程，有点像这样： for j=2:(Nx-1) %%% simple diffusion请注意，为了清楚起见，我在这里使用了简单的扩散方程。问题是，

浏览 0提问于2015-05-27得票数 0

1回答

Lapack程序dstev中的故障分割

c++、vector、numerical-methods、lapack

我目前正在编写一个小程序，用C++中的有限差分法求解微分方程。问题是，我正在使用lapack来解决我的三对角矩阵，但是在例程被调用之后，我得到了一个分段错误。这个例程仍然存在于info=0中，所以我不知道问题出在哪里。这是节目。ev.begin()->begin()), &vdim, &*

浏览 0提问于2017-03-05得票数 1

回答已采纳

1回答

python中任意网格偏导数的有限差分估计

python、performance、mathematics、numerical-methods

我目前正在用有限差分测试python上的一些代码来求解微分方程，但我正在处理的一些问题和方程如下：因此，为了避免错误和避免穷尽的手工选择方程，我决定建立一个函数，

浏览 0提问于2023-03-06得票数 2

回答已采纳

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回答

pyplot生成两个不同的图形？

numpy、multidimensional-array、matplotlib

我正在尝试用中心有限差分法求解波动方程。我从一个基本的二维高斯函数开始作为我的初始条件u[x,y,t] = u[:,:,0] = init_fn(x,y)def init_fn(x,y):下

浏览 1提问于2014-11-30得票数 0

1回答

如何使用我在python中使用swig生成的包装器函数？

python、c、swig、numerical

好吧，我是个新手。我终于用swig和numpy.i成功地包装了我的python程序中最昂贵的部分。程序是二维波偏微分方程组的有限差分格式。我的问题是，我现在如何使用它？在IPython中导入后我可以看到它。：$ gcc -c -fpic wave2.c wave

浏览 4提问于2014-11-16得票数 0

2回答

如何利用SymPy加速长函数的符号导数？

python、optimization、sympy、integral、derivative

我正在用Python编写一个程序，用自由ICI方法求解Schr dinger方程(现在，SICI方法.但自由ICI是它将变成的)。如果这听起来不太熟悉，那是因为关于这个主题的信息很少，而且绝对没有可以使用的示例代码。这个过程需要迭代地得到一个偏微分方程的解。在这样做的过程中，有许多象征性的衍生工具需要执行。问题是，当程序运行时

浏览 0提问于2019-05-16得票数 9

回答已采纳

1回答

使用Mathematica的Poisson求解器

wolfram-mathematica、solver、poisson

我正在寻找一些帮助的泊松解决方案，我正在编写的数学。插入数组的代码相当长，但是详细信息可以在上找到。我用从泊松方程导出的中心差分法计算给定电荷密度的电压。在计算了电压后，我测试了数据集的收敛性。我将收敛阈值设置为10^-1000+。我将循环设置为在

浏览 2提问于2012-10-09得票数 2

1回答

我如何修改我的图形，以便它显示正确的信息为轴？

python、matplotlib、plot

我编写了一个用有限差分法数值求解热方程(u_t = k * u_xx)的程序。对于我的问题，u是x和t的函数，其中0 < x < L和t > 0。我已经为我的问题指定了L = 1 (杆的长度)和终端时间T = 10秒，所以我希望在域(x,t) \in {(0,1) x (0, 10)}上显示该图形。不过，我<

浏览 3提问于2016-07-12得票数 0

回答已采纳

6回答

Python:多维数组的迭代是不是太慢了？

python、performance、numpy

我必须在二维整数数组中迭代所有项，并更改值(根据某些规则，不重要)。a = numpy.ndarray((5000,5000), dtype = numpy.int32) x = 123

浏览 3提问于2012-04-11得票数 6

回答已采纳

点击加载更多