python,克鲁格曼产业同构

求python计算克鲁格曼产业结构差异度计算代码
[size=10.4637pt]其中[size=10.4637pt], [size=10.4637pt]n [size=10.4637pt]、[size=10.4637pt]m [size=10.4637pt]表示两个地区[size=10.4637pt], [size=10.4637pt]t [size=10.4637pt]表示时间[size=10.4637pt], [size=10.4637pt]i [size=10.4637pt]表示行业 [size=10.4637pt]。[size=10.4637pt]Y[size=5.97924pt]nm[size=5.97924pt], [size=5.97924pt]t[size=10.4637pt]表示的是 [size=10.4637pt]:[size=10.4637pt]在 [size=10.4637pt]t [size=10.4637pt]期[size=10.4637pt], [size=10.4637pt]地区 [size=10.4637pt]n [size=10.4637pt]、[size=10.4637pt]m [size=10.4637pt]之间的[size=10.4637pt]产业趋同状况[size=10.4637pt]。[size=10.4637pt]绝对值符号内是两个地区某个行业产量占地区总产量的比率之差[size=10.4637pt], [size=10.4637pt]然后再对 [size=10.4637pt]I [size=10.4637pt]个[size=10.4637pt]行业加总 [size=10.4637pt]。[size=10.4637pt]如果地区间产业结构差异越大[size=10.4637pt], [size=10.4637pt]那么这个指标的值应该越大[size=10.4637pt], [size=10.4637pt]反之[size=10.4637pt], [size=10.4637pt]如果两地产业结构[size=10.4637pt]趋同[size=10.4637pt], [size=10.4637pt]这个指标将越小[size=10.4637pt]。

谢谢分享！