值错误：使用 numpy 和 Bessel 函数将数组元素设置为序列

Question

我一直在尝试解决值错误（将数组元素设置为序列），但我尝试的任何方法都无法纠正问题。代码围绕使用第一类贝塞尔函数并将它们与jv(v,z) 函数进行比较，但a[k] = B1(theta[k]) 的行导致了我的错误。我只使用 Python 3 编码大约 6 周，所以我对此还是很陌生。

编辑缩进

import numpy as np
import scipy
import pylab as plt
from numpy import pi
from scipy import integrate
from scipy import special
from numpy import cos
from numpy import sin

x = 1.0
def B1(theta):
  j1 = cos(theta-x*sin(theta))
  return (j1)

theta = np.linspace(0.0, 20.0, 1000)
a = np.linspace(0.0, 20.0, 1000)
b = a
r = np.linspace(0.0e-6, 1.0e-6, 1000)
x = np.linspace(-1.0, 1.0, 1000)
y = x
Ij = (1.0/pi)*integrate.simps(B1(theta), theta)

for k in range(0, len(theta)):
  if k>0: 
    a[k] = B1(theta[k])
    b[k] = scipy.special.jv(1, theta[k])
  else:
    a[k] = 1.0
    b[k] = 1.0
plt.plot(theta, a,'r-')
plt.plot(theta, b,'k-')
plt.show() 

任何人的帮助将不胜感激，因为我已经为此苦苦挣扎了几个小时。

编辑：作为参考，我的跟踪日志如下 文件“Ag1P1.py”，第 32 行，在 a[k] = B1(theta[k]) ValueError: 使用序列设置数组元素。

我认为现在我的错误在于我使用以下内容的地方

     x = 1.0
    def B1(theta):
      j1 = cos(theta-x*sin(theta))
      return(j1)

我不确定为什么，但我输入的任何值都是这样；

      a[k] = B1(theta[k])

产生这个错误。但是我不知道有什么解决方案可以纠正这个问题。

Answer 1

我记得，贝塞尔函数只为正参数定义。也许您正在传递一个导致ValueError 的负参数。我建议添加这样的断言作为调试手段。

for k in range(0, len(theta)):
  if k>0: 
    assert theta[k] > 0
    a[k] = B1(theta[k])
    b[k] = scipy.special.jv(1, theta[k])
  else:
    a[k] = 1.0
    b[k] = 1.0

此外，您应该发布您的堆栈跟踪（错误说明了什么），这将有助于我和其他人调试。