使用MATLAB软件中函数impseq.m及stepseq.m产生单位脉冲序列及单位阶跃序列，并分别绘制出序列时移m=0、m= -1及m=2时单位脉冲序列δ(n-m)及单位阶跃序列u(n-m)波形的代码

在MATLAB中，`impseq.m` 和 `stepseq.m` 都是内置函数，用于生成离散时间的单位脉冲序列（Impulse sequence，也称作Dirac delta sequence）和单位阶跃序列（Step sequence）。这两个函数通常在数字信号处理（DSP）实验或理论分析中使用。

以下是创建单位脉冲序列（δ(n-m)）和单位阶跃序列（u(n-m)）并展示它们在时移 `m` 分别为0、-1和2的情况下的代码示例：

% 导入必要的库
clear all;
close all;
clc;

% 定义序列长度和时延
n = 100; % 序列长度
m = [0, -1, 2]; % 时延向量

% 单位脉冲序列（Impulse sequence）
imp_seq = impseq(n);
for i = 1:length(m)
    shifted_impulse = imp_seq .* (delta(n-i)); % 时移操作
    stem(shifted_impulse, 'b', 'filled');
    title(sprintf('Impulse Sequence at Shift = %d', m(i)));
    xlabel('Sample Index');
    ylabel('Amplitude');
    hold on;
end
hold off;

% 单位阶跃序列（Step sequence）
step_seq = stepseq(n);
for i = 1:length(m)
    shifted_step = step_seq .* (Heaviside(n-i)); % 时移操作，Heaviside 函数在MATLAB中不是内置，需要自定义
    plot(shifted_step, 'r', 'LineWidth', 2);
    title(sprintf('Step Sequence at Shift = %d', m(i)));
   xlabel('Sample Index');
    ylabel('Amplitude');
    legend('Impulse', 'Shifted Impulse', 'Step', 'Shifted Step');
    hold on;
end
hold off;

注意，上面的代码中有一个假设，即你需要自己编写一个 `Heaviside` 函数，因为MATLAB的`heaviside`函数并不是内置的。这个函数通常定义为 `Heaviside(x) = x > 0;`。确保在运行此代码之前添加这个函数。