啁啾光纤光栅matlab仿真程序

% chirp grating
clc
L=.02;R=1;S=0;
% power reflectoin coefficient
for k=1:1000
   %fringe visibility of induced index change
   v=1;
	% “ac“ index change
	delta=0.0001;
	% efficient induced index
	n=1.45;
	% design wavelength
	lambdaD=（1550+k/1000）*10e-10;
	% wavelength
	lambda=[1549.8:.001:1555.4]*10e-10;
	% “AC“ coupling coefficient
	kai=pi*v*delta./lambdaD;
	% “dc“ coupling coefficient