倾斜界面求波的反透射系数matlab

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