1、实验三 蒙特卡 罗方法论仿真和 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);s
2、ubplot(2,1,1)plot(t,s);for i=1:length(a)if a(i)=1t=i-1:0.01:i;s=sin(2*pi*t);endif a(i)=0t=i-1:0.01:is=-sin(4*pi*t); endsubplot(2,1,2);plot(t,s);hold onend2. 利用 MATLAB 编程，对 蒙特卡罗仿真的差错概率与理论差错概率进行比较。要求大家理解蒙特卡罗仿真的思路，能读懂书本上的源程序，对每行程序进行解释。回答下面问题：N=input(Enter number of symbols ); %设置仿真系统数据符号的个数snrdB_min=-5
3、;snrdB_max=8; %设置信噪比取值的上下限：dB为单位snrdB=snrdB_min:1:snrdB_max; %信噪比取值向量，dB 为单位snr=10.(snrdB/10); %计算相应的信噪比值len_snr=length(snrdB); %信噪比取 值的个数%对每个 给定信噪比的通信系统，计算 BER 的估计值for j=1:len_snrsigma=sqrt(1/(2*snr(j); %由给定的信噪比计算加性白噪声标准差Ne=0; %差错计数器初始化%传输 N 个数据符号，统计差错符号（即比特）for k=1:Nd=round(rand(1); %产 生一个数据符号：可能取
4、值为0，1 round 代表 4 舍 5 入x_d=2*d-1; %得到发送器发送的数据符号：可能取值为-1，+1n_d=sigma*randn(1); %产生信道的加性白噪声y_d=x_d+n_d; %接收器接收的数据%检测器判决接收的数据if y_d0d_est=1; %判决 发送器发送的数据为 1elsed_est=0; %判决 发送器发送的数据为 0end%检测器判决结果与发送器发送的数据符号进行比较if (d_est=d)Ne=Ne+1; %差错计数器累加 1endend%计算差 错符号个数errors(j)=Ne;%计算 误符号率（即误比特率 BER）的估计值ber_sim(j)=
5、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， 仿真结果的数据点数有什么变化？为什么？数据点减小，因为误码概率是确定的，所以 N 越大，误码数目也增加6、所有代码和参数都保持不变，运行多次，每次结 果是否相同？为什么？不同，因为发送器发送的数据符号是随机产生的。
