oqpsk matlab仿真,课程设计OQPSK通信系统的matlab仿真分析

课程设计OQPSK通信系统的matlab仿真分析

说 明1. 课程设计任务书由指导教师填写，并经专业学科组审定，下达到学生。2. 学生根据指导教师下达的任务书独立完成课程设计。3. 本任务书在课程设计完成后，与论文一起交指导教师，作为论文评阅和课程设计答辩的主要档案资料。一、课程设计的主要内容和基本要求1 产生等概率且相互独立的二进制序列，画出时域和频域的波形；2 产生均值为 0，方差为 1 的加性高斯随机噪声；3 进行 OQPSK 调制，画出波形；4 进行误码率分析，并与理论值比较；5 解调 OQPSK，画出眼图。6 画出星座图二、课程设计图纸内容及张数本实验没有规定的图纸内容，在实验结果中附有本次实验的结果图三、课程设计应完成的软硬件的名称、内容及主要技术指标MATLAB 7.0四、主要参考资料通信原理基础 (北京邮电大学出版社)通信原理 (国防工业出版社)樊昌信 曹丽娜 编著信号与系统——MATLAB 综合实验(高等教育出版社)MATLAB7 辅助信号处理技术与应用(电子工业出版社)飞思科技产品研发中心编著OQPSK 通信系统的 matlab 仿真分析1 设计目的和意义1. 对 oqpsk 进行调制和解调，通过 MATLAB 编程，掌握 MATLAB 的使用，熟练掌握 OQPSK 的调制原理，解调原理。2. 对 OQPSK 通信系统进行 matlab 仿真分析，分析起信噪比和差错率。为现实中通信系统的调制，解调，及信道传输进行理论指导。2 设计原理1. OPSK 的调制它和 QPSK 有着同样的相位关系，也是把输入码流分成两路，然后进行正交调制。随着数字通信技术的发展和广泛应用，人们对系统的带宽、频谱利用率和抗干扰性能要求越来高。而与普通的 QPSK 比较，交错正交相移键控的同相与正交两支路的数据流在时问上相互错开了半个码元周期，而不像 QPSK 那样I、Q 两个数据流在时间上是一致的(即码元的沿是对齐的)。由于 OQPSK 信号中的 I(同相)和 Q(正交)两个数据流，每次只有其中一个可能发生极性转换，所以，每当一个新的输入比特进入调制器的 I 或 Q 信道时，其输出的 OQPSK 信号中只有 0°、+90°三个相位跳变值，而根本不可能出现 180°相位跳变。所以频带受限的 OQPSK 信号包络起伏比频带受限的 QPSK 信号要小，而经限幅放大后的频带展宽也少。2．OQPSK 的基本原理OQPSK 信号的数学公式可以表示为：⑴OQPSK 的调制方法与 QPSK 类似，仅在一条正交支路上引入了一个比特的延时，以使得两支路的数据不会同时发生变化，降低最大相位跳变。其中电平映射关系为:1→1,0→-1.图 3：OQPSK 调制原理框图经 OQPSK 调制后，调制点的星相图以及状态转移图如图 4 所示。⑵OQPSK 解调OQPSK 信号可采用正交相干解调方式解调，其解调原理如图 a 所示。由图 a 可以看出，OQPSK 与 QPSK 信号的解调原理基本相同，其差别仅在于对 Q 支路信号抽样判决后要延迟 Tb／2，这是因为在调制时，Q 支路信号在时间上偏移了Tb／2，所以抽样判决时刻也相应偏移了 Tb／2，以保证对两支路的交错抽样。⑶眼图是信号由垂直扫描进入与同周期的水平扫描锯齿波叠加到示波器上时到得图案。眼图能够反映信号在传输过程中受到的信道噪声影响的强度，眼图越模糊，眼睛越闭合，则说明噪声越强，反之，则说明噪声强度弱，也能说明信道性能更优良。3 设计结果及分析⑴输入的二进制序列以及串并转换成的上下两路信号，下支路信号已延时，由产生的随机序列可以知道前十个码元为 1001001101，映射电平为 1-1-11-1-111-11，下面的抽样脉冲与 ts=1 的矩形脉冲相卷即为输入信号⑵上下支路分别经过调制后的信号⑶调制后的信号以及加了高斯白噪声的信号⑷调制后的信号乘以相干载波后的信号，通过低通滤波器以及经过抽样判决后的信号(前面上下支路的信号的与判决后的信号对比，可以看出都有延时)⑸恢复成的最终信号与原始信号的对比⑹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=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 小于 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:Lt));axis([tx1,tx2,ty1,ty2]);xlabel( 时间(S) );ylabel( 幅度 );title( 基带信