M-QAM调制在衰落信道下的误比特率及峰均比性能分析，信道可以采用单抽头衰落信道（即无频率选择性衰落）+AWGN信道，分析不同调制阶数的QAM在不同信噪比下的误比特率。具体matlab代码

以下是使用 MATLAB 实现 M-QAM 调制在衰落信道下的误比特率及峰均比性能分析的代码。其中，假设调制阶数 M = 16，信道采用瑞利衰落信道，信噪比范围为 -10dB 到 20dB：

%% 参数设置
M = 16; % 调制阶数
k = log2(M); % 每个符号的比特数
n = 10000; % 发送的符号数
EbN0dB = -10:2:20; % 信噪比范围

%% 生成M-QAM调制的星座图
x = [-(sqrt(M)-1):2:(sqrt(M)-1)];
y = [-(sqrt(M)-1):2:(sqrt(M)-1)];
[X,Y] = meshgrid(x,y);
map = (X + 1j*Y)/sqrt(10);

%% 发送端
% 生成随机发送符号
tx = randi([0 M-1],1,n);
% 星座映射
st = map(tx+1);

%% 信道
% 随机生成瑞利衰落信道
h = (randn(1,n) + 1j*randn(1,n))/sqrt(2);
% 瑞利衰落信道
rx = h.*st;
% 添加AWGN噪声
for i = 1:length(EbN0dB)
    % 计算噪声方差
    Es = mean(abs(st).^2);
    Eb = Es/k;
    N0 = Eb/(10^(EbN0dB(i)/10));
    noise = sqrt(N0/2)*(randn(1,n) + 1j*randn(1,n));
    % 加噪声
    rxn = rx + noise;
    % 信道估计
    h_hat = rxn./st;
    % 接收信号均衡
    rxe = rxn./h_hat;
    % 信号检测
    [dist,idx] = min(abs(rxe(:)-map(:)).^2,[],1);
    rx_sym = idx-1;
    % 统计误比特率
    err(i) = sum(tx~=rx_sym)/n;
    % 统计峰均比
    peak(i) = max(abs(rxn).^2)/mean(abs(rxn).^2);
end

%% 绘制误比特率和峰均比性能曲线
figure;
semilogy(EbN0dB,err,'bo-','LineWidth',2);
hold on;
semilogy(EbN0dB,peak,'r*-','LineWidth',2);
grid on;
legend('误比特率','峰均比');
xlabel('信噪比(dB)');
ylabel('性能');
title(['M-QAM 调制在瑞利衰落信道下的性能分析(M = ' num2str(M) ')']);

以上代码主要分为三个部分：

1. 参数设置：设置调制阶数、每个符号的比特数、发送的符号数和信噪比范围。
2. 发送端：随机生成发送符号，并进行星座映射得到发送的调制符号。
3. 信道：首先随机生成瑞利衰落信道，然后对于每个信噪比值，计算噪声方差并添加AWGN噪声，接着进行信道估计和接收信号均衡，最后进行信号检测并统计误比特率和峰均比性能。

最后，绘制误比特率和峰均比性能曲线，分析不同信噪比下的调制性能。