matlab heaviside函数生成矩形脉冲

意思就是下面这张图

其“微分”是狄拉克 δ 函数，事实上，x=0 的值在函数应用上并不重要，可以任意取。在matlab中，x=0取值为1/2。

1. 产生一个信号 u(t)

syms x
fplot(heaviside(x));

之所以看到在0的地方有竖线，就是因为fplot函数把断开部分连接起来了。

2 产生一个 u(t-1) 的信号

3. 产生一个矩形脉冲

clear all;
syms t;
x = heaviside(t + 1/4) - heaviside(t - 1/4);
fplot(x);

4. 产生一个矩形脉冲序列

clear all;%一定要加这行代码，否则可能你之前的代码能运行生成正确图形，改错之后，运行结果还是之前的图形，每次运行前必须先做清理
syms t;
T = 1;%矩形脉冲的周期
N = 5;%每个for循环生成的矩形脉冲数

%生成[T-1/4*T, N*T - 1/4*T]区间内的信号
for k = 1: N,
    x(k) = heaviside(t - T * k + T/4) - heaviside(t - T * k - T/4);
end

%生成[-N*T-1/4*T， T-1/4*T]区间内的信号
for k = N: -1:1,
    x2(k) = heaviside(t + T * k + T/4) - heaviside(t + T * k - T/4);
end

%生成[T-1/4*T, T-1/4*T]区间内的信号
x3 = heaviside(t  + T/4) - heaviside(t - T/4);

%把以上信号都整合起来，然后幅值乘以2
f = 2 * [x x2 x3];

fplot(f, [-10,10]);