EMD模态分解，怎样计算分解后的IMF正交性指标？

Question

如图，为EMD原始论文中计算正交性指标的公式，根据此公式编写代码如下：

IO=[]
for i in range(8928):
    sum_io=0
    for j in range(12):
        a=imf[j]
        a1=a[i]
        for k in range(j+1,12):
            b=imf[k]
            b1=b[i]
            chenji=a1*b1
            sum_io=sum_io+chenji
    if(oringal[i]==0):
        shang=0
    else:
        shang=sum_io/(oringal[i]*oringal[i])
        shang=('%.4f'%(shang))
    IO.append(shang)

代码中oringal[i]对应原始数据，对应公式中的x(t)。
根据这段代码计算出来的结果很大，如图

请各位看一看，代码哪块编写有错误呢?