function?[RTp?RTs]=visoelasticRT（kxkpzkszmiu）
%?反射p波???反射s波?透射p波?透射s波
%有局部坐標系中的慢度及參數，推導出全局坐標系中的反透射系數
%kx:x方向的局部坐標系中水平慢度
%kpz：z方向分界面上下的P波局部垂直慢度
%KSZ：z方向分界面上下s波的局部垂直慢度
????dp1=kpz（1）;?dp2=kpz（2）;
????%dp1：p波在上界面局部坐標系中的垂直慢度
????%dp2：p波在下界面局部坐標系中的垂直慢度
????ds1=ksz（1）;?ds2=ksz（2）;
????%ds1：s波在上界面局部坐標系中的垂直慢度
????%ds2：s波在下界面局部坐標系中的垂直慢度
????miu1=miu（1）;miu2=miu（2）;
????%上下界面的miu
????A=zeros（44）;
????A（1:）=[kx?-ds1