基于MATLAB的PAM通信系统仿真实现

一、系统架构设计

PAM（脉冲幅度调制）通信系统仿真包含以下核心模块：

二、关键代码实现

1. 参数设置与信号生成

%% 系统参数设置
nSym = 1e5;          % 符号数
M = 4;               % 调制阶数（2/4/8/16）
EbN0dB = 0:2:20;     % 信噪比范围(dB)
k = log2(M);         % 每符号比特数

% 生成随机信息序列（Gray编码）
msg = randi([0 M-1],1,nSym);
grayCode = pam_gray_map(M);  % Gray映射表
txSym = grayCode(msg+1);     % 映射为符号

2. PAM调制实现

%% PAM调制（使用MATLAB内置函数）
s = pammod(txSym, M, 0, 'gray');  % 生成复基带信号
t = (0:length(s)-1)/1e3;           % 时间向量（1kHz采样率）

% 基带脉冲成形（矩形脉冲）
pulse = ones(1, 10);             % 脉冲宽度10个符号周期
txPulse = conv(s, pulse, 'same');% 卷积成型

3. AWGN信道仿真

%% 信道传输（添加高斯噪声）
snrLinear = 10.^(EbN0dB/10);     % 线性信噪比
EsN0dB = 10*log10(k) + EbN0dB;   % 转换为Es/N0
rxSignal = awgn(real(txPulse), EsN0dB, 'measured');  % 添加噪声

4. PAM解调实现

%% PAM解调
rxSym = pamdemod(real(rxSignal), M, 0, 'gray');  % 解调为符号
ber = sum(msg ~= rxSym)/nSym;                    % 计算误码率

% 理论误码率计算（16-PAM示例）
Eb = 10^(EbN0dB(1)/10);  % 参考Eb值
Es = k*Eb;               % Es值
ser_theory = 2*(1-1/M)*0.5*erfc(sqrt(3*Es/(M^2-1)));  % 理论SER公式

5. 性能可视化

%% 结果可视化
figure;
subplot(2,1,1);
semilogy(EbN0dB, ber, 'bo-', 'LineWidth',1.5);
hold on;
semilogy(EbN0dB, ser_theory, 'r--', 'LineWidth',1.5);
grid on;
xlabel('Eb/N0 (dB)');
ylabel('误码率 (BER)');
legend('仿真结果', '理论曲线');
title('4-PAM误码率性能');

subplot(2,1,2);
scatterplot(s);
title('PAM星座图');
xlabel('实部');
ylabel('虚部');

三、关键算法详解

1. Gray编码映射

function grayCode = pam_gray_map(M)
    % 生成Gray编码表
    if M == 2
        grayCode = [0 1];
    elseif M == 4
        grayCode = [0 1 3 2];
    elseif M == 8
        grayCode = [0 1 3 2 6 7 5 4];
    elseif M == 16
        grayCode = [0 1 3 2 6 7 5 4 12 13 15 14 10 11 9 8];
    end
end

2. 匹配滤波接收

%% 匹配滤波优化（提高信噪比）
pulse = fir1(31, 0.5);  % 设计升余弦滤波器
rxMatched = filter(pulse, 1, rxSignal);  % 匹配滤波

参考代码 利用matlab实现PAM通信系统的仿真 www.youwenfan.com/contentalh/99941.html

四、仿真结果分析

五、扩展应用场景

1. 高速光纤通信

   实现16-PAM高速传输系统：

   M = 16;
   upsampleFactor = 8;  % 提升符号率
   s_up = upsample(s, upsampleFactor);
   

2. 电力线通信

   添加脉冲噪声干扰：

   noise = 0.1*randn(size(s));  % 高斯噪声
   impulsiveNoise = 5*randi([0 1],1,length(s));  // 二进制脉冲噪声
   rxSignal = s + noise + impulsiveNoise;
   

3. 无人机通信

   实现移动信道仿真：

   h = rayleighchan(1/1e3, 30);  % 瑞利信道模型
   rxSignal = filter(h, s);      % 信道衰落