1.读入dem文件( dem_local.grd )，得到高程矩阵zmat 2.构造同尺寸的高斯噪声矩阵，标准差为1.5，叠加至zmat3用pcolor函数生成伪彩色地形图 4.编写一个局域子窗滤波函数(边缘位置数值不变，其它每个位置，计算周围9个点高程的中位数，作为该位置的输出值) 5.以zmat为函数输入，得到输出zmat1，成图对比处理前后效果

好的，让我来逐步回答您的问题。

1. 读入dem文件( dem_local.grd )，得到高程矩阵zmat

可以使用Matlab中的`gridread()`函数来读取.grd文件，将其转化为高程矩阵zmat。示例代码如下：

[zmat, refvec] = gridread('dem_local.grd');

其中，`zmat`为高程矩阵，`refvec`为参考向量。

2. 构造同尺寸的高斯噪声矩阵，标准差为1.5，叠加至zmat

可以使用Matlab中的`randn()`函数生成同尺寸的高斯噪声矩阵，标准差为1.5，然后将其加到zmat中即可。示例代码如下：

noise = 1.5 * randn(size(zmat));
zmat = zmat + noise;

3. 用pcolor函数生成伪彩色地形图

可以使用Matlab中的`pcolor()`函数生成伪彩色地形图。示例代码如下：

figure;
pcolor(zmat);
shading interp;
colorbar;

4. 编写一个局域子窗滤波函数

可以编写一个局域子窗滤波函数，以处理高程矩阵zmat。示例代码如下：

function zmat1 = local_filter(zmat)
    zmat1 = zmat;
    [m, n] = size(zmat);
    for i = 2:m-1
        for j = 2:n-1
            window = zmat(i-1:i+1, j-1:j+1);
            zmat1(i, j) = median(window(:));
        end
    end
end

5. 以zmat为函数输入，得到输出zmat1，成图对比处理前后效果

可以调用局域子窗滤波函数，对高程矩阵zmat进行滤波处理，得到输出矩阵zmat1。然后，可以使用`subplot()`函数将处理前后的伪彩色地形图进行对比显示。示例代码如下：

zmat1 = local_filter(zmat);

figure;
subplot(1, 2, 1);
pcolor(zmat);
shading interp;
colorbar;
title('Original');

subplot(1, 2, 2);
pcolor(zmat1);
shading interp;
colorbar;
title('Filtered');

希望这些回答可以帮助到您！