2．用Matlab 程序对DSB调制、相干解调系统进行仿真建模。 其中：解调器所用载波要求用两种方法得到：① 理想载波，即从发送端直接传送给接收端；②用载波提取电路恢复得到（原理参考教材3.6节）选做。 要求：测试调制前后信号谱频，测试各输出点的波形。

要使用MATLAB对DSB（双边带）调制和相干解调系统进行仿真建模，你可以按照以下步骤进行：

1. **导入必要的库**：

   % 加载基本图形和通信工具箱
   addpath('toolbox'); % 假设toolbox目录包含通信相关的工具函数

2. **定义信号参数**：
- 选择信号类型（模拟音频或数字数据）
- 设定采样频率（fs）
- 定义信息比特率（Baud rate）

   signal = 'analog'; % 'analog' 或 'digital'
   fs = 44100; % 例如，如果用于音频信号，则取常见的采样频率
   baud_rate = 1000; % 比特率

3. **生成载波**：
- 理想载波（理想情况下无失真地传输）

   if strcmp(signal, 'analog')
       carrier_freq = 1000; % 频率（Hz）
       ideal_carrier = sin(2*pi*carrier_freq/fs * (0:length(fs)-1));
   else
       % 对于数字信号，可能需要创建正交相位载波
       [carrier, ~] = mseq(baud_rate, fs);
       ideal_carrier = cos(carrier);
   end

- 使用载波提取电路恢复载波（参阅教材3.6节，可能涉及到滤波和同步操作）

   extracted_carrier = some_recovery_circuit(ideal_carrier, noise_signal); % 假设noise_signal代表噪声或其他影响因素

4. **DSB调制**：

   if strcmp(signal, 'analog')
       message = randi([0 1], length(fs)/baud_rate, 1); % 产生随机二进制序列
       dsb_modulated = ideal_carrier .* message;
   else
       dsb_modulated = qammod(message, 4, 'PSK', 'BitOutput', true); % 对于数字信号，用QAM调制
   end

5. **显示信号谱分析**：

   subplot(2,1,1)
   plot(freqz(dsb_modulated, 1, fs))
   title('Signal Spectrum After DSB Modulation')

   subplot(2,1,2)
   plot(t, abs(fftshift(fft(dsb_modulated))'))
   xlabel('Frequency (Hz)')
   ylabel('Magnitude')

6. **相干解调**：

   demodulated_message = corrcoeff(dsb_modulated, extracted_carrier); % 用相关系数解调

7. **输出波形和测试结果**：

   time = (0:length(fs)-1)'/fs;
   subplot(2,1,1)
   hold on
   plot(time, ideal_carrier, 'r', time, dsb_modulated, 'b')
   legend('Ideal Carrier', 'Modulated Signal')
   hold off

   subplot(2,1,2)
   plot(time, demodulated_message)
   title('Demodulated Waveform')

8. **评估性能**：
- 分析误码率（BER），如果信号是数字的

   if strcmp(signal, 'digital')
       ber = berawgn(demodulated_message, EbNodB, 'qam'); % EbNodB表示信噪比
       fprintf('BER: %f\n', ber)
   end

9. **清理工作**：

   clearvars -except signal fs baud_rate carrier ideal_carrier extracted_carrier dsmodulated ...
               demodulated_message t time ber EbNodB % 清除内存中的变量