传输矩阵法算缺陷层吸收的matlab编程，下面是算TE波的，怎么改成算TM波的？据说要改矩阵的一个地方还有q0，可是小白不懂啊！
%TE
%q0=1/c/mu0*cos(theta);epsilon0=8.85*10^-12;
%transmission=2*q0/(totalT(1,1)*q0+totalT(1,2)*q0*q0+totalT(2,1)+totalT(2,2)*q0);
%tt(p)=abs(transmission)^2;
%reflection=(totalT(1,1)*q0+totalT(1,2)*q0*q0-totalT(2,1)-totalT(2,2)*q0)/(total(1,1)*q0+totalT(1,2)*q0*q0+totalT(2,1)+totalT(2,2)*q0);
%rr(p)=abs(reflection)^2;
%absorption(p)=1-tt(p)-rr(p);