实验三 蒙特卡罗方法论仿真和2PSK调制.doc

实验三 蒙特卡罗方法论仿真和2PSK调制（综合性实验）1. （考试）要求编写2PSK调制程序，任意给定一组二进制数，计算经过这种调制方式的输出信号。程序书写要规范，加必要的注释；经过程序运行的调制波形要与理论计算出的波形一致。1） 熟悉2PSK调制原理。2） 编写2PSK调制程序。3） 画出原信号和调制信号的波形图。实验报告要求如下内容：1）2PSK调制原理；对给定信号画出理论调制波形；2) 程序设计思想，画出流程图;3) 源程序代码4) 测试结果（打印）和理论计算结果对比是否一致a=1,1,0,0,1,0,1; t=0:0.01:10; s=sin(2*pi*t); subplot(2,1,1) plot(t,s); for i=1:length(a) if a(i)=1 t=i-1:0.01:i; s=sin(2*pi*t);end if a(i)=0 t=i-1:0.01:i s=-sin(4*pi*t); end subplot(2,1,2); plot(t,s); hold onend2. 利用MATLAB编程，对蒙特卡罗仿真的差错概率与理论差错概率进行比较。要求大家理解蒙特卡罗仿真的思路，能读懂书本上的源程序，对每行程序进行解释。回答下面问题：N=input(Enter number of symbols ); %设置仿真系统数据符号的个数snrdB_min=-5;snrdB_max=8; %设置信噪比取值的上下限：dB为单位snrdB=snrdB_min:1:snrdB_max; %信噪比取值向量，dB为单位snr=10.(snrdB/10); %计算相应的信噪比值len_snr=length(snrdB); %信噪比取值的个数%对每个给定信噪比的通信系统，计算BER的估计值for j=1:len_snr sigma=sqrt(1/(2*snr(j); %由给定的信噪比计算加性白噪声标准差Ne=0; %差错计数器初始化%传输N个数据符号，统计差错符号（即比特） for k=1:N d=round(rand(1); %产生一个数据符号：可能取值为0，1 round 代表4舍5入 x_d=2*d-1; %得到发送器发送的数据符号：可能取值为-1，+1 n_d=sigma*randn(1); %产生信道的加性白噪声 y_d=x_d+n_d; %接收器接收的数据 %检测器判决接收的数据 if y_d0 d_est=1; %判决发送器发送的数据为1 else d_est=0; %判决发送器发送的数据为0 end %检测器判决结果与发送器发送的数据符号进行比较 if (d_est=d) Ne=Ne+1; %差错计数器累加1 endend%计算差错符号个数errors(j)=Ne;%计算误符号率（即误比特率BER）的估计值 ber_sim(j)=errors(j)/N;end%对每个给定信噪比的通信系统，计算BER的理论值ber_theor=qfunc(sqrt(2*snr); %可以看书上对2psk的误码率分析%BER曲线：理论值和估计值对比图semilogy(snrdB,ber_theor,snrdB,ber_sim,o)axis(snrdB_min snrdB_max 0.0001 1)xlabel(SNR in dB)ylabel(BER)legend(Theoretical,Simulation)1、程序中哪行代码相当于BPSK? 12行2、程序中哪行代码表示信号经过了噪声信道？噪声强度由哪个参数来控制？14至15行；snrdB_min,snrdB_max;3、改变参数N,观察结果随N如何变化？并分析原因。随着N的增大，误码率随之减小4、为什么曲线是单调降的？信噪比越大，受干扰能力越强，误码率越小5、当N太少时，比如50， 仿真结果的数据点数有什么变化？为什么？