CN104079524B - 一种基于ofdm的qam调制下的已失真通信信号识别方法 - Google Patents
一种基于ofdm的qam调制下的已失真通信信号识别方法 Download PDFInfo
- Publication number
- CN104079524B CN104079524B CN201410355612.2A CN201410355612A CN104079524B CN 104079524 B CN104079524 B CN 104079524B CN 201410355612 A CN201410355612 A CN 201410355612A CN 104079524 B CN104079524 B CN 104079524B
- Authority
- CN
- China
- Prior art keywords
- signal
- cluster
- barycenter
- training sequence
- judgement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种基于OFDM的QAM调制下的已失真通信信号识别方法,它包括发送端的训练符号添加和接收端已失真信号的判别步骤。本发明简化了平坦瑞利衰落信道下接收机的结构,接收机无需使用信道估计和信道均衡模块,使用导频信号对现有k‑means聚类算法进行了初始聚类质心的改进,并将其应用于已失真通信信号的判决和识别,具有迭代次数少,计算量小等特点,对信号相位失真和幅度失真的判决和识别性能更优。
Description
技术领域
本发明属于移动通信技术领域,具体涉及一种在多径衰落信道下基于正交频分复用(OFDM)的正交幅度调制(QAM)系统的信号聚类识别方法。
背景技术
正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)是一种多载波复用调制技术,它将高比特率的串行数据流进行分流,实现低速率多状态的数据符号来调制对应的子载波,从而实现数据的并行传输。由于它的各个子载波是相互正交的,允许相互之间的频谱发生重叠现象,这在很大程度上实现对频谱利用率的提高,也减小了子载波间的相互干扰,因此在现代通信中得到广泛应用。
对于无线通信而言,消息信号到达接收端,可能会经过多种障碍物,以至于在接收端收到的信号会来自直射路径、反射路径、绕射路径和散射路径等多种路径的综合。由于不同路径所经历的路程基本不同,导致到达的时间有所不同,因此每条路径的到达相位也不相同,产生多径效应现象,再加上信号在传输过程中,也会受到一定程度的衰减影响,即频率选择性衰落、时间选择性衰落和空间选择性衰落。由此会在接收端的综合信号某些部分会发生消减,而在另外一部分会发生加强现象。如果发送端与接收端之间再有相对的移动,也会产生远近效应和多普勒效应,对于接收端的相对移动,有时甚至会使接收端处于深衰落的位置,导致在接收端不能正常接收信号。
对基于OFDM的QAM调制后的发送信号进行建模,调制后产生的数据经过串/并的变换,从而产生并行的N个频域符号X(k),使其通过IFFF变换,进行子载波调制,产生的时域信号为x(n),则时域符号可以表示为
假设信道为一个多径时变衰落信道,可用线性抽头延迟冲激响应来表示
式中,hl(t)为t时刻时域多径信道冲激响应的第l条路径的复增益,fD(t)是多普勒频移,为第l条路径的时变相移,πl(t)为第l条路径的时变等效延迟函数。若信道是慢衰落信道,即在一帧OFDM符号内信号基本不变,则hl(t)与πl(t)在一个符号周期Ts内近似为恒值,同理fD(t)和θl(t)也近似为恒值,故把它们记为hl、πl、fD、θl,并设θl=0,因此公式(2)可以表示为
再将接收时由于不同步引起的载波频偏fD和高斯白噪声考虑进去,则得到接收的时域信号为
式中,w(n)为加性高斯白噪声,ε=ΔfTs=(fD+fG)Ts为由多普勒频偏和载波频偏共同决定的归一化频偏,把公式(1)代入公式(4)可得
式中,表示第k个子载波信道的频率冲激响应函数,在接收端,将添加的循环前缀去除掉,然后对y(n)进行FFT变换,可以得到它的频域形式为
式中,0≤k≤N-1,W(k)是对高斯白噪声w(n)的FFT变换,从公式(6)中可以看出,多径信道对信号造成了相位和幅度失真,载波相位偏移将使接收星座图的样点发生旋转,采用基本的信道估计(频域最小二乘法,LS)或均衡算法能在一定程度上补偿信号的失真,并不能获得更好的效果。为了解决上述问题,现有方法主要是改进信道估计算法或者改进信道均衡算法,然而这些方法都将在一定程度上增加接收端电路的复杂度。
发明内容
本发明的目的在于针对平坦瑞利衰落信道下采用OFDM技术和QAM调制的信号,提出使用一种改进的聚类算法用于失真信号的识别。在无需使用信道估计和信道均衡情况下,仍能获得较好的识别性能。
本发明的目的是通过以下技术方案来实现的:一种基于OFDM的QAM调制下的已失真通信信号识别方法,其特征在于:它包括发送端的训练符号添加和接收端的失真信号识别步骤,其中,所述的发送端的训练符号添加包括以下子步骤:
S11,对拟发送数据依次进行QAM调制、串/并转换和IFFT快速傅立叶逆变换;
S12,在经步骤S11处理后得到的数据中,加入循环前缀和训练序列:
S13,对加入循环前缀和训练序列后的数据依次进行并/串转换、I/Q调制和数模转换,并通过无线信道发射;
所述的接收端已失真信号判别包括以下子步骤:
S21,接收发送端发射的含有循环前缀和训练序列的射频信号,并依次进行模数转换、I/Q解调和串/并转换,得到含有循环前缀和训练序列的基带信号;
S22,去除基带信号中的循环前缀,提取训练序列,得到无循环前缀和训练序列的基带信号;
S23,对步骤S22得到的基带信号进行FFT快速傅立叶变换和并/串转换处理;
S24,对步骤S23处理得到的数据进行聚类判决,正确识别出接收信号。
发送端添加训练符号时采用均匀分布,对于K阶QAM调制,Ns为1帧中频域训练序列的长度,则训练序列的组成为每个星座个数据点,剩余的个数据点随机产生,其中,为向下取整。
所述的聚类判决包括以下子步骤:
S241:初始化聚类质心,从解调的基带信号中分离出导频信号,通过导频信号计算出初始聚类质心,并对每个质心进行编号;
S242:信号聚类,分别计算信号点与各聚类质心的欧式距离,为每个信号点选取最短距离所对应的聚类质心的编号作为自己的簇编号;
S243:更新聚类质心,根据簇编号对每一簇的信号点进行求平均,计算出下次迭代使用的聚类质心;
S244:迭代停止判断,计算新质心与原质心的差值,并判断迭代是否结束,若需继续迭代,则返回信号聚类步骤S242,若不需继续迭代,则结束。
本发明的有益效果是:(1)本发明简化了平坦瑞利衰落信道下接收机的结构,接收机无需使用信道估计和信道均衡模块。(2)使用导频信号对现有k-means聚类算法进行了初始聚类质心的改进,并将其应用于失真通信信号的判决和识别,具有迭代次数少,计算量小等特点,对信号相位失真和幅度失真的判决和识别性能更优。
附图说明
图1为本发明系统结构图;
图2为本发明与采用LS信道估计和ZF信道均衡后在瑞利平坦衰落信道下信号判决性能比较示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案:如图1所示,一种基于OFDM的QAM调制下的已失真通信信号识别方法,其特征在于:它包括发送端的训练符号添加和接收端的失真信号识别步骤,其中,所述的发送端的训练符号添加包括以下子步骤:
S11,对拟发送数据依次进行QAM调制、串/并转换和IFFT快速傅立叶逆变换;
S12,在经步骤S11处理后得到的数据中,加入循环前缀和训练序列:
S13,对加入循环前缀和训练序列后的数据依次进行并/串转换、I/Q调制和数模转换,并通过无线信道发射;
所述的接收端已失真信号判别包括以下子步骤:
S21,接收发送端发射的含有循环前缀和训练序列的射频信号,并依次进行模数转换、I/Q解调和串/并转换,得到含有循环前缀和训练序列的基带信号;
S22,去除基带信号中的循环前缀,提取训练序列,得到无循环前缀和训练序列的基带信号;
S23,对步骤S22得到的基带信号进行FFT快速傅立叶变换和并/串转换处理;
S24,对步骤S23处理得到的数据进行聚类判决,正确识别出接收信号。
发送端添加训练符号时采用均匀分布,对于K阶QAM调制,Ns为1帧中频域训练序列的长度,则训练序列的组成为每个星座个数据点,剩余的个数据点随机产生,其中,为向下取整。
所述的聚类判决包括以下子步骤:
S241:初始化聚类质心,从解调的基带信号中分离出导频信号,通过导频信号计算出初始聚类质心,并对每个质心进行编号;
S242:信号聚类,分别计算信号点与各聚类质心的欧式距离,为每个信号点选取最短距离所对应的聚类质心的编号作为自己的簇编号;
S243:更新聚类质心,根据簇编号对每一簇的信号点进行求平均,计算出下次迭代使用的聚类质心;
S244:迭代停止判断,计算新质心与原质心的差值,并判断迭代是否结束,若需继续迭代,则返回聚类判决步骤S242,若不需继续迭代,则结束。
在接收端先后需进过IQ解调、去除循环前缀、FFT、并/串转换,然后通过聚类判决解调出接收符号,聚类判决步骤的具体描述为:
(1)计算初始聚类质心:发送的训练序列已知为且T中的元素可分为K类,即ti∈Uk,i=0,1,…,NS,k=1,2,…,K,其中Uk为信号聚类的第k类的集合。从接收信号中分离出训练序列为对接收到的训练序列进行分类,分类原则为{tj'∈Uk|ti∈Uk,i=0,1,…,NS,j=1,2,…,K},再根据分类结果计算初始聚类质心,假设Nk为集合Uk中的元素个数,计算质心的公式如下
(2)根据聚类质心对信号进行聚类:分别计算信号点xi与各聚类质心ck的欧式距离d(*),为每个信号点选取最短距离所对应的质心编号为聚类编号,即
(3)更新聚类质心:分别对属于每一类的数据点xi求平均,计算出新的聚类质心,计算公式如下
(4)判断迭代是否结束:计算新旧质心之间的误差和判断迭代是否结束,若E≤ε,则算法收敛迭代结束并继续下一步,若E>ε,则返回执行步骤(2)(3)继续迭代,直至算法收敛,ε为可设置的迭代停止阈值。
(5)聚类判决结束,对xi∈Uk,有xi星座编号判决为k。
在本实例中,发送端OFDM系统的参数如下:系统子载波数量为52,FFT长度为64,循环前缀长度为16,每帧中包含的OFDM符号数为6,其中训练符号数为1,调制方式为16QAM,多普勒频移为100Hz。图2为本专利方法的判决性能与采用LS信道估计和ZF信道均衡后判决性能的对比图。
Claims (2)
1.一种基于OFDM的QAM调制下的已失真通信信号识别方法,其特征在于:它包括发送端的训练符号添加和接收端已失真信号的判别步骤,其中,所述的发送端的训练符号添加包括以下子步骤:
S11,对拟发送数据依次进行QAM调制、串/并转换和IFFT快速傅立叶逆变换;
S12,在经步骤S11处理后得到的数据中,加入循环前缀和训练序列:
S13,对加入循环前缀和训练序列后的数据依次进行并/串转换、I/Q调制和数模转换,并通过无线信道发射;
所述的接收端已失真信号判别包括以下子步骤:
S21,接收发送端发射的含有循环前缀和训练序列的射频信号,并依次进行模数转换、I/Q解调和串/并转换,得到含有循环前缀和训练序列的基带信号;
S22,去除基带信号中的循环前缀,提取训练序列,得到无循环前缀和训练序列的基带信号;
S23,对步骤S22得到的基带信号进行FFT快速傅立叶变换和并/串转换处理;
S24,对步骤S23处理得到的数据进行聚类判决,正确识别出接收信号;
所述的聚类判决包括以下子步骤:
S241:初始化聚类质心,从解调的基带信号中分离出导频信号,通过导频信号计算出初始聚类质心,并对每个质心进行编号;具体地,发送的训练序列已知为且T中的元素可分为K类,即ti∈Uk,i=0,1,…,NS,k=1,2,…,K,其中Uk为信号聚类的第k类的集合;从接收信号中分离出训练序列为对接收到的训练序列进行分类,分类原则为{tj'∈Uk|ti∈Uk,i=0,1,…,NS,j=1,2,…,K},再根据分类结果计算初始聚类质心,假设Nk为集合Uk中的元素个数,计算质心的公式如下:
S242:信号聚类,分别计算信号点xi与各聚类质心ck的欧式距离d(*),为每个信号点选取最短距离所对应的聚类质心的编号作为自己的簇编号,即:
S243:更新聚类质心,根据簇编号对每一簇的信号点xi进行求平均,计算出下次迭代使用的聚类质心,计算公式如下:
S244:迭代停止判断,计算新旧质心之间的误差和判断迭代是否结束,若E≤ε,则算法收敛迭代结束并继续下一步,若E>ε,则返回执行步骤S242和步骤S243继续迭代,直至算法收敛,ε为可设置的迭代停止阈值;
S245:聚类判决结束,对xi∈Uk,有xi星座编号判决为k。
2.根据权利要求1所述的一种基于OFDM的QAM调制下的已失真通信信号识别方法,其特征在于:发送端添加训练符号时采用均匀分布,对于K阶QAM调制,Ns为1帧中频域训练序列的长度,则训练序列的组成为每个星座个数据点,剩余的个数据点随机产生,其中,为向下取整。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410355612.2A CN104079524B (zh) | 2014-07-24 | 2014-07-24 | 一种基于ofdm的qam调制下的已失真通信信号识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410355612.2A CN104079524B (zh) | 2014-07-24 | 2014-07-24 | 一种基于ofdm的qam调制下的已失真通信信号识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104079524A CN104079524A (zh) | 2014-10-01 |
CN104079524B true CN104079524B (zh) | 2017-07-28 |
Family
ID=51600573
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410355612.2A Expired - Fee Related CN104079524B (zh) | 2014-07-24 | 2014-07-24 | 一种基于ofdm的qam调制下的已失真通信信号识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104079524B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107017929B (zh) * | 2017-04-28 | 2020-02-18 | 电子科技大学 | Mimo系统信号发送和接收方法 |
CN107566039A (zh) * | 2017-09-04 | 2018-01-09 | 复旦大学 | 一种基于聚类判决的可见光系统非线性补偿方法 |
CN107707494B (zh) * | 2017-10-10 | 2020-02-11 | 苏州大学 | 用于64-qam相干光通信系统的光纤非线性均衡方法 |
CN108900460B (zh) * | 2018-06-12 | 2020-11-13 | 南京邮电大学 | 一种基于k均值聚类的抗相位噪声的鲁棒符号检测方法 |
CN109088837B (zh) * | 2018-08-06 | 2021-04-09 | 河海大学 | 一种基于聚类分析的多种雷达与通信信号分类方法 |
CN111131108B (zh) * | 2019-12-30 | 2022-05-20 | 哈尔滨工程大学 | 一种非合作水声ofdm子载波调制方式识别方法 |
CN111181880B (zh) * | 2020-01-03 | 2022-09-02 | 深圳大学 | 基于集成聚类的无线通信系统的相位噪声补偿方法及系统 |
CN115913868A (zh) * | 2022-11-25 | 2023-04-04 | 中科南京移动通信与计算创新研究院 | 一种基于聚类的矢量调制信号分析方法及其系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6314146B1 (en) * | 1998-06-05 | 2001-11-06 | The Board Of Trustees Of The Leland Stanford Junior University | Peak to average power ratio reduction |
CN101835167A (zh) * | 2010-05-01 | 2010-09-15 | 华中科技大学 | 用于非连续频谱ofdm动态频谱接入的频谱同步方法 |
CN101867425A (zh) * | 2010-07-08 | 2010-10-20 | 华中科技大学 | 用于非连续ofdm动态频谱接入的可用子载波检测同步方法 |
CN103873406A (zh) * | 2014-04-09 | 2014-06-18 | 厦门大学 | 水声正交频分复用通信系统帧间干扰消除方法 |
-
2014
- 2014-07-24 CN CN201410355612.2A patent/CN104079524B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6314146B1 (en) * | 1998-06-05 | 2001-11-06 | The Board Of Trustees Of The Leland Stanford Junior University | Peak to average power ratio reduction |
CN101835167A (zh) * | 2010-05-01 | 2010-09-15 | 华中科技大学 | 用于非连续频谱ofdm动态频谱接入的频谱同步方法 |
CN101867425A (zh) * | 2010-07-08 | 2010-10-20 | 华中科技大学 | 用于非连续ofdm动态频谱接入的可用子载波检测同步方法 |
CN103873406A (zh) * | 2014-04-09 | 2014-06-18 | 厦门大学 | 水声正交频分复用通信系统帧间干扰消除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104079524A (zh) | 2014-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104079524B (zh) | 一种基于ofdm的qam调制下的已失真通信信号识别方法 | |
CN1777166B (zh) | 在正交频分复用的移动通信系统中获取同步的装置及方法 | |
CN102404268B (zh) | 高速移动环境下莱斯信道中多普勒频偏估计与补偿方法 | |
CN106161328B (zh) | 基于载波索引调制的mimo-ofdm系统的检测方法 | |
CN101778069A (zh) | 一种新型ofdm信道估计联合ici自消除方法 | |
CN103051578B (zh) | 带有ici消除的迭代差错扩散判决ofdm信道估计方法 | |
CN104735014B (zh) | 一种基于前导符号差分相关的定时同步方法 | |
US10334605B2 (en) | Time domain pilot of single-carrier MIMO system and synchronization method thereof | |
CN104253772B (zh) | 正交频分复用系统的信道估计方法 | |
CN100493056C (zh) | 带时域包络加权的正交频分复用系统的频域信道估计方法 | |
CN101895309B (zh) | 一种基于信噪比估计的频偏调整方法 | |
CN107124194A (zh) | 海上数字中频发射设备 | |
KR20060009434A (ko) | 빠른 레일리 페이딩 채널 환경의 이동통신 시스템에서최우 주파수 오프셋 추정 방법 | |
CN102113285A (zh) | 用于多载波系统中的分布式资源分配的简化均衡方案 | |
CN101442512A (zh) | Ofdm信号的识别方法 | |
CN104022980B (zh) | 一种ofdm系统信干噪比盲估计方法及系统 | |
CN106789774B (zh) | 用于多载波系统的信道估计方法 | |
CN102664858A (zh) | Ofdm系统降低峰均比与载波频率跟踪联合方法 | |
KR101341611B1 (ko) | 최대 우도 판정 기반 복호 장치 및 방법, 및 무선 통신 시스템에서 최대 우도 판정을 기반으로 하는 수신 장치 및 방법 | |
CN102223336B (zh) | 无线通信方法和设备 | |
CN106302279A (zh) | 基于干扰方差统计的fbmc系统均衡方法 | |
Liu et al. | Study on channel estimation technology in OFDM system | |
CN111817990A (zh) | Ofdm系统中基于最小均方误差的一种信道估计改进算法 | |
CN104348765B (zh) | 信道估计方法 | |
Charrada et al. | Complex support vector machines regression for robust channel estimation in LTE downlink system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170728 Termination date: 20200724 |