能否提供一个用于实现马赫曾德尔干涉仪(Mach-Zehnder Interferometer, MZI)功能的MATLAB代码示例？

Mach-Zehnder干涉仪是一种光学干涉实验装置，常用于测量光的相位差。下面是一个简单的MATLAB代码示例，它模拟了双缝干涉的基本原理，并展示了如何创建一个MZI的仿真。这个例子假设你知道基本的光学波函数操作。

% 导入必要的库
clear all; close all; clc;
import matlab.graphics.*;
import optics.*

% 设定光源
wavelength = 532e-9; % 波长，这里以纳米计
source = Source('laser', 'Wavelength', wavelength);

% 创建光栅
grating = Grating('diffraction', 'GratingPeriod', 0.6e-6, 'Direction', [0, -1]);

% 创建第一个分束镜
BS1 = BeamSplitter('bs', 'R', 0.5); % 反射率为50%

% 创建双缝
slits = Slit('doubleSlit', 'Separation', 0.2e-6);

% 创建第二个分束镜
BS2 = BeamSplitter('bs', 'T', 1); % 穿透率为100%

% 构建干涉路径
path1 = BS1:slits:BS2; % 直通路径
path2 = BS1^' : grating : slits^' : BS2^'; % 干涉路径

% 合并两条路径
interferometer = path1 + path2;

% 计算干涉图案
pattern = propagate(interferometer, source);
image(pattern, 'ColorMethod', 'phase');

% 显示结果
figure;
title('Mach-Zehnder Interferometer Pattern');
xlabel('Distance (m)');
ylabel('Distance (m)');
colorbar;