bv_4.m

通过调节入射光的波矢和偏振的配置

syms lamda l1 l2 l3 l4 m1 m2 m3 m4 n1 n2 n3 n4  % a b1 b2 c1 c2 c3 
a=1;b1=0.33;b2=1.5;c1=3;c2=1.2;c3=0.8;

l12=l1-l2;l13=l1-l3;l14=l1-l4;
m12=m1-m2;m13=m1-m3;m14=m1-m4;
n12=n1-n2;n13=n1-n3;n14=n1-n4;

eq1=l12-lamda/a;
eq2=l13;
eq3=l14;
eq4=m12+b1*lamda/(a*b2);
eq5=m13+lamda/b2;
eq6=m14;
eq7=n12-lamda/(a*c3)*(c2*b1/b2-c1);
eq8=n13-c2*lamda/(c3*b2);
eq9=n14-lamda/c3;
eq10=l1^2+m1^2+n1^2-1;
eq11=l2^2+m2^2+n2^2-1;
eq12=l3^2+m3^2+n3^2-1;
eq13=l4^2+m4^2+n4^2-1;
s=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12,eq13,'l1,l2,l3,l4,m1,m2,m3,m4,n1,n2,n3,n4,lamda');

lamda=subs(s.lamda(3));

ek=subs([s.l1(3),s.m1(3),s.n1(3);s.l2(3),s.m2(3),s.n2(3);s.l3(3),s.m3(3),s.n3(3);s.l4(3),s.m4(3),s.n4(3)]);
k=2*pi/lamda*ek;
phi=zeros(4,1);
theta=zeros(4,1);
for i=1:4
    phi(i,1)=acos(ek(i,3));
    if(phi(i,1)==0)
        theta(i,1)=0;
    else
        if(ek(i,2)>0)
            if(ek(i,1)>0)
                theta(i,1)=asin(ek(i,2)/sin(phi(i,1)));
            else
                theta(i,1)=pi/2+asin(ek(i,2)/sin(phi(i,1)));
            end
        elseif(ek(i,2)<0)
            if(ek(i,1)>0)
                theta(i,1)=asin(ek(i,2)/sin(phi(i,1)));
            else
                theta(i,1)=asin(ek(i,2)/sin(phi(i,1)))-pi/2;
            end
        else
            theta(i,1)=0;
        end
    end
end