用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 。

首先，我们需要定义啁啾光栅的传输矩阵。啁啾光栅的传输矩阵可以表示为： [T_{grating} = \begin{bmatrix} \cos(\beta L) & -i\sin(\beta L)/\beta \ -i\beta\sin(\beta L) & \cos(\beta L) \end{bmatrix}] 其中，(\beta)为啁啾系数，(L)为啁啾区域的长度。

然后，我们可以使用MATLAB的plot函数绘制反射率随波长变化的曲线图。具体步骤如下：

1. 定义啁啾光栅的参数：

beta = [7e-9, 10e-9, 13e-9, 16e-9, 19e-9]; % 啁啾系数
L = 1; % 啁啾区域的长度

2. 定义波长范围和步长：

wavelength = linspace(1549, 1551, 1000); % 波长范围

3. 计算反射率：

reflection = zeros(length(beta), length(wavelength)); % 反射率矩阵，每一行对应一个啁啾系数
for i = 1:length(beta)
    T_grating = [cos(beta(i)*L), -1i*sin(beta(i)*L)/beta(i); -1i*beta(i)*sin(beta(i)*L), cos(beta(i)*L)]; % 计算传输矩阵
    for j = 1:length(wavelength)
        k = 2*pi/wavelength(j); % 波矢
        T = [exp(1i*k*L), 0; 0, exp(-1i*k*L)]; % 光纤段的传输矩阵
        M = T_grating * T; % 总的传输矩阵
        reflection(i, j) = abs(M(2, 1)/M(1, 1))^2; % 计算反射率
    end
end

4. 绘制反射率随波长变化的曲线图：

figure;
plot(wavelength, reflection);
xlabel('波长 (nm)');
ylabel('反射率');
legend(cellstr(num2str(beta', '%.1e')));
title('啁啾光纤布拉格光栅的反射谱');

运行以上代码，即可得到反射率随波长变化的曲线图，图中每条曲线对应一个啁啾系数。图像的横坐标为波长，纵坐标为反射率，图例表示不同的啁啾系数。

注意：以上代码中的参数和计算方式仅为示例，实际应用中需要根据具体情况进行调整