CN113904713A - 一种卫星通信系统低速ofdm突发信号解调方法 - Google Patents
一种卫星通信系统低速ofdm突发信号解调方法 Download PDFInfo
- Publication number
- CN113904713A CN113904713A CN202111136463.7A CN202111136463A CN113904713A CN 113904713 A CN113904713 A CN 113904713A CN 202111136463 A CN202111136463 A CN 202111136463A CN 113904713 A CN113904713 A CN 113904713A
- Authority
- CN
- China
- Prior art keywords
- value
- subcarrier
- timing
- estimation
- burst signal
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004891 communication Methods 0.000 title claims abstract description 27
- 239000000969 carrier Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 5
- 238000007619 statistical method Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000001502 supplementing effect Effects 0.000 claims description 4
- 238000011084 recovery Methods 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 11
- 230000000875 corresponding effect Effects 0.000 description 11
- 238000009825 accumulation Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 102100026758 Serine/threonine-protein kinase 16 Human genes 0.000 description 1
- 101710184778 Serine/threonine-protein kinase 16 Proteins 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
- H04L27/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Radio Relay Systems (AREA)
Abstract
本发明公开了一种卫星通信系统低速OFDM突发信号解调方法,它涉及卫星通信系统体制设计、信关站和终端设计技术。本发明针对卫星通信系统在卫星上下行链路中存在频谱占用,特别是复杂干扰或其他因素造成信道资源受限、信道质量较差时,根据干扰情况在有限的信道资源中,利用单个或数个非连续OFDM子载波进行信号传输和解调处理,并在解调中采用能够适应非连续极少个数子载波的频偏估计和定时估计方法,从而准确的完成载波恢复和数据解调。本发明具有信道灵活、适应性强、误码性能高的特点,适合信道资源受限情况下卫星通信系统的低速OFDM突发信号解调和信道设计。
Description
技术领域
本发明属于卫星通信技术领域,涉及一种低速OFDM突发信号解调方法。特别适用于卫星通信系统信道资源深度受限情况下的信关站及终端低速OFDM突发信号的解调。
背景技术
随着卫星通信系统的发展,系统传输在实际运行中的信道条件和信号带宽差别较大,有的传输需要大带宽高速率,有的传输只需要较少带宽、较低速率,而有的传输需要保证低速下的基本通信需求。常规的OFDM解调方法一般适用于宽带的连续频率信号,在信道资源深度受限、可用带宽极低、可用频谱离散和不固定的情况下不能适用,而传统的单载波MPSK调制方法,难以在频域上进行非连续扩展(单个用户难以利用不连续的频域资源)。因此本发明根据信道条件,提出结合单载波和多载波系统优点,能够适应个数极少的离散子载波,包括只有单个子载波时的低速OFDM突发信号解调方法。同时由于信道估计对解调性能影响很大,传统的OFDM信道估计方法通常利用频偏和时偏产生的相位偏移,基于宽带连续多子载波进行估计和抵抗噪声影响,但卫星通信系统传输距离远,信噪比低,特别是信道资源深度受限情况下,可利用传输的子载波个数极少且位置离散、随机,传统方法不能适用。为此本发明又提出了基于单子载波频域符号的低速OFDM突发信号频率估计方法和基于统计和置信概率的极少个数子载波低速OFDM突发信号定时估计方法,在低信噪比下仍可较好的消除噪声影响,提高误码性能。
因此,设计一种基于极少子载波的低速OFDM突发信号解调方法,对卫星通信系统在复杂条件下应用,特别是信道资源深度受限情况下的卫星通信传输具有重要意义。
发明内容
本发明针对卫星通信系统在信道资源深度受限情况下,只有个数极少的离散子载波甚至单个子载波可利用的特点,根据信道情况,利用非连续的数个或单个子载波,将频偏估计、定时估计和时频资源利用、多用户复用等问题统一考虑,提出能够适应信道资源深度受限情况的OFDM解调方法,使用基于单子载波频域符号的OFDM突发信号频率估计方法进行频率估计使用基于极少数子载波的能量峰值方法进行定时估计,并使用统计方法进行多子载波和多帧的二维统计平滑和基于置信概率的定时调整值判断,从而准确完成信道估计和数据解调。
为了实现上述目的,本发明所采用的技术方案是:
一种卫星通信系统低速OFDM突发信号解调方法,包括如下步骤:
(1)对接收的突发信号进行解映射和子载波提取,得到各子载波的频域符号数据;
(2)对各子载波的频域符号数据进行独特字符号提取处理;
(3)对各子载波的独特字频域符号进行去调制运算;
(4)对各子载波的去除调制的独特字频域符号,分别进行基于单子载波频域符号的频率估计;
(5)对各单子载波频域符号的频率估计结果,进行多子载波和多帧的二维统计平滑,得到准确的突发信号频率;
(6)利用准确的突发信号频率,对突发信号进行去频差运算;
(7)对去除调制和频差的独特字频域符号,进行基于统计和置信概率的极少个数子载波定时估计,得到定时调整值;
(8)判断定时调整值是否超出设定的阈值范围,若超出则进行定时调整,当本端为终端侧时,利用定时调整值调整接收突发信号的定时误差,本端为基站侧时,将定时调整值反馈给终端,终端调整发送突发信号的定时误差,本端对当次突发的解调继续转入步骤(9),并返回步骤(1)接收调整后的突发信号;若没有超出则直接转入步骤(9);
(9)利用去除调制和频差的独特字频域符号,进行基于单子载波的信噪比自适应的信道估计和均衡,得到数据字的解调结果;
完成卫星通信系统低速OFDM突发信号解调。
其中,所述步骤(4)中基于单子载波频域符号的频率估计,具体包括以下步骤:
(401)对全部M个子载波按位置排列并编号,编号从0到M-1,同时子载波计数器NSC清零;其中,当只有一个子载波时,编号为0;
(402)判断去调制后独特字频域符号是否等间隔排列,如果是则转入步骤(404),否则转入步骤(403);
(403)根据独特字的位置特征,内插虚拟独特字符号,使内插后的独特字符号和虚拟独特字符号等间隔排列;
(404)将等间隔排列的独特字符号和虚拟独特字符号序列补零,使序列长度等于2^NFFT,其中NFFT为阶数;
(405)对步骤(404)的输出序列,进行快速傅里叶变换,得到长度为2^NFFT的输出序列;
(406)对步骤(405)的输出序列,进行功率值寻峰,并记录功率峰值位置;
(407)对功率峰值位置进行数值变换,得到当前子载波的频率估计值;
(408)子载波计数器NSC加1,判断NSC是否大于M-1,如果否,则转入步骤(402),如果是,则结束。
其中,所述步骤(7)具体包括以下步骤:
(701)对去调制和去频差的独特字频域符号,每个子载波分别进行时域上的累加,得到各子载波上的累加值;
(702)将计数器POS值清零;
(703)以计数器POS值作为定时估计值,根据各个子载波的位置,计算定时偏移在各子载波上产生的相位;
(704)对步骤(701)中各子载波上的累加值,根据步骤(703)中得到的相位,去除累加值中的定时偏移相位;
(705)对各个子载波的去除定时偏移相位的累加值相加,得到的结果记为SUM,记录计数器POS值和对应的SUM值;
(706)将计数器POS值加1,并判断POS值是否大于定时偏移最大值PMAX,如果否,则转入步骤(703),如果是,则转入步骤(707);
(707)对所有POS值对应的SUM值进行功率值寻峰,并记录峰值对应的POS值;
(708)对步骤(707)的输出值进行数值变换,得到当次定时估计值,返回步骤(1)重复进行N次;其中N为设定数值;
(709)对步骤(708)得到的N次定时估计值进行基于置信概率的统计分析,得到定时调整值。
其中,所述步骤(709)具体包括以下步骤:
(7091)对N次的定时估计值依次排列,作为统计序列输入;
(7092)计数器K清零;
(7093)对输入序列中序号为K的点,统计落在以该点为中心,前后F点范围内的个数,记为集中度C(K);其中,F为设定值;
(7094)计数器K加1,判断K是否大于N-1,如果否,则转入步骤(7093),如果是,则转入步骤(7095);
(7095)将集中度C(K)从高到低排列,得到集中度最大值CMAX;
(7096)判断CMAX是否大于置信阈值,如果是,则定时估计有效,否则定时估计无效;
(7097)如果定时估计无效,则将统计序列末位值调出,补充最新一次的定时估计值,然后循环步骤(7092)至(7096),直至定时估计有效;如果定时估计有效,对集中度最大值对应的点,计算以该点为中心,前后F点范围内的定时估计值的平均数,作为定时调整值。
本发明相比现有技术具有如下优点:
本发明能够针对卫星通信系统中只有个数极少的离散子载波或单个子载波的OFDM突发信号进行处理,从而准确完成信道估计和数据解调。
本发明在存在严重干扰和信道资源深度受限情况下,仍可保证卫星通信系统的传输,得到准确的载波恢复和解调结果,满足卫星通信系统多用户实时通信的需求。
附图说明
图1为本发明的流程图;
图2为本发明突发信号结构示意图;
图3为本发明基于单子载波频域符号的频率估计方法流程图;
图4为本发明非连续极少个数子载波的定时估计方法流程图;
图5为基于统计和置信概率的的定时估计值处理流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用于解释本发明,并不用于限定本发明。
图2显示了卫星通信系统在信道资源深度受限情况下,只有个数极少的离散子载波甚至单个子载波可利用时的突发信号结构示意图。图中上部显示了突发信号的频域特征,其中深色表示突发信号所占的传输子载波,其利用有限可用的信道资源穿插离散分布在宽带深度受限的信道之间,在深度受限的信道频域位置分布空子载波,以浅色表示,传输子载波和空子载波的位置可随时根据信道受限情况实时调整。图中下部显示了对于某个传输子载波N,突发信号的时域特征。其中DATA为传输的数据字符号,UW1,2..K为传输的独特字符号,突发信号在时域上可长可短,可连续也可不连续,独特字符号可以等间隔排列,也可以非等间隔排列。
参照图1流程框图,一种卫星通信系统低速OFDM突发信号解调方法,适用于信道资源深度受限情况下的卫星通信系统,其包括以下步骤:
(1)对接收的突发信号进行解映射和子载波提取,得到各子载波的频域符号数据;
(2)对各子载波的频域符号数据进行独特字符号提取处理;
(3)对各子载波的独特字频域符号进行去调制运算;
记接收信号提取的独特字频域符号序列为UWrecv(n),本地已知的独特字符号序列为UWlocal(n),n∈{1,2,...L},L为独特字个数,那么相关运算公式为:
UWdemod(n)=UWrecv(n)×UWlocal *(n)
其中,UWlocal *(n)是UWlocal(n)的复共轭,UWdemod(n)为去调制后的独特字符号。由于接收信号的独特字符号序列在经过信道之前与本地已知的独特字符号序列相同,其与自身的复共轭序列求相关,即去除了发送调制信息。
(4)对各子载波的去除调制的独特字频域符号,分别进行基于单子载波频域符号的频率估计;
(5)对各单子载波频域符号的频率估计结果,进行多子载波和多帧的二维统计平滑,得到准确的突发信号频率;
(6)利用准确的突发信号频率,对突发信号进行去频差运算;
(7)对去除调制和频差的独特字频域符号,进行基于统计和置信概率的极少个数子载波定时估计,得到定时调整值;
(8)判断定时调整值是否超出设定的阈值范围,若超出则进行定时调整,当本端为终端侧时,利用定时调整值调整接收突发信号的定时误差,本端为基站侧时,将定时调整值反馈给终端,终端调整发送突发信号的定时误差,本端对当次突发的解调继续转入步骤(9),并返回步骤(1)接收调整后的突发信号;若没有超出则直接转入步骤(9);
(9)利用去除调制和频差的独特字频域符号,进行基于单子载波的信噪比自适应的信道估计和均衡,得到数据字的解调结果;
参照图3,所述步骤(4)中基于单子载波频域符号的频率估计具体包括以下步骤:
(401)对全部M个子载波按位置排列并编号,编号从0到M-1,同时子载波计数器NSC清零;其中,当只有一个子载波时,编号为0;
(402)判断去调制后独特字频域符号是否等间隔排列,如果是则转入步骤(404),否则转入步骤(403);
不妨假设条件一:突发信号每帧7个OFDM符号,共N帧,突发信号符号序列记为C(n),n∈{1,2,...N},独特字占每帧第一个符号,记为d1(i),i∈{1,2,...N},其中d1(i)=C(7×(i-1)+1)。则此时独特字符号为等间隔排列,转入步骤(404);
假设条件二:突发信号每帧7个OFDM符号,共N帧,突发信号符号序列记为C(n),n∈{1,2,...N},独特字占每帧第二个符号和第五个符号,记为d2(i),i∈{1,2,...2N},其中d2(i)=C(7×k+2),当i=2×k;d2(i)=C(7×k+5),当i=2×k+1;k∈{0,1,2,...N}。此时独特字符号为非等间隔排列,转入步骤(403);
(403)根据独特字的位置特征,内插虚拟独特字符号,使内插后的独特字符号和虚拟独特字符号等间隔排列;
根据假设条件二,内插的虚拟独特字符号为:
d3(i)=UWdemod(n),当i=7×k+1or i=7×k+3;n=2k+1
d3(i)=UWdemod(n),当i=2×k+4or i=7×k+6;n=2k+2
当i=2×k+7;n1=2k+2,n2=2(k+1)+1,
其中k∈{0,1,2,...N};
内插后等间隔排列的独特字符号/虚拟独特字符号记为duwdemod(i):
duwdemod(i)=UWdemod(i),当i=7×k+2or i=7×k+5;
duwdemod(i)=d3(i),当i=else;其中k∈{0,1,2,...N};
(404)将等间隔排列的独特字符号和虚拟独特字符号序列补零,使序列长度等于2^NFFT,其中NFFT为阶数;
记补零后的序列为fi(n),当由假设条件一转入步骤(406):
fi(n)=UWdemod(n),n∈{1,2,...N}
fi(n)=0,n∈{N+1,N+2,...NFFT}
当由假设条件二转入步骤(406):
fi(n)=duwdemod(n),n∈{1,2,...7N}
fi(n)=0,n∈{7N+1,7N+2,...NFFT}
(405)对步骤(404)的输出序列,进行快速傅里叶变换,得到长度为2^NFFT的输出序列;
对输入信号fi(n)进行NFFT点的离散傅里叶变换,处理公式为:
其中Y(k)为输出序列,k∈{1,2,...NFFT};
(406)对步骤(405)的输出序列,进行功率值寻峰,并记录功率峰值位置;
处理公式为:Pmax=Pindex{max(|Y(k)|)},其中Pmax为输出的峰值位置,max(|Y(k)|)为Y(k)功率值序列的最大值,Pindex{.}函数表示功率最大值对应的位置索引值;
(407)对功率峰值位置进行数值变换,得到当前子载波的频率估计值;
当由假设条件一转入步骤(403),则数值变换公式为:
fest=fsc+Δf
其中Ts为突发信号的符号周期,fest为当前子载波的频率估计值;
当由假设条件二转入步骤(403),则数值变换公式为:
fest=fsc+Δf
其中Ts为突发信号的符号周期,fest为当前子载波的频率估计值;
(408)子载波计数器NSC加1,判断NSC是否大于M-1,如果否,则转入步骤(402),如果是,则结束。
参照图4,所述步骤(7)中基于统计和置信概率的极少个数子载波定时估计,其具体方式为:
(701)对去调制的独特字符号,根据步骤(5)得到的估计频率,进行去频差处理,然后对其在每个子载波分别进行时域上的累加,得到各子载波上的累加值;
例如满足假设条件一时,记去调制和频偏补偿后的独特字符号为UWcomp(n),则去频差处理公式为: 其中UWdemod(n)为去调制的独特字符号, n∈{1,2,...L},L为每个子载波上的独特字个数,Ts为突发信号的符号周期;
对UWcomp(n)进行每个子载波上的累加,公式为:
其中Ssc(Nsc)为第Nsc个子载波上的累加值,Nsc∈{1,2,...M}
(702)计数器POS值清零;
(703)以计数器POS值作为定时估计值,根据各个子载波的位置,计算定时偏移在各子载波上产生的附加相位;
不妨假设子载波位置序列为{p1,p2,...pM},其中M为子载波个数,pnsc∈{1,2,...NFFT};
此时由定时偏移产生的各子载波上的附加相位为:
(704)对(701)步中各子载波上的累加值,根据(703)步中得到的相位,去除累加值中的定时偏移相位;
(705)对各个子载波的去除定时偏移相位的累加值相加,得到的结果记为Psum,记录计数器POS值和对应的Psum值;
处理公式为
(706)计数器POS值加1,判断POS是否大于定时偏移最大值PCP,如果否则转入步骤(703),如果是则转入步骤(707);
(707)对所有POS的Psum值进行功率值寻峰,并记录峰值对应的POS值;
处理公式为:Pmax=Pindex{max(|Psum(POS)|)},其中Pmax为输出的峰值位置,max(|Psum(POS)|)为Psum(POS)功率值序列的最大值,Pindex{.}函数表示功率最大值对应的位置索引值;
(708)对步骤(707)的输出进行数值变换,得到当次定时估计值,返回步骤(1)重复进行N次;其中N为设定数值;
数值变换公式为:TA(n)=Pmax-PCP
(709)对步骤(708)的定时估计值进行基于置信概率的统计分析,得到定时调整值。
基于置信概率的统计分析步骤和公式详见下节。
可选的,所述步骤(709)具体包括以下步骤:
(7091)对最近N次的定时估计值依次排列,作为统计序列输入;
记统计序列输入为:{TA(1),TA(2),...TA(N)},其中TA(1)为最近1次的定时估计值,TA(N)为前第N次的定时估计值;
(7092)计数器K清零;
(7093)对输入序列中序号为K的点,统计落在以该点为中心,前后F点范围内的个数,记为集中度C(K);其中,F为设定值;
(7094)计数器K加1,判断K是否大于N-1,如果否则转入(7093)步,如果是则转入步骤(7095);
(7095)集中度C(K)从高到低排列,得到集中度最大值CMAX;
(7096)判断CMAX是否大于置信阈值,如果是,则定时估计有效,否则定时估计无效;
置信阈值与突发长度L和突发子载波个数N有关,记为Th(N)L:其中Th(N)L由事先仿真,程序中查表得出。
(7097)如果定时估计无效,则将统计序列末位值调出,补充最新一次的定时估计值,然后循环(7092)至(7096)步,直至定时估计有效;
记最新一次的定时估计值为TAnew(1),新一次循环的输入序列为{TAnew(1),TAnew(2),...TAnew(N)},则TAnew(n)=TA(n-1),n∈{2,3...N}
如果定时估计有效,对集中度最大值对应的点,计算以该点为中心,前后F点范围内的定时估计值的平均数,作为TA调整值;
TA调整值公式为:
其中TACMAX(n)为以CMAX对应的点为中心,前后F点范围内的定时估计值;
参照图5,所述步骤(709)具体包括以下步骤:
(7098)对N次的定时估计值依次排列,作为统计序列输入;
(7099)计数器K清零;
(7100)对输入序列中序号为K的点,统计落在以该点为中心,前后F点范围内的个数,记为集中度C(K);其中,F为设定值;
(7101)计数器K加1,判断K是否大于N-1,如果否,则转入步骤(7093),如果是,则转入步骤(7095);
(7102)将集中度C(K)从高到低排列,得到集中度最大值CMAX;
(7103)判断CMAX是否大于置信阈值,如果是,则定时估计有效,否则定时估计无效;
(7104)如果定时估计无效,则将统计序列末位值调出,补充最新一次的定时估计值,然后循环步骤(7092)至(7096),直至定时估计有效;如果定时估计有效,对集中度最大值对应的点,计算以该点为中心,前后F点范围内的定时估计值的平均数,作为定时调整值。
所述步骤(9)具体包括以下步骤:
(901)利用全部N个子载波的独特字频域符号,进行噪声功率估计和信号功率估计,得到信噪比估计值;
σs为信号功率均方差,
NUW为一个子载波上的UW字符号数;
(902)根据信噪比估计值,按照信噪比和平衡长度分段关系,得到平衡长度L;
处理公式为L=β×rSN,其中β为平衡系数,由仿真得到,在程序中查表得出。
(903)对全部N个子载波按位置排列并编号,编号从0到N-1(当只有一个子载波时,编号为0),同时子载波计数器NSC清零;
(904)对编号NSC的子载波的去调制和频偏补偿后的独特字符号依次排列,进行平衡长度L的滑动平均滤波处理,得到独特字位置的信道估计值;
处理公式为:
(905)根据独特字位置的信道估计值和位置特征,进行内插和平滑处理,得到编号NSC的子载波的数据字位置的信道估计值;
例如满足假设条件一时,处理公式为:Hdata(n)=γ×HUW(i)+(1-γ)×HUW(i+1)
(906)根据步骤(5)得到的突发信号频率,对编号NSC的子载波的接收数据字符号进行去频差处理;
例如满足假设条件一时,处理公式为: 其中d(n)NSC为接收数据字符号,n=6×i+k, i∈{1,2,...N},N为帧数,k∈{1,2,...6},Ts为突发信号的符号周期,dcomp(n)NSC为去频差的数据字符号;
(907)根据步骤(905)的信道估计值,对去频差的数据字符号进行信道均衡处理,得到编号NSC的子载波的解调数据;
例如LS均衡的处理公式为:ddemod(n)NSC=dcomp(n)NSC×Hdata *(n)
其中Hdata *(n)为Hdata(n)的复共轭,ddemod(n)NSC为编号NSC的子载波的解调数据,也可采用MMSE均衡等方式;
(908)子载波计数器NSC加1,判断NSC是否大于N-1,如果否则转入步骤(904),如果是则完成解调。
Claims (4)
1.一种卫星通信系统低速OFDM突发信号解调方法,其特征在于,包括如下步骤:
(1)对接收的突发信号进行解映射和子载波提取,得到各子载波的频域符号数据;
(2)对各子载波的频域符号数据进行独特字符号提取处理;
(3)对各子载波的独特字频域符号进行去调制运算;
(4)对各子载波的去除调制的独特字频域符号,分别进行基于单子载波频域符号的频率估计;
(5)对各单子载波频域符号的频率估计结果,进行多子载波和多帧的二维统计平滑,得到准确的突发信号频率;
(6)利用准确的突发信号频率,对突发信号进行去频差运算;
(7)对去除调制和频差的独特字频域符号,进行基于统计和置信概率的极少个数子载波定时估计,得到定时调整值;
(8)判断定时调整值是否超出设定的阈值范围,若超出则进行定时调整,当本端为终端侧时,利用定时调整值调整接收突发信号的定时误差,本端为基站侧时,将定时调整值反馈给终端,终端调整发送突发信号的定时误差,本端对当次突发的解调继续转入步骤(9),并返回步骤(1)接收调整后的突发信号;若没有超出则直接转入步骤(9);
(9)利用去除调制和频差的独特字频域符号,进行基于单子载波的信噪比自适应的信道估计和均衡,得到数据字的解调结果;
完成卫星通信系统低速OFDM突发信号解调。
2.根据权利要求1所述的一种卫星通信系统低速OFDM突发信号解调方法,其特征在于,所述步骤(4)中基于单子载波频域符号的频率估计,具体包括以下步骤:
(401)对全部M个子载波按位置排列并编号,编号从0到M-1,同时子载波计数器NSC清零;其中,当只有一个子载波时,编号为0;
(402)判断去调制后独特字频域符号是否等间隔排列,如果是则转入步骤(404),否则转入步骤(403);
(403)根据独特字的位置特征,内插虚拟独特字符号,使内插后的独特字符号和虚拟独特字符号等间隔排列;
(404)将等间隔排列的独特字符号和虚拟独特字符号序列补零,使序列长度等于2^NFFT,其中NFFT为阶数;
(405)对步骤(404)的输出序列,进行快速傅里叶变换,得到长度为2^NFFT的输出序列;
(406)对步骤(405)的输出序列,进行功率值寻峰,并记录功率峰值位置;
(407)对功率峰值位置进行数值变换,得到当前子载波的频率估计值;
(408)子载波计数器NSC加1,判断NSC是否大于M-1,如果否,则转入步骤(402),如果是,则结束。
3.根据权利要求1所述的一种卫星通信系统低速OFDM突发信号解调方法,其特征在于,所述步骤(7)具体包括以下步骤:
(701)对去调制和去频差的独特字频域符号,每个子载波分别进行时域上的累加,得到各子载波上的累加值;
(702)将计数器POS值清零;
(703)以计数器POS值作为定时估计值,根据各个子载波的位置,计算定时偏移在各子载波上产生的相位;
(704)对步骤(701)中各子载波上的累加值,根据步骤(703)中得到的相位,去除累加值中的定时偏移相位;
(705)对各个子载波的去除定时偏移相位的累加值相加,得到的结果记为SUM,记录计数器POS值和对应的SUM值;
(706)将计数器POS值加1,并判断POS值是否大于定时偏移最大值PMAX,如果否,则转入步骤(703),如果是,则转入步骤(707);
(707)对所有POS值对应的SUM值进行功率值寻峰,并记录峰值对应的POS值;
(708)对步骤(707)的输出值进行数值变换,得到当次定时估计值,返回步骤(1)重复进行N次;其中N为设定数值;
(709)对步骤(708)得到的N次定时估计值进行基于置信概率的统计分析,得到定时调整值。
4.根据权利要求3所述的一种卫星通信系统低速OFDM突发信号解调方法,其特征在于,所述步骤(709)具体包括以下步骤:
(7091)对N次的定时估计值依次排列,作为统计序列输入;
(7092)计数器K清零;
(7093)对输入序列中序号为K的点,统计落在以该点为中心,前后F点范围内的个数,记为集中度C(K);其中,F为设定值;
(7094)计数器K加1,判断K是否大于N-1,如果否,则转入步骤(7093),如果是,则转入步骤(7095);
(7095)将集中度C(K)从高到低排列,得到集中度最大值CMAX;
(7096)判断CMAX是否大于置信阈值,如果是,则定时估计有效,否则定时估计无效;
(7097)如果定时估计无效,则将统计序列末位值调出,补充最新一次的定时估计值,然后循环步骤(7092)至(7096),直至定时估计有效;如果定时估计有效,对集中度最大值对应的点,计算以该点为中心,前后F点范围内的定时估计值的平均数,作为定时调整值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111136463.7A CN113904713A (zh) | 2021-09-27 | 2021-09-27 | 一种卫星通信系统低速ofdm突发信号解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111136463.7A CN113904713A (zh) | 2021-09-27 | 2021-09-27 | 一种卫星通信系统低速ofdm突发信号解调方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113904713A true CN113904713A (zh) | 2022-01-07 |
Family
ID=79029842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111136463.7A Pending CN113904713A (zh) | 2021-09-27 | 2021-09-27 | 一种卫星通信系统低速ofdm突发信号解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113904713A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115379550A (zh) * | 2022-07-29 | 2022-11-22 | 西安空间无线电技术研究所 | 一种基于离散导频的突发帧同步方法、装置及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018054236A1 (zh) * | 2016-09-26 | 2018-03-29 | 珠海全志科技股份有限公司 | 一种ofdm系统同步跟踪方法和装置 |
CN110932770A (zh) * | 2019-11-13 | 2020-03-27 | 北京邮电大学 | 一种低轨卫星星地链路同步序列设计及频偏估计方法 |
CN112803992A (zh) * | 2021-04-08 | 2021-05-14 | 成都星联芯通科技有限公司 | 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法 |
-
2021
- 2021-09-27 CN CN202111136463.7A patent/CN113904713A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018054236A1 (zh) * | 2016-09-26 | 2018-03-29 | 珠海全志科技股份有限公司 | 一种ofdm系统同步跟踪方法和装置 |
CN110932770A (zh) * | 2019-11-13 | 2020-03-27 | 北京邮电大学 | 一种低轨卫星星地链路同步序列设计及频偏估计方法 |
CN112803992A (zh) * | 2021-04-08 | 2021-05-14 | 成都星联芯通科技有限公司 | 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法 |
Non-Patent Citations (1)
Title |
---|
张庆业 等: "卫星通信中一种突发信号的解调算法", 数字技术与应用, no. 05, 15 May 2016 (2016-05-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115379550A (zh) * | 2022-07-29 | 2022-11-22 | 西安空间无线电技术研究所 | 一种基于离散导频的突发帧同步方法、装置及设备 |
CN115379550B (zh) * | 2022-07-29 | 2024-05-31 | 西安空间无线电技术研究所 | 一种基于离散导频的突发帧同步方法、装置及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4900447B2 (ja) | 送信装置及び送信方法 | |
CN1777166B (zh) | 在正交频分复用的移动通信系统中获取同步的装置及方法 | |
CN110493156B (zh) | 5g移动通信系统中基于星座点分集的频偏估计方法 | |
KR100635534B1 (ko) | 고속 이동 환경을 위한 하이브리드 채널 추정 방법 및시스템 | |
CN102812679B (zh) | 用于无线通讯系统中准确时间同步的方法及装置 | |
CN102185822B (zh) | 一种ofdm/oqam系统及其时频同步方法 | |
CN110868369B (zh) | 基于5g nr系统的上行信道估计方法及装置 | |
US8416868B2 (en) | Method and system for diversity and mask matching in channel estimation in OFDM communication networks using circular convolution | |
CN102215205B (zh) | 一种ofdm/oqam系统及其时频同步方法 | |
CN111327551B (zh) | 数据与导频频域复用的超奈奎斯特传输方法及传输装置 | |
CN109861939B (zh) | 一种oqpsk频域均衡无线数据传输方法 | |
CN102215204A (zh) | 基于反馈迭代的ofdm/oqam系统及其时频同步方法 | |
US20070133393A1 (en) | Multi-carrier receiving method and multi-carrier receiving apparatus | |
CN110581813A (zh) | 一种多载波系统导频信号的传输方法 | |
CN113973031B (zh) | 一种ofdm系统的信道均衡方法 | |
CN113904713A (zh) | 一种卫星通信系统低速ofdm突发信号解调方法 | |
CN109412987A (zh) | 一种ofdm系统信道跟踪方法 | |
Li et al. | Subspace based blind channel estimation for OFDM by exploiting virtual carrier | |
KR20060001646A (ko) | 직교주파수분할다중 통신 시스템의 채널 추정 방법 및 장치 | |
Chen et al. | Pilot power allocation for OFDM systems | |
CN107204953B (zh) | 一种cp-fbmc通信系统中的盲频偏估计方法 | |
CN102281224B (zh) | Ofdm接收机及其实现方法 | |
CN1937602A (zh) | 一种拓宽多载波通信相干带宽的方法及其装置 | |
Ge et al. | Efficient method for carrier offset correction in OFDM system | |
CN109861938A (zh) | 高数据传输效率的抑制信号中相位旋转误差的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |