给出一个图像，用MATLAB编写实现图像全局直方图、局部直方图统计程序、基于直方图的全局均值及方差计算程序、局部方差及方差计算程序、选择性局部图像增强程序

在MATLAB中，要实现上述图像处理功能，你可以按照以下步骤编写相应的函数：

1. **图像全局直方图**:
使用`imhist()`函数可以获取整个图像的直方图，例如：

   img = imread('your_image.jpg'); % 加载图像
   global_hist = imhist(img, 'Normalization', 'probability');

2. **局部直方图统计**:
可以通过循环遍历图像的窗口或感兴趣区域，使用`imlocalmin`等工具提取局部像素，再分别计算直方图：

   window_size = [50, 50]; % 窗口大小
   local_hists = zeros(size(img));
   for i = 1:size(img, 1) - window_size(1)
       for j = 1:size(img, 2) - window_size(2)
           local_window = img(i:i+window_size(1)-1, j:j+window_size(2)-1);
           local_hists(i:i+window_size(1)-1, j:j+window_size(2)-1) = histcounts(local_window);
       end
   end

3. **全局均值及方差计算**:
使用`mean2()`和`var2()`函数：

   global_mean = mean2(img(:));
   global_var = var2(img(:));

4. **局部方差及方差计算**:
类似地，从局部窗口计算：

   local_vars = std2(local_windows(:), 'all')^2;

5. **选择性局部图像增强**:
使用`imadjust()`函数可以根据局部直方图调整对比度，例如：

   enhanced_local_img = imadjust(local_windows, [0.2; 0.8], []);

这里 `[0.2; 0.8]` 表示调整后的灰度范围，`[]` 表示使用线性变换。

以上代码仅为概述，实际编写时需要根据具体需求调整参数，并可能需要对结果进行可视化。