用matlab实现4-PSK通过AWGN信道

以下是matlab的代码实现：

% 4-PSK调制
M = 4;  % 符号数
k = log2(M);  % 每个符号的比特数
Fs = 1000;  % 采样率
fc = 10;  % 载波频率
T = 1/Fs;  % 采样间隔
t = 0:T:1-T;  % 时间序列
fcos = cos(2*pi*fc*t);  % 载波

% 生成符号序列
data = randi([0 M-1], 1, 100);  % 生成100个随机符号
symbols = pskmod(data, M);  % 对符号进行4-PSK调制

% 将符号序列转换为数字信号
bits = de2bi(data, k);  % 将符号转换为比特
tx = reshape(bits', 1, []);  % 将比特序列转换为行向量

% 通过AWGN信道发送数字信号
SNR = 10;  % 信噪比
rx = awgn(tx, SNR, 'measured');  % 添加高斯噪声

% 将数字信号转换为符号序列
rx_bits = reshape(rx, k, [])';  % 将比特序列转换为列向量
rx_data = bi2de(rx_bits);  % 将比特转换为符号

% 对接收到的符号进行解调
rx_symbols = pskdemod(rx_data, M);  % 进行4-PSK解调

% 绘制调制前后的星座图
subplot(2,1,1);
plot(symbols, 'o');
title('调制前');
subplot(2,1,2);
plot(rx_symbols, 'o');
title('调制后');

% 计算误码率
[num, ber] = biterr(data, rx_data);
fprintf('误码率 = %f\n', ber);

这段代码实现了4-PSK调制，发送符号序列，通过AWGN信道添加高斯噪声，接收数字信号并解调成符号序列，最后计算误码率。其中，载波频率为10Hz，每个符号周期内采样100个点，符号周期为1s。