OFDM系统的PAPR曲线的matlab代码

以下是一个简单的 OFDM 系统 PAPR 曲线的 Matlab 代码示例，其中使用了 IFFT 和 FFT 进行信号处理：

clear all;
close all;
clc;

N = 64;                 % 子载波数量
M = 16;                 % 星座大小
L = 10000;              % 仿真帧数
SNR = [0:2:20];         % 信噪比范围
PAPR = zeros(1, length(SNR));

for i = 1:length(SNR)
    for j = 1:L
        bits = randi([0, 1], 1, N*log2(M));      % 生成随机比特序列
        qam = qammod(bits, M);                  % QAM调制
        ifft_sig = ifft(qam)*sqrt(N);           % IFFT
        PAPR(j) = 10*log10(max(abs(ifft_sig).^2)/mean(abs(ifft_sig).^2));  % 计算PAPR
        noise = sqrt(N/(2*10^(SNR(i)/10)))*(randn(1, N) + 1j*randn(1, N)); % 添加噪声
        ofdm_sig(i, j, :) = ifft_sig + noise;   % 添加噪声后的OFDM信号
    end
    PAPR_avg(i) = mean(PAPR);                   % 平均PAPR
end

% 绘制PAPR曲线
figure;
plot(SNR, PAPR_avg, 'b-o', 'LineWidth', 2);
grid on;
xlabel('SNR (dB)');
ylabel('PAPR (dB)');
title('OFDM系统PAPR曲线');