2PSK数字仿真MATLAB代码与详解（误码率分析）

2PSK数字仿真MATLAB代码与详解

本篇博客将介绍如何使用MATLAB 2020编写一个简单的2PSK数字仿真程序，并对代码进行详细解释。

代码示例

% 设置仿真参数
fs = 1000; % 采样率
T = 1/fs; % 采样周期
numBits = 1000; % 传输比特数
EbNo = 10; % 信噪比（以分贝为单位）

% 生成随机二进制数据
bits = randi([0 1], 1, numBits);

% 将二进制数据转换为PSK调制信号
symbols = 2*bits - 1; % 映射为-1和+1
t = 0:T:(numBits*T)-T; % 时间向量
carrierFreq = 10; % 载波频率
phase = 2*pi*carrierFreq*t; % 载波相位
pskSignal = symbols .* exp(1j*phase); % PSK信号

% 添加高斯噪声
snr = 10^(EbNo/10); % 信噪比转换为线性比例
noiseVar = 1/(2*snr); % 噪声方差
noise = sqrt(noiseVar)*(randn(1, numBits) + 1j*randn(1, numBits)); % 高斯噪声
receivedSignal = pskSignal + noise; % 接收信号

% 解调接收信号
receivedSymbols = real(receivedSignal); % 只提取实部
decodedBits = receivedSymbols > 0; % 查找两个可能的符号之间的边界

% 计算误码率
numErrors = sum(bits ~= decodedBits);
ber = numErrors/numBits;

% 显示结果
disp(['信噪比（Eb/No）：', num2str(EbNo), ' dB']);
disp(['误码率：', num2str(ber)]);

代码详解

设置仿真参数
- fs: 采样率，决定了信号的采样频率。
- T: 采样周期，是采样率的倒数。
- numBits: 传输比特数，表示要传输的二进制数据的位数。
- EbNo: 信噪比，以分贝（dB）为单位，表示信号能量与噪声能量的比值。
生成随机二进制数据
- randi([0 1], 1, numBits): 生成一个大小为 1 x numBits 的矩阵，其中包含0和1的随机数，代表要传输的二进制数据。
将二进制数据转换为PSK调制信号
- symbols = 2*bits - 1: 将二进制数据（0和1）映射为PSK符号（-1和+1）。
- t = 0:T:(numBits*T)-T: 生成时间向量，用于表示信号的时间轴。
- carrierFreq: 载波频率，决定了PSK信号的载波频率。
- phase = 2*pi*carrierFreq*t: 计算载波相位。
- pskSignal = symbols .* exp(1j*phase): 生成2PSK信号，将符号与载波相乘。
添加高斯噪声
- snr = 10^(EbNo/10): 将信噪比从分贝转换为线性比例。
- noiseVar = 1/(2*snr): 计算噪声方差。
- noise = sqrt(noiseVar)*(randn(1, numBits) + 1j*randn(1, numBits)): 生成高斯白噪声，其均值为0，方差为noiseVar。
- receivedSignal = pskSignal + noise: 将噪声添加到2PSK信号中，模拟接收到的信号。
解调接收信号
- receivedSymbols = real(receivedSignal): 提取接收信号的实部，因为PSK信号的相位信息包含在实部中。
- decodedBits = receivedSymbols > 0: 通过比较接收符号与0的大小来解调接收到的数据，大于0则解码为1，否则解码为0。
计算误码率
- numErrors = sum(bits ~= decodedBits): 计算传输比特与解码比特之间的误码数。
- ber = numErrors/numBits: 计算误码率，即误码数与总传输比特数之比。
显示结果
- 显示信噪比和计算得到的误码率。

总结

这段代码演示了如何使用MATLAB进行2PSK数字仿真的基本步骤，包括信号生成、调制、添加噪声、解调和误码率计算。你可以根据需要修改代码中的参数，例如信噪比、载波频率等，以观察它们对误码率的影响。