资源描述:
《matlab作业》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、多幅度电平信号差错率仿真设计一、设计任务完成多幅度电平信号差错率的仿真,运用MATLAB语言,仿真它们的误码性能,并与理论误码曲线相比较。先产生多组符号序列,再将符号序列映射到对应的幅度电平上,以四电平为例,先产生一个四组符号的序列,将该序列映射到对应的幅度电平上。随机变量信道中的加性噪声可以认为是统计独立,均值为0的高斯分量,分别作出二、四、八、十六电平仿真得到的误码率曲线与和理论上的误码率曲线。二、物理模型多幅度电平信号仿真系统模型如下图所示:图1三、预期目标通过以上的物理模型,作出二、四、八、十六电平信号仿真得到的误码率曲线与和理论上的误码率曲线。四
2、、系统中所涉及的公式如下:1、平均比特SNR定义为:即可得:二电平时平均比特SNR为:四电平时平均比特SNR为:八电平时平均比特SNR为:十六电平时平均比特SNR为:2、M电平PAM系统最佳检测器的差错概率为:其中:五、仿真结果图2六、源程序代码main_plot.m%主程序源代码,调用其他5个函数得出仿真结果图clear;clc;echooff;closeall;%二电平单极性SNR_in_dB1=0:1:17;SNR_in_dB2=0:0.1:16;fori=1:length(SNR_in_dB1)smld_err_prb21(i)=smldPe21(
3、SNR_in_dB1(i));%实测差错概率end;fori=1:length(SNR_in_dB2)SNR=exp(SNR_in_dB2(i)*log(10)/10);%信噪比换算成真值theo_err_prb21(i)=1/2*erfc(sqrt(SNR/4));%理论差错概率计算endsemilogy(SNR_in_dB1,smld_err_prb21,'blueO');%实测差错概率结果图holdon;semilogy(SNR_in_dB2,theo_err_prb21,'blue');%理论差错概率结果图holdon;%二电平双极性SNR_in_
4、dB3=0:1:11;SNR_in_dB4=0:0.1:10;fori=1:length(SNR_in_dB3)smld_err_prb22(i)=smldPe22(SNR_in_dB3(i));%实测差错概率end;fori=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,'blackX');%实测差错概率结
5、果图holdon;semilogy(SNR_in_dB4,theo_err_prb22,'black');%理论差错概率结果图holdon;%四电平双极性SNR_in_dB5=0:1:16;SNR_in_dB6=0:0.1:14;fori=1:length(SNR_in_dB5)smld_err_prb42(i)=smldPe42(SNR_in_dB5(i));%实测差错概率end;fori=1:length(SNR_in_dB6)SNR_per_bit=exp(SNR_in_dB6(i)*log(10)/10);theo_err_prb42(i)=3/4
6、*erfc(sqrt(2/5*SNR_per_bit));%理论差错概率计算end;semilogy(SNR_in_dB5,smld_err_prb42,'red.');holdon;semilogy(SNR_in_dB6,theo_err_prb42,'red');holdon;%八电平双极性SNR_in_dB7=0:1:19;SNR_in_dB8=0:0.1:19;fori=1:length(SNR_in_dB7)smld_err_prb82(i)=smldPe82(SNR_in_dB7(i));%实测差错概率end;fori=1:length(SNR
7、_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*');holdon;semilogy(SNR_in_dB8,theo_err_prb82,'green');holdon;%十六电平双极性SNR_in_dB9=0:1:19;SNR_in_dB10=0:0.1:19;fori=1:length(SNR
8、_in_dB9)smld_err_prb162(i)=smldPe