Matlab代码的Python实现--有限差分法

文章/答案/技术大牛

发布

1回答

python、matlab、numpy、numerical-methods

根据我的老师创建的这段Matlab代码： function [] = explicitWave(T,L,N,J)% T: LengthPython中实现此功能： import numpy as np plt.plot(x, Une

浏览 25提问于2020-11-08得票数 1

回答已采纳

1回答

Matlab有限差分法

matlab

我在matlab中的代码有个小问题。keep the above graphhold off; 我需要编辑代码，而不是固定的"x^3+x^2-1.25x-0.75“，我可以自己在matlab的主窗口中定义参数。

浏览 1提问于2015-06-09得票数 0

回答已采纳

1回答

用matlab实现有限差分法

matlab、differential-equations

我试图在matlab中实现有限差分法。我做了一些计算，得到y(i)是y(i-1)和y(i+1)的函数，当我知道y(1)和y(n+1)时。但是，我不知道如何实现这一点，因此y的值是以正确的方式更新的。

浏览 1提问于2015-01-12得票数 0

回答已采纳

1回答

基于python边界条件的二阶ODE

python、odeint

我试图用y'' + ay' + by + c = 0 (二阶微分方程)的形式用odeint来求解python中的方程。据我理解，odeint只适用于初始条件为y(0) = y1，y'(0) = y2。而我的条件是: y'( 0 ) = 0，y'(pi/4) =0。在这种情况下，有什么方法可以使用odeint吗？下面是我的代码：import nump

浏览 1提问于2018-02-26得票数 0

回答已采纳

2回答

用有限差分计算给定函数的导数，Python

python、function、reference、derivative

我可能首先应该说，我对python还比较陌生，但我以前用java和Matlab编写过代码。在python中，代码 return fprint(g(0))>1.0def derivator(f, h = 1e-8):g

浏览 1提问于2015-04-09得票数 2

回答已采纳

1回答

有限差分法的实现

code-golf、math、number

有限差分法是一种在数列中寻找下一项的技巧，因为该序列是由某一多项式的连续值决定的。给定一个n项的列表，该方法将能够确定多项式的程度，这是n+1的最大值。然后，任何数字与其右边的邻居之间的区别(假设它有一个右边的邻居)写在下面的行中。注意，条目必须在前一列最长条目的右侧写入一个空格(ASC

浏览 0提问于2015-12-01得票数 8

1回答

Hessian在Python中的有限差分方法

python、numerical-methods、differentiation

现在，我想在python中实现这个函数z的hessian。到目前为止，我已经用有限差分法计算了导数，如下所示- def derivative(x, y, f, h): return [(f(x + h, y) - f(x - h, y )) / (2*

浏览 28提问于2019-09-18得票数 0

1回答

梯度计算/装配非常慢

openmdao

我正在使用OpenMDAO在Python中实现现有的MATLAB优化代码。这基本上是与MATLAB代码相同的结构。在MATLAB方面，我正在计算问题级的解析梯度，但是对于该部分的有限差分梯度，我将大组调用14次。在Python中，我使用手工计算的梯度和该大组及其子组中组件的自动微分工具的混合。在MATL

浏览 8提问于2020-11-09得票数 0

回答已采纳

2回答

“pdepe”(MATLAB函数的)使用的是哪种数值方法？

matlab、function、numerical-methods、pde

我用MATLAB的函数'pdepe‘来解一个偏微分方程的问题，一个抛物方程。谢谢你帮我。

浏览 4提问于2017-04-02得票数 0

回答已采纳

1回答

求解圆形区域上的泊松方程

physics、poisson

假设存在一个具有圆区域的泊松方程，且边界条件为零。我想用fininte差分法进行数值求解。然而，据我所知，有限差分在角区域有着广泛的应用。那么，如果我仍然使用有限差分在圆区域，我应该做什么？

浏览 7提问于2022-04-13得票数 0

1回答

用Matlab求解差分方程

matlab、math、calculus、differentiation

作为一个编程练习，我编写了一个Matlab函数，它用有限差分法找到函数的导数。在一个脚本中，我调用了函数，并希望使用内置函数来检查它，但我在实现这个函数时遇到了困难。为了得到一个检查，我使用diff(eqn)来得到微分方程，但是我不知道如何用这个方程来求解一个特定的值，例如。X= 2。这是我的代码： %set value for the derivative to be evaluated

浏览 2提问于2014-08-09得票数 0

1回答

离散周期数据的导数计算

methods、derivative、spectral

我有一个描述周期信号的数组yx，x=0,1,2，...,10^6，y(10^6)=y(0)，我想用一种快速的方法计算它的导数dy/dx。我尝试了光谱差分法，即所得结果与(yx+1-yx-1)/2 (有限差分法)的结果不同。以下是了解结果差异的努力：如

浏览 1提问于2012-04-18得票数 1

回答已采纳

2回答

用matlab/mathematica数值计算它不能符号化的不定积分

matlab、wolfram-mathematica

我试图在Matlab和Mathematica中计算一个函数的积分，这是软件无法象征性地完成的。a = @(t) sin(t);F = int(asin(sin(t在数学堆栈交换中，我找到了一个，它在一点上导出了一个像积分一样的有限</e

浏览 2提问于2018-09-12得票数 1

回答已采纳

2回答

Fick第一扩散定律在模拟细胞膜水运动中的应用

python

我想用fick定律做上面的模拟，用有限差分法。我的问题是整个问题的开始。物理特性：边界条件： C = sc.zeros(5)F = -

浏览 4提问于2012-01-14得票数 1

回答已采纳

1回答

如果无法获得梯度，是否有类似于梯度下降的方法？

gradient-descent

函数f的输入是一组参数，输出是精度。f的梯度不能得到，有没有类似梯度下降的方法，比如有限差分法？

浏览 0提问于2017-10-03得票数 0

1回答

如何改进Levenberg-Marquardt多项式曲线拟合方法？

machine-learning、regression、numerical-methods、non-linear-regression、levenberg-marquardt

几周前，我开始在Matlab中从头开始编写Levenberg-Marquardt算法。我对数据的多项式拟合感兴趣，但是我还没有达到我想要的精度水平。在尝试其他多项式之后，我使用了一个五阶多项式，这似乎是最好的选择。无论我尝试实现什么改进，该算法总是收敛到相同的函数最小化。到目前为止，我没有成功地添加以下功能：作为二阶correctionDelayed满足项的测地加速度项，用于根据有限差分法的ite

浏览 2提问于2020-06-06得票数 0

回答已采纳

1回答

欧拉-拉格朗日方程数值解的算法建议

python、numerical-methods、pde

在我的课程作业中，我被告知要用Python编写一个算法来求解具有Dirichlet边界条件的欧拉-拉格朗日方程。非常感谢!

浏览 0提问于2018-11-17得票数 0

1回答

我正在尝试用有限差分法在python中编写2D Schrödinger方程的曲线图

python、physics

我已经开始编写一个程序，用有限差分法计算二维Schrödingers方程的解。我想要使用等高线图或一些其他图形显示，通过接受用户输入的尺寸和网格点数量，以图形方式显示解决方案。我将hbar^2/2m设置为1，并将电位(V)设置为0，从而简化了Schrödinger方程，从而得出： -(dψ^2/dx^2 + dψ^2/dy^2) = E*ψ 使用有限差分法，方程的左侧变成如下形式的矩阵：

浏览 16提问于2020-04-27得票数 0

2回答

需要帮助使用Numpy计算导数和积分

python、methods、difference、integral、derivative

我需要帮助来计算函数的导数和积分使用有限差分法和Numpy而不使用循环。整个任务是:对m=0和s=0.5,5，在-10,10的区间上，将高斯函数f(x) = (1./(sqrt(2.*pi)*s))*e**(-0.5*((x-m)/s)**2)列成表，用有限差分法计算函数的导数和积分，不需要循环。创建函数及其导数的曲线图。使用Numpy和Matplotlib。下面是程序的开始部

浏览 0提问于2014-12-21得票数 0

1回答

求解微分方程的Python有限差分法

python、numpy、differential-equations

我必须解Euler Bernoulli微分光束方程，它是：(x是梁轴点的坐标)u(0)=0, u’(0)=0, u’’(1)=0,u’’’(1)=a U’k = (1/2*h)(Uk+1 - Uk-1) U’’’’k = (1/h4)(Uk+2 - 4 U

浏览 4提问于2014-11-16得票数 3

点击加载更多