课程设计OQPSK通信系统的matlab仿真分析Word格式文档下载.doc
《课程设计OQPSK通信系统的matlab仿真分析Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《课程设计OQPSK通信系统的matlab仿真分析Word格式文档下载.doc(22页珍藏版)》请在冰豆网上搜索。
通信原理(国防工业出版社)樊昌信曹丽娜编著
信号与系统——MATLAB综合实验(高等教育出版社)
MATLAB7辅助信号处理技术与应用(电子工业出版社)飞思科技产品研发中心编著
窗体顶端
窗体底端
OQPSK通信系统的matlab仿真分析
㈠设计目的和意义
1.对oqpsk进行调制和解调,通过MATLAB编程,掌握MATLAB的使用,熟练掌握OQPSK的调制原理,解调原理。
2.对OQPSK通信系统进行matlab仿真分析,分析起信噪比和差错率。
为现实中通信系统的调制,解调,及信道传输进行理论指导。
㈡设计原理
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,以保证对两支路的交错抽样。
⑶眼图是信号由垂直扫描进入与同周期的水平扫描锯齿波叠加到示波器上时到得图案。
眼图能够反映信号在传输过程中受到的信道噪声影响的强度,眼图越模糊,眼睛越闭合,则说明噪声越强,反之,则说明噪声强度弱,也能说明信道性能更优良。
㈢设计结果及分析
⑴输入的二进制序列以及串并转换成的上下两路信号,下支路信号已延时,由产生的随机序列可以知道前十个码元为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/Ts
B1=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的变为-1
dd=sigexpand(d,fc*N_sample);
%将序列d的周期变为Ts
gt=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,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=[];
fori=1:
N
ifrem(i,2)==1
d1((i+1)/2)=d(i);
else
d2(i/2)=d(i);
end
end
dd1=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))