用matlab写代码实现的信号包含了一个20Hz的简谐波，并在0.5s和0.52s 处叠加了两个50Hz的雷克子波，在0.11s和0.29s之间叠加了两个30Hz的子波， 并在0.81s出叠加了一个100Hz的子波

以下是实现该信号的 MATLAB 代码：

% 设置采样率和时间向量
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量

% 生成20Hz的简谐波
x = sin(2*pi*20*t);

% 叠加两个50Hz的雷克子波
r1 = rectpuls(t-0.5, 0.02); % 第一个雷克子波
r2 = rectpuls(t-0.52, 0.02); % 第二个雷克子波
x = x + 3*sin(2*pi*50*t).*r1 + 2*sin(2*pi*50*t).*r2;

% 叠加两个30Hz的子波
s1 = sin(2*pi*30*t); % 第一个子波
s2 = sin(2*pi*30*t); % 第二个子波
s1(0.11*fs:0.13*fs) = s1(0.11*fs:0.13*fs) + 2; % 在0.11s处叠加
s2(0.29*fs:0.31*fs) = s2(0.29*fs:0.31*fs) + 1; % 在0.29s处叠加
x = x + s1 + s2;

% 在0.81s处叠加一个100Hz的子波
x(0.81*fs:end) = x(0.81*fs:end) + 5*sin(2*pi*100*t(0.81*fs:end));

% 绘制信号波形图
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Signal with Harmonic and Subharmonic Components');

执行代码后，将会生成一个包含20Hz简谐波和多个子波的信号波形图。