利用Simulink的模块，搭建MPSK/MQAM经过AWGN信道的传输模型，推导系统的理论性能，画出比特信噪比EbNo与误码率BER（包括仿真与理论值）的关系曲线，其中M取2,4,8,16；请提供相应代码

利用Simulink搭建MPSK (Modified Phase Shift Keying) 和 MQAM (M-QAM) 经过AWGN (Additive White Gaussian Noise) 信道的传输模型，通常涉及以下几个步骤：

1. **创建系统模块**:
- 创建MPSK或MQAM调制器模块，如`CommModulator`或自定义函数。
- 创建AWGN噪声生成模块，如`awgn`。

2. **设置参数**:
- 设置MPSK的阶数（例如，M=2对应BPSK，M=4对应QPSK等）。
- 设定信噪比（SNR或EbNo）。

3. **构建通信链路**:
- 将调制器模块连接到AWGN模块，然后连接到解调器模块，如`CommDemodulator`。

4. **计算理论性能**:
- 使用香农公式推导在高斯白噪声环境下的信息传输速率和误码率（BER）理论值。

   EbNo = ...; % 信噪比
   theoretical_BER = 10^(-EbNo * shannon_theorem( modulation_order, snr ) / 10); % 假设shannon_theorem是自定义函数

5. **仿真与绘制结果**:
- 进行多次仿真实验，改变信噪比并记录BER。

   EbNos = ...; % 一系列信噪比点
   BERs_simulated = zeros(size(EbNos)); % 初始化存储仿真结果
   for i = 1:length(EbNos)
       EbNo = EbNos(i);
       % 仿真过程...
       BERs_simulated(i) = your_simulation_function(EbNo, M);
   end
   
   plot(EbNos, BERs_simulated, 'o', EbNos, theoretical_BERs, '-');
   xlabel('EbNo (dB)');
   ylabel('BER');
   legend('Simulation', 'Theory');

请注意，由于这是一个较为复杂的Simulink模型设计，具体的Matlab代码需要你自己编写，并且上述步骤可能需要结合MPSK和MQAM的具体数学模型（比如相位旋转、星座图等）。如果你需要进一步的帮助，可以提供更详细的需求或者具体的问题，以便我能给出更精确的指导。