帮我写一份matlab的代码 要求内容灰度图像膨胀、腐蚀、开运算、闭运算 （1）读取一张灰度图像，对灰度图像进行膨胀、腐蚀、开运算、闭运算。 （2）验证开和闭运算具有幂等性。 将灰度图像和四个运算结果图像在同一个窗口带标题显示；将一次开运算、二次开运算、一次闭运算、二次闭运算的图像在另外一个窗口带标题显示。

以下是您需要的 MATLAB 代码：

% 读取灰度图像
img = imread('your_image.jpg');
gray_img = rgb2gray(img);

% 定义结构元素
se = strel('square', 3);

% 膨胀
dilated_img = imdilate(gray_img, se);

% 腐蚀
eroded_img = imerode(gray_img, se);

% 开运算
opened_img = imopen(gray_img, se);

% 闭运算
closed_img = imclose(gray_img, se);

% 将原始图像和四个运算结果图像在同一个窗口带标题显示
figure(1);
subplot(2,3,1); imshow(gray_img); title('Gray Image');
subplot(2,3,2); imshow(dilated_img); title('Dilated Image');
subplot(2,3,3); imshow(eroded_img); title('Eroded Image');
subplot(2,3,4); imshow(opened_img); title('Opened Image');
subplot(2,3,5); imshow(closed_img); title('Closed Image');

% 验证开和闭运算具有幂等性
opened_twice = imopen(opened_img, se);
closed_twice = imclose(closed_img, se);

% 将一次开运算、二次开运算、一次闭运算、二次闭运算的图像在另外一个窗口带标题显示
figure(2);
subplot(2,2,1); imshow(opened_img); title('One-time Opened Image');
subplot(2,2,2); imshow(opened_twice); title('Twice Opened Image');
subplot(2,2,3); imshow(closed_img); title('One-time Closed Image');
subplot(2,2,4); imshow(closed_twice); title('Twice Closed Image');

希望这份代码可以帮助到您。