应用背景:求解光在单轴晶体中,各传播方向的光波与折射率的关系
将(48)、(49)代入光的波法线菲涅尔方程
得到
利用Matlab求解符号代数方程求解(50)
>> syms n n0 ne a
y=(n^2-n0^2)*(n^2*(n0^2*(sin(a))^2+ne^2*(cos(a))^2)-n0^2*ne^2)
>> solve(y,n)
ans =
n0
-n0
(n0*ne)/(n0^2*sin(a)^2 + ne^2*cos(a)^2)^(1/2)
-(n0*ne)/(n0^2*sin(a)^2 + ne^2*cos(a)^2)^(1/2)