经过借鉴前人成果，结合自己实践，终于完成了，废话不多说：
首先由Geoda生成距离权重矩阵——保存——生成Wd.gwt

打开stata
*安装软件包 ssc install sppack，包含-spmat-, -spreg-, and -spivreg-等命令
*使用-spmat import- 读取.GWT文件
spmat import M using Wd.gwt, geoda

*如果还是想生成一个 Stata 格式的权重矩阵文件，需要从spmat对象释出权重矩阵并将其保存为一般的Mata matrix
*假设其名称为W

*该命令会创建一个名叫W的spmat对象，存储从Geoda导入的权重矩阵信息.
spmat getmatrix M W

*进一步保存为 Stata 格式，数据集中的变量会被命名为x1, x2, x3, ……
getmata (x*)=W


之后会在stata数据编辑器中看到空间距离权重矩阵，但是很多小伙伴发现x1, x2, x3, ……和自己定义的编码是不对应的（换句话说：是矩阵乱的）
对应方法：
Arcgis打开生成空间权重矩阵的.shp面打开属性表，属性表中的FID列的数字是和stata生成的权重矩阵中的x1, x2, x3, ……相对应。之后在同一替换，就会得到你想要的矩阵！