更多pythonic方法来计算numpy中广播加法的导数？

文章/答案/技术大牛

发布

1回答

python、numpy、array-broadcasting

假设是c = a + b，但是a和b是ndarrays，它们的形状不一定相同。也就是说，它们可以是后面的任意两个数组。我有一些输出dl/dc的导数，我想要计算dl/da。如果a和b的形状相同，则为dl/da = dl/db = dl/dc。通常，dl/da是在a中广播的所有轴上的dl/dc的总和。为了在一般情况下计算a和b的链式

浏览 1提问于2017-08-01得票数 0

1回答

我想知道，如果有一个得分函数可以描述，或者评级可以与其他曲线进行比较，如何评估曲线的平滑度。我最初使用标准差来评估曲线的平滑度，当曲线在平坦的条件下是平均的时，它效果很好。然而，如果曲线的斜率使它变得更高，它就不能给出预期的结果，因为它可能会受到最大值或最小值的影响，但我只想考虑平滑度。下面是我的想法的一个例子： ? ?(只看到绿线)第一张图片中绿线的分数或比率应该更高，因为它看起来比第二张更平滑。我在想，我是

浏览 1441提问于2021-06-18得票数 0

1回答

Python中数据点曲线的弧长

python、numpy、scipy、sympy、robotics

我正在进行机器人模拟，并试图计算机器人沿着计划的轨道距离目标的距离。轨迹曲线是为了避免障碍物，它是由一个坐标列表给出的。要找到进展，我需要找到从当前位置到目标的弧长。我熟悉函数的弧长方程：我计划使用的方法是用NumPy的polynomial.polyfit从数据点创建轨迹函数的多项式近似，然后找到它的导数，乘以1，取平方根，最后积分。然而，多项式的平方根并不总是存在的，所以

浏览 7提问于2022-05-26得票数 0

2回答

如何确定是否需要numpy.vectorize()？

python、numpy

作为背景，我问这个问题，是因为我在一个程序中使用Numpy从文献中可用的热力学函数(基于CALPHAD)计算相图。在给定的温度下，计算自由能函数，确定接触凹面的共切曲线(二阶导数> 0)，确定相共存的组成范围。对此，直接定义二阶导数函数是很好的。在我尝试用一个简单的抛物线自由能来测试它之前，所有的一切都很好，它有一个恒定的二阶导数。这破坏了我的</e

浏览 0提问于2014-01-30得票数 3

回答已采纳

2回答

在python中寻找带Abs函数的导数

python-2.7、numpy、sympy、derivative

我想用python计算y=Abs(0.5-0.5*sqrt(1-4*x))在0.1中的导数。这是我的密码：y=Abs(0.5-0.5*sqrt(1-4*x))d=lambdify(x,deriv,'numpy') print d(的新手，所以我想我用了错误的方法来确定导数。d=lambdify(x,deriv,&#

浏览 5提问于2013-12-28得票数 5

回答已采纳

1回答

在numpy数组中扫描相同值的组

python、numpy

我有一个numpy数组，其中0表示空格，1表示位置已填充。我正在尝试寻找一种快速的方法来扫描numpy数组，查找彼此相邻的多个零值，并返回中心零的位置。例如，如果我有以下数组 [0 1 0 1][0 1 0 1] [1 1 1 1] 我想返回在中心零的两边都有相邻零的位置。e.g [1,1] 因为这是3个零的中心，也就是说，在这个位置上，零的两边都有一个零我知道这可以使用if语句来计算

浏览 17提问于2021-01-04得票数 1

回答已采纳

2回答

是否有更多的Pythonic方法可以通过函数的字符串名访问函数？

python、function、variables

我有几个函数(硬编码的，我不希望Python计算导数，实际上函数更复杂)：def g(x): return x ** 2def g_prime(x): return 2 * xmyfunctions = [f, f, g, f] 如何获得相关“素数”函数的列表？0x00000000022CF048>, <function g_prime at 0x00000000022CF0B8>, <f

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

回答已采纳

6回答

“扩展”一个麻木ndarray的好方法？

python、numpy、multidimensional-array

有没有好的方法来“扩展”一个numpy ndarray？假设我有一个这样的ndarray： [3 4]][[1 2 0 0 0]我知道肯定有一些暴力的方法来做到这一点(比如用零构造一个更大的数组，然后从旧的更小的数组中复制元素)，只是想知道有没有pythonic的</

浏览 1提问于2012-10-01得票数 53

回答已采纳

1回答

功能、广播和元素智能操作的矢量化评价

python、numpy、vectorization、elementwise-operations

我必须解释这段代码做什么，知道它执行矢量化的计算F，使用广播和元素的操作概念. D = x_pos.reshape(1,-1) - x_pos.reshape/alpha) * (alpha.reshape(1,-1) * R(D)).sum(axis=1)函数F的第一行接收x_pos和α作为参数(都是numpy数组)，在第二行中，矩阵D是通过广播计算

浏览 1提问于2016-09-20得票数 1

1回答

嵌套for循环中的Python优化整形操作

python、python-3.x、numpy、optimization、reshape

我正在寻找一种更多的pythonic/广播方法来优化以下两种阵列整形功能： """ Reshaping input任何建议，以提高我对numpy整形技术的理解，并避免使用for循环，将不胜感激。谢谢。

浏览 2提问于2019-07-26得票数 0

回答已采纳

2回答

如何从图像中提取指定坐标的像素

python、numpy、opencv、image-processing

我正在使用以下内容从图像中提取一些像素：import numpy as nptargets = np.array(, [333, 333]], dtype=np.uint8)for x, y in targets:有没有更多的pythonic或矢量化的方法来代替这个for循环在一行

浏览 0提问于2021-02-03得票数 0

1回答

潘达斯块经理是如何提高性能的？

python、pandas、numpy、performance

Pandas文档指出: BlockManager的主要好处是提高了某些操作的性能(从2D数组构造、二进制操作、跨列减少)，特别是对于wide DataFrames。我以为我理解了BlockManager是如何通过一篇伟大的文章()来提高性能的，但是我意识到这个例子中有一个小错误。如果我纠正了示例中的错误：a2 = np.arange(128 * 1024 * 1024) a_both = np.empty(. dev

浏览 7提问于2022-06-21得票数 1

回答已采纳

1回答

python中的最快多项式求值

python、algorithm、performance、numpy

我使用的代码如下：import numpy as npfrom numpy import polyval 问题是，我很确定numpy.polyval()会是最快的，但在Linux上，不管x0的

浏览 2提问于2020-11-11得票数 2

回答已采纳

1回答

交换在python类中定义的数组中的元素

python、class、swap

我正在模拟两个不同物种的粒子碰撞。我使用一个类来定义每个物种，并使用一个类来定义保存粒子位置和速度的数组。class Particle(object): self.x=numpy.zeros((num,1))#Position self.v=numpy.zeros((num,3)) #

浏览 0提问于2016-11-01得票数 0

3回答

Numpy数组变量避免计算错误

python、arrays、numpy

我在用幂律方程来设定一个矮小的阵列。问题是，当x是负且n不是整数时，我的部分域尝试执行numpy.power(x，n)。在域的这一部分中，我希望值为0.0。下面是一个具有正确行为的代码，但是是否有更多的Pythonic方法来做到这一点呢？# note mesh.x is a numpy array of length nxpara = [5.8

浏览 3提问于2015-12-17得票数 2

回答已采纳

3回答

平均超过2d的变片numpy阵列

performance、numpy、indexing、slice

我需要计算二维numpy数组中每列切片变化的平均值。例如，我有一个数组每列的结束索引定义为然后，第一列的平均值将是有什么合适的(即Pythonic +高效)的方法来做到这一点？我的

浏览 5提问于2016-06-14得票数 3

回答已采纳

1回答

python中对数图的导数

python、numpy、derivative

我们有x和y值，我用logx = np.log10(x)和logy = np.log10(y)记录他们的日志。我试着计算出w.r.t logx的导数，所以dlogy/dlogx。我曾经成功地使用numpy梯度，更准确地说但是由于一些奇怪的原因，它似乎不再起作用了，产生了错误："Tracebackin <module> grady = np.gradient(

浏览 5提问于2018-01-09得票数 0

回答已采纳

2回答

根据列值过滤numpy* ndarray (矩阵)*

python、matrix、numpy

这个问题是关于根据某些列值过滤NumPy ndarray。我有一个相当大的NumPy ndarray (300000，50)，我正在根据某些特定列中的值对其进行过滤。第一列名为category_code，我需要过滤矩阵，只返回category_code在("A", "B", "C")中的行。结果将需要另一个NumPy ndarray，它的列仍然可以通过dtype名称访问。我

浏览 2提问于2012-08-23得票数 7

回答已采纳

4回答

Numpy纯函数用于性能、缓存

python、numpy、optimization、memoization、blas

我正在用numpy编写一些性能适中的关键代码。这段代码将在内部最循环中，计算的运行时间以小时为单位。快速计算表明，在计算的某些变体中，该代码将执行大约10^12次。因此，函数是计算sigmoid(X)，另一个计算它的导数(梯度)。Sigmoid具有用于在python中，numpy如

浏览 8提问于2014-01-14得票数 13

回答已采纳

1回答

预分配多个相同形状的numpy数组

python、arrays、matlab、python-2.7、numpy

我正在为相同大小的numpy数组(模型网格数据)计算许多属性。我更熟悉matlab语法，在该语法中，可以使用类似于以下命令的命令预先分配多个大小相同的数组：目前，我的python2.7.7 (numpy1.9.0)语法看起来非常糟糕、重复，而且非常不像pythonic语法：ar

浏览 2提问于2014-09-17得票数 1

回答已采纳

点击加载更多