0

我想从 python 中的数组中减去一个数字,但我想保留原始数组。

这是我的代码:

for X0 in [X0 * 0.01 for X0 in range(-2,2)]:
    for Y0 in [Y0 * 0.01 for Y0 in range(6,10)]:
        RAm = RA
        Decm = Dec
        RAm[:] = [x - X0 for x in RAm]
        Decm[:] = [x - Y0 for x in Decm]
        RAs = np.sum(RAm)
        Decs = np.sum(Decm)
        #a = np.sqrt(((RAs**2*Decs**2)**2 - RAs**4*Decs**4)/(RAs**2*Decs**2*Decs**2 -Decs**2*RAs**4))
        #b = np.sqrt(((RAs**2*Decs**2)**2 - RAs**4*Decs**4)/(RAs**2*Decs**2*RAs**2 - RAs**2*Decs**4))
        print RA
        print RAm

这让我知道 RA 正在发生变化,即使它应该保持不变,因为我想在每个循环中从 RA 中减去一个不同的数字。我怎样才能解决这个问题?

4

3 回答 3

2

如果您仍然在使用 NumPy,您可以利用它同时支持对所有元素的操作这一事实。

RAm = np.array(RA) - X0
Decm = np.array(Dec) - Y0

或者甚至跳过这些作业,如果您没有将结果用于除总和之外的任何内容:

RAs = np.sum(np.array(RA) - X0)
Decs = np.sum(np.array(Dec) - Y0)

我的印象是你是一个(Python)初学者,所以也许坚持使用基本的 Python,直到你觉得舒服为止。当然,所有这些都可以在没有 NumPy 的情况下完成:

RAs = sum(x - X0 for x in RA)
Decs = sum(y - Y0 for y in Dec)

但是,如果您确实想使用 NumPy,请正确使用它。这意味着:使用矩阵计算而不是循环。例如,您可以执行以下操作* 来获得array不同Y0值的总和:

Y0_values = np.array([Y0 * 0.01 for Y0 in range(6,10)])
y_matrix = np.tile(RA, [len(Y0_values), 1])
Y0_matrix = np.tile(Y0_values[np.newaxis].T, [1, len(RA)])
sums = np.sum(y_matrix - Y0_matrix, axis=1)

*仅用于演示目的。我远不是 NumPy 专家,所以可能有更好的方法来做到这一点。出于同样的原因,我不会解释它是如何工作的。我只是鼓励你更多地了解 NumPy,这对你有利。

于 2013-04-25T18:01:17.020 回答
1

而不是遍历列表并对其进行更改,请执行以下操作:

for element in list(your_list_here):
    if element.something():
        your_list_here.remove(element)

另外..分配Ram=Ra意味着对对象的每个方法调用,实际上是指同一个对象。Ram = list(Ra)如果您想安全,请使用。

于 2013-04-25T17:51:40.350 回答
1

RAm = RA为;创建RAm别名 修改列表。RARAm[:] = ...

相反,只做

RAm = [x - X0 for x in RA]
于 2013-04-25T17:55:35.743 回答