Question

MATLAB用for循环解个方程

syms n i q
H=1.8;lanmuda=1.4;n1=1.55;n2=1.67;n3=1.50;
k0=2*pi/lanmuda;
for q=0:2
i=q;
Si=solve(sqrt(n2*n2-n^2)*k0*H-atan(sqrt(n^2-n1*n1)/sqrt(n2*n2-n^2))-atan(sqrt(n^2-n3*n3)/sqrt(n2*n2-n^2))-q*pi==0,n);
end
主要是是想求在q=0,1,2,3.....下的不同的n值，循环程序不对，求帮忙修改。

Answer 1

把“solve”那行里最后的==0去掉即可。

追问

请问下如何让输出结果是S0，S1,S2......

追答

最经济的是用数组S(i)=...，当然你也可以用eval函数。例如
eval(['S', num2str(i), '=solve(sqrt(n2*n2-n^2)*k0*H-atan(sqrt(n^2-n1*n1)/sqrt(n2*n2-n^2))-atan(sqrt(n^2-n3*n3)/sqrt(n2*n2-n^2))-q*pi,n);']);

追问

for q=0:1
i=q;
eval(['S', num2str(i), '=solve(sqrt(n2*n2-n^2)*k0*H-atan(sqrt(n^2-n1*n1)/sqrt(n2*n2-n^2))-atan(sqrt(n^2-n3*n3)/sqrt(n2*n2-n^2))-q*pi,n);']);
end
没有输出啊！回车后什么也没有，

追答

你是说输入S0或S1回车后什么都没有？

追问

好了，主要是不太懂，所以问的比较繁琐，非常感谢你的耐心解答。

追答

不客气。