通过数值逼近(与符号计算相反)解决数学问题的算法。
我为此方法生成的代码使用x的行空间中的N + 1,以便包含初始数据。绘制时,它给出的解决方案接近于实际情况,但是,绘制错误时,...
我正在使用Mupad计算机代数系统来尝试计算称为A的3 X 3矩阵的特征向量。A:= matrix([[3,2,1],[1,3,2],[0 ,2,6]]);似乎有两种方式...
我一直在尝试构造矩阵Dij,如该图片所示,Dmatrix我想将其绘制为间隔[-1,1上xi = -cos [pi *(2i + 1)/ 2N]上的点]以因此采用...
QuyRk in numpy-使用反射器在QR分解中将所有归一化元素分配为0
[数组中所有归一化的元素都设置为0。我正在尝试为数组创建反射器。当我通过第一个元素(x [1:] / = x [0],x [0] = 1)归一化数组中的所有元素时...
我选择了Julia来做一些数值分析工作,并试图实现完整的枢轴LU分解(例如,试图获得尽可能稳定的LU分解)。我以为...
我正在使用interpolate.interp2d()将2D样条拟合到函数上。我如何获得花键w.r.t.的一阶导数每个因变量?到目前为止,这是我的代码,Z是...
我们在我的数值分析类中使用python进行编码,对此我还很陌生,所以我不太确定我做对了吗,感谢您的帮助。使用python,我需要计算:i)0.01 ^ 2和...
我正在尝试使用for循环在python中编写Simpson规则,但我不断收到断言错误,无法找到原因。 def integration_numeric(xmin,xmax,N):'''f的数值积分...
这对我来说似乎是一个显而易见的问题,但是我在SO的任何地方都找不到。我有一个三次多项式,我需要找到该函数的实根。这样做的方式是什么?我发现了...
R问题:寻找最快的方法以数值方式求解一堆任意立方,这些立方具有有效系数和三个真实根。据报道,R中的多根函数使用Jenkins-Traub的...
错误来自尝试从scipy.integrate导入Radau方法(之所以需要,是因为Oregonator模型是一个刚性系统)。我正在尝试以数字方式整合Oregonator模型...
我正在尝试开发Matlab程序来平衡化学方程式。我可以通过求解线性方程组来平衡它们。目前,我的输出是带有...
我想知道迭代求解器是否是求解线性系统(非稀疏,对称,正定)的更快方法。我尝试了R包Rlinsolve和cPCG中的共轭梯度方法,但是两者都...
我正在尝试通过家用反射器实施QR分解。在一个非常简单的数组上进行尝试时,我得到了奇怪的数字。任何可以告诉我的人,为什么还要使用@ vs * ...
带有Runge Kutta的Lotka Volterra达不到所需的精度
人口随时间推移(应该是每个峰的高度都相同,我已经编写了一个代码,使用Runge-Kutta 4阶模拟老鼠和狐狸的种群。但是结果却不尽人意。.每个...] >
使用python包scipy,可以使用“柯西”加权方法找到函数的原理值(假设极点为低阶,请参阅scipy.integrate.quad(例如考虑...
我是MPI的新手,并将其作为大学课程学习。任务是使用MPI_Send()和MPI_Recv()从数字上找到const e的值。我找到的唯一合适的方法是在2、3 ...