移动通信实验报告精讲.doc

序号（学号）： 178840056 沈阳农业大学学生实验报告书 实验类别 移动通信实验 学 院 信息工程学院 专 业 通信工程 班 级 13级电信二班 姓 名 恽星彤 指导教师 张大鹏老师 2015 年 12 月 10 日 实验一：AWGN信道中BPSK调制系统的?BER仿真计算 一、实验目的 掌握二相BPSK调制的工作原理。 掌握利用MATLAB进行误比特率测试BER的方法。 掌握AWGN信道中BPSK调制系统的BER仿真计算方法。 二、实验原理 1. 仿真概述及原理 在数字领域进行的最多的仿真任务是进行调制解调器的误比特率测试，在相同的条件下进行比较的话，接收器的误比特率性能是一个十分重要的指标。误比特率的测试需要一个发送器、一个接收器和一条信道。首先需要产生一个长的随机比特序列作为发送器的输入，发送器将这些比特调制成某种形式的信号以便传送到仿真信道，我们在传输信道上加上一定的可调制噪声，这些噪声信号会变成接收器的输入，接收器解调信号然后恢复比特序列，最后比较接收到的比特和传送的比特并计算错误。 误比特率性能常能描述成二维图像。纵坐标是归一化的信噪比，即每个比特的能量除以噪声的单边功率谱密度，单位为分贝。横坐标为误比特率，没有量纲。 2. 仿真过程及计算 运行发生器：通过发送器将伪随机序列变成数字化的调制信号。 设定信噪比：假定SNR为m dB,则Eb/N0=10^(m/10),用MATLAB假设SNR单位为分贝。 确定Eb。 计算N0。 计算噪声的方差σn。 产生噪声：因为噪声具有零均值，所以其功率和方差相等。我们产生一个和信号长度相同的噪声向量，且该向量方差为σn。 加上噪声，运行接收器。 确定时间延迟。 产生误差向量。 统计错误比特：误差向量“err”中的每一个非零元素对应着一个错误的比特。 最后计算误比特率BER：每运行一次误比特率仿真，就需要传输和接收固定数量的比特，然后确定接收到的比特中有多少错误的。使用MATLAB计算BER: ber=te/length(tx)。 三、实验结果 MATLAB仿真程序如下： % Simulation of BPSK AWGN Max_SNR=10; N_trials=1000; N=200; Eb=1; ber_m=0; for trial=1:1:N_trials; trial msg=round(rand(1,N)); % 1,0 sequence s=1-msg.*2; %0--1,1--1 n=randn(1,N)+j.*randn(1,N); %generate Guassian white noise ber_v=[]; for snr_dB=1:2:Max_SNR snr=10.^(snr_dB./10); %snr(db)--snr(decimal) N0=Eb./snr; sgma=sqrt(N0./2); y=sqrt(Eb).*s+sgma.*n; y1=sign(real(y)); y2=(1-y1)./2; %1, 0 sequence error=sum(abs(msg-y2)); %error bits ber_snr=error./N; %ber ber_v=[ber_v,ber_snr]; end %for snr ber_m=ber_m+ber_v; end ber=ber_m./N_trials; ber_theory=[]; for snr_db=1:2:Max_SNR snr=10.^(snr_db./10); snr_1=qfunc(sqrt(2*snr)); ber_theory=[ber_theory,snr_1]; end i=1:2:Max_SNR; semilogy(i,ber,-r,i,ber_theory,*b); xlabel(E_b/N_0(dB)) ylabel(BER) legend(Monte Carlo,Theoretic) 仿真结果如下： 四、实验小结 用蒙特卡罗分析方法得到的误码率曲线与理论公式 得到的误码率曲线吻合，可见蒙特卡罗分析是一种很好的分析