function[x,n]=delta(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0]; %定义函数1 function[x,n]=step(n0,n1,n2) n=[n1:n2]; x=[(n-n0)>=0]; %定义函数2 [x1,n]=delta(3,1,10); [x2,n]=step(-2,1,10); [x3,n]=delta(-3,1,10); [x4,n]=step(2,1,10); subplot(2,3,1);stem(n,(x1+x2));title('f1(t)');grid on subplot(2,3,2);stem(n,(x3+x4));title('f2(t)');grid on; subplot(2,3,3);stem(n,(x1+x2)+(x3+x4));title('f1+f2');grid on; subplot(2,3,4);stem(n,(x1+x2)-(x3+x4));title('f1-f2');grid on; subplot(2,3,5);stem(n,(x1+x2).*(x3+x4));title('f1*f2');grid on;

这段代码定义了两个函数delta和step，并利用这两个函数绘制了5个离散时间信号的图像。其中，delta函数表示一个以n0为中心的单位脉冲信号，step函数表示一个从n0开始的单位阶跃信号。代码中的5个信号分别为：

1. f1(t) = delta[n-3] + u[n+2]，其中u[n]表示从n=0开始的单位阶跃信号。
2. f2(t) = delta[n+3] - u[n-2]。
3. f1(t) + f2(t)
4. f1(t) - f2(t)
5. f1(t) * f2(t)

subplot函数用于将这5个信号的图像绘制在同一张图中，方便对比和分析。每个图像的标题表示对应的信号名称，grid on命令用于显示图像的网格线。注意，这段代码中的变量名和函数名不能与MATLAB的内置函数重名，否则会出现错误。