gusucode.com > PPML_v1.2 > general/smpropag_fw.m

    function [S1r,S2r] = smpropag_fw(S1,S2,p1,p2,A1,A2,f1,f2)

% Forward propagation of scattering matrix. 

X = p1\p2; Y = A1\A2;
I = (X+Y)/2;
J = (X-Y)/2;

S1r = (I - diag(f1)*S2*J)\(diag(f1)*S1);
S2r = (I - diag(f1)*S2*J)\(diag(f1)*S2*I - J)*diag(f2);

end