numerical-methods 相关问题

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

Python-用LU分解求解线性方程组的数值问题

信息操作系统:Manjaro Linux开发环境:Spyder4(Python 3.7)库:笨拙的问题您好,我按照以下三种方法编写了一些函数来求解线性方程组:LU ...

回答 1 投票 0

将具有两个输入的数值函数映射为一个具有一个输入的数值函数

我在编程方面很烂,所以请多多包涵。我什至不确定我现在需要的概念是什么,所以我不知道要用谷歌搜索或写这篇文章的标题。我的问题是,...

回答 1 投票 0

如何在我的代码中产生此溢出?

因此,我试图使用Bulirsch-Stoer自适应方法来解决EDO,但是却收到了溢出警告。如何解决呢?从数学导入sin,cos,pi从numpy导入范围,数组,空,...

回答 1 投票 0

在Matlab中加速程序

我有2个函数:ccexpan-用于计算函数f的插值多项 式与第一种基础的Chebyshew多项式中的N个节点的系数。 csum-为...

回答 2 投票 1

计算增量均值/平均值对精度是否有益?

关于“计算平均值的最佳数字方法是什么”的建议是,计算滚动平均值,即平均值= a [n] / n +(n-1)/ n *平均值可能在数值上更稳定...

回答 1 投票 0

如何将这个情节变成动画?我有一个情节,但想以某种方式将其变成动画

在最近的几个晚上,我一直在努力如何在每个时间步长或每个x步长后将此处的波形转换为某种动画。如何修改和编写代码,以便...

回答 1 投票 0

四阶非线性微分方程的稳定解

我已经在python中使用bvp求解器解决了以下bvp问题。 d4y / dx4 = 0.00033 * V /(0.000001-y)^(2),y(0)= y'(0)= y(1)= y'(1)= 0在上面的等式'V'是一个使用...

回答 1 投票 0

IndexError是什么原因:列表索引超出范围?

因此,我的教授要求我创建一个程序,该程序将使用Runge-Kutta方法(二阶)来解决该问题。但是我得到一个IndexError:列表索引超出范围,这可以回到我的第一个...

回答 2 投票 0

您如何在Python中实现基于Householder的QR分解?

我目前正在尝试对矩形矩阵实施基于Householder的QR分解,如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf(第3、4、5页)所述。显然...

回答 1 投票 2

如何在Python中为矩形矩阵实现基于Householder的QR分解?

我目前正在尝试对矩形矩阵实施基于Householder的QR分解,如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf(第3、4、5页)所述。显然...

回答 1 投票 2

PDE解决:'NoneType'没有属性'toarray'错误

我正在学习如何在Python中采用有限差分方法,以泊松方程为例。但是,当我运行时,一切都很好,直到A = ASparse.toarray()逆向代码行...

回答 1 投票 0

单面Jacobi SVD和Divide&Conquer SVD稳定性和误差

我正在MATLAB中研究SVD,尤其是Jacobi SVD和Divide&Conquer SVD算法。在这些方法的稳定性和错误分析上,我找不到任何东西。有人可以显示...

回答 1 投票 0

为隐式Runge-Kutta方法编写函数(第四阶)

我正在尝试使用第4阶隐式Runge-Kutta方法(IRK)编写一个能够解决ODES系统的函数,但是我无法正确定义循环。在这里,我们通过...

回答 1 投票 2

我如何求解带有空气阻力的弹丸运动的二阶微分方程?

等式为:d ^ 2 r / dt ^ 2 = -c / m(dr / dt)+ g,其中r是弹丸的位置,c是阻力系数,m是弹丸的质量,g是重力引起的加速度。假设...

回答 2 投票 -1

方程组-scipy的解决方法给出的错误答案

我正在尝试求解以下耦合方程:x = 1; y-0.5 * y-0.7 * v = 0; w-0.7 * x-0.5 * x = 0; v =1。(我知道方程= 1似乎是不必要的,但以后我需要它们...

回答 1 投票 0

我如何在PYTHON中求解带有空气阻力的弹丸运动的二阶微分方程?

等式为:d ^ 2 r / dt ^ 2 = -c / m(dr / dt)+ g,其中r是弹丸的位置,c是阻力系数,m是弹丸的质量,g是重力引起的加速度。假设...

回答 1 投票 -1

高精度计算导数时除以零误差

我正在使用斜率方法计算函数的导数。我需要放大图表以查看实际情况,但是数据点太低而无法给出平滑的曲线。所以当我...

回答 1 投票 0

Matlab中的函数fzero不收敛

我正在解决我的宏观经济学课上的问题。考虑以下方程式:在此,k是固定的,而c(k)是通过Matlab中的``interp1''函数定义的。这是我的代码:beta = 0.98; ...

回答 1 投票 1

公差标准布伦特方法

Brent方法的停止标准是,如果abs(m)<= tol或fb == 0.0,则//找到根(间隔足够小)发现:= true;但是,如果abs(m)达到上述容差以下,但...

回答 1 投票 0

使用fsolve和Euler 3BDF的SIR模型

嗨,我被要求使用MATLAB中的fsolve命令来求解SIR模型,而Euler 3向后指向。我对如何继续感到很困惑,请帮忙。这是我到目前为止所拥有的。我创建了一个函数...

回答 1 投票 1

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