CN108649968A - 一种基带信号的提取方法 - Google Patents
一种基带信号的提取方法 Download PDFInfo
- Publication number
- CN108649968A CN108649968A CN201810449706.4A CN201810449706A CN108649968A CN 108649968 A CN108649968 A CN 108649968A CN 201810449706 A CN201810449706 A CN 201810449706A CN 108649968 A CN108649968 A CN 108649968A
- Authority
- CN
- China
- Prior art keywords
- signal
- baseband signal
- error
- radiofrequency signal
- formula
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000000605 extraction Methods 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims description 15
- 239000000284 extract Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/0003—Software-defined radio [SDR] systems, i.e. systems wherein components typically implemented in hardware, e.g. filters or modulators/demodulators, are implented using software, e.g. by involving an AD or DA conversion stage such that at least part of the signal processing is performed in the digital domain
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
一种基带信号的提取方法,包括以下步骤:步骤1:获取探测器的下行射频信号s(t);所述射频信号s(t)包括至少两个基带信号,基带信号的个数为N;步骤2:对射频信号s(t)进行处理;步骤3:从射频信号s(t)中提取基带信号sb(t);步骤4:对得到的基带信号sb(t)进行误差分析;本发明能够直接把基带信号从射频信号中提取出来,该方法简单、误差小,提高了工作效率。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及一种基带信号的提取方法。
背景技术
在深空探测、大地测量过程中,接收机接收的目标探测器发射的下行信号为已调制的射频信号(S波段载波频率2GHz,X波段载波频率8GHz),地面接收站需要对该信号进行解调、滤波等处理,将射频信号变换到基带,然后进行记录、传输及相关处理。由于接收机接收到探测器下行信号为射频频信号,如果用软件设计一个该发射接收系统,所需要的数字采样率高,产生的数据量大,且其中的大部分为噪声信号,有用(有效)带宽只有几兆到数十兆;如果设计数字滤波器提取基带信号,所需要设计的滤波器采样率很高(S波段至少4GHz,X波段至少16GHz),根据分析滤波器的阶数需要上千阶才能达到设计要求,运算量大,计算缓慢,效率低下。
目前把射频信号提取到基带进行处理的方法是根据载波频率,构建一个余弦函数,对射频信号进行解调,然后滤波,该方法同样有数据量大,设计的滤波器阶数高,运算缓慢,效率低下的缺点。
为解决以上问题,本发明提供了一种基带信号的提取方法,通过对接收机接收的典型探测器下行信号的频域特点进行分析,给出了一种在不运用滤波器的情况下,直接把基带信号从射频信号中提取出来的简单有效的方法;通过该方法,产生的数据量与直接对射频信号进行分析产生的数据量以及采用滤波器提取基带信号产生的数据量相比,至少小了2到3个数量级,运算速度和效率大大提升,而信号的信息并没有丢失;同时通过对对该方法做误差分析,使得可以根据所需要的精度,更加灵活的提取所需的基带信号。
发明内容
针对现有技术中存在的不足,本发明的目的在于提供一种基带信号的提取方法,能够直接把基带信号从射频信号中提取出来,该方法简单、误差小,提高了工作效率。
一种基带信号的提取方法,包括以下步骤:
步骤1:获取探测器的下行射频信号s(t);所述射频信号s(t)包括至少两个基带信号,基带信号的个数为N;
步骤2:对射频信号s(t)进行处理;
步骤3:从射频信号s(t)中提取基带信号sb(t);
步骤4:对得到的基带信号sb(t)进行误差分析;在提取基带信号过程中产生的总误差为RMT,所述总误差包括误差RM和误差RM′两部分。
优选地,所述射频信号s(t)为:
其中,PT为探测器射频信号的总功率,fc为载波频率,d(t)或sin(2πft)为副载波,θ为其调制系数;mTM为射频信号的调制系数,fTM为射频信号的频率;f1,f2,fN分别为基带信号1、基带信号2和基带信号N的频率,m1,m2,mN分别为基带信号1、基带信号2和基带信号N的调制系数。
优选地,所述步骤2具体包括:
步骤2.1:对射频信号s(t)进行化简;
其中,fi为各个基带信号的频率,mi为各个基带信号的调制系数,i的取值范围为1~N;
步骤2.2:根据欧拉公式式(2)中的射频信号s(t)转换为:
由于ω=2πf,对式(3)进行整理,得到射频信号s(t)的泰勒展开式(4):
步骤2.3:取射频信号s(t)的泰勒展开式(4)的正频部分,忽略负频部分,得到下式:
优选地,所述步骤3具体包括:
步骤3.1:假设射频信号s(t)信号中只有两个基带信号,即取N=2,根据式(5)在正频域的表现形式,其载波信号和各个基带信号的频域幅值强度分别为PC,P1,P2;
步骤3.2:根据该射频信号s(t)的频谱特点,构建一个只含有基带信号的表达式sb(t):
其中,式(6)中的aik为:
优选地,所述步骤4具体包括:
步骤4.1:得到所述误差RM;在得到式(5)泰勒展开式时,对k取有限的自然数M,产生误差RM,误差RM为:
步骤4.2:得到所述误差RM′;根据需要的精度,对式(6)中的k可取有限的自然数M,取值时产生误差RM′,误差RM′为:
其中,RMi为各个基带信号展开项的误差,由下式表示:
步骤4.3:在提取基带信号过程中产生的总误差RMT为误差RM与误差RM′的乘积,如下式:
本发明的有益效果在于:
本发明一种基带信号的提取方法,通过对接收机接收的典型探测器下行射频信号的频域特点进行分析,在不运用滤波器的情况下,直接把基带信号从射频信号中提取出来,该方法简单有效;通过该方法,产生的数据量与直接对射频信号进行分析产生的数据量以及采用滤波器提取基带信号产生的数据量相比,至少小了2到3个数量级,运算速度和效率大大提升,而信号的信息并没有丢失;同时通过对该方法得出的基带信号做了误差分析,使得可以根据所需要的精度,更加灵活的提取所需的基带信号。
附图说明
图1是本发明一种基带信号的提取方法的流程图。
图2是本发明的实施例中射频信号的频谱示意图。
图3是本发明的实施例中根据射频信号频谱构建的基带信号的频谱示意图。
图4是本发明的实施例中射频信号的频谱图。
图5是本发明的实施例中两个基带信号的频谱图。
图6是本发明的实施例中射频信号降频之后与基带信号的频谱比较图。
具体实施例
下面结合附图和实施例对本发明做进一步阐述。
实施例1:
如图1所示,本发明的一种基带信号的提取方法,包括以下步骤:
步骤1:获取探测器的下行射频信号s(t);所述射频信号s(t)包括至少两个基带信号,基带信号的个数为N;
步骤2:对射频信号s(t)进行处理;
步骤3:从射频信号s(t)中提取基带信号sb(t);
步骤4:对得到的基带信号sb(t)进行误差分析;在提取基带信号过程中产生的总误差为RMT,所述总误差包括误差RM和误差RM′两部分。
实施例2:
本发明的另一种基带信号的提取方法,包括以下步骤:
步骤1:获取探测器的下行射频信号s(t);所述射频信号s(t)包括至少两个基带信号,基带信号的个数为N;
所述射频信号s(t)为:
其中,PT为探测器射频信号的总功率,fc为载波频率,d(t)或sin(2πft)为副载波,θ为其调制系数;mTM为射频信号的调制系数,fTM为射频信号的频率;f1,f2,fN分别为基带信号1、基带信号2和基带信号N的频率,m1,m2,mN分别为基带信号1、基带信号2和基带信号N的调制系数。
步骤2:对射频信号s(t)进行处理;
所述步骤2具体包括:
步骤2.1:对射频信号s(t)进行化简;
其中,fi为各个基带信号的频率,mi为各个基带信号的调制系数,i的取值范围为1~N;
步骤2.2:根据欧拉公式式(2)中的射频信号s(t)转换为:
由于ω=2πf,对式(3)进行整理,得到射频信号s(t)的泰勒展开式(4):
步骤2.3:取射频信号s(t)的泰勒展开式(4)的正频部分,忽略负频部分,得到下式:
步骤3:从射频信号s(t)中提取基带信号sb(t);
所述步骤3具体包括:
步骤3.1:假设射频信号s(t)信号中只有两个基带信号,即取N=2,根据式(5)在正频域的表现形式,其载波信号和各个基带信号的频域幅值强度分别为PC,P1,P2;
步骤3.2:根据该射频信号s(t)的频谱特点,构建一个只含有基带信号的表达式sb(t):
其中,式(6)中的aik为:
步骤4:对得到的基带信号sb(t)进行误差分析;在提取基带信号过程中产生的总误差为RMT,所述总误差包括误差RM和误差RM′两部分。
所述步骤4具体包括:
步骤4.1:得到所述误差RM;在得到式(5)泰勒展开式开式时,对k取有限的自然数M,产生误差RM,误差RM为:
步骤4.2:得到所述误差RM′;根据需要的精度,对式(6)中的k取有限的自然数M,取值时产生误差RM′,误差RM′为:
其中,RMi为各个基带信号展开项的误差,由下式表示:
步骤4.3:在提取基带信号过程中产生的总误差RMT为误差RM与误差RM′的乘积,如下式:
作为一种可实施的方式,为了验证本方法的可行性及有效性,才用Matlab仿真软件构建了一个含有两个基带信号的射频信号s(t)的函数,根据公式(5)画出其在正频域的表现形式,如图2所示为射频信号的频谱示意图,图中取N=2;如图3所示为根据射频信号频谱构建的基带的频谱示意图。
假设射频信号的发射总功率PT为1,则式(1)转化为:
其中,fc,f1,f2分别取900kHz、20kHz和50kHz,并假设载波和各个基带信号的频域幅值强度为PC,P1,P2;根据该信号的频谱特点,并对式(5)做泰勒七次展开式,同时构建一个只含有两个基带信号的函数sb(t):
sb(t)=a1sin(2πf1t)+a2sin(2πf2t),
其中,系数a1,a2由式(9)得出;
图4为得到的射频信号s(t)的频谱图,图5为得到的两个基带信号的sb(t)的频谱图,其中,f1,f2分别为20kHz和50kHz;为了方便分析,将射频信号降频之后并和基带信号的频谱进行比较,如图6所示为射频信号降频后与基带信号的频谱比较图;从图中可以看出,降频后的射频信号的频谱与基带信号的频谱吻合的很好,误差低于0.1%;证明了该基带信号提取方法的可行性及正确性。
以上,是为了本领域技术人员理解本发明,而对本发明所进行的详细描述,但可以想到,在不脱离本发明的权利要求所涵盖的范围内还可以做出其它的变化和修改,这些变化和修改均在本发明的保护范围内。
Claims (5)
1.一种基带信号的提取方法,其特征在于,包括以下步骤:
步骤1:获取探测器的下行射频信号s(t);所述射频信号s(t)包括至少两个基带信号,基带信号的个数为N;
步骤2:对射频信号s(t)进行处理;
步骤3:从射频信号s(t)中提取基带信号sb(t);
步骤4:对得到的基带信号sb(t)进行误差分析;在提取基带信号过程中产生的总误差为RMT,所述总误差包括误差RM和误差RM′两部分。
2.根据权利要求1所述的一种基带信号的提取方法,其特征在于,所述射频信号s(t)为:
其中,PT为探测器射频信号的总功率,fc为载波频率,d(t)或sin(2πft)为副载波,θ为其调制系数;mTM为射频信号的调制系数,fTM为射频信号的频率;f1,f2,fN分别为基带信号1、基带信号2和基带信号N的频率,m1,m2,mN分别为基带信号1、基带信号2和基带信号N的调制系数。
3.根据权利要求1所述的一种基带信号的提取方法,其特征在于,所述步骤2具体包括:步骤2.1:对射频信号s(t)进行化简;
其中,fi为各个基带信号的频率,mi为各个基带信号的调制系数,i的取值范围为1~N;
步骤2.2:根据欧拉公式式(2)中的射频信号s(t)转换为:
由于ω=2πf,对式(3)进行整理,得到射频信号s(t)的泰勒展开式(4):
步骤2.3:取射频信号s(t)的泰勒展开式(4)的正频部分,忽略负频部分,得到下式:
4.根据权利要求1所述的一种基带信号的提取方法,其特征在于,所述步骤3具体包括:
步骤3.1:假设射频信号s(t)信号中只有两个基带信号,即取N=2,根据式(5)在正频域的表现形式,其载波信号和各个基带信号的频域幅值强度分别为PC,P1,P2;
步骤3.2:根据该射频信号s(t)的频谱特点,构建一个只含有基带信号的表达式sb(t):
其中,式(6)中的aik为:
5.根据权利要求1所述的一种基带信号的提取方法,其特征在于,所述步骤4具体包括:
步骤4.1:得到所述误差RM;在得到式(5)泰勒展开式时,对k取有限的自然数M,产生误差RM,误差RM为:
步骤4.2:得到所述误差RM′;根据需要的精度,对式(6)中的k取有限的自然数M,取值时产生误差RM′,误差RM′为:
其中,RMi为各个基带信号展开项的误差,由下式表示:
步骤4.3:在提取基带信号过程中产生的总误差RMT为误差RM与误差RM′的乘积,如下式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810449706.4A CN108649968B (zh) | 2018-05-11 | 2018-05-11 | 一种基带信号的提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810449706.4A CN108649968B (zh) | 2018-05-11 | 2018-05-11 | 一种基带信号的提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108649968A true CN108649968A (zh) | 2018-10-12 |
CN108649968B CN108649968B (zh) | 2021-04-27 |
Family
ID=63754713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810449706.4A Expired - Fee Related CN108649968B (zh) | 2018-05-11 | 2018-05-11 | 一种基带信号的提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108649968B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770660A (zh) * | 2004-10-06 | 2006-05-10 | 美国博通公司 | 接收器内处理信号的方法和系统 |
CN1933377A (zh) * | 2006-09-28 | 2007-03-21 | 上海大学 | 基于插入导频法毫米波光纤传输系统的双向传输结构及信号传输方法 |
KR20080067166A (ko) * | 2007-01-15 | 2008-07-18 | 엘지전자 주식회사 | 광대역 수신기 |
CN103546400A (zh) * | 2013-05-13 | 2014-01-29 | 上海数字电视国家工程研究中心有限公司 | 单多载波信号识别装置、方法及接收机 |
US20140119485A1 (en) * | 2010-09-24 | 2014-05-01 | Electronics And Telecommunications Research Institute | Apparatus and method for receiving dual band rf signals simultaneously |
CN103973353A (zh) * | 2009-03-31 | 2014-08-06 | 马维尔国际贸易有限公司 | 用于无线通信的探测和引导协议 |
US9473181B1 (en) * | 2015-12-10 | 2016-10-18 | International Business Machines Corporation | Digital combiner for phased-array RF receivers |
-
2018
- 2018-05-11 CN CN201810449706.4A patent/CN108649968B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770660A (zh) * | 2004-10-06 | 2006-05-10 | 美国博通公司 | 接收器内处理信号的方法和系统 |
CN1933377A (zh) * | 2006-09-28 | 2007-03-21 | 上海大学 | 基于插入导频法毫米波光纤传输系统的双向传输结构及信号传输方法 |
KR20080067166A (ko) * | 2007-01-15 | 2008-07-18 | 엘지전자 주식회사 | 광대역 수신기 |
CN103973353A (zh) * | 2009-03-31 | 2014-08-06 | 马维尔国际贸易有限公司 | 用于无线通信的探测和引导协议 |
US20140119485A1 (en) * | 2010-09-24 | 2014-05-01 | Electronics And Telecommunications Research Institute | Apparatus and method for receiving dual band rf signals simultaneously |
CN103546400A (zh) * | 2013-05-13 | 2014-01-29 | 上海数字电视国家工程研究中心有限公司 | 单多载波信号识别装置、方法及接收机 |
US9473181B1 (en) * | 2015-12-10 | 2016-10-18 | International Business Machines Corporation | Digital combiner for phased-array RF receivers |
Also Published As
Publication number | Publication date |
---|---|
CN108649968B (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104104400B (zh) | 高灵敏度星载ads-b信号接收机 | |
CN109886116B (zh) | 基于ads-b的辐射源个体识别方法及装置 | |
CN109307862A (zh) | 一种目标辐射源个体识别方法 | |
CN103051458B (zh) | 一种从bpsk信号中获取稳健射频指纹的方法 | |
CN106130583B (zh) | 一种基于fpga的数字中频信号解调方法 | |
CN107896203A (zh) | Bpsk解调 | |
CN110113288A (zh) | 一种基于机器学习的ofdm解调器的设计和解调方法 | |
CN108226921A (zh) | 基于cpu和gpu架构的二次雷达信息处理系统 | |
Gong et al. | A generative adversarial network based framework for specific emitter characterization and identification | |
CN106169070A (zh) | 基于协作表示的通信辐射源个体识别方法及系统 | |
CN108123723A (zh) | 一种基于无线电技术的定位系统及方法 | |
CN105515711B (zh) | 联合压缩感知与接收分集的测距仪干扰抑制系统 | |
CN107656254A (zh) | 非正交被动mimo雷达模糊函数分析方法 | |
CN108134639B (zh) | 一种基于频谱监测数据的定频信号整编方法 | |
CN108649968A (zh) | 一种基带信号的提取方法 | |
CN103763234A (zh) | 一种基于多路时间相位调制信号的信息传输方法 | |
CN104483682B (zh) | S波段统一测控系统及其侧音提取模块和测距音转发方法 | |
CN106059708B (zh) | 一种多码率数据无线传输系统 | |
CN109309541A (zh) | 一种基于无人机数据链抗干扰的半实物验证系统及验证方法 | |
CN103345648A (zh) | 一种rfid标签设备指纹认证系统 | |
CN103487814B (zh) | 基于场强定向/卫星定位的无线目标识别搜寻方法与装置 | |
CN101283293A (zh) | 闪电探测 | |
CN106130944A (zh) | 脉冲调制信号接收处理系统和方法 | |
CN205176303U (zh) | 卫星三频信标接收装置 | |
CN108770082A (zh) | 基于平滑伪Winger-Ville分布和最优时频分布的通信基站 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210427 |
|
CF01 | Termination of patent right due to non-payment of annual fee |