1．试编写程序，仿真4PAM调制信号在高斯信道下的性能，并与理论分析结果相比。(1）画出两条性能曲线，一条是根据理论平均错误概率画出，另一条是仿真曲线；(2）程序的基本流程：信源产生信息比特、调制、将调制信号送入信道（产生高斯白噪声的程序）、接收端检测、将检测结果与信源原始信息比较计算误符号率和误比特率；(3）在给定信噪比下，第二步需多次重复，以得到一个平均错误概率；(4）信噪比范围：4PAM(0dB-14dB)，间隔是1dB;也可在 BER =106左右终止。(5)信噪比计算 SNR =10log( Es /N0)=10log( REb /N0)。注意调制不能运用MATLAB内置函数pammod。注意标注中文注释。注意检测方法使用多进制调制信号软输出检测。注意一定用MATLAB编写。编写程序时注意矩阵维度要一致不要出错。注意索引值不能超出数组边界，不要出错。检测多进制调制信号软输出检测，注意不要出现矩阵维度必须一致的错误

以下是一个可能的程序实现，其中包括了编码、调制、通信信道、接收端检测和误码率计算等步骤。程序中的注释是中文的，尽可能详细地解释了各个步骤的实现方法和意义。

%% 4PAM调制信号在高斯信道下的性能仿真
clear all;
close all;
clc;

%% 参数定义
M = 4;  % 调制阶数
k = log2(M);  % 每个符号所需的比特数
Eb = 1;  % 平均比特能量
Es = Eb * k;  % 平均符号能量
N = 10^6;  % 发送的符号数
SNR_dB = 0:1:14;  % 信噪比范围
noise_var = Es ./ (2 * 10.^(SNR_dB / 10));  % 计算噪声方差
err_symbol_theory = zeros(length(SNR_dB), 1);  % 理论平均错误概率
err_symbol_sim = zeros(length(SNR_dB), 1);  % 仿真平均错误概率
err_bit_theory = zeros(length(SNR_dB), 1);  % 理论平均误比特率
err_bit_sim = zeros(length(SNR_dB), 1);  % 仿真平均误比特率

%% 信源产生信息比特
source_bits = randi([0, 1], N, k);  % 产生随机比特流

%% 编码
% 4PAM调制使用格雷码编码
gray_code = [0 1 3 2];  % 4PAM调制的格雷码表
source_gray = bi2de(source_bits, 'left-msb') + 1;  % 将二进制比特流转化为十进制符号
mod_signal = gray_code(source_gray)';  % 格雷码映射后得到的调制信号

%% 通信信道
for i = 1:length(SNR_dB)
    % 产生高斯白噪声
    noise = sqrt(noise_var(i)) * randn(1, N);
    % 加入高斯白噪声
    channel_out = mod_signal + noise;
    
    %% 接收端检测
    % 多进制调制信号软输出检测，此处采用基于后验概率的判决方法
    % 这里为了简化，假设接收端已知噪声方差
    P = zeros(M, N);  % 初始化后验概率矩阵
    for j = 1:M
        % 计算每个符号的后验概率
        P(j, :) = exp(-abs(channel_out - (2 * j - M - 1) * sqrt(Es) / sqrt(10^(SNR_dB(i) / 10))) .^ 2 / (2 * noise_var(i))) / sqrt(2 * pi * noise_var(i));
    end
    [~, detected] = max(P);  % 选择后验概率最大的符号作为检测结果
    detected = detected' - 1;  % 恢复到十进制符号
    
    %% 计算误码率
    % 计算平均错误概率和平均误比特率
    err_symbol_theory(i) = 3 / 2 * erfc(sqrt(0.5 * 10^(SNR_dB(i) / 10)));  % 理论平均错误概率
    err_symbol_sim(i) = sum(detected ~= source_gray) / N;  % 仿真平均错误概率
    err_bit_theory(i) = err_symbol_theory(i) * k;  % 理论平均误比特率
    err_bit_sim(i) = sum(sum(dec2bin(abs(detected - source_gray)')) ~= '0') / (N * k);  % 仿真平均误比特率
end

%% 画图
figure;
semilogy(SNR_dB, err_symbol_theory, 'bo-', 'LineWidth', 1.5, 'MarkerSize', 8);
hold on;
semilogy(SNR_dB, err_symbol_sim, 'r+-', 'LineWidth', 1.5, 'MarkerSize', 8);
grid on;
legend('理论曲线', '仿真曲线');
xlabel('信噪比（dB）');
ylabel('平均错误概率');
title('4PAM调制在高斯信道下的性能');

figure;
semilogy(SNR_dB, err_bit_theory, 'bo-', 'LineWidth', 1.5, 'MarkerSize', 8);
hold on;
semilogy(SNR_dB, err_bit_sim, 'r+-', 'LineWidth', 1.5, 'MarkerSize', 8);
grid on;
legend('理论曲线', '仿真曲线');
xlabel('信噪比（dB）');
ylabel('平均误比特率');
title('4PAM调制在高斯信道下的性能');

注意：由于这个程序使用了随机数生成器，因此每次运行的结果可能会略有不同。此外，由于噪声的影响，仿真曲线会有一定的抖动。