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

实验三蒙特卡罗方法论仿真和2PSK调制(综合性实验)(考试)要求编写2PSK调制程序，任意给定一组二进制数，计算经过这种调制方式的输出信号。程序书写要规范，加必要的注释；经过程序运行的调制波形要与理论计算出的波形一致。熟悉2PSK调制原理。编写2PSK调制程序。画出原信号和调制信号的波形图。实验报告要求如下内容：2PSK调制原理；对给定信号画出理论调制波形；程序设计思想，画出流程图；源程序代码测试结果(打印)和理论计算结果对比是否一致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);fori=1:length(a)ifa(i)==1t=i-1:0.01:i;s=sin(2*pi*t);endifa(i)==0t=i-1:0.01:is=-sin(4*pi*t);endsubplot(2,1,2);plot(t,s);

holdonend2.利用MATLAB2.利用MATLAB编程，对蒙特卡罗仿真的差错概率与理论差错概率进行比较。要求大家理解蒙特卡罗仿真的思路，能读懂书本上的源程序，对每行程序进行解设置仿真系统数据符号的设置仿真系统数据符号的%设置信噪比取值的上下限：dB%信噪比取值向量，dB为单位%计算相应的信噪比值%信噪比取值的个数释。回答下面问题：N=input('Enternumberofsymbols>');个数snrdB_min=-5;snrdB_max=8;为单位snrdB=snrdB_min:1:snrdB_max;snr=10.A(snrdB/10);len_snr=length(snrdB);%对每个给定信噪比的通信系统，计算BER的估计值forj=1:len_snrsigma=sqrt(1/(2*snr(j))); %由给定的信噪比计算加性白噪声标准差Ne=0; %差错计数器初始化%传输N个数据符号，统计差错符号(即比特)fork=1:N%产生一个数据符号：可能取值d=round(rand(1));

%产生一个数据符号：可能取值%得到发送器发送的数据符号：可能取值为%产生信道的加性白噪声%接收器接收的数据%判决发送器发送的数据为1%判决发送器发送的数据为0为0,1round%得到发送器发送的数据符号：可能取值为%产生信道的加性白噪声%接收器接收的数据%判决发送器发送的数据为1%判决发送器发送的数据为0-1,+1n_d=sigma*randn(1);y_d=x_d+n_d;%检测器判决接收的数据ify_d>0d_est=1;elsed_est=0;end%检测器判决结果与发送器发送的数据符号进行比较if(d_est~=d)Ne=Ne+1;%差错计数器累加1Ne=Ne+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_minsnrdB_max0.00011])xlabel('SNRindB')ylabel('BER')legend('Theoretical','Simulation')1、 程序中哪行代码相当于BPSK?12行2、 程序中哪行代码表示信号经过了噪声信道？噪声强度由哪个参数来控制?14至15彳亍；snrdB_min,snrdB_max;3、 改变参数虬观察结果随N如何变化？并分析原因。随着N的增大，误码率随之减小4、 为什么曲线是单调降的？信噪比越大，受干扰能力越强，误码率越小5、 当N太少时，比如50,仿真结果的数据点数有什么变化？为什么？