帮我用MATLAB编程利用传输矩阵法画出啁啾光栅在设置不同的啁啾系数7e-9、10e-9、13e-9 、16e-9、19e-9时的反射光谱图要求绘制的图的横坐标波长在1549 nm 1551 nm范围纵坐标反射率设置在01范围。要求：加图名啁啾光纤布拉格光栅的反射谱、横坐标注释波长 nm、纵坐标注释反射率 、图例 7e-9、10e-9、13e-9 、16e-9 、19e-9 。

下面是一个使用MATLAB编程绘制啁啾光栅反射光谱图的示例代码：

% 设置波长范围和步长
wavelength = linspace(1549, 1551, 1000);

% 设置啁啾系数
chirp_coefficients = [7e-9, 10e-9, 13e-9, 16e-9, 19e-9];

% 初始化反射率矩阵
reflection_spectrum = zeros(length(wavelength), length(chirp_coefficients));

% 计算反射率
for i = 1:length(chirp_coefficients)
    chirp = chirp_coefficients(i);
    
    % 利用传输矩阵法计算反射率
    % 在这里进行你的计算，将结果存储在reflection_spectrum矩阵中
    
    % 这里只是一个示例，你需要根据你的具体计算方法进行修改
    reflection_spectrum(:, i) = sin(chirp * wavelength);
end

% 绘制反射光谱图
figure;
plot(wavelength, reflection_spectrum);

% 设置图名、横纵坐标注释和图例
title('啁啾光纤布拉格光栅的反射谱');
xlabel('波长 (nm)');
ylabel('反射率');
legend('7e-9', '10e-9', '13e-9', '16e-9', '19e-9');

请根据你的具体需求进行修改计算反射率的部分