qam调制

M = 16; % 调制阶数
fc = 10; % 载波频率（Hz）
fs = 1000; % 采样频率（Hz）
T = 1 / fs; % 采样周期（秒）
Eb = 1; % 符号能量
N = 1000;    % 信号长度
A = 1;          % 振幅
phi = 0;        % 初始相位
t=(0:N-1)*(1/fs);
num_symbols = randi([0 M-1], N, 1); % 随机符号
% QAM调制
qam_symbols = qammod(num_symbols, M);
% 生成调制信号
carrier = A * cos(2*pi*fc*t + phi);
modulated_signal = real(qam_symbols .* carrier.');


% 添加高斯白噪声
SNR_dB = 20; % 信噪比
modulated_signal_noisy = awgn(modulated_signal, SNR_dB, 'measured');
% 信号解调
demodulated_signal = demodulate_qam(modulated_signal_noisy, A, fc, phi, fs, M);
% 计算欧式距离和能量
distance = sqrt(sum(abs(qam_symbols - demodulated_signal).^2) / N);
energy = sum(abs(modulated_signal).^2) / N;

disp(['欧式距离：', num2str(distance)]);
disp(['平均能量：', num2str(energy)]);


% 绘制星座图
figure;
plot(qam_symbols,'o');
title('QAM星座图');
xlabel('I路');
ylabel('Q路');

% 画出调制信号波形
figure;
plot(t, modulated_signal);
title('QAM调制信号波形');
xlabel('时间 (秒)');
ylabel('幅度');
figure;
plot(t, modulated_signal_noisy);
title('已调信号（含噪声）');
xlabel('时间');
ylabel('幅度');

% QAM信号解调函数
function demodulated_signal = demodulate_qam(signal, A, f_c, phi, fs, M)
    t = (0:length(signal)-1) * (1/fs);
    carrier = A * cos(2*pi*f_c*t + phi);
    demodulated_signal = signal .* carrier.';
    demodulated_signal = downsample(demodulated_signal, fs/(100*f_c));
    demodulated_signal = qamdemod(demodulated_signal, M);
end