怎样得出夫琅和费多封衍射图形，并能表示缺级现象。

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('衍射光强空间分布'); 

原代码如上
是否能通过哪种方法得到下图

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答

本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。

​​​​因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。

您好，我是有问必答小助手，您的问题已经有小伙伴帮您解答，感谢您对有问必答的支持与关注！
PS：问答VIP年卡 【限时加赠：IT技术图书免费领】，了解详情>>> https://vip.csdn.net/askvip?utm_source=1146287632