CN112803992B - 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法 - Google Patents

一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法 Download PDF

Info

Publication number
CN112803992B
CN112803992B CN202110376239.9A CN202110376239A CN112803992B CN 112803992 B CN112803992 B CN 112803992B CN 202110376239 A CN202110376239 A CN 202110376239A CN 112803992 B CN112803992 B CN 112803992B
Authority
CN
China
Prior art keywords
frequency offset
user
frequency
signal
pilot
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.)
Active
Application number
CN202110376239.9A
Other languages
English (en)
Other versions
CN112803992A (zh
Inventor
钱斌
刘波
干永刚
荣元政
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Xinglian Xintong Technology Co ltd
Original Assignee
Chengdu Xinglian Xintong Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Xinglian Xintong Technology Co ltd filed Critical Chengdu Xinglian Xintong Technology Co ltd
Priority to CN202110376239.9A priority Critical patent/CN112803992B/zh
Publication of CN112803992A publication Critical patent/CN112803992A/zh
Application granted granted Critical
Publication of CN112803992B publication Critical patent/CN112803992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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

本发明涉及一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，该方法使用单一导频即可估计出频偏值与时延值，能够及时对用户信号进行校正，而不必等待多个导频；同时，该方法可以一步地同时估计出时延与频偏值，一次计算即可估计出两个量，比两步估计的方法更为简单；再次，该方法可以估计任意地整数倍的频偏，没有频偏估计范围的限制，不受导频在时域上分配间隔的限制。最后，该方法在整数倍频偏估计后进行定时，有效避免了传统ZC序列定时受整数倍频偏影响的问题。

Description

一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法
技术领域
本发明涉及卫星通信领域，具体涉及一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法。
背景技术
OFDM是一种宽带通信中普遍采用的技术，OFDM利用大量的子载波来承载数据，将不同的子载波分配给不同用户，用户间的子载波保持正交，以实现多址接入。
当前，低轨宽带卫星系统多采用OFDM技术实现多用户接入，多个用户被分配到不同的子载波集上，各自信号单独OFDM调制，经过卫星信道后发送给基站，基站侧接收到的是多个用户信号的时频叠加。由于低轨卫星相对于地面的移动速度较高，用户信号到达基站存在较大的多普勒频偏，多用户各自不同的频偏将导致接收端子载波的正交性遭到破坏，产生载波间干扰（ICI），从而严重影响系统性能。因此基站需要能够分别估计多用户的频偏，并进行补偿，如果频偏比较严重，更需要通知用户调整。
低轨宽带卫星OFDM系统与地面OFDM系统的另外一个不区别是：卫星系统中的信号延迟较大，而上行多用户OFDM信号能够成功解调的前提是多个用户的信号到达基站的时间延迟均不能超过循环前缀（CP）长度，否则超出CP长度的信号会造成符号间干扰（ISI）。为避免ISI，CP的长度要保证大于多径的时延，而传播导致的时延不会被考虑在CP范围之内，否则会导致CP长度过长，物理资源浪费严重。但是，低轨卫星系统的小区覆盖范围大，不同用户到卫星的时延差异也较大，为了避免这种情况导致的ISI，一种比较正常的处理方式是在基站侧对用户时延进行估计，然后通过时间调整命令告知用户调整发送时间，从而保证多个用户信号的到达时间保持在一个较小的误差范围之内，从而提升系统效率。所以对低轨卫星系统，为避免ISI，基站对多用户时延的准确估计十分重要。综上所述，不同与地面OFDM通信系统中，频率同步和时延同步的问题在低轨宽带OFDM系统中更加重要。
OFDM系统中的频偏估计一般有三种方法：1）导频方法是利用两个已知导频间的相位差来估计频偏，这种方法的前提是需要两个导频，而且该方法能够估计的最大频偏取决与两个导频的间隔，间隔越大可估范围越小，一般最大可估范围是在两个导频连续排列时，此时最大估计范围为子载波间隔的正负0.5倍。但在实际OFDM系统中导频很少这样排列，一般情况是间隔几个符号插入一个导频，因此实际可估的频偏范围要小很多，在低轨卫星系统中的应用受到一定限制；2）基于CP的方法是利用CP的重复数据，使用CP数据和符号后部的某些数据的相位差来估计频偏，这种方法同样有估计范围限制，最大可估频偏为子载波间隔的正负0.5倍，同时，由于多径等影响，可利用的CP数据较少，也会导致低信噪比下影响估计精度；3）还有一类盲估计的方法，比如基于ESPRIT旋转不变性质的方法，这类盲方法的原理是利用信号的统计特征，所以需要相关矩阵的精确计算，这是建立在有大量的样本数据的基础之上，因此，在少量观测的情况下效果有限，同时后续处理中一般要进行SVD分解和求伪逆等复杂运算，在工程实现上比较困难。
因此，如何保证低轨宽带OFDM系统中的频率同步实现精准的频偏估计，是现阶段继续解决的问题。
发明内容
本发明的目的在于克服现有技术的不足，提供一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，该方法使用单一导频即可估计出频偏值与时延值，能够及时对用户信号进行校正，而不必等待多个导频；再次，该方法可以估计多倍子载波间隔的频偏，没有频偏估计范围的限制，不受导频在时域上分配间隔的限制。
本发明的目的是通过以下技术方案来实现的：
一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，包括：
1）定义发射信号模型：
假设小区内共有M个用户，定义用户m的导频参考信号为
Figure 430214DEST_PATH_IMAGE001
，则导频基带发送信号
Figure 843877DEST_PATH_IMAGE002
可以表示为：
Figure 632580DEST_PATH_IMAGE003
……（1）
其中，m是用户索引，p代表是导频信号，上角标T代表发射，
Figure 247232DEST_PATH_IMAGE004
为序列长度，
Figure 357270DEST_PATH_IMAGE005
为根值为
Figure 449991DEST_PATH_IMAGE006
的ZC序列，
Figure 719298DEST_PATH_IMAGE007
；
2）定义接收信号模型：
在基站，接收到的基带信号
Figure 3387DEST_PATH_IMAGE008
，是M个用户信号的叠加：
Figure 600721DEST_PATH_IMAGE009
（2）
其中，
Figure 762712DEST_PATH_IMAGE010
为噪声，
Figure 355368DEST_PATH_IMAGE011
为接收到的第m个用户的导频信号，上角标R代表接收，导频信号接收模型
Figure 577402DEST_PATH_IMAGE011
为：
Figure 178146DEST_PATH_IMAGE012
其中，
Figure 878248DEST_PATH_IMAGE013
为第m用户的离散信道冲击响应（CIR）向量，
Figure 731935DEST_PATH_IMAGE014
为用户m的多径数量，
Figure 718346DEST_PATH_IMAGE015
为频率偏差，
Figure 555852DEST_PATH_IMAGE016
为采样周期整数倍的时间偏移，
Figure 558180DEST_PATH_IMAGE017
为采样时间。
Figure 797532DEST_PATH_IMAGE018
是信道冲击响应
Figure 95789DEST_PATH_IMAGE019
的DFT，
Figure 155012DEST_PATH_IMAGE020
；
3）频偏估计：
定义用户m的频偏超出[-0.5，+0.5]倍子载波间隔的范围内，将用户频偏划分为整数倍频偏
Figure 961032DEST_PATH_IMAGE021
和小数倍频偏
Figure 789310DEST_PATH_IMAGE022
两部分，则接收到的用户m导频信号
Figure 258469DEST_PATH_IMAGE011
为：
Figure 398463DEST_PATH_IMAGE023
对接收时域信号做傅里叶变换到频域得到
Figure 509639DEST_PATH_IMAGE024
后，与本地导频
Figure 690959DEST_PATH_IMAGE025
做相关，
Figure 331019DEST_PATH_IMAGE026
搜索相关峰
Figure 99255DEST_PATH_IMAGE027
的最大值的位置，即可得到整数倍频偏
Figure 342017DEST_PATH_IMAGE021
的估计值；
得到整数倍频偏
Figure 879309DEST_PATH_IMAGE021
估计值后，找到峰值位置，并在峰值位置上取相应的N个频域值，与本地参考频域导频的共轭对应一一相乘，然后对共轭相乘后的数据进行傅里叶逆变换，搜索峰值，利用主峰值与次峰值进行小数倍频偏
Figure 194665DEST_PATH_IMAGE022
估计值，小数倍频偏
Figure 574831DEST_PATH_IMAGE022
估计值为：
Figure 496650DEST_PATH_IMAGE028
4）时延估计：
经过整数倍频偏的校正以后，整数倍的频偏不再影响相关后的峰值位置，主峰值的位置既为信号的时间发送位置，由此可以判断信号的时延是否大于CP（循环前缀），并借此为用户TA调整的依据。
进一步的，因为已经经过整数倍频偏校正，所述[-K，+K]中，K的取值为0.5倍子载波间隔。
进一步的，所述2）中第m个用户的导频信号：
Figure 154028DEST_PATH_IMAGE029
在接收信号中分离出信道、频偏与时延参数，则得到所述导频信号接收模型，其中，
Figure 135890DEST_PATH_IMAGE030
为频偏，
Figure 642833DEST_PATH_IMAGE031
为多径索引，
Figure 492977DEST_PATH_IMAGE014
为多径数量，
Figure 4861DEST_PATH_IMAGE032
为信道系数，
Figure 157625DEST_PATH_IMAGE033
为发射信号，
Figure 122170DEST_PATH_IMAGE034
为用户时延，
Figure 776005DEST_PATH_IMAGE035
为第l径的时延。
本发明的有益效果是：本方案以低轨宽带卫星OFDM系统为基础，推导了基于单导频的时延-频偏联合估计算法。该方法可以对抗大的频偏和时延，并且频偏估计范围不受导频配置方式的限制，更加适用于低轨卫星系统频偏和频偏变化率均较大的情况。同时，频偏和时延的估计只依赖于导频相关值的计算，无需多余步骤，减少了联合估计的计算量。最后，仿真时延验证了算法的有效性。
附图说明
图1是本发明信噪比-10dB下的相关峰仿真图；
图2是采用CP数据相关的相关峰仿真图；
图3是小数倍频偏估计的性能仿真均值图；
图4是小数倍频偏估计的性能仿真方差图。
具体实施方式
下面结合具体实施例进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，应用于一种对抗低轨宽带卫星系统大频偏大时延的联合估计方法中的整数倍子载波间隔频偏估算，其中，一种对抗低轨宽带卫星系统大频偏大时延的联合估计方法包括如下：
发射信号模型：
考虑低轨卫星系统中上行采用DFT-S-OFDM情况，在一个小区内共有M个用户，用户数据调制后做DFT（离散傅里叶变换），并承载到分配给用户的子载波上，然后按照规定格式进行组帧，做IDFT（离散傅里叶逆变换）后添加CP形成基带时域信号。其中，根据上层配置，导频和数据在不同的符号上发送，而且，不同于数据符号，导频或SRS信号不做第一步的DFT运算，而是直接映射到子载波上做IDFT运算后生成基带信号。假设用户m的导频参考信号为
Figure 375351DEST_PATH_IMAGE036
，则导频基带发送信号
Figure 699016DEST_PATH_IMAGE037
可以表示为：
Figure 541070DEST_PATH_IMAGE038
其中, m是用户索引，p代表是导频信号，上角标T代表发射，
Figure 608383DEST_PATH_IMAGE004
为序列长度，
Figure 829280DEST_PATH_IMAGE005
为根值为
Figure 87961DEST_PATH_IMAGE006
的ZC序列，定义如下：
Figure 151732DEST_PATH_IMAGE039
接收信号模型：
在基站，接收到的
Figure 22736DEST_PATH_IMAGE008
是M个用户信号的叠加：
Figure 98139DEST_PATH_IMAGE040
其中，
Figure 29186DEST_PATH_IMAGE041
为噪声，
Figure 580253DEST_PATH_IMAGE011
为接收到的第m个用户的导频信号，上角标R代表接收，
Figure 482045DEST_PATH_IMAGE042
可以表示为：
Figure 677534DEST_PATH_IMAGE043
其中，
Figure 638537DEST_PATH_IMAGE044
为第m用户的离散信道冲击响应（CIR）向量，
Figure 552266DEST_PATH_IMAGE045
为多径数量，
Figure 765073DEST_PATH_IMAGE046
为频率偏差，
Figure 844762DEST_PATH_IMAGE047
为采样周期整数倍的时间偏移，
Figure 586453DEST_PATH_IMAGE017
为采样时间。在基站端，对导频符号的利用，第一步是先估计每个用户的
Figure 112112DEST_PATH_IMAGE046
与
Figure 863031DEST_PATH_IMAGE048
，然后，做频偏补偿，再进行信道估计均衡等处理。
一般估计这两个参数需要利用已知的导频信息，本文假设多用户上行采用ZC序列为导频（例如LTE中的DMRS/SRS或NR中的SRS等），然后利用单一导频实现时间偏移与频率偏移的联合估计。
考虑用户m的
Figure 33112DEST_PATH_IMAGE045
个信道扩展在OFDM符号的CP扩展之内，则接收信号
Figure 709819DEST_PATH_IMAGE042
可以表示为：
Figure 722774DEST_PATH_IMAGE049
Figure 277383DEST_PATH_IMAGE050
其中，
Figure 301971DEST_PATH_IMAGE051
是信道冲击响应
Figure 916623DEST_PATH_IMAGE052
的DFT：
Figure 151295DEST_PATH_IMAGE053
传统情况下，将时延
Figure 476972DEST_PATH_IMAGE054
导致的影响合并到信道以内：
Figure 887225DEST_PATH_IMAGE055
然后在均衡中消除
Figure 266254DEST_PATH_IMAGE054
的影响，所以不需要进行时间恢复。这种方法要求UE到基站的传播时延和信道扩展之和在CP的范围以内，如果超出CP范围，则会造成符号间干扰。因此限制了所能忍受的最大传播时延，也既限制了用户与基站间的最大距离。
在低轨卫星通信中，地面设备通过卫星联络，通常具有较大的传播时延，要保证CP长与传播时延，会导致CP长度过长，物理资源浪费严重。一种比较正常的处理方式是在基站侧估计用户时延，然后通过时间调整命令告知用户调整发送时间，从而保证多个用户信号的到达时间保持在一个较小的CP范围之内，从而提升系统效率。
本方案要解决问题是利用多个用户的单一导频信号对每一个用户分别进行频偏和时延的联合估计，从而解决基站侧多用户的时间频率同步问题。在接收信号中分离出信道、频偏与时延参数如下
Figure 863588DEST_PATH_IMAGE056
频偏估计：
整数倍子载波间隔频偏估计：
假设用户m的频偏超出[-0.5，+0.5]倍子载波间隔的范围，可将用户频偏划分为整数倍频偏
Figure 760000DEST_PATH_IMAGE057
和小数倍频偏
Figure 618235DEST_PATH_IMAGE058
两部分，则接收到的用户m导频信号为：
Figure 344663DEST_PATH_IMAGE059
在接收端对时域信号
Figure 429294DEST_PATH_IMAGE042
进行傅里叶变换运算，则子载波
Figure 863817DEST_PATH_IMAGE060
上的接收信号为
Figure 717504DEST_PATH_IMAGE061
：
Figure 969494DEST_PATH_IMAGE062
Figure 39956DEST_PATH_IMAGE063
其中，
Figure 278170DEST_PATH_IMAGE064
，上式中第一项为主分量，后一项是因为频偏导致载波间不正交带来的干扰。经过简化后
Figure 845418DEST_PATH_IMAGE065
可以表示如下：
Figure 143675DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE067
Figure 639116DEST_PATH_IMAGE068
其中，第一项为距离频偏
Figure 415442DEST_PATH_IMAGE046
最近的子载波
Figure 368355DEST_PATH_IMAGE069
上的值，后一项是因为频偏导致载波间不正交带来的干扰。
对接收信号做傅里叶变换得到
Figure 837513DEST_PATH_IMAGE070
后，与本地导频序列做相关运算；
Figure 384032DEST_PATH_IMAGE071
忽略后一项频偏导致的载波间不正交干扰；
Figure 728164DEST_PATH_IMAGE072
Figure 270003DEST_PATH_IMAGE073
由于导频序列的自相关特性，在
Figure 175643DEST_PATH_IMAGE074
时，
Figure 943878DEST_PATH_IMAGE075
有最大值。因此对接收时域信号做傅里叶变换到频域后，与本地导频做相关，计算相关
Figure 62007DEST_PATH_IMAGE075
， 搜索大最值的位置，即可得到整数倍频偏的估计值。
小数倍子载波间隔频偏估计：
接收时域信号的傅里叶变换为：
Figure 693491DEST_PATH_IMAGE076
将
Figure 629086DEST_PATH_IMAGE065
乘以本地频域参考信号的共轭
Figure DEST_PATH_IMAGE077
；
然后做傅里叶逆变换
Figure 556722DEST_PATH_IMAGE078
；
可以证明上式等效为时域相关：
Figure 478542DEST_PATH_IMAGE079
所以接收处理流程为：在接收端对接
Figure 260553DEST_PATH_IMAGE080
进行FFT后，搜索频域序列的最大相关值，定出频域整数倍偏移，然后在相应位置取出频域接收序列与本地频域参考序列的共轭相乘，相乘后做逆傅里叶变换。则相关值
Figure 740951DEST_PATH_IMAGE081
如下，其中
Figure 749358DEST_PATH_IMAGE031
为搜索值索引。
Figure 209289DEST_PATH_IMAGE082
Figure 845807DEST_PATH_IMAGE083
其中，
Figure 998571DEST_PATH_IMAGE084
为小数倍频偏，
Figure 461651DEST_PATH_IMAGE085
为子载波间隔。然后计算
Figure 990852DEST_PATH_IMAGE086
处的
Figure 481877DEST_PATH_IMAGE087
值：
Figure 805542DEST_PATH_IMAGE088
Figure 788541DEST_PATH_IMAGE089
假设：
Figure 354389DEST_PATH_IMAGE090
则：
Figure 434341DEST_PATH_IMAGE091
然后，计算较大副峰处的
Figure 194486DEST_PATH_IMAGE092
值，假设频偏为正，则MOD(
Figure 133624DEST_PATH_IMAGE093
：
Figure 739048DEST_PATH_IMAGE094
Figure 204665DEST_PATH_IMAGE095
则：
Figure 640106DEST_PATH_IMAGE096
假设频偏为正切在0.5倍子载波间隔内，则
Figure 66540DEST_PATH_IMAGE097
，当N较大时，取如下近似：
Figure 741235DEST_PATH_IMAGE098
最终，小数倍频偏的估计为：
Figure 936724DEST_PATH_IMAGE099
时延估计：
经过整数倍频偏的校正以后，整数倍的频偏不再影响相关后的峰值位置，主峰值的位置既为信号的时间发送位置，由此可以判断信号的时延是否大于CP，并借此为用户TA调整的依据。证明如下：
由上节可知，发送信号为
Figure 632147DEST_PATH_IMAGE100
，假设接收信号经过延迟
Figure 44412DEST_PATH_IMAGE101
到达接收机，则接收信号为：
Figure 257218DEST_PATH_IMAGE102
其中，
Figure 431848DEST_PATH_IMAGE103
。则相关值为：
Figure 704697DEST_PATH_IMAGE104
当
Figure 840144DEST_PATH_IMAGE105
时，相关值最大。所以，搜索到的相关峰最大峰值的位置既为信号的实际发送位置。
通过上面分析可知，本方法在频偏与时延估计时，均依赖于相关值
Figure 715696DEST_PATH_IMAGE106
，无需再进行多余的计算，优于两步的联合估计方法，因为两步方法的时延估计也是依赖于相关值，而频偏估计需要额外的运算。
仿真结果：
考虑低轨卫星系统采用ka波段，为避免相位噪声影响，采用较大的子载波间隔，在本仿真中，子载波间隔为常用的120kHz。低轨系统中的频偏相对较大，用户的频偏校正分为多个步骤进行，首先，UE根据星历估算卫星轨道能够进行一部分多普勒频偏的补偿，其次，在UE接入网络的时候，使用PRACH信道，基站侧可以进一步估算UE的频偏并进行补偿，本发明的频偏补偿方案是用在用户接入以后，随着时间、位置和相对运动的变化，用户频偏产生变化，基站在与用户通信的时候利用用户的导频进行频偏估计与调整。仿真中设定用户的频偏为子载波间隔的1.5倍，时延为20个采样点偏差。给用户分配的子载波个数为839个，用户导频采用根植为300的ZC序列，长度为839。发送中，UE生成839的ZC序列，然后映射到分配的子载波上，然后IDFT生成基带信号。在基站接收端，对小区带宽进行一次FFT计算，然后取出用户分配的子载波上的数据进行整数倍频偏估计、小数倍频偏估计和时延估计。
整数倍频偏的估计，如前述，在接收端用本地参考ZC序列与接收机FFT后的数据做相关值，值最大对应的偏移，即为整数倍的频偏。去掉整数倍的频偏后，用对应的数据与本地参考ZC序列做共轭相乘后IDFT得到相关序列，在相关序列内搜索最大峰值点，做时延估计。相关峰如图1所示。
图1中是信噪比-10dB下的相关峰，图1中可见，相关峰在索引值21的位置上，对应发射信号的20个采样点的偏移。仿真实验证明，在大于-10dB的信噪比下，时延估计均能得到正确值。与之对比，图2是采用CP数据相关的相关峰，其中符号长度为2048，CP长度为144，信噪比为0Db，在-10dB情况下已经无法得到明显的峰值。由图1与图2对比可见，本方案的方法具有更准确的定位精度和可以对抗更低的信噪比。
图3和图4为小数倍频偏估计的性能仿真图，根据系统参数设置，去除整数倍频偏后，小数倍频偏的残留值为36000Hz。图3为频偏估计的均值，由图3可见，随着信噪比的增大，估计性能越来越好，能够正确估计出0.3倍的频偏。图4为估计方差图，由图4可知，随着信噪比增加，估计方差逐渐变小。与本发明的方法对比，利用2个导频的方法估计频偏的最大值受到导频间隔的影响，比如两个导频符号之间插入两个数据符号的配置，其频偏估计的最大范围为1/6的子载波间隔，此时，已经无法估计出0.3倍的频偏。而导频的密度不可能再大，太大的导频密度，将导致系统效率变低。所以本发明方法再对抗大频偏和大时延方面具有很好的优势，更适用于低轨卫星的应用场景。
本方案以低轨宽带卫星OFDM系统为基础，推导了基于单导频的时延-频偏联合估计算法。该方法可以对抗大的频偏和时延，并且频偏估计范围不受导频配置方式的限制，更加适用于低轨卫星系统频偏和频偏变化率均较大的情况。同时，频偏和时延的估计只依赖于导频相关值的计算，无需多余步骤，减少了联合估计的计算量。最后，仿真时延验证了算法的有效性。
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (3)

1.一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，其特征在于，该方法包括：
频偏估算步骤：判断用户频偏是否满足预设的子载波间隔条件，将用户频偏划分为整数倍频偏和小数倍频偏，对接收时域信号做傅里叶变换到频域后，与本地导频做相关运算，搜索相关峰最大值的位置，得到整数倍频偏的估计值；
得到整数倍频偏估计值后，找到峰值位置，在峰值位置上取N个频域值与本地参考频域导频的共轭 一一相乘，然后对共轭相乘后的数据进行傅里叶逆变换，搜索峰值，利用主峰值与次峰值得到小数倍频偏的估计值；
该方法还包括构建发射信号模型步骤和构建接收信号模型步骤；在执行所述频偏估算步骤之前需要依次完成所述构建发射信号模型步骤和构建接收信号模型步骤；
所述构建发射信号模型步骤包括：
B1、将用户m的导频参考信号
Figure DEST_PATH_IMAGE002
调制后做DFT，并承载到分配给用户的子载波上；
B2、按照规定格式进行组帧，在进行IDFT后添加CP得到导频基带发送信号
Figure DEST_PATH_IMAGE004
，其中, m是用户索引，p代表是导频信号，j表示虚数，上角标T代表发射，N为序列长度，
Figure DEST_PATH_IMAGE006
为用户索引m的序列长度，k为子载波索引，n为时域索引值；
所述构建接收信号模型步骤包括：
C1、从接收到的M个用户信号叠加的信号中分离出第m个用户的导频信号
Figure DEST_PATH_IMAGE008
；
其中，
Figure DEST_PATH_IMAGE010
为第m个用户的多径数量，
Figure DEST_PATH_IMAGE012
为发射信号，l为多条路径，N为序列长度，
Figure DEST_PATH_IMAGE014
为第m用户的离散信道冲击响应向量，
Figure DEST_PATH_IMAGE016
为频率偏差，
Figure DEST_PATH_IMAGE018
为采样周期整数倍的时间偏移，int表示整数数据类型的数据，即取括号里面公式的整数倍，τm表示时延，
Figure DEST_PATH_IMAGE020
为采样时间，n为时域索引值；
C2、当信道扩展满足
Figure DEST_PATH_IMAGE022
在OFDM符号的CP扩展之内条件时，将导频信号转换为
Figure DEST_PATH_IMAGE024
，
Figure DEST_PATH_IMAGE026
表示用户m的导频参考信号；
其中，
Figure DEST_PATH_IMAGE028
表示信道冲击响应
Figure DEST_PATH_IMAGE030
的DFT，
Figure DEST_PATH_IMAGE032
表示频偏的影响，进而从中将信道
Figure 851635DEST_PATH_IMAGE028
与频偏
Figure 34355DEST_PATH_IMAGE032
进行了分离。
2.根据权利要求1所述的一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，其特征在于，设置子载波间隔范围为[-K，+K]，当用户频偏超出子载波间隔范围[-K，+K]时，则表示用户频偏满足预设的子载波间隔条件，否则表示不满足。
3.根据权利要求2所述的一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法，其特征在于，所述子载波间隔范围[-K，+K]中，K的取值为大于等于0.5倍的子载波间隔。
CN202110376239.9A 2021-04-08 2021-04-08 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法 Active CN112803992B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110376239.9A CN112803992B (zh) 2021-04-08 2021-04-08 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110376239.9A CN112803992B (zh) 2021-04-08 2021-04-08 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法

Publications (2)

Publication Number Publication Date
CN112803992A CN112803992A (zh) 2021-05-14
CN112803992B true CN112803992B (zh) 2021-07-13

Family

ID=75816571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110376239.9A Active CN112803992B (zh) 2021-04-08 2021-04-08 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法

Country Status (1)

Country Link
CN (1) CN112803992B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112887241B (zh) 2021-04-29 2021-07-16 成都星联芯通科技有限公司 频偏估计方法、装置、通信设备及存储介质
CN113904713A (zh) * 2021-09-27 2022-01-07 中国电子科技集团公司第五十四研究所 一种卫星通信系统低速ofdm突发信号解调方法
CN114389930B (zh) * 2021-12-29 2024-03-15 中电科思仪科技(安徽)有限公司 一种终端上行信号基带解调同步与频偏补偿装置及方法
CN117014278B (zh) * 2023-09-25 2023-12-22 为准（北京）电子科技有限公司 整数倍频偏估计方法、装置及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1480401A2 (en) * 2003-05-20 2004-11-24 Samsung Electronics Co., Ltd. Frequency-offset estimator for multi-carrier reception
KR20070056881A (ko) * 2005-11-30 2007-06-04 삼성전자주식회사 직교 주파수 다중 접속 시스템에서의 주파수 복원 장치 및방법
CN101009513A (zh) * 2006-01-26 2007-08-01 上海原动力通信科技有限公司 宽带时分双工蜂窝系统的小区同步方法及小区初搜方法
CN101106557A (zh) * 2006-07-14 2008-01-16 华为技术有限公司 一种组帧同步方法
CN101437005A (zh) * 2008-12-12 2009-05-20 中国人民解放军理工大学通信工程学院 通信同步过程中存在定时误差的整数倍频偏估计方法
CN101640650A (zh) * 2008-07-31 2010-02-03 三星电子株式会社 整数倍频偏和符号定时的联合估计方法和设备
CN101667990A (zh) * 2009-07-23 2010-03-10 中国科学院微电子研究所 Ofdm频率偏移联合估计方法
KR20100136867A (ko) * 2009-06-19 2010-12-29 유정시스템(주) 와이브로 시스템에서 정수배 주파수 오차 추정 장치 및 그 방법
CN110545247A (zh) * 2019-07-31 2019-12-06 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法
CN112584538A (zh) * 2019-09-30 2021-03-30 华为技术有限公司 卫星通信方法和相关通信设备

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104320363B (zh) * 2014-10-22 2017-09-05 西安电子科技大学 单载波频域均衡系统时频二维联合同步方法
US10257776B2 (en) * 2016-09-14 2019-04-09 Indian Institute Of Technology, Hyderabad Method of performing cell search by user equipment (UE)
CN110932770B (zh) * 2019-11-13 2020-10-16 北京邮电大学 一种低轨卫星星地链路同步序列设计及频偏估计方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1480401A2 (en) * 2003-05-20 2004-11-24 Samsung Electronics Co., Ltd. Frequency-offset estimator for multi-carrier reception
KR20070056881A (ko) * 2005-11-30 2007-06-04 삼성전자주식회사 직교 주파수 다중 접속 시스템에서의 주파수 복원 장치 및방법
CN101009513A (zh) * 2006-01-26 2007-08-01 上海原动力通信科技有限公司 宽带时分双工蜂窝系统的小区同步方法及小区初搜方法
CN101106557A (zh) * 2006-07-14 2008-01-16 华为技术有限公司 一种组帧同步方法
CN101640650A (zh) * 2008-07-31 2010-02-03 三星电子株式会社 整数倍频偏和符号定时的联合估计方法和设备
CN101437005A (zh) * 2008-12-12 2009-05-20 中国人民解放军理工大学通信工程学院 通信同步过程中存在定时误差的整数倍频偏估计方法
KR20100136867A (ko) * 2009-06-19 2010-12-29 유정시스템(주) 와이브로 시스템에서 정수배 주파수 오차 추정 장치 및 그 방법
CN101667990A (zh) * 2009-07-23 2010-03-10 中国科学院微电子研究所 Ofdm频率偏移联合估计方法
CN110545247A (zh) * 2019-07-31 2019-12-06 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路载波频偏估计与补偿方法
CN112584538A (zh) * 2019-09-30 2021-03-30 华为技术有限公司 卫星通信方法和相关通信设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单/多载波通信系统中载波同步技术的研究与实现;任文龙;《中国优秀硕士学位论文全文数据库 信息科技辑》;20170315 *

Also Published As

Publication number Publication date
CN112803992A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
CN112803992B (zh) 一种低轨宽带卫星系统整数倍子载波间隔频偏估算方法
US7907593B2 (en) Staggered pilot transmission for channel estimation and time tracking
US7526020B2 (en) Circuit for improving channel impulse response estimation and compensating for remnant frequency offset in the orthogonal frequency division multiplexing (OFDM) baseband receiver for IEEE 802.11a/g wireless LAN standard
US8553822B2 (en) Time filtering for excess delay mitigation in OFDM systems
US8761274B2 (en) Least squares channel identification for OFDM systems
US7554902B2 (en) Preamble design for frequency offset estimation and channel equalization in burst OFDM transmission system
JP5497800B2 (ja) 直交周波数分割多重無線通信システムにおけるキャリア周波数オフセットを補償するための方法及び装置
US7844006B2 (en) Method of non-uniform doppler compensation for wideband orthogonal frequency division multiplexed signals
JP4164364B2 (ja) 複雑さが低減したチャンネル応答推定を有するマルチキャリヤ伝送システム
US20040005010A1 (en) Channel estimator and equalizer for OFDM systems
US20070086328A1 (en) Method and circuit for frequency offset estimation in frequency domain in the orthogonal frequency division multiplexing baseband receiver for IEEE 802.11A/G wireless LAN standard
US9762414B2 (en) Least squares channel identification for OFDM Systems
US8223862B2 (en) OFDM inter-carrier interference cancellation method
US10644904B2 (en) System and method for channel estimation
US20080137723A1 (en) Equalization method and apparatus for time domain synchronous orthogonal frequency division multiplexing receiver
KR101485785B1 (ko) 무선 통신 시스템에서 주파수 추정 방법 및 장치
CN113098819B (zh) 一种低轨宽带卫星系统小数倍子载波间隔频偏估算方法
US8199845B2 (en) Up-link SDMA receiver for WiMAX
CN113179233B (zh) 一种对抗低轨宽带卫星系统大频偏大时延的联合估计方法
Takahashi et al. A novel symbol synchronization algorithm with reduced influence of ISI for OFDM systems
CN110943949B (zh) 一种基于sc-fde传输体制的宽带角跟踪方法
Wang et al. A low-complexity joint time synchronization and channel estimation scheme for orthogonal frequency division multiplexing systems
Gunther et al. A new approach for symbol frame synchronization and carrier frequency estimation in OFDM communications
Witschnig et al. The advantages of a unique word for synchronisation and channel estimation in a SC/FDE system
US9083596B1 (en) Non-data-aided joint time and frequency offset estimation method for OFDM systems using channel order based regression

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
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A method for estimating integer multiple subcarrier interval frequency offset in low orbit broadband satellite systems

Granted publication date: 20210713

Pledgee: Chengdu Rural Commercial Bank Co.,Ltd. Jinquan Branch

Pledgor: Chengdu Xinglian Xintong Technology Co.,Ltd.

Registration number: Y2024980019990