求解符号代数方程实例（1）

应用背景：求解光在单轴晶体中，各传播方向的光波与折射率的关系

将（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)