MATLAB环境下多幅度电平信号仿真程序的说明 .doc

MATLAB环境下多幅度电平信号仿真程序说明 一、多幅度信号仿真系统模型 二、系统中涉及的公式 1、平均比特SNR定义为： 即可得： 二电平时平均比特SNR为： 四电平时平均比特SNR为： 八电平时平均比特SNR为： 2、M电平PAM系统最佳检测器的差错概率为： 其中： 三、仿真结果图 附源程序代码： main_plot.m %主程序源代码，调用其他5个函数得出仿真结果图 clear;clc;echo off;close all; % 二电平单极性 SNR_in_dB1=0:1:17; SNR_in_dB2=0:0.1:16; for i=1:length(SNR_in_dB1) smld_err_prb21(i)=smldPe21(SNR_in_dB1(i)); %实测差错概率 end; for i=1:length(SNR_in_dB2) SNR=exp(SNR_in_dB2(i)*log(10)/10); %信噪比换算成真值 theo_err_prb21(i)=1/2*erfc(sqrt(SNR/4)); %理论差错概率计算 end; semilogy(SNR_in_dB1,smld_err_prb21,o); %实测差错概率结果图 hold on; semilogy(SNR_in_dB2,theo_err_prb21); %理论差错概率结果图 hold on; % 二电平双极性 SNR_in_dB3=0:1:11; SNR_in_dB4=0:0.1:10; for i=1:length(SNR_in_dB3) smld_err_prb22(i)=smldPe22(SNR_in_dB3(i)); %实测差错概率 end; for i=1:length(SNR_in_dB4) SNR=exp(SNR_in_dB4(i)*log(10)/10); %信噪比换算成真值 theo_err_prb22(i)=1/2*erfc(sqrt(SNR)); %理论差错概率计算 end; semilogy(SNR_in_dB3,smld_err_prb22,black.); %实测差错概率结果图 hold on; semilogy(SNR_in_dB4,theo_err_prb22,black); %理论差错概率结果图 hold on; % 四电平双极性 SNR_in_dB5=0:1:16; SNR_in_dB6=0:0.1:14; for i=1:length(SNR_in_dB5) smld_err_prb42(i)=smldPe42(SNR_in_dB5(i)); %实测差错概率 end; for i=1:length(SNR_in_dB6) SNR_per_bit=exp(SNR_in_dB6(i)*log(10)/10); theo_err_prb42(i)=3/4*erfc(sqrt(2/5*SNR_per_bit)); %理论差错概率计算 end; semilogy(SNR_in_dB5,smld_err_prb42,red+); hold on; semilogy(SNR_in_dB6,theo_err_prb42,red); hold on; % 八电平双极性 SNR_in_dB7=0:1:19; SNR_in_dB8=0:0.1:19; for i=1:length(SNR_in_dB7) smld_err_prb82(i)=smldPe82(SNR_in_dB7(i)); %实测差错概率 end; for i=1:length(SNR_in_dB8) SNR_per_bit=exp(SNR_in_dB8(i)*log(10)/10); theo_err_prb82(i)=7/8*erfc(sqrt(9/63*SNR_per_bit)); %理论差错概率计算 end; semilogy(SNR_in_dB7,smld_err_prb82,green*); hold on; semilogy(SNR_in_dB8,theo_err_prb82,green); hold on; xlabel(SNR in dB); ylabel(Pe); title(多幅度电平信号仿真结果与理论值比较); legend(二电平单极性实测值,二电平单极性理论值,二电平双极性实测值,二电平双极性理论值,... 四电平实测值,四电平理论值,八电平实测值,八电平理论值); gngaus