对单极性信号进行16PSK调制解调

上篇文章介绍了数字基带信号的产生，这里给出16PSK的调制方法，本文基于MATLAB，使用了较多MATLAB现有函数。

Ts = 1; % 码元周期
N_sample =100; % 单个码元抽样点数
dt = Ts / N_sample; % 抽样时间间隔
N = 100; % 码元数
t = 0 : dt : (N * N_sample - 1) * dt; % 序列传输时间
gt1 = ones(1, N_sample); % NRZ
%% 生成随机序列
RAN = round(rand(1, N)); % 随机0 1序列 N长
NRZ= [];
for i = 1 : N % 生成序列
   if RAN(i)==1
       NRZ = [NRZ gt1];
   else
       NRZ = [NRZ 0*gt1];
   end
end

这里已经产生了信号NRZ

接下来对其进行调制

产生调制信号

使用pskmod函数，定义调制信号，其中pskmod表示进行进行PSK调制，输入参数NRZ为基带信号，16是进行16PSK。

psk_symbols = pskmod(NRZ, 16);

绘制已调信号

figure
plot(t,real(psk_symbols));
title('16PSK调制后的时域图（实部）');
xlabel('时间/s');
ylabel('幅度');

分析已调信号的频谱特性

figure
fft_modulated = fftshift(fft(psk_symbols));
PE1 = 10 * log10(abs(fft_modulated) .^ 2 / (N * Ts)); %功率：能量除以时间，公式法求概率谱密度，单位db
PEL1 = (-length(fft_modulated) / 2 : length(fft_modulated) / 2 - 1) /10; % 求区间长度
plot(PEL1,PE1);
title('16PSK调制后的功率谱');
xlabel('频率/Hz');
ylabel('幅度/dB');
axis([-50 50 -50 50]); 

绘制已调信号的星座图

scatterplot(psk_symbols)				%绘星座图 FIGURE 4
title('16PSK调制后的星座图');

在SNR=10情况下，存在加性高斯白噪声，接收已调信号

snr = 10; % 信噪比
received_signal = awgn(psk_symbols, snr, 'measured');
figure;
plot(t,real(received_signal));
title(['叠加噪声的已调信号 (SNR = ' num2str(snr) ' dB, 实部)']);
xlabel('时间');
ylabel('幅度');

解调信号

demodulated_signal = pskdemod(received_signal, 16);
figure(8);
plot(t,real(demodulated_signal));
title('解调后基带信号的时域图');
xlabel('时间');
ylabel('幅度');
figure(9);
fft_demodulated = fftshift(fft(demodulated_signal));
PE1 = 10 * log10(abs(fft_demodulated) .^ 2 / (N * Ts)); %功率：能量除以时间，公式法求概率谱密度，单位db
PEL1 = (-length(fft_demodulated) / 2 : length(fft_demodulated) / 2 - 1) /10; % 求区间长度
plot(PEL1,PE1);
title('加噪后16PSK解调的功率谱');
xlabel('频率/Hz');
ylabel('幅度/dB');
axis([-50 50 -50 50]); 

绘制误码率

snr_range = -10:20; % 信噪比范围
ber_values = zeros(size(snr_range));
 
for i = 1:length(snr_range)
    received_signal = awgn(psk_symbols, snr_range(i), 'measured');
    demodulated_signal = pskdemod(received_signal, 16);
    num_errors = sum(NRZ ~= demodulated_signal);
    ber_values(i) = num_errors / N;
end
 
% 画出不同信噪比下的系统误码率曲线
figure;
semilogy(snr_range, ber_values, '-o');
title('系统误码率曲线');
xlabel('信噪比 (SNR, dB)');
ylabel('误码率 (BER)');
grid on;