CN105954586B - 基于变采样率的空间电磁信号频谱获取方法 - Google Patents
基于变采样率的空间电磁信号频谱获取方法 Download PDFInfo
- Publication number
- CN105954586B CN105954586B CN201610561491.6A CN201610561491A CN105954586B CN 105954586 B CN105954586 B CN 105954586B CN 201610561491 A CN201610561491 A CN 201610561491A CN 105954586 B CN105954586 B CN 105954586B
- Authority
- CN
- China
- Prior art keywords
- transmitter
- transmitting signal
- communication
- radar
- spectrum
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于变采样率的空间电磁信号频谱获取方法,主要解决现有技术运算速度缓慢,造成电子设备无法正常工作的问题。其实现步骤是:1)设置雷达发射机、通信发射机的数量及参数;2)计算雷达发射机和通信发射机的发射信号采样频率及采样点数;3)计算观测点处雷达发射机和通信发射机的发射信号幅值;4)根据发射信号幅值计算观测点处雷达发射机和通信发射机的发射信号频谱;5)根据发射信号频谱计算观测点处雷达发射机和通信发射机的发射信号叠加频谱;6)对整个观测区域重复步骤3至步骤5,获取空间电磁信号频谱。本发明大幅度减小了运算量,提高了电子设备反应速度,可用于避免各种电磁设备相互干扰。
Description
技术领域
本发明属于空间电磁态势感知技术领域,特别涉及一种空间电磁信号频谱获取方法。可用于分析电磁环境,避免各种电磁设备相互干扰。
背景技术
随着信息化的发展,现实环境中的电磁设备越来越多,造成电磁环境异常复杂,如何得到空间电磁信号频谱对研究电磁环境具有重要价值:一是有助于避免各电磁设备间的相互干扰;二是可以引导飞行器在欲动设备中避免电磁干扰,以使自身安全起飞和降落。
针对该问题的研究中,传统获取空间电磁信号频谱的方法是:在观测点处先将所有经过该点的时域信号叠加,再对叠加信号进行频谱变换获得空间电磁信号频谱。该方法需要在每个观测点处均进行一次运算量较大的频谱变换计算,在有大量观测点的情况下,运算速度十分缓慢,将造成电子设备无法正常工作或飞行器无法及时躲避干扰等情况。
近年来,针对空间电磁信号频谱的获取提出的新思路是:在观测点处先将所有经过该点的时域信号进行频谱变换,再对所有变换后的信号进行叠加获得空间电磁信号频谱。西安电子科技大学在其申请的专利“获取空间电磁强度数据的方法”(专利申请号201410103112.x)中公开了一种获取空间电磁信号频谱方法,该方法需要对每个发射信号进行一次全采样,很大程度上增加了系统的负担,同时由于信号中心频率的不同,在对变换后的信号进行叠加时,增加了运算量,将造成空间电磁信号频谱获取出现误差,进而导致电子设备无法正常工作。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于变采样率的空间电磁信号频谱获取方法,以在降低运算量的情况下,大大减轻系统的负担,提高电子设备反应速度,避免各种电磁设备的相互干扰。
为实现上述目的,本发明技术方案包括如下:
(1)在空间放置M台雷达发射机和N台通信发射机,设每台雷达发射机的功率Prm和天线增益Grm均大于零,每台雷达发射机的发射信号为srm(t);设每台通信发射机的功率Pcn和天线增益Gcn均大于零,每台通信发射机的发射信号为scn(t),其中,M≥1,N≥1,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N;
(2)根据各雷达发射机的发射信号带宽Brm和各通信发射机的发射信号带宽Bcn,分别计算各雷达发射机的发射信号采样频率frm、采样点数Krm和各通信发射机的发射信号采样频率fcn、采样频率Kcn;
(3)设置一个观测点,根据各雷达发射机的雷达发射机的功率Prm、天线增益Grm和各通信发射机的功率Pcn、天线增益Gcn,分别计算该观测点处各雷达发射机的发射信号幅值Arm和各通信发射机的发射信号幅值Acn;
(4)根据各雷达发射机的发射信号采样频率frm和各雷达发射机的发射信号幅值Arm、各通信发射机的发射信号采样频率fcn和各通信发射机的发射信号幅值Acn,分别计算观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f);
(5)根据观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f),分别计算观测点处雷达发射机的发射信号叠加频谱Sr(f)和通信发射机的发射信号叠加频谱Sc(f):
(5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f):
(5a1)以雷达发射机的发射信号采样频率范围作为x轴坐标,观测点处雷达发射机的发射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,其横坐标轴单位长度为:
(5a2)找到各雷达发射机的发射信号频谱Srm(f)中心频率点的值在雷达频谱坐标系的横坐标轴上所对应的点yrm;
(5a3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中心频率点与在雷达频谱坐标系的横坐标轴上所对应的点yrm重合;
(5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Srm(f)重叠部分进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f);
(5b)计算观测点处通信发射机的发射信号叠加频谱Sc(f):
(5b1)以通信发射机的发射信号采样频率范围作为x轴坐标,观测点处通信发射机的发射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,其横坐标轴单位长度为
(5b2)找到各通信发射机的发射信号频谱Scn(f)中心频率点的值在通信频谱坐标系的横坐标轴上所对应的点zcn;
(5b3)将各通信发射机的发射信号频谱Scn(f)放置在通信频谱坐标系中,使其中心频率点与在通信频谱坐标系的横坐标轴上所对应的点zcn重合;
(5b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Scn(f)重叠部分进行相加,得到观测点处的通信发射机的发射信号叠加频谱Sc(f);
(6)对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域的空间电磁信号频谱。
本发明与现有技术相比具有以下优点:
1)本发明设置各发射信号时间长度一致,仅仅取其各自带宽内的信号进行采样,灵活构建采样系数,避免了对每个发射信号进行一次全采样,因而很大程度上减轻了系统的负担。
2)本发明通过构建频谱坐标系,在每个观测点处,只需将各发射信号频谱搬移至频谱坐标系内并将重叠部分叠加,即可获得空间电磁信号频谱,避免了在每个观测点处均进行一次运算量较大的频域变换,大大提高了电子设备的运算速率,节省了运算时间,使电子装备更具灵活性。
附图说明
图1是本发明的实现流程图。
具体实施方式
以下参照附图,对本发明的技术方案和效果作进一步的详细说明。
参照图1,本发明的实现步骤如下:
步骤1:设置雷达发射机和通信发射机的数量及相关参数。
在空间放置M台雷达发射机和N台通信发射机;
设每台雷达发射机的功率Prm和天线增益Grm均大于零,每台雷达发射机的发射信号为srm(t);
设每台通信发射机的功率Pcn和天线增益Gcn均大于零,每台通信发射机的发射信号为scn(t),
其中,M≥1,N≥1,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N。
步骤2:分别计算各雷达发射机的发射信号采样频率frm、采样点数Krm和各通信发射机的发射信号采样频率fcn、采样频率Kcn。
(2a)计算各雷达发射机的发射信号采样频率frm、采样点数Krm:
(2a1)根据各雷达发射机的发射信号带宽Brm,计算各雷达发射机的发射信号采样频率frm:
frm=ξBrm
其中,ξ表示采样系数,ξ的取值范围是1.2<ξ<1.3;
(2a2)根据各雷达发射机的发射信号时间长度Trm和各雷达发射机的发射信号采样频率frm,计算雷达发射信号采样点数Krm:
Krm=Trmfrm;
(2b)计算各通信发射机的发射信号采样频率fcn、采样频率Kcn:
(2b1)根据各通信发射机发射信号的带宽Bcn,计算各通信发射机发射信号的采样频率fcn:
fcn=ξBcn;
(2b2)根据各通信发射信号时间长度Tcn和各通信发射信号的采样频率fcn,计算通信发射信号采样点数Kcn:
Kcn=Tcnfcn。
步骤3:计算观测点处各雷达发射机的发射信号幅值Arm和各通信发射机的发射信号幅值Acn。
(3a)根据雷达发射机的功率Prm和天线增益Grm,计算观测点处各雷达发射机的发射信号幅值Arm:
其中,Rrm表示雷达发射机观测点的距离;
(3b)根据通信发射机的功率Pcn和天线增益Gcn,计算观测点处各通信发射机的发射信号幅值Acn:
其中,Rcn表示通信发射机到观测点的距离。
步骤4:分别计算观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f)。
(4a)计算观测点处各雷达发射机的发射信号频谱Srm(f):
(4a1)对各雷达发射机的发射信号srm(t)以采样频率frm进行采样,并对采样信号进行傅里叶变换,得到雷达发射机的发射信号采样频谱
(4a2)根据雷达发射机的发射信号采样频谱计算观测点处各雷达发射机的发射信号频谱Srm(f):
其中,τrm=Rrm/c表示观测点与雷达发射机的时间延迟,Rrm表示雷达发射机到观测点的距离,c表示光速;
(4b)计算观测点处各通信发射机的发射信号频谱Scn(f):
(4b1)对各通信发射机的发射信号scn(t)以采样频率fcn进行采样,并对采样信号进行傅里叶变换,得到通信发射机的发射信号采样频谱
(4b2)根据通信发射机的发射信号采样频谱计算各通信发射机的发射信号频谱Scn(f):
其中,τcn=Rcn/c表示观测点与通信发射机的时间延迟,Rcn表示通信发射机到观测点的距离,c表示光速。
步骤5:分别计算观测点处雷达发射机的发射信号叠加频谱Sr(f)和通信发射机的发射信号叠加频谱Sc(f)。
(5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f):
(5a1)以雷达发射机的发射信号采样频率范围作为x轴坐标,观测点处雷达发射机的发射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,该频谱坐标系的横坐标轴单位长度为:
(5a2)找到各雷达发射机的发射信号频谱Srm(f)中心频率点的值在雷达频谱坐标系横坐标轴上所对应的点yrm;
(5a3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中心频率点与在雷达频谱坐标系的横坐标轴上所对应的点yrm重合;
(5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Srm(f)重叠部分进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f);
(5b)计算观测点处通信发射机的发射信号叠加频谱Sc(f):
(5b1)以通信发射机的发射信号采样频率范围作为x轴坐标,观测点处通信发射机的发射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,该频谱坐标系的横坐标轴单位长度为
(5b2)找到各通信发射机的发射信号频谱Scn(f)中心频率点的值在通信频谱坐标系的横坐标轴上所对应的点zcn;
(5b3)将各通信发射机的发射信号频谱Scn(f)放置在通信频谱坐标系中,使其中心频率点与在通信频谱坐标系的横坐标轴上所对应的点zcn重合;
(5b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Scn(f)重叠部分进行相加,得到观测点处的通信发射机的发射信号叠加频谱Sc(f)。
步骤6:对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域的空间电磁信号频谱。
以上描述仅是本发明的一个具体实例,显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修改和改变,但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。
Claims (4)
1.一种基于变采样率的空间电磁信号频谱获取方法,包括:
(1)在空间放置M台雷达发射机和N台通信发射机,设每台雷达发射机的功率Prm和天线增益Grm均大于零,每台雷达发射机的发射信号为srm(t);设每台通信发射机的功率Pcn和天线增益Gcn均大于零,每台通信发射机的发射信号为scn(t),其中,M≥1,N≥1,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N;
(2)根据各雷达发射机的发射信号带宽Brm和各通信发射机的发射信号带宽Bcn,分别计算各雷达发射机的发射信号采样频率frm、采样点数Krm和各通信发射机的发射信号采样频率fcn、采样点数Kcn;
(3)设置一个观测点,根据各雷达发射机的雷达发射机的功率Prm、天线增益Grm和各通信发射机的功率Pcn、天线增益Gcn,分别计算该观测点处各雷达发射机的发射信号幅值Arm和各通信发射机的发射信号幅值Acn;
(4)根据各雷达发射机的发射信号采样频率frm和各雷达发射机的发射信号幅值Arm、各通信发射机的发射信号采样频率fcn和各通信发射机的发射信号幅值Acn,分别计算观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f);
(5)根据观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f),分别计算观测点处雷达发射机的发射信号叠加频谱Sr(f)和通信发射机的发射信号叠加频谱Sc(f):
(5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f):
(5a1)以雷达发射机的发射信号采样频率范围作为x轴坐标,观测点处雷达发射机的发射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,其横坐标轴单位长度为:
(5a2)找到各雷达发射机的发射信号频谱Srm(f)中心频率点的值在雷达频谱坐标系的横坐标轴上所对应的点yrm;
(5a3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中心频率点与在雷达频谱坐标系的横坐标轴上所对应的点yrm重合;
(5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Srm(f)重叠部分进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f);
(5b)计算观测点处通信发射机的发射信号叠加频谱Sc(f):
(5b1)以通信发射机的发射信号采样频率范围作为x轴坐标,观测点处通信发射机的发射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,其横坐标轴单位长度为
(5b2)找到各通信发射机的发射信号频谱Scn(f)中心频率点的值在通信频谱坐标系的横坐标轴上所对应的点zcn;
(5b3)将各通信发射机的发射信号频谱Scn(f)放置在通信频谱坐标系中,使其中心频率点与在通信频谱坐标系的横坐标轴上所对应的点zcn重合;
(5b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Scn(f)重叠部分进行相加,得到观测点处的通信发射机的发射信号叠加频谱Sc(f);
(6)对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域的空间电磁信号频谱。
2.根据权利要求1所述的方法,其中步骤(2)中分别计算各雷达发射机的发射信号采样频率frm、采样点数Krm和各通信发射机的发射信号采样频率fcn、采样点数Kcn,按如下步骤进行:
(2a)计算各雷达发射机的发射信号采样频率frm、采样点数Krm:
(2a1)根据各雷达发射机的发射信号带宽Brm,计算各雷达发射机的发射信号采样频率frm:
frm=ξBrm
其中,ξ表示采样系数,对复信号采样时ξ的取值范围是1.2<ξ<1.3,对实信号采样时ξ的取值范围是2.4<ξ<2.6;
(2a2)根据各雷达发射机的发射信号时间长度Trm和各雷达发射机的发射信号采样频率frm,计算雷达发射信号采样点数Krm:
Krm=Trmfrm;
(2b)计算各通信发射机的发射信号采样频率fcn、采样点数Kcn:
(2b1)根据各通信发射机发射信号的带宽Bcn,计算各通信发射机发射信号的采样频率fcn:
fcn=ξBcn;
(2b2)根据各通信发射信号时间长度Tcn和各通信发射信号的采样频率fcn,计算通信发射信号采样点数Kcn:
Kcn=Tcnfcn。
3.根据权利要求1所述的方法,其中步骤(3)中分别计算观测点处各雷达发射机的雷达信号幅值Arm和各通信发射机的通信信号幅值Acn,按如下公式计算:
其中,Rrm和Rcn分别表示雷达发射机和通信发射机到观测点的距离。
4.根据权利要求1所述的方法,其中步骤(4)中分别计算观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频谱Scn(f),按如下步骤进行:
(4a)计算观测点处各雷达发射机的发射信号频谱Srm(f):
(4a1)对各雷达发射机的发射信号srm(t)以采样频率frm进行采样,并对采样信号进行傅里叶变换,得到雷达发射机的发射信号采样频谱
(4a2)根据雷达发射机的发射信号采样频谱计算观测点处各雷达发射机的发射信号频谱Srm(f):
其中,τrm=Rrm/c表示观测点与雷达发射机的时间延迟,Rrm表示雷达发射机到观测点的距离,c表示光速;
(4b)计算观测点处各通信发射机的发射信号频谱Scn(f):
(4b1)对各通信发射机的发射信号scn(t)以采样频率fcn进行采样,并对采样信号进行傅里叶变换,得到通信发射机的发射信号采样频谱
(4b2)根据通信发射机的发射信号采样频谱计算各通信发射机的发射信号频谱Scn(f):
其中,τcn=Rcn/c表示观测点与通信发射机的时间延迟,Rcn表示通信发射机到观测点的距离,c表示光速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610561491.6A CN105954586B (zh) | 2016-07-15 | 2016-07-15 | 基于变采样率的空间电磁信号频谱获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610561491.6A CN105954586B (zh) | 2016-07-15 | 2016-07-15 | 基于变采样率的空间电磁信号频谱获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105954586A CN105954586A (zh) | 2016-09-21 |
CN105954586B true CN105954586B (zh) | 2018-11-30 |
Family
ID=56900832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610561491.6A Active CN105954586B (zh) | 2016-07-15 | 2016-07-15 | 基于变采样率的空间电磁信号频谱获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105954586B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112711032B (zh) * | 2020-12-07 | 2022-05-13 | 中国人民解放军海军航空大学 | 一种基于图数据和gcn的雷达目标检测方法及系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3338370B2 (ja) * | 1998-05-14 | 2002-10-28 | 株式会社アドバンテスト | 周波数分析方法及びこの方法を用いた掃引型スペクトラム・アナライザ |
CN100538394C (zh) * | 2007-04-06 | 2009-09-09 | 清华大学 | 一种采用多发多收频分信号的宽带雷达及其成像方法 |
CN101452073B (zh) * | 2007-11-30 | 2011-12-28 | 清华大学 | 基于多发多收频分雷达的宽带信号合成方法 |
CN103852648B (zh) * | 2014-03-20 | 2016-07-06 | 西安电子科技大学 | 获取空间电磁强度数据的方法 |
JP6263143B2 (ja) * | 2015-03-20 | 2018-01-17 | 大井電気株式会社 | 周波数解析装置 |
CN104898107B (zh) * | 2015-06-17 | 2017-06-20 | 西安电子科技大学 | 一种多发多收合成孔径激光雷达信号处理方法 |
CN205228470U (zh) * | 2015-11-26 | 2016-05-11 | 陈昊 | 一种无人飞行器预警定位装置 |
-
2016
- 2016-07-15 CN CN201610561491.6A patent/CN105954586B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105954586A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10909864B2 (en) | Drone obstacle avoidance using real-time wind estimation | |
CN105867115B (zh) | 一种非平稳随机振动试验控制方法 | |
CN113204015B (zh) | 用于分布式雷达目标模拟的多无人机航迹规划方法及系统 | |
CN111507008B (zh) | 一种基于参数化的通用雷达模型建模方法 | |
CN109471068B (zh) | 基于射频阵列信号doa估计的无人机定位系统和方法 | |
CN105572473B (zh) | 高分辨率线性时频分析方法 | |
CN105954586B (zh) | 基于变采样率的空间电磁信号频谱获取方法 | |
CN110632553A (zh) | 定位方法、装置、系统、计算机可读存储介质及无人机 | |
CN103852648B (zh) | 获取空间电磁强度数据的方法 | |
CN103675927B (zh) | 固定翼航空电磁系统接收吊舱摆动角度的校正方法 | |
Nasution et al. | UAV telemetry communications using ZigBee protocol | |
CN106291126B (zh) | 基于变采样率的复杂电磁环境时域信号模拟方法 | |
CN109188353B (zh) | 基于多普勒频率差和压缩感知的单站无源定位方法 | |
CN112710984B (zh) | 一种基于同伦延拓的无源定位方法及系统 | |
Xia et al. | Over-the-air radar emitter signal classification based on SDR | |
CN112422214B (zh) | 一种适用于航空信道的通信效果演示验证系统 | |
CN106371072B (zh) | 一种基于单个脉冲频域采样的空间信号频谱普查方法 | |
CN109828276A (zh) | 一种基于干涉算法的风廓线雷达 | |
WO2017150323A1 (ja) | 位置推定システム | |
CN208399498U (zh) | 一种应用于智慧城市的土壤检测装置 | |
Shaheen et al. | Robust Air-to-Air Channel Model for Swarms of Drones in Search and Rescue Missions | |
CN106291528A (zh) | 一种用于多旋翼飞行器的微波雷达测速测距系统及方法 | |
CN205404801U (zh) | 一种快速电离层探测系统 | |
CN109030773A (zh) | 一种应用于智慧城市的土壤检测装置、系统及方法 | |
Fu et al. | A DNN-Based Channel State and Scenario Identifications for Millimeter-Wave Communications |
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 |