请编写一段matlab程序，运用最大方差阈值法, 实现对jpg格式图像的分割

以下是一段用 Matlab 实现最大方差阈值法对 JPG 格式图像进行分割的程序：

% 读取 JPG 图像
img = imread('example.jpg');

% 将彩色图像转为灰度图像
gray_img = rgb2gray(img);

% 计算图像的方差
img_var = var(double(gray_img(:)));

% 设置阈值范围
threshold_range = 1:255;

% 初始化最大方差和最佳阈值
max_variance = 0;
best_threshold = 0;

% 遍历阈值范围，找到最大方差和最佳阈值
for threshold = threshold_range
    % 将图像分为两部分
    img_part1 = gray_img(gray_img <= threshold);
    img_part2 = gray_img(gray_img > threshold);
    
    % 计算两部分的方差
    part1_var = var(double(img_part1(:)));
    part2_var = var(double(img_part2(:)));
    
    % 计算总方差
    total_var = length(img_part1) / numel(gray_img) * part1_var ...
                + length(img_part2) / numel(gray_img) * part2_var;
            
    % 更新最大方差和最佳阈值
    if total_var > max_variance
        max_variance = total_var;
        best_threshold = threshold;
    end
end

% 根据最佳阈值将图像分割
binary_img = gray_img > best_threshold;

% 显示原图和分割后的二值图像
subplot(1, 2, 1), imshow(img), title('Original Image');
subplot(1, 2, 2), imshow(binary_img), title('Binary Image');

上述程序首先读取一个 JPG 格式的图像，将其转换为灰度图像，并计算其方差。然后设置阈值范围并遍历每个阈值，将图像分为两部分并计算它们的方差，最后计算总方差并更新最大方差和最佳阈值。根据最佳阈值将图像分割为二值图像，并显示分割前后的图像。