通过数值逼近（与符号计算相反）解决数学问题的算法。
我试图创建一个计算Givens旋转QR分解的函数，按照这个伪代码。 function [Q,R] = givens(A) [m,n] = size(A); indexI = zeros(m,n); indexJ = zeros(m,n); C ...
计算(sin(x)-x)*x^{-3}的优化算法(matlab中)
我的任务是写出计算矩阵Y的最优程序，给定矩阵X，其中：y=(sin(x)-x)*x^{-3}。这是我目前写的代码： n = size(X, 1); m = size(X, 2); Y = zeros(n, m); d ...
不同步长大小的欧拉法 。如何改变算法的代码，以考虑步长大小的不同值？
我有一个数值微分方程问题的算法，叫做欧拉法。本质上，欧拉法是对微分方程的近似解。我的函数工作在一个...
如何实现log_softmax()，以更好的速度和数值稳定性计算其值（和梯度）？
MXNet 和 PyTorch 都提供了计算 log(softmax())的特殊实现，速度更快，数值更稳定。但是，我找不到这个Python的实际实现 ...
旨在求解这个耦合微分方程系统：$ frac{dx}{dt} = -y $ $/frac{dy}{dt} = x $ 按照下面的隐式演化方案：$ y(t_{n+1}) = y(t_{n}) + \frac{/Delta t}......
我试图应用scipy的solve_bvp来解决以下问题 T'''''(z) = -k^4 * T(z) 在大小为l的域和一些常数A的边界条件下： T(0) = T'''(0) = T''''(l) = 0 T'(l) = A 所以 ...
我试图用欧拉方法找到微分方程在某点的解，但我得到的答案是错误的。请看我下面的代码。#初始条件x0，y0=0，1 h=0. 1 #步长... ...
我试着比较正弦函数的解析导数和近似导数，用后向、前向、中心、3阶后向差分法，它们之间有巨大的差异。它...
我想在C++中实现第一种和第二种复数贝塞尔函数的描述。现在我正在寻找在我的源代码中引入它们的可能性。...
在这个特殊的问题上，我有一个想象中的城市被划分为若干个方块--基本上是一个覆盖城市的MxN方格。M和N可以比较大，所以我有超过40,000平方的案例......。
如何改进Levenberg-Marquardt的多项式曲线拟合方法？
几周前，我开始在Matlab中从头开始编码Levenberg-Marquardt算法。我对数据的多项式拟合很感兴趣，但我还不能达到准确度的水平......
我正在尝试从本文中实现方程式16：上面的方程式16是粘性波动方程式，应该从大开始并随着时间的流逝逐渐消失，最终接近0。但是...
我正在尝试在Python中创建泰勒级数，但我不知道如何消除此错误：x = Symbol（“ x”）def f（x）：return（（math.e）** x） * sin（x）y = f（x）print（y.diff（x））def Taylor（f，x，m，a）：y = ...
[如何在python中创建函数以使它们返回x作为变量？ （插值）
[如何在Python中创建函数，例如：从数学import *从numpy import *从sympy import *从numpy.linalg import * n = int（input（“结数：”））xsolmed = []对于i in ...
我有一个RKF7（8）积分器，我已经用几个简单的测试功能验证了其输出。但是，当我在感兴趣的运动方程式上使用它时，局部截断误差为...
如何使用scipy.brentq函数，使其等效于pascal代码？
我正在尝试进行一些软件考古。尝试用python中的现代脚本替换旧的pascal文件。在某些时候，pascal脚本正在使用如下所示的例程：FUNCTION ...
我正在实现RKF4（5）集成器，但我无法弄清楚我的代码是否有效，并且我不了解本地截断错误，或者我的代码是否无效。我为代码的大小表示歉意...
供参考，我正在使用此页面。我了解原始的rankrank公式，但是我无法理解为什么稀疏矩阵实现是正确的。下面是它们的代码复制：def ...
为什么我的Runge-Kutta解算器无法随着时间步长的减小而收敛？
出于某种原因，我需要在PyTorch中实现Runge-Kutta4方法（因此，不，我将不使用scipy.odeint）。我尝试过，但在最简单的测试用例上却得到了怪异的结果，用x（0）= 1（...
我正在尝试解决问题。运行此代码时，我实现了想要的功能。（。8372e-03 1.8372e-06 1.8372e-09）truncation.error