ValueError: 操作数不能与形状 (7,) (6,) (7,) 一起广播答案

【问题标题】：ValueError: operands could not be broadcast together with shapes (7,) (6,) (7,)ValueError: 操作数不能与形状 (7,) (6,) (7,) 一起广播
【发布时间】：2019-12-20 07:03:40
【问题描述】：

我收到一条错误消息，尝试在 Python 中运行脚本以重复进行金箔辐照冷却时间。这是错误：

Traceback (most recent call last):
  File "C:/Users/wolfk/PycharmProjects/Reactor Phys In-House/iterations.py", line 26, in <module>
    cr_tot /= math.e ** (-2.98256E-6 * out_times[l])
ValueError: operands could not be broadcast together with shapes (7,) (6,) (7,)

代码非常混乱，因为我大部分使用数值而不是变量，但数学并不重要，我可以稍后验证。我需要做的是让这个东西运行起来！所以你知道，进出时间的向量在设计上是不同的长度，我试着让它们的长度相同，所以这不是问题。代码如下：

import numpy as np
import math
import matplotlib.pyplot as plt

in_times = np.array([[697800, 603840, 74460, 161520, 428880, 75960, 62760]])
out_times = np.array([[248880, 3960, 9240, 5400, 5880, 8760]])

cr_tot = 0
cr_fast = 0
i = 1
j = 1
k = 0
l = 0

while i < 13:
    if j == 1:
        cr_tot /= math.e ** (-2.98256E-6 * in_times[k])
        cr_fast /= math.e ** (-2.98256E-6 * in_times[k])

        cr_tot += ((2.98256E-6 * 4950377) / 0.021554) / (1 - math.e ** (-2.98256E-6 * in_times[k]))
        cr_fast += ((2.98256E-6 * 2141114) / 0.043109) / (1 - math.e ** (-2.98256E-6 * in_times[k]))

        j -= 1
        k += 1
    else:
        cr_tot /= math.e ** (-2.98256E-6 * out_times[l])
        cr_fast /= math.e ** (-2.98256E-6 * out_times[l])

        j += 1
        l += 1

    i += 1

xi_tot = (2.98256E-6 * 4950377) / (0.021554 * cr_tot)
xi_fast = (2.98256E-6 * 2141114) / (0.043109 * cr_fast)

cr_tot *= (2.98256E-6 * 4950377) / (0.021554 * xi_tot * math.e ** (-2.98256E-6 * 660) * (1 - math.e ** (-2.98256E-6 * 1000)))
cr_fast *= (2.98256E-6 * 2141114) / (0.043109 * xi_fast * math.e ** (-2.98256E-6 * 660) * (1 - math.e ** (-2.98256E-6 * 2000)))

print(cr_tot)
print(cr_fast)

如果能得到任何帮助，我将不胜感激。 :)

【问题讨论】：

标签： python arrays numpy vector

【解决方案1】：

从初始的 numpy 数组中删除双括号实际上会有所帮助。

...
in_times = np.array([697800, 603840, 74460, 161520, 428880, 75960, 62760])
out_times = np.array([248880, 3960, 9240, 5400, 5880, 8760])
...

输出：

0.476018036910266
0.9882585377315589
0.9728174262826453
0.9840231812760558
0.9826154332766018
0.9742111371275044
7480153046556.081
175166298775.0072

【讨论】：

【解决方案2】：

我不明白你的逻辑。在您的代码中，开始时 k=0 和 in_times[k]in if 分支将返回数组（[697800, 603840, 74460, 161520, 428880, 75960, 62760]）。因此，cr_tot 的长度为 7。然后 j -= 1，j 等于 0。下一次迭代转到 else 分支。 out_times[l] 表示 out_times[0] 并返回长度为 6 的数组 ([248880, 3960, 9240, 5400, 5880, 8760])。但是您的代码是 cr_tot /= math.e ** (-2.98256E-6 * out_times[l])，使用长度为 7 的数组来划分 cr_tot(length is 6 return from if branch)。这就是你得到错误的原因。

【讨论】：