pam_main.m

simulink框图来实现多进制PAM调制

% 2PAM Monte Carlo Simulation
clear all
clc
tic
% parameters:
BitRate=10000;

x=0:2.5:25;
for i=1:length(x)
    M=2;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y1(i)=mean(BitErrorRate);
end   

for i=1:length(x)
    M=4;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y2(i)=mean(BitErrorRate);
end

for i=1:length(x)
    M=8;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y3(i)=mean(BitErrorRate);
end

for i=1:length(x)
    M=16;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y4(i)=mean(BitErrorRate);    
end

for i=1:length(x)
    M=32;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y5(i)=mean(BitErrorRate);    
end

for i=1:length(x)
    M=64;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y6(i)=mean(BitErrorRate);
end

for i=1:length(x)
    M=128;
    d=sqrt(3*log2(M)/(M^2-1));
    SNR=x(i);
    sim('PAM_sim');
    y7(i)=mean(BitErrorRate);  
end
%plot
figure
semilogy(x,y1,'-b.');hold on
semilogy(x,y2,'--r.');hold on
semilogy(x,y3,'-.k.');hold on
semilogy(x,y4,':m.');hold on
semilogy(x,y5,'--k.');hold on
semilogy(x,y6,'-.y.');hold on
semilogy(x,y7,':k.');
axis([0 25 10^(-3) 1]);grid
xlabel('SNR per bit expressed in dB');
ylabel('Bit Error Probability');
legend('2-PAM','4-PAM','8-PAM','16-PAM','32-PAM','64-PAM','128-PAM');
hold off
toc