怎样得出夫琅禾费多缝衍射缺级现象图形？

clear,clc;
Lambda = 632.8e-9;
a = 0.1e-3;                 %光栅狭缝宽度单位转为：m
d = 0.3e-3;                 %光栅缝距长度单位转为：m
N = 3;      %缝数
L = 1;      %光栅到屏的距离
I0 = 1;

Screen_Xmax = 2*Lambda*L/a;                  %接收屏x方向位置(边界)--单位为：m
Screen_x = linspace(-Screen_Xmax,Screen_Xmax,501);  %屏x坐标取样501份

sin_angle = Screen_x./L;
alpha = pi*a/Lambda.*sin_angle;      %公式里的其中一个相位差
deta = 2*pi*a/Lambda.*sin_angle;     %公式里的另一个相位差
I = I0* (sin(alpha)./alpha).^2 .* (sin(N/2.*deta)./sin(deta./2)).^2;
I = I/max(I);

% 衍射光强函数分布图*************************
subplot(2,1,2);
plot(Screen_x,I);
axis([-Screen_Xmax,Screen_Xmax,0,1]);
xlabel('空间坐标x'), ylabel('衍射光强度');
title('多缝衍射光强空间分布函数');
grid on;

% 衍射光强灰度分布图*************************
I = I*1024;                          %光强归一,扩大显示
subplot(2,1,1);
image(Screen_x,Screen_Xmax,I);       %设置x和y的像素，显示数值
colormap gray;                       %妆点色彩
xlabel('空间坐标x');
title('衍射光强空间分布'); 

原代码如上

能否根据原代码改成下图

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，目前超出我们的服务范围，暂时无法为您解答。

首次提问人员可免费体验一次有问必答服务。目前首次提问的问题服务范围为：编程语言、Java开发、python、数据库、前端开发 领域专业技术问题，为您提供问题的解决思路和指导。不提供源码代写、项目文档代写、论文代写、作业代写、安装包资源发送或安装、软件使用指导等服务。
我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。