查了一些资料，在论坛上看了一下，好像没有这个问题的解决方法，所以把自己试探性找到的解决方案放在下面，如果哪里不恰当，欢迎交流。
eststo:  xtivreg2 y c year (x1 x2=l.x1 l.x2),fe cluster(id) first savefirst savefprefix(alpha)
esttab alphax1 alphax2 est1 using tset1.rtf
这样就能导出多个工具变量的2SLS第一阶段回归结果了。