作业帮MATLAB方程求解问题
来源：学生作业帮 编辑：拍题作业网作业帮 分类：数学作业 时间：2024/06/03 01:42:16
MATLAB方程求解问题
clear
syms n
H=1.60;
lanmuda=1.55;
n1=1.55;
n2=1.67;
n3=1.5;
k0=2*pi/lanmuda;
c1=n2^2/n1^2;
c2=n2^2/n3^2;
r1=sqrt(n.^2-n1^2)*k0;
r2=sqrt(n2^2-n.^2)*k0;
r3=sqrt(n.^2-n3^2)*k0;
fun=r2*H-atan(1.1608*r1/r2)-atan(1.2395*r3/r2)
solve(fun==0,n)
这个方程解不出来,但我画图显示是有解的,请问问题出在哪?而且如果变成fun=r2*H-atan（r1/r2)-atan（r3/r2)就可以求解.
clear
syms n
H=1.60;
lanmuda=1.55;
n1=1.55;
n2=1.67;
n3=1.5;
k0=2*pi/lanmuda;
c1=n2^2/n1^2;
c2=n2^2/n3^2;
r1=sqrt(n.^2-n1^2)*k0;
r2=sqrt(n2^2-n.^2)*k0;
r3=sqrt(n.^2-n3^2)*k0;
fun=r2*H-atan(1.1608*r1/r2)-atan(1.2395*r3/r2)
solve(fun==0,n)
这个方程解不出来,但我画图显示是有解的,请问问题出在哪?而且如果变成fun=r2*H-atan（r1/r2)-atan（r3/r2)就可以求解.
fun =
(64*pi*(27889/10000 - n^2)^(1/2))/31 - atan((2479*(n^2 - 9/4)^(1/2))/(2000*(27889/10000 - n^2)^(1/2))) - atan((1451*(n^2 - 961/400)^(1/2))/(1250*(27889/10000 - n^2)^(1/2)))
ans =
-1.6344863105805439966305079859321
这是我运行出来的结果
(64*pi*(27889/10000 - n^2)^(1/2))/31 - atan((2479*(n^2 - 9/4)^(1/2))/(2000*(27889/10000 - n^2)^(1/2))) - atan((1451*(n^2 - 961/400)^(1/2))/(1250*(27889/10000 - n^2)^(1/2)))
ans =
-1.6344863105805439966305079859321
这是我运行出来的结果