load ('e1.mat')
syms y
d=1.8;%input('请输入芯层厚度:');
for i=0:N
q=i;
n=teneff1{i+1,1}
r1=sqrt(n^2-n1^2)*k0;
r2=sqrt(n2^2-n^2)*k0;
r3=sqrt(n^2-n3^2)*k0;
a=-atan(r3/r2)+q*pi;
if y>d
EY1=cos(r2*d+a)*exp(-r1*(y-d))
else if 0<y<d
EY1=cos(r2*y+a);
else y<0
EY1=cos(a)*exp(r3*y)
end
end
end
这个程序是分段函数,运行时d有问题:
无法从 sym 转换为 logical。
出错 EYy (line 12)
if y>d
拉莫斯之舞