课程设计OQPSK通信系统的matlab仿真分析Word格式文档下载.doc

1、通信原理 （国防工业出版社）樊昌信 曹丽娜 编著信号与系统MATLAB综合实验（高等教育出版社）MATLAB7辅助信号处理技术与应用（电子工业出版社）飞思科技产品研发中心编著窗体顶端窗体底端 OQPSK通信系统的matlab仿真分析 设计目的和意义1. 对oqpsk进行调制和解调，通过MATLAB编程，掌握MATLAB的使用，熟练掌握OQPSK的调制原理，解调原理。2. 对OQPSK通信系统进行matlab仿真分析，分析起信噪比和差错率。为现实中通信系统的调制，解调，及信道传输进行理论指导。 设计原理1. OPSK的调制它和QPSK有着同样的相位关系，也是把输入码流分成两路，然后进行正交调制。

2、随着数字通信技术的发展和广泛应用，人们对系统的带宽、频谱利用率和抗干扰性能要求越来高。而与普通的QPSK比较，交错正交相移键控的同相与正交两支路的数据流在时问上相互错开了半个码元周期，而不像QPSK那样I、Q两个数据流在时间上是一致的（即码元的沿是对齐的）。由于OQPSK信号中的I（同相）和Q（正交）两个数据流，每次只有其中一个可能发生极性转换，所以，每当一个新的输入比特进入调制器的I或Q信道时，其输出的OQPSK信号中只有0、+90三个相位跳变值，而根本不可能出现180相位跳变。所以频带受限的OQPSK信号包络起伏比频带受限的QPSK信号要小，而经限幅放大后的频带展宽也少。2OQPSK的基本

3、原理OQPSK信号的数学公式可以表示为：OQPSK的调制方法与QPSK类似，仅在一条正交支路上引入了一个比特的延时，以使得两支路的数据不会同时发生变化，降低最大相位跳变。其中电平映射关系为:11,0-1.图3：OQPSK调制原理框图经OQPSK调制后，调制点的星相图以及状态转移图如图4所示。OQPSK解调OQPSK信号可采用正交相干解调方式解调，其解调原理如图a所示。由图a可以看出，OQPSK与QPSK信号的解调原理基本相同，其差别仅在于对Q支路信号抽样判决后要延迟Tb2，这是因为在调制时，Q支路信号在时间上偏移了Tb2，所以抽样判决时刻也相应偏移了Tb2，以保证对两支路的交错抽样。眼图是信号

4、由垂直扫描进入与同周期的水平扫描锯齿波叠加到示波器上时到得图案。眼图能够反映信号在传输过程中受到的信道噪声影响的强度，眼图越模糊，眼睛越闭合，则说明噪声越强，反之，则说明噪声强度弱，也能说明信道性能更优良。 设计结果及分析输入的二进制序列以及串并转换成的上下两路信号，下支路信号已延时，由产生的随机序列可以知道前十个码元为1001001101，映射电平为1-1-11-1-111-11，下面的抽样脉冲与ts=1的矩形脉冲相卷即为输入信号上下支路分别经过调制后的信号调制后的信号以及加了高斯白噪声的信号调制后的信号乘以相干载波后的信号，通过低通滤波器以及经过抽样判决后的信号（前面上下支路的信号的与判决

5、后的信号对比，可以看出都有延时）恢复成的最终信号与原始信号的对比oqpsk的星座图oqpsk的眼图oqpsk系统的理论误码率此次实验中实际误码率加入不同的噪声进行循环 rt=1.8%（四）matlab程序及其功能主程序以及注释clc;A=1; % 载波幅度fc=2; % 载波频率Ts=1; % 码元宽度fs=1/TsB1=fs; %低通滤波器的宽度N_sample=32; % 基带码元抽样点数N=500; % 码元数 dt=Ts/fc/N_sample; % 抽样时间间隔T=N*Ts; % 信号持续时间长度t=0:dt:T-dt; % 时间向量Lt=length（t）; % 时间向量长度tx1

6、=0; % 时域波形图横坐标起点tx2=10; % 时域波形图横坐标终点ty1=-2; % 时域波形图纵坐标起点ty2=2; % 时域波形图纵坐标终点fx1=-10; % 功率谱图横坐标起点fx2=10; % 功率谱图横坐标终点fy1=-40; % 功率谱图纵坐标起点fy2=25; % 功率谱图纵坐标终点EsN0dB = 3:0.5:10 ; %设定EbNo范围EsN0 = 10.（EsN0dB/10）;rt=zeros（1,length（EsN0dB）; %初始化误码率向量M=4;%产生二进制信源m=randn（1,N）; % 产生1到n的随机数d=sign（m）; % 将大于0的变为1小于

7、0的变为-1dd=sigexpand（d,fc*N_sample）; %将序列d的周期变为Tsgt=ones（1,fc*N_sample）; % 产生宽度为Ts的矩形窗d_NRZ=conv（dd,gt）; % 卷积产生基带信号figure（5）;subplot（2,2,1）;plot（t,dd（1:Lt）;axis（tx1,tx2,ty1,ty2）;figure（1）;subplot（2,4,1）;plot（t,d_NRZ（1:xlabel（时间（S）;ylabel（幅度title（基带信号时域波形图grid;f,d_NRZf=T2F（t,d_NRZ（1:% 进行傅里叶变换subplot（2,

8、4,5）;plot（f,10*log10（abs（d_NRZf）.2/T）;axis（fx1,fx2,fy1,fy2）;频率（Hz）功率谱密度（dB/Hz）基带信号功率谱图figure（4）;subplot（2,2,2）;% 串并转换d1=;d2=;for i=1:N if rem（i,2）=1 d1（i+1）/2）=d（i）; else d2（i/2）=d（i）; endenddd1=sigexpand（d1,2*fc*N_sample）; %功能同上 gt1=ones（1,2*fc*N_sample）; d_NRZ1=conv（dd1,gt1）; subplot（2,4,2）;plot（t,d_NRZ1（1:上支路基带信号时域波形图f1,d_NRZ1f=T2F（t,d_NRZ1（1:subplot（2,4,6）;plot（f1,10*log10（abs（d_NRZ1f）.2/T）;上支路基带信号功率谱图figure（3）;Lt）