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

以下是一个可能的实现，仅供参考。在实际应用中，可能需要进行更多的测试和优化。

clc;
clear all;
close all;
%% 参数设置
M = 4; % 4PAM
Eb = 1; % 符号能量
T = 1; % 符号周期
nBits = 1000000; % 发送比特数
nSnr = 15; % 信噪比范围
nErrs = zeros(1,nSnr); % 错误比特数
nRuns = 20; % 每个信噪比运行次数
nSymPerFrame = 100; % 每帧符号数

%% 信源产生信息比特
txBits = randi([0 M-1],1,nBits);

%% 调制
txSym = pammod(txBits,M);

%% 信道
for iSnr = 1:nSnr
    EbN0dB(iSnr) = iSnr - 1; % 信噪比，单位dB
    EbN0(iSnr) = 10^(EbN0dB(iSnr)/10); % 转换为线性值
    n0(iSnr) = Eb/(2*EbN0(iSnr)); % 噪声功率
    for iRun = 1:nRuns
        %% 加入高斯白噪声
        noise = sqrt(n0(iSnr)/2)*(randn(1,nSymPerFrame)+1j*randn(1,nSymPerFrame));
        rxSym = txSym + noise;
        %% 接收端检测，软输出检测
        rxSymDemod = real(rxSym) > 0;
        %% 计算误符号率和误比特率
        nErrs(iSnr) = nErrs(iSnr) + sum(rxDemod ~= txBits(1:nSymPerFrame))/nSymPerFrame;
    end
end

%% 计算平均错误概率
pErrSim = nErrs/(nBits*nRuns);

%% 理论分析
pErrTheory = 2*(1-1/sqrt(M))*qfunc(sqrt(3*EbN0/(M-1)));

%% 画图
figure;
semilogy(EbN0dB,pErrSim,'o-');
hold on;
semilogy(EbN0dB,pErrTheory,'r-');
grid on;
xlabel('SNR(dB)');
ylabel('P_e');
legend('仿真','理论分析');