在上面整形信号上加上它的二次谐波，谐波的幅度a1=0.2a，正弦波信号和谐波的初始相位都为0，持续时间为1s，画出波形，并播放出这个信号；做fft，画出幅度谱。

下面是添加二次谐波、绘制波形和播放信号的修改代码：

% 参数设置
fs = 8000;    % 采样频率
f = 440;      % 信号频率
a = 1;        % 幅度
theta = 0;    % 初始相位
duration = 1; % 持续时间

% 生成时间轴
t = 0:1/fs:duration;

% 生成正弦波信号
x = a * sin(2*pi*f*t + deg2rad(theta));

% 添加二次谐波
a1 = 0.2 * a; % 谐波幅度
x_harmonic = a1 * sin(2*pi*2*f*t);

% 合并信号
x_total = x + x_harmonic;

% 绘制波形
subplot(2, 1, 1);
plot(t, x_total);
xlabel('Time (s)');
ylabel('Amplitude');
title('Sinusoidal Waveform with Harmonic');
grid on;

% 播放信号
sound(x_total, fs);

% 做FFT，并绘制幅度谱
subplot(2, 1, 2);
N = length(x_total);
frequencies = (0:fs/N:fs-fs/N);
X = abs(fft(x_total))/N;
plot(frequencies, X);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Amplitude Spectrum');
grid on;

这段代码会在 MATLAB 中生成一个带有二次谐波的正弦波信号的波形图和幅度谱，并通过扬声器播放出来。