通过MATLAB程序来仿真一串随机消息代码的基带信号波形，首先产生1000个随机信号序列，用双相码的编码规则进行编码，并且求平均功率谱密度。

在MATLAB中，你可以按照以下步骤来模拟一串随机消息的基带信号波形：

1. **生成随机消息序列**:
使用`randn`函数生成一个包含1000个元素的正态分布随机数列，代表消息的原始数据。

message = randn(1, 1000);

2. **双相编码**:
双相码是一种二进制脉冲编码调制(Pulse Amplitude Modulation, PAM)的一种变体，其中每个比特由两个周期长度的一系列+1和-1组成。你可以创建一个双相码映射，然后对随机消息进行映射。

dual_phase_mapping = [-1 1; -1 -1]; % 例如对于二电平双相码
encoded_message = sign(message) .* dual_phase_mapping;

3. **基带信号波形**:
创建一个时间轴，比如`t = linspace(0, 1, length(encoded_message))`，然后计算对应的基带信号。假设每个符号持续的时间为1秒，频率采样率为每秒多少次（如44100Hz），可以使用`squeeze`函数压缩维度，然后乘以相应的频率信号。

fs = 44100; % 假设采样率
time_axis = t';
baseband_signal = squeeze(time_axis) * encoded_message;

4. **求平均功率谱密度**:
对基带信号进行窗函数（如汉明窗或矩形窗）加窗处理，然后利用`fft`函数计算其频域表示，再取实部和平方，最后除以窗口长度和总的点数，得到平均功率谱密度。

window = hamming(length(encoded_message)); % 使用汉明窗
windowed_signal = baseband_signal .* window;
psd = abs(fft(windowed_signal)) ./ (length(windowed_signal) / fs);

5. **可视化**:
可以绘制基带信号波形以及平均功率谱密度图，以便观察编码效果。

plot(time_axis, baseband_signal);
title('Baseband Signal Waveform');
xlabel('Time (s)');
ylabel('Amplitude');

figure;
spectrogram(psd, [], fs, [], 'yaxis');
title('Average Power Spectral Density');
xlabel('Frequency (Hz)');
ylabel('Power');