通过数值逼近(与符号计算相反)解决数学问题的算法。
我正在尝试使用最初用C语言编写的Python复制DSP算法。诀窍是,我还需要保留C版本中32位定点变量的相同行为,...
我尝试了以下代码,使用增量搜索技术找到方程式的根。当我将精度设置为5或7位数字时,代码可以正常工作。但是当我将精度设置为6 ...
用scipy Solve_bvp解决边界值问题(扩散反应方程)
我正在努力解决以下二阶边值问题:y''+ 2 / x * y'+ k ** 2.0 * F(y)= 0 y(x = 1)= 1,y'( x = 0)= 0 F(y)= -y或F(y)= -y * exp(AB *(1-y)/(1 + B(1-y)))我以某种方式无法设置。 。
因此该方法本身可以正常工作:double trap_method(double a,double b,double n){int i = 1;两倍的面积= 0;双倍h =(b-a)/(n-1);双x = a + h; while(i <= n-2){area = ...
我正在尝试实现中点公式y [n + 1/2] = y [n] + h / 2 f(x [n],y [n])和y [n + 1] = y [n] + h * f(x [n] + h / 2,y [n + 1/2]),因此可以使用中点法求解ODE。我的函数是函数[...
我正在尝试使用角谱方法在数值上传播给定的(电场)。为此,我遵循“傅立叶光学的原理和应用”(Robert K. Tyson)第3章,...
我正在尝试实现光谱聚类算法,以解决图问题中的社区。我有一个非常大的矩阵来计算其特征向量,矩阵> 1Mx1M。 Numpy和Scipy需要...
我目前正在使用Richard L Burden的《数值分析》第10版作为近似积分技术的参考。此处描述了输入...
我需要以下问题的帮助。我需要为此的MATLAB代码。我的数据集是:x = [0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5]; y = [1 1.127626 1.543081 2.352409615 3.762195691 6.13228948 10.067662 16 ....
我在R和Python中执行数值方法,我在python中应用了jumpfrog方法,效果很好,但是我想在R中做更苗条的事情。在这里您可以看到我的代码,可能我已经尝试过...]]] >>
[尝试做解决分配根的方法。但是,我的代码仅使用每个列表的第一个索引输出我的根v。误差= 1温度= [273,313,353,373]压力= [1,2,3] v = ...
在python中实现Backward Euler方法来解决摆锤
我正在尝试为摆F设置一个隐式求解器,而dF的定义为:def func(y,t):###简化了消除摩擦的函数,因为它消除了x,v = y [:3 ],y [3:6] ...
因此,我试图在python中实现一些数值方法,并且遇到一些问题,即我的所有函数或多或少都输出与常规euler方法相同的东西。我认为这是...
嵌入在tkinter中的matplotlib动画:从未调用过的更新函数
我一直试图将matplotlib动画嵌入到tkinter中。该应用程序的目标是使用rk4方法仿真一些微分方程,并在仿真过程中显示实时图形。在...
我试图了解如何使用MATLAB进行FFT的积分或微分。但是,我认为我在某处做错了事,想知道我在想什么......'
在Python中实现Adams Bashforth Moulton方法
我目前正在致力于实施亚当斯·巴什福斯·莫尔顿方法来解决摆的问题。我当前的代码如下:将numpy导入为np ###################################### ## ...
我正在尝试使我的runge-kutta第四阶代码模块化。我不想每次使用时都必须编写和声明代码,而是在.hpp和.cpp文件中声明以分别使用它。但是我是...
我正在尝试使用math.net库来计算矩阵的伪逆,但是到目前为止,当我尝试使用以下代码时,我得到一个错误,告诉我:B = A.Inverse();矩阵...
嘿,我正在尝试通过使用sympy.solvers从几个方程式中得出值'b'。从excel工作表中导入了一些此处未给出的值R = 0.287 cp = 1 cv = 0.714 ...
我想构建一个包含公差级别,上下限和函数的函数,以计算自适应梯形正交,并绘制一个图形,例如这样的图形:...