numerical-methods 相关问题

通过数值逼近(与符号计算相反)解决数学问题的算法。

mesh_size的含义

下面的“mesh_size”是指每个单元的最大周长还是最大面积? 使用 pygmsh.geo.Geometry() 作为 geom: geom.add_polygon( [ [0.0, 0....

回答 1 投票 0

使用 CUDA 计算数百个小矩阵的特征值/特征向量

我有一个关于使用 CUDA 对数百个小矩阵进行特征分解的问题。 我需要计算数百个(例如 500 个)小(64×64)实对称的特征值和特征向量

回答 2 投票 0

为什么我们不能使用大界来替代 mosek 中的无穷大?

在文档中(https://docs.mosek.com/latest/pythonfusion/debugging-numerical.html) 我们看: 切勿使用非常大的数字来代替无穷大 。而是定义变量或

回答 1 投票 0

使用 NumPy 向量化 python 循环

由于我的数学建模的原因,我有一个 vector_U 函数 将 numpy 导入为 np def 矢量_U(U_0, t, func, dt): res = np.empty((len(t), 4)) 资源[0] = U_0 对于范围内的 i(1, len(t)): ...

回答 1 投票 0

使用一些 Python 库对正半定矩阵进行 Cholesky 分解

我正在寻找一个内置函数,可以在Python中实现正半定矩阵的Cholesky分解。 存在 NumPy (numpy.linalg.cholesky) 和 SciPy (sc...

回答 1 投票 0

方程解的错误图,通过Python中的RK4算法求解

我正在尝试制作一个火箭从太阳系的一个行星飞往另一个行星的二维模型。在模拟太空飞行之前,我必须制作行星本身的模型。所以,对于每个星球我...

回答 1 投票 0

编码二维阵列中多个谐振子的四阶龙格-库塔方法。陷入单一振荡器

我之前在 Python 中使用四阶龙格库塔方法对一般(非线性)摆摆进行了数值求解。这是与此相关的问题,尝试求解第二个非线性摆......

回答 1 投票 0

如何使用 python 脚本对一组两个二阶双变量偏微分方程进行数值近似解(边界值问题)?

我有一组两个偏微分方程,如下所示(伪代码): r'' = (s')^2 * r - G*M*r^-2 r' = -(r*s'')/(2*s') s 和 r 是时间的函数,' 的导数,G 和 M 是常数。 它们是从...

回答 1 投票 0

尝试使用四阶龙格-库塔法求解非线性摆二阶微分方程,未得到预期结果

我是Python新手,所以我只将过程限制在matplotlib上,而不是进入NumPy。 我按照 Abhijit Kar Gupta 的书《Python 中的科学计算》,编写了 Python

回答 1 投票 0

有没有什么方法或者公式可以计算欧拉法和龙格-库塔4的误差

使用简单欧拉方法对步骤h1的系统进行建模时,得到误差a1,而4阶龙格-库塔方法得到误差b1。 我怎样才能确定误差值...

回答 1 投票 0

使用python的牛顿法

p= 拉姆达 x: 3 * x**3 - 5.03 * x**2 - 1.95 * x + 0.02 dp = 拉姆达 x: 9 * x**2 - 10.06 * x - 1.95 x0 = -1 公差 = 10**-4 最大迭代数 = 100 def newton_method(f, df, x0, tol, max_iter=100): 对于我...

回答 1 投票 0

在 R 中迭代处理数值积分误差的最佳方法?

考虑 R 中的以下数值积分: 被积函数 = 函数(u){log(u) * u ^ (2 * (H - 1)) * (exp(a * u) - exp(- a * u))} 积分(f = 被积函数,下限 = 0,上限 = 100, 相对于...

回答 1 投票 0

如何使用solve_ivp用谱法求解偏微分方程?

我想用谱法求解偏微分方程。这样的方程,公式,初始条件为u(t=0,x)=(a^2)*sech(x),u'_t (t=0)=0。 为了解决这个问题,我使用了 pyt...

回答 2 投票 0

尝试进行verlet集成

我尝试过调整欧拉方法代码,并使用欧拉来计算第一个值,所以我可以在verlet中使用两个值,但是当我绘制图形时,我只得到两条垂直的直线...

回答 1 投票 0

Maple 2023plot() 函数对于复合函数的工作非常奇怪。我是不是做错了什么?

我正在尝试制作一个在 1 和 0 之间振荡的周期函数。我制作了一个在间隔上按照我需要的方式振荡的函数 <0,7), when i tried to make it periodical i wrote a mod...

回答 1 投票 0

Python 2.7 的用户输入功能存在问题,无法让这些输入与程序一起使用

我对 python 世界很陌生,尽管我已经用 php 做了很多工作......这是我的案例...... 我正在用 python 2.7 为我的小程序编写一些代码。 在该程序中,我需要接受 2 个用户输入,机器人...

回答 4 投票 0

求解矩阵的最小二乘而不是向量

问题是找到 Z 使得 epsilon(E)(方程 21)最小化。 Z 是一个 MxN 矩阵,这就是我们需要找到的。 Zx和Zy也是已知的MxN矩阵。 Dx 和 Dy 是...

回答 1 投票 0

Haskell 中缺乏有限元方法实现 - 有什么具体原因吗?

我很好奇为什么 Haskell 或任何函数式语言中似乎缺乏有限元方法 (FEM) 实现。鉴于 Haskell 的纯函数性质,我预计......

回答 1 投票 0

我如何对以下表达式进行数值积分?

我目前正在尝试在Python中计算以下表达式 该表达式来自 t 分布,v 是表示自由度的标量,mu 和 sigma 是向量和矩阵

回答 1 投票 0

修复Python代码以实现雅可比方法

我正在尝试编写代码来在线性方程组上实现雅可比方法。 到目前为止我已经写过这个: #雅可比方法 我=0 z = np.array([[-1.,3.,2.,11.],[1.,1.,1.,6.],[0...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.