%---------------sinc切趾光纖光柵反射譜和透射譜----------------

%纖芯折射率1.452纖芯折射率平均增加值7.5*10^（-4）光柵周期535nm切趾函數(shù)參數(shù)1mmk0表示入射光波長
clear;
nn=15000;
a=0.014/nn;
di=a;
i=1;
for?z=0:a:1.4*10^（-2）;
????n（i）=1.452+7.5*10^（-4）*（1+（1+0.18）*cos（2*pi*z/0.014）+0.18*cos（4*pi*z/0.014））/（2+2*0.18）*（1+cos（2*pi*z/（535*10^（-9））））;???
????i=i+1;
end
i=1;
t=1;
for?k0=1550*10^（-9）:0.02*10^（-9）:1557*10^（-9）
????M=[1?0;0?1];
????for?i=1:1:nn
????????M1=[n（i）+n（i+1）n（i）-n（i+1）;n（i）-n（i+1）n（i）+n（i+1）]/（2*n（i））;
????????M2=[exp（j*（（2*pi）/k0*n（i）*di））0;0exp（j*（-1）*（2*pi）/k0*n（i）*di）];
????????M=M*M2*M1;
????end
????r=M（21）/M（11）;
????R（t）=（abs（r））^2;??????%反射率
????e=1/M（11）;
????T（t）=（abs（e））^2;??????%折射率
????t=t+1;
end

k0=1550*10^（-9）:0.02*10^（-9）:1557*10^（-9）;
plot（k0R）
%figure;
%k0=1550*10^（-9）:0.02*10^（-9）:1557*10^（-9）;
%plot（k0T）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????806??2012-12-14?19:53??even_FBG.m

?????文件????????856??2012-12-14?19:06??Gauss_apodized_FBG.m

?????文件????????872??2012-12-14?19:18??Hamming_apodized_FBG.m

?????文件????????857??2012-12-14?19:02??sinc_apodized_FBG.m

?????文件????????870??2012-12-14?19:23??tanh_apodized_FBG.m

?????文件????????901??2012-12-14?19:37??Blackman_apodized_FBG.m

?????文件????????875??2012-12-14?19:36??Cauthy_apodized_FBG.m

-----------?---------??----------?-----??----

?????????????????6037????????????????????7