1

首先，我想为我的英语道歉。这不是我的第一语言。

这是我正在开发的一个程序，用于对重复负载三轴测试产生的数据进行排序。第一列是一个段。三段形成一个循环。加载、卸载、暂停。大约 50 个数据点用于加载，同样用于卸载，100 个用于暂停。我从 121 开始加载。对于本节，我需要第四列（索引 3）的加载和卸载阶段的最大值。

我的数组有：

[[  1.21000000e+02   1.00313720e+02   2.00015190e-02 ...,   2.25933480e-01
2.95645450e-01  -3.33373370e-01]
[  1.21000000e+02   1.00318600e+02   2.00071220e-02 ...,   2.25933600e-01
2.95629110e-01  -3.33358880e-01]
[  1.21000000e+02   1.00323490e+02   2.00045150e-02 ...,   2.25932690e-01
2.95642500e-01  -3.33374260e-01]
..., 
[  1.50000000e+02   1.10347900e+02   2.00072340e-02 ...,   2.24460500e-01
2.94727620e-01  -3.38975370e-01]
[  1.50000000e+02   1.10352780e+02   1.99971700e-02 ...,   2.24458930e-01
2.94705120e-01  -3.38966550e-01]
[  1.50000000e+02   1.10357670e+02   2.00063640e-02 ...,   2.24455860e-01
2.94704710e-01  -3.38963510e-01]]

在我通过那个循环运行它之后：

g = np.loadtxt('test.txt')

Sigmad = []


DataCol = np.hsplit(g, g.shape[1])
DataCharge = DataCol[3] #Charge mean load
DataCycle = DataCol[0]


ld = 0 #ligne de début, start line of the load segment
fc = 0 #Fin de chargement, end line of the unload segment
seg1 = DataCycle[0] #Segment initial, numerical value of the first segment

chargemax = []
i = 0


while i < len(DataCycle):

    if DataCycle[i] == seg1 and DataCycle[i+1] == seg1 and DataCycle[i-1] == seg1 - 1:
        ld = i        

    elif DataCycle[i] == seg1 + 2 and DataCycle[i+1] ==  seg1 + 2 and DataCycle[i-1] == seg1 + 1 :

        fc = i
        print seg1
        chargemax.append(np.max(DataCharge[ld:fc]))
        ld = i + 1
        seg1 += 3

    i+=1

Sigmad.append(1000*np.mean(chargemax[:])/aire)

g 被修改了，我不知道为什么， g 现在等于：

[[  1.51000000e+02,   1.00313720e+02,   2.00020000e-02, ...,
      2.25933000e-01,   2.95645000e-01,  -3.33373000e-01],
   [  1.21000000e+02,   1.00318600e+02,   2.00070000e-02, ...,
      2.25934000e-01,   2.95629000e-01,  -3.33359000e-01],
   [  1.21000000e+02,   1.00323490e+02,   2.00050000e-02, ...,
      2.25933000e-01,   2.95642000e-01,  -3.33374000e-01],
   ..., 
   [  1.50000000e+02,   1.10347900e+02,   2.00070000e-02, ...,
      2.24460000e-01,   2.94728000e-01,  -3.38975000e-01],
   [  1.50000000e+02,   1.10352780e+02,   1.99970000e-02, ...,
      2.24459000e-01,   2.94705000e-01,  -3.38967000e-01],
   [  1.50000000e+02,   1.10357670e+02,   2.00060000e-02, ...,
      2.24456000e-01,   2.94705000e-01,  -3.38964000e-01]]

有人可以向我解释为什么会这样。第 4 列中的一些值也会被修改。

谢谢！

4

1 回答 1

2

在这个部分：

DataCol = np.hsplit(g, g.shape[1])
DataCharge = DataCol[3] #Charge mean load
DataCycle = DataCol[0]

[...]

seg1 = DataCycle[0]

您不是在制作副本，您只是为原始数组的视图赋予新名称。所以当你这样做时

seg1 += 3

在循环内部，您正在修改g自己。例如：

>>> a = np.array([[1,2,3], [4,5,6]])
>>> b = a[0]
>>> b
array([1, 2, 3])
>>> b += 100
>>> b
array([101, 102, 103])
>>> a
array([[101, 102, 103],
       [  4,   5,   6]])

如果你想要一个副本，你可以.copy()显式调用：

>>> a = np.array([[1,2,3], [4,5,6]])
>>> b = a[0].copy()
>>> b
array([1, 2, 3])
>>> b += 1000
>>> b
array([1001, 1002, 1003])
>>> a
array([[1, 2, 3],
       [4, 5, 6]])
于 2013-01-15T02:49:46.303 回答