在matlab中表示输出的参考信号是振幅为±1、周期为4秒的方波信号，而参考信号为零。在半个模拟时间进入模拟的控制信号中添加一个单位步长输入干扰信号

在MATLAB中，你可以使用`square`函数生成一个振幅为±1的方波信号，其周期可通过`period`确定，比如一个周期为4秒的方波可以这样创建：

Ts = 0.5; % 模拟采样时间
T = 4;   % 方波周期
f = 1/T; % 频率

% 创建一个长度等于模拟时间（半个周期）的方波信号
time_vec = 0:Ts:(T-Ts); % 时间向量
ref_signal = square(2*pi*f*time_vec) - 1; % ±1幅度的方波，中心为0

这里假设模拟时间是半个周期，即从0到2秒。

对于单位步长输入干扰信号，MATLAB有一个内置函数`step`，它通常用于离散系统的阶跃响应分析。如果你想在半个周期的时间点加入一个单位步长干扰，可以在相应位置添加：

扰动_time = T/2; % 单位步长干扰的位置
unit_step = step(time_vec >=扰动_time); % 对应时间点为1，其他时间为0

% 如果希望在整个时间向量上平滑地引入干扰，可以用Heaviside函数
heaviside = heaviside(time_vec -扰动_time);

% 将干扰信号与原始信号相加得到带干扰的信号
total_signal = ref_signal + unit_step * 1; % 或者 total_signal = ref_signal + (1 - ref_signal).*heaviside;

这样就得到了一个初始为0的方波信号，然后在半个周期处添加了单位步长的干扰。