我通过双for循环计算元素如下。

N,l=20,10
a=np.random.rand(N,l)
b=np.random.rand(N,l)
r=np.zeros((N,N,l))

for i in range(N):

for j in range(N):
    r[i,j]=a[i]*a[j]*(b[i]-b[j])-a[i]/a[j]

题：

如何对阵列进行矢量化并通过广播进行计算？

我还想设置索引i不等于j，这意味着将对角元素保留为零。我也可以通过矢量化来做到这一点吗？