AWGN信道中BPSK调制系统BER仿真计算

/序号（学号）：学生实验报告书2014年4月27日实验一：信道中调制系统的仿真计算一、实验目的掌握二相调制的工作原理2．掌握利用进行误比特率测试的方法3．掌握信道中调制系统的仿真计算方法二．实验内容利用仿真程序在环境下完成信道中调制系统的仿真计算，得到仿真结果，写出实验小结，完成实验报告。三．实验仪器：计算机软件四、实验原理在数字领域进行的最多的仿真任务是进行调制解调器的误比特率测试，在相同的条件下进行比较的话，接收器的误比特率性能是一个十分重要的指标。误比特率的测试需要一个发送器、一个接收器和一条信道。首先需要产生一个长的随机比特序列作为发送器的输入，发送器将这些比特调制成某种形式的信号以便传送到仿真信道，我们在传输信道上加上一定的可调制噪声，这些噪声信号会变成接收器的输入，接收器解调信号然后恢复比特序列，最后比较接收到的比特和传送的比特并计算错误。误比特率性能常能描述成二维图像。纵坐标是归一化的信噪比，即每个比特的能量除以噪声的单边功率谱密度，单位为分贝。横坐标为误比特率，没有量纲。五．实验步骤①运行发生器：通过发送器将伪随机序列变成数字化的调制信号。②设定信噪比：假定为m,则0=10,用假设单位为分贝。③确定④计算N0⑤计算噪声的方差σn⑥产生噪声：因为噪声具有零均值，所以其功率和方差相等。我们产生一个和信号长度相同的噪声向量，且该向量方差为σn。⑦加上噪声，运行接收器⑧确定时间延迟⑨产生误差向量⑩统计错误比特：误差向量“”中的每一个非零元素对应着一个错误的比特。最后计算误比特率：每运行一次误比特率仿真，就需要传输和接收固定数量的比特，然后确定接收到的比特中有多少错误的。使用计算:()。实验结果及分析程序：10;1000;200;1;0;1:1((1));%1,01.*2;%0>1,1>1(1).*(1);[];1:210.^(10);()>()N0;(N02);().*.*n;y1((y));y2=(11)2;%1,0((y2));;[];;;[];1:210.^(10);1((2*));[1];1:2;(,'','*b');('0()')('')('','')程序分析：做1000次试验，每次试验取200个抽样点，求出每次试验的误比特率，然后对1000次试验的误比特率取平均值，即得仿真误比特率，然后将此误比特率及理论值进行比较。程序运行结果如图3所示。图1误比特率仿真结果及理论值比较结果分析：由图3可知，仿真结果及理论值基本相符，只是在信噪比较大时，仿真误码率及理论值存在误差，这是因为在高信噪比情况下，误比特率会呈下降趋势，若要仿真实际的误比特率，就需要进行更多次试验，而且每次试验取的抽样点数也应该增加。1000;100;10;1;0;0;1;s10((1));(s10*2-1)(2);S1(1:2);S2(2:2);1.*S2;(12).*(12);[];[];0:1;(1/2)*(10.^(10));.*;((Y))(2);((Y))(2);[];1();[,[(k)(k)]];;*;;;0;1(k)0;1;;;0;12;(k)0;1;;;[];[(2)];;;%;;[];[];0:1;10.^(10);((2.*));1-(1-(1/2).*(())).^2;[];[];;0:1;(,'','*b');('','');('0()');('');0:1;(,'','*y');('','');('0()');('');程序分析：及类似，做1000次试验，每次试验取100个点，然后每两个点形成一个符号，首先计算仿真误比特率和误码率，然后及理论值进行比较。 程序运行结果如图4和图5所示，其中图4是图，图5是图。图2仿真结果及理论值结果分析：由图4和图5可知，在低信噪比的情况下，和的仿真结果及理论值完全相符，在高信噪比情况下，两者存在差异，这是因为随着信噪比的增加，误比特率和误码率也会增加，这就需要更多的试验次数和更多的抽样点数。将和仿真结果进行比较可知，在相同信噪比的情况下，比小，这是因为只要一个比特错，相应的符号就会出错，而一个要正确接收，就必须保证它的两个比特全部正确接收。图3仿真结果及理论值实验二：无线通信信道建模一、实验目的1.掌握无线移动通信信道的特点，根据相关模型进行仿真。2.掌握语言对上述参数进行仿真。二．实验内容利用仿真程序在环境下完成移动信道建模的仿真分析，得到仿真结果，写出实验小结，完成实验报告。三．实验仪器计算机软件四、实验原理当移动台在一个较小的范围（小于20个工作波长）运动时，引起接收信号的幅度、相位和到达角等的快速变化，这种变化称为小尺度衰落。典型的小尺度衰落有、衰落，因为当信号在传播过程中经过许多反射路径后，接收到的信号幅度可以用或概率密度函数来描述。在接受信号有直达信号的情况下，幅度的衰落呈现分布，而当在接收端没有直达信号的情况下，幅度的衰落呈现分布。采用小尺度衰落模型的信道，衰落幅度是服从或分布的随机变量，这些变量将会影响到接收信号的幅度和功率。五．实验步骤1．选择路径数2．按均匀分布产生各条路径的延迟3．按功率时延谱确定对应的各径的功率4．按模型产生各径的瑞利衰落系数5．对瑞利衰落系数进行统计分析并及理论值相比较说明：1．路径数目2-4自己确定，或采用某个国际标准2．每条路径时间延迟满足（0，）范围内均匀分布，为自己选择的最大采样步长数200-600间比较合适，或采用国际标准3．功率可以按时延迟谱求得，也可用国际标准测量值。功率延迟谱：①若采用等功率分配产生功率：Pit；②采用指数分布的功率延迟谱产生功率：1/6*(6)六．实验结果及分析程序：%;=30;M=8;%#-1N=4*2;1.02404;=2(N);=1(2);=0;%=0;t0=0.001;t=0:0.5%=+1g()=0;n=1:1,n<=M()=2**(*);%()=2**(*);%()=*(2**);%()=*(2**);%()=(2)*(4);()=(2)*(4);()=;()=;;%()=()*(2**()*(0)+);%()=()*(2**()*(0)+);%;n()=*((,11))1()=*((,11));%n(.^21.^2);((.^2));[]=(,'');%[]=(1,'');%();[22]=(1,'')(1)/2:1:(1)/2(2001)2...1:120011000:1:1000;*;2.*.**;0=1;0.^2.*(0);;((),'-','*')('t()');('');('');((),'-','*')('t()');('');('')1=1:1000;(1*(1:1000));('t()');('.');();(1,501);0;(100.*);1(k)<=5001;(1(k)+1)(1(k)+1)+1;0.01;2=0.5;[0:0.01:5];(2).*(-1.*x.^2(2.*2));(,'-','*');('','');('r');('r');程序分析：首先选定模型得到参数8、4*2，然后计算各路的同相及正交分量的增益和多普勒频率，以及它们的表达式，该程序假设各路的延时相等，均为t0=0.001，各路相位均为0。然后分别求同相分量和正交分量的自相关和互相关，分别作图，并分别及理论值相比较。接着作出瑞利衰落系数图形，最后作出仿真的和理论的瑞利分布的概率密度函数，并进行比较。程序运行结果如图6-图9所示。图4.同相分量的自相关仿真图形和理论图形图5.正交分量的自相关仿真图形和理论图形图6.瑞利衰落系数图形图7.瑞利分布概率密度函数仿真图形和理论图形结果分析：图4和图5的自相关和互相关函数图形仿真结果及理论结果大致相同，说明模型很好地分离出了接收信号同相分量和正交分量，图6是瑞利衰落系数的时域图形，图7中瑞利分布概率密度函数的仿真结果的大致变化趋势及理论值相同，只不过在概率密度较大的区域两者有较大差异，这是因为抽样点数太少，无法满足理论计算的要求，而且模型的M、N取值也会影响概率密度函数图形的仿真结果。实验三：通信仿真实验目的1.通信具有很多通信特点，不仅被95移动通信系统使用，目前已成为3G的主要技术。2.通过实验：①掌握直接序列扩频发射机及接收机的组成及仿真；②仿真验证信道下单用户直接序列扩频系统的性能；③仿真验证平坦瑞利信道下单用户直接序列扩频系统的性能；④观察存在干扰用户时的系统性能变化。实验内容设计一个系统，用进行仿真，统计或随信噪比的关系，绘出曲线。对统计试验的结果及单用户的理论值进行比较，对仿真结果进行分析。三．实验仪器计算机软件四．实验原理仿真基带直接序列扩频系统：采用或映射。扩频序列可以是随机产生，可以是m序列，也可以是码，长度自选。最后对或随信噪比变化画图及理论单用户的结果比较，并对仿真结果进行分析。五．实验步骤确定用户数目、信道特征以及调制方式。确定基带扩频仿真系统的原理结构图，按照框图设计一个系统，并进行仿真。信源卷积编码信源卷积编码交织编码加扰码扩频调制信道调制解扩信宿解卷积解交织去扰码用进行仿真，统计或随信噪比的关系，绘出曲线。对统计试验的结果及单用户的理论值进行比较。对仿真结果进行分析。六．实验结果及分析程序：95%此函数用于95前向链路系统的仿真，包括扩%频调制，匹配滤波，接收等相关通信模块。%仿真环境:加性高斯白噪声信道.%数据速率=9600%('');Rj;=0;%控制程序运行中的显示=0;%选择软/硬判决接收主要的仿真参数设置=9600;%比特率=1228800;%码片速率N=184;%源数据数=1;%匹配滤波器类型升余弦R=5;生成多项式=[111101011;101110001]生成多项式矩阵K=(,2);%列数L=(,1);%行数矩阵=64;码的长度=([1;0]*(1,2),,1);%32个10行=(,1);扩频调制码的生成多项式=[1010001110100001]';=[1001110001111001]';=[15,13,9,8,7,5,0]';路码生成多项式参数=[15,12,11,10,6,5,4,3,0]';路码生成多项式参数=(16,1);%16×1的0矩阵(16)=(())根据配置i路码生成多项式=[(()-1,1);1];%I路信道码生成器的初始状态=(16,1);%16×1的0矩阵(16)=(());%根据配置q路码生成多项式=[(()-1,1);1];%Q路信道码生成器的初始状态扰码生成多项式=[42,35,33,31,27,26,25,22,21,19,18,17,16,10,7,6,5,3,2,1,0]';=(43,1);%43×1的0矩阵(43)=(());%根据配置扰码生成多项式=[(()-1,1);1];%长序列生成器的初始状态信道=10*10()处理增益=10*10(L);=[-1:0.5:1];%仿真信噪比范围()=[-2:0.5:-1.5];主程序=[];=[];=[];(1)%判断软/硬判决接收('\n\n\n');('\n\n\n');1()%根据多次运行('\%1.1f()',(i))输出当前值=0;=0;=0;(<300)&(<150);发射机=((N,1)>0)生成源数据%速率为19.2[,]=(,,);%产生95前向链路系统的发送数据包%速率为1.2288[x]=(,,)实现95前向链路系统的数据调制信道=1(2)*(2)*(((x))+j*((x)))*10^(-((i)-)/20)生成噪声序列r=加入噪声接收机=(r,,,);%软判决，速率为19.2=(>0);%定义接收码片的硬判决()[]=(,,);%软判决[]=(,,);%硬判决()(311);(,'');('');%软判决结果图(312);((,),'');('')接收机输入符号及