产生13位巴克码信号，子脉冲长1us，载频2MHz，采样率32MHz。绘制载波、巴克码包络以及调制后信号时域波形;绘制巴克码包络及频谱、相同长度的矩形信号包络及频谱对比图。

在MATLAB中，生成13位巴克码信号（Barker Code）的过程包括以下几个步骤：

1. **定义巴克码序列**:
巴克码是一种特殊的数字序列，它具有低自相关特性。对于13位的Barker码，你需要查找预先存在的代码，比如MATLAB自带的BarkerCodeGenerator函数。

barker_code = barkercode(13);

2. **载波和子脉冲**:
载波频率是2MHz，子脉冲宽度是1us，因此可以创建一个时间序列来表示每个子脉冲的位置。

subpulse_width = 1e-6; % 子脉冲长度 (s)
carrier_freq = 2e6; % 载波频率 (Hz)
t = 0:subpulse_width:1/sample_rate; % 从0到1秒，采样率为32MHz
carrier = cos(2*pi*carrier_freq*t); % 载波信号

3. **巴克码调制**:
将Barker码应用到载波上，生成调制后的信号。

modulated_signal = carrier .* sign(barker_code);

4. **绘制时域波形**:
使用`plot`函数绘制载波、巴克码包络以及调制后信号。

figure;
subplot(2,1,1);
plot(t, carrier, 'r', t, modulated_signal, 'b');
legend('Carrier', 'Modulated Signal');
xlabel('Time (s)');
ylabel('Amplitude');
title('Barker Code Modulation (Time Domain)');

5. **包络和频谱**:
   对调制后的信号和矩形信号分别计算包络和频谱。

% 对比矩形信号（假设是1s持续的矩形，全占空比）
rectangular_signal = ones(size(t));
rectangular_signal envelope_rect = max(abs(rectangular_signal));
rect_envelope = abs(fftshift(fft(rectangular_signal)));
rect_freq = (-sample_rate/2:sample_rate/length(rect_freq):sample_rate/2 - 1)';

   
   绘制两个包络的对比图和它们的频谱图。

subplot(2,1,2);
plot(t, abs(modulated_signal), 'b', t, envelope_rect, 'g');
legend('Barker Code Envelope', 'Rectangular Envelope');
xlabel('Time (s)');
ylabel('Amplitude');

figure;
plot(rect_freq, rect_envelope, 'g', freq_range, abs(fft_domain), 'b');
xlabel('Frequency (Hz)');
ylabel('Magnitude');
legend('Rectangular Signal Spectrum', 'Barker Code Spectrum');