L=100000;
n=1.46;
b1=1550;
b0=1549:0.001:1551;
a=(2*n*pi)./b0-(2*n*pi)/b1;
k=(atanh(0.9^0.5))/L;
s=(k*k-a.*a).^0.5;
c=(k./s).*sinh(s*L);
d=(cosh(s*L)).*(cosh(s*L));
e=(a./s).*sinh(s*L)
j=e.*e;
f=abs(d+e);
l=abs(c.*c);
R=l/f;
plot(b0,R);
大家能帮我看下那有错误吗?给我注释下也可以,是关于光栅的程序!~
參考答案:把R=l/f改成R=l./f即可，否则R是一个常数，结果当然是一条直线了。
