matlab实现一个函数，对灰度级图像按系数n进行子采样，n是2的倍数。该函数应该能够在水平和垂直方向独立地进行子采样，或者同时在两个方向进行子采样。 显示使用以下因子对图像Lena.ppm进行子采样的结果： 2个水平、2个垂直、2个垂直和8个水平、4个垂直和4个水平。在提交的文件中包括结果的文件

好的，我可以帮你实现这个函数。以下是MATLAB代码：

function [subsampled] = subsample_image(image, n, dim)
% INPUTS:
%   image: 输入的灰度级图像
%   n: 子采样系数，必须是2的倍数
%   dim: 子采样方向，可以是 'h'、'v' 或 'hv'
% OUTPUTS:
%   subsampled: 子采样后的图像

% 检查输入参数
if mod(n,2) ~= 0
    error('子采样系数必须是2的倍数！');
end

if ~(strcmp(dim,'h') || strcmp(dim,'v') || strcmp(dim,'hv'))
    error('无效的子采样方向！');
end

% 获取图像尺寸
[height, width] = size(image);

% 计算子采样后的图像尺寸
new_height = height / n;
new_width = width / n;

% 初始化子采样后的图像
subsampled = zeros(new_height, new_width);

% 按指定方向进行子采样
if strcmp(dim,'h')
    % 只在水平方向进行子采样
    for i = 1:new_height
        for j = 1:new_width
            subsampled(i,j) = image(n*i, n*j);
        end
    end
    
elseif strcmp(dim,'v')
    % 只在垂直方向进行子采样
    for i = 1:new_height
        for j = 1:new_width
            subsampled(i,j) = image(n*i, n*j);
        end
    end
    
else
    % 同时在水平和垂直方向进行子采样
    for i = 1:new_height
        for j = 1:new_width
            subsampled(i,j) = image(n*i, n*j);
        end
    end
    
end

end

你可以使用以下代码来测试该函数，并生成所需的图像：

% 读取图像
image = imread('Lena.ppm');
image = rgb2gray(image);

% 进行子采样并显示图像
subsampled1 = subsample_image(image, 2, 'h');
subsampled2 = subsample_image(image, 2, 'v');
subsampled3 = subsample_image(image, 2, 'hv');
subsampled4 = subsample_image(image, 8, 'h');
subsampled5 = subsample_image(image, 4, 'v');
subsampled6 = subsample_image(image, 4, 'hv');

figure;
subplot(2,3,1);
imshow(subsampled1);
title('2个水平');

subplot(2,3,2);
imshow(subsampled2);
title('2个垂直');

subplot(2,3,3);
imshow(subsampled3);
title('2个水平和2个垂直');

subplot(2,3,4);
imshow(subsampled4);
title('8个水平');

subplot(2,3,5);
imshow(subsampled5);
title('4个垂直');

subplot(2,3,6);
imshow(subsampled6);
title('4个水平和4个垂直');

这将生成一个包含6个子图的图像，每个子图显示使用不同的子采样因子进行子采样后的图像。你可以将其保存为一个文件。