导数计算中的双标量除零我正在学习使用python进行数值计算。我想用中心差分法计算导数。当我试图设置我的间隔dx时,python将其设置为0,即使实际值是(1/6)?有什么办法摆脱这个吗? 代码如下: import num ...2024-05-14 已阅读: n次
Python有限差分样条我试着做一个程序,用有限差分法在点和点之间插值。它必须能够返回xy坐标,以便可以绘制到屏幕上 我的样条线类: class Spline(): def __init__(self): ...2024-05-14 已阅读: n次
数字梯度指定轴(中心差分法)我目前使用的是numpy的1.92版本,我相信这是最新的公开版本。我希望在n维数组上快速使用中心差分法,并指定执行计算的轴。在 我的第一个想法是数字差异允许使用axis规范,但是这返回一个右差而不是中 ...2024-05-14 已阅读: n次
如何使用有限差分法计算阵列的导数?我用这个代码创建了一个合成数据- # x variable x_value = np.linspace(0,2*np.pi,1000) # y variable for i in x_value: ...2024-05-14 已阅读: n次
如何使用NumPy数组算术Python对有限差分方法进行矢量化?尝试使用Numpy算法将有限差分法矢量化为数值数组(不包括边值)。应该取周围值(上、下、左、右)的平均值来计算它。在 例如,应该能够使用一个numpy数组,比如: u = np.array([[100 ...2024-05-14 已阅读: n次
需要帮助用Numpy计算导数和积分吗我需要帮助用有限差分法和Numpy来计算函数的导数和积分,而不使用循环。在 整个任务:在区间[-10,10]上对m=0和s=[0.5,5]的高斯函数f(x)=(1./(sqrt(2.*pi)*s))* ...2024-05-14 已阅读: n次
Python中的行顺序与笛卡尔平面的抽象 如何将二维笛卡尔平面(x,y)的坐标系转换为索引为a[i,j]的矩阵?在 假设你必须用有限差分法来模拟二维平面中的热传导。系统的边界是左下角的一个点pos_inf,右上角是pos_sup。在 结果形 ...2024-05-14 已阅读: n次
使用python的有限差分法的矢量化版本我需要用python实现有限差分法。但是,我想使用矢量化版本来实现它。我想解方程‘dA(z,t)/dz=-alpha*A(z,t)’。其中z,t是函数A的空间坐标和时间坐标,该函数被视为向量 我已经初 ...2024-05-14 已阅读: n次
在python中使用有限差分我尝试使用Python和Numpy来使用有限差分法来求解一个基本方程。代码给了我正确的第一个值,也就是说它给了我一个[1];但是,之后的其他值都是零? 我不知道我做错了什么,因为它显然适用于第一个值, ...2024-05-14 已阅读: n次
如何有效地lambdify稀疏矩阵?我正在使用Symphy构建一个稀疏(N^2 x N^2)矩阵,并尝试将其转换为稀疏scipy矩阵。它用于有限差分法中,根据k向量k=(kx,ky)在二维网格上求解薛定谔方程 稠密矩阵到SciPy 我目 ...2024-05-14 已阅读: n次
python中高阶偏微分方程的有限差分法我一直在尝试求解一些偏微分方程,虽然我已经能够成功地用有限差分法求解扩散方程或波动方程,但对于高阶偏微分方程,我很难做到。在 例如,假设我想解下列方程 U_{xxxx}-x-2x*x_{tx}-2x_ ...2024-05-14 已阅读: n次
尝试理解python矩阵中索引切片的语法我一直在查看和解释一些示例代码,在理解一些切片语法时遇到了一些问题。在某些情况下,我正在编写一个计算程序,该程序使用有限差分法近似正方形区域的温度 以下是一些代码: maxIter = 500 len ...2024-05-14 已阅读: n次
mimp===mimpy==图片:https://api.travis-ci.org/ohinai/mimpy.png?branch=master:target:https://travis ci.org/o ...2024-05-14 已阅读: n次
fdmFDM: Finite Difference Methods fdm用有限差分估计导数。 另请参见FDM.jl。 请参阅docs。 Installation Multivariate Deri ...2024-05-14 已阅读: n次
findiffp 用最陡下降法、共轭梯度法和预处理共轭下降法求解椭圆型方程的有限差分法 此包Python名称:findiffp 目前版本: findiffpy 1.0 ...2024-05-14 已阅读: n次