橫波時差計算公式

if (ISAC == 1)//用AC代替DTS,用AC代替DTC

{

////使用SH和POR1曲線計算的

if (POR1 >= 1.F)

RMSC = 5.F;

else

{

float fPOR = POR1; //原來的

float VSAD = 1 - fPOR - SH1;

float VDSH = 1 - fPOR;

float RCOM = VPSM * VSAD / VDSH * SH1 / VDSH;

//VPSM 引數

RCOM = VPSM * VSAD / VDSH + VPSH*SH1 / VDSH;

float fTempR = 2.F - fPOR;

RMSC = RCOM * (fPOR * (CLOG / log10(fTempR) - 1.F) + 1);

}

if (SH1 > 1)

RMSC *= SH1;

else if (SH1 < 0)

RMSC /= log10(10 - SH1);

DTS = RMSC * DTC;

}

else if (ISAC == 3)//使用SH引數計算的DTS

{

float ASH = (DTC - DTCM) / (DTCL - DTCM);

ASH = (float)((powf(2, ASH * 2) - 1.F) / 3.F);

float VSAD = 1 - POR1 - ASH;

float VDSH = 1 - POR1;

float RCOM = VPSM * VSAD / VDSH + VPSH * ASH / VDSH;

RMSC = RCOM * (POR1 * (CLOG / log10(2 - POR1) - 1.F) + 1);

DTS = RMSC * DTC;

}