基于变采样率的空间电磁信号频谱获取方法

文档序号:10592749阅读:299来源:国知局
基于变采样率的空间电磁信号频谱获取方法
【专利摘要】本发明公开了一种基于变采样率的空间电磁信号频谱获取方法,主要解决现有技术运算速度缓慢,造成电子设备无法正常工作的问题。其实现步骤是:1)设置雷达发射机、通信发射机的数量及参数;2)计算雷达发射机和通信发射机的发射信号采样频率及采样点数;3)计算观测点处雷达发射机和通信发射机的发射信号幅值;4)根据发射信号幅值计算观测点处雷达发射机和通信发射机的发射信号频谱;5)根据发射信号频谱计算观测点处雷达发射机和通信发射机的发射信号叠加频谱;6)对整个观测区域重复步骤3至步骤5,获取空间电磁信号频谱。本发明大幅度减小了运算量,提高了电子设备反应速度,可用于避免各种电磁设备相互干扰。
【专利说明】
基于变采样率的空间电磁信号频谱获取方法
技术领域
[0001] 本发明属于空间电磁态势感知技术领域,特别设及一种空间电磁信号频谱获取方 法。可用于分析电磁环境,避免各种电磁设备相互干扰。
【背景技术】
[0002] 随着信息化的发展,现实环境中的电磁设备越来越多,造成电磁环境异常复杂,如 何得到空间电磁信号频谱对研究电磁环境具有重要价值:一是有助于避免各电磁设备间的 相互干扰;二是可W引导飞行器在欲动设备中避免电磁干扰,W使自身安全起飞和降落。
[0003] 针对该问题的研究中,传统获取空间电磁信号频谱的方法是:在观测点处先将所 有经过该点的时域信号叠加,再对叠加信号进行频谱变换获得空间电磁信号频谱。该方法 需要在每个观测点处均进行一次运算量较大的频谱变换计算,在有大量观测点的情况下, 运算速度十分缓慢,将造成电子设备无法正常工作或飞行器无法及时躲避干扰等情况。
[0004] 近年来,针对空间电磁信号频谱的获取提出的新思路是:在观测点处先将所有经 过该点的时域信号进行频谱变换,再对所有变换后的信号进行叠加获得空间电磁信号频 谱。西安电子科技大学在其申请的专利"获取空间电磁强度数据的方法"(专利申请号 201410103112.X)中公开了一种获取空间电磁信号频谱方法,该方法需要对每个发射信号 进行一次全采样,很大程度上增加了系统的负担,同时由于信号中屯、频率的不同,在对变换 后的信号进行叠加时,增加了运算量,将造成空间电磁信号频谱获取出现误差,进而导致电 子设备无法正常工作。

【发明内容】

[0005] 本发明的目的在于针对上述现有技术存在的不足,提出一种基于变采样率的空间 电磁信号频谱获取方法,W在降低运算量的情况下,大大减轻系统的负担,提高电子设备反 应速度,避免各种电磁设备的相互干扰。
[0006] 为实现上述目的,本发明技术方案包括如下:
[0007] (1)在空间放置M台雷达发射机和N台通信发射机,设每台雷达发射机的功率Prm和 天线增益Grm均大于零,每台雷达发射机的发射信号为Sr"(t);设每台通信发射机的功率P。。 和天线增益Gcn均大于零,每台通信发射机的发射信号为Scn(t),其中,1,修1,m表示雷 达发射机序号,m=l,2,. . .,M,n表示通信发射机序号,n = l,2,. . .,N;
[0008] (2)根据各雷达发射机的发射信号带宽Brm和各通信发射机的发射信号带宽B。。,分 别计算各雷达发射机的发射信号采样频率fr"、采样点数Krm和各通信发射机的发射信号采 样频率f。。、采样频率K。。;
[0009] (3)设置一个观测点,根据各雷达发射机的雷达发射机的功率Pr"、天线增益G?和 各通信发射机的功率P。。、天线增益G。。,分别计算该观测点处各雷达发射机的发射信号幅值 Arm和各通信发射机的发射信号幅值A。。;
[0010] (4)根据各雷达发射机的发射信号采样频率和各雷达发射机的发射信号幅值 Arm、各通信发射机的发射信号采样频率f。。和各通信发射机的发射信号幅值A。。,分别计算观 测点处各雷达发射机的发射信号频谱Sr"(f)和各通信发射机的发射信号频谱Stn(f);
[0011] (5)根据观测点处各雷达发射机的发射信号频谱Sr"(f)和各通信发射机的发射信 号频谱Scn(f),分别计算观测点处雷达发射机的发射信号叠加频谱Sr(f)和通信发射机的发 射信号叠加频谱Sc(f):
[0012] (5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f):
[0013] (5al) W雷达发射机的发射信号采样频率范围作为X轴坐标,观测点处雷达发射机 的发射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,其横坐标轴单位长度为:
[0014] (5a2)找到各雷达发射机的发射信号频谱Srm(f)中屯、频率点的值在雷达频谱坐标 系的横坐标轴上所对应的点yrm;
[001引(5曰3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中屯、 频率点与在雷达频谱坐标系的横坐标轴上所对应的点yrm重合;
[0016] (5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Sr"(f)重叠 部分进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f);
[0017] 师)计算观测点处通信发射机的发射信号叠加频谱Sc(f):
[0018] 巧bl) W通信发射机的发射信号采样频率范围作为X轴坐标,观测点处通信发射机 的发射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,其横坐标轴单位长度为
[0019] (5b2)找到各通信发射机的发射信号频谱Scn(f)中屯、频率点的值在通信频谱坐标 系的横坐标轴上所对应的点Zcn ;
[0020] 巧b3)将各通信发射机的发射信号频谱Scn(f)放置在通信频谱坐标系中,使其中屯、 频率点与在通信频谱坐标系的横坐标轴上所对应的点Z。。重合;
[0021] 巧b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Stn(f)重叠部分进行 相加,得到观测点处的通信发射机的发射信号叠加频谱ScXf );
[0022] (6)对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域的 空间电磁信号频谱。
[0023] 本发明与现有技术相比具有W下优点:
[0024] 1)本发明设置各发射信号时间长度一致,仅仅取其各自带宽内的信号进行采样, 灵活构建采样系数,避免了对每个发射信号进行一次全采样,因而很大程度上减轻了系统 的负担。
[0025] 2)本发明通过构建频谱坐标系,在每个观测点处,只需将各发射信号频谱搬移至 频谱坐标系内并将重叠部分叠加,即可获得空间电磁信号频谱,避免了在每个观测点处均 进行一次运算量较大的频域变换,大大提高了电子设备的运算速率,节省了运算时间,使电 子装备更具灵活性。
【附图说明】
[0026] 图I是本发明的实现流程图。
【具体实施方式】
[0027] W下参照附图,对本发明的技术方案和效果作进一步的详细说明。
[00%]参照图1,本发明的实现步骤如下:
[0029] 步骤1:设置雷达发射机和通信发射机的数量及相关参数。
[0030] 在空间放置M台雷达发射机和N台通信发射机;
[0031] 设每台雷达发射机的功率Prm和天线增益Grm均大于零,每台雷达发射机的发射信 号为Srm(t);
[0032] 设每台通信发射机的功率Pen和天线增益Gcn均大于零,每台通信发射机的发射信 号为Sc;n(t),
[0033] 其中,1,N> 1,m表示雷达发射机序号,m= 1,2,. . .,M,n表示通信发射机序号,n =1,2,... ,No
[0034] 步骤2:分别计算各雷达发射机的发射信号采样频率frm、采样点数Krm和各通信发 射机的发射信号采样频率f。。、采样频率K。。。
[0035] (2a)计算各雷达发射机的发射信号采样频率frm、采样点数Krm:
[0036] (2al)根据各雷达发射机的发射信号带宽Brm,计算各雷达发射机的发射信号采样 频率frm:
[0037] frm=巧 rm
[003引其中,C表示采样系数,C的取值范围是1.2作<1.3;
[0039] (2a2)根据各雷达发射机的发射信号时间长度Trm和各雷达发射机的发射信号采样 频率frm,计算雷达发射信号采样点数Krm:
[0040] Krm = Trmfrm;
[0041 ] (2b)计算各通信发射机的发射信号采样频率f。。、采样频率Kcn:
[0042] (化1)根据各通信发射机发射信号的带宽B。。,计算各通信发射机发射信号的采样 频率fen:
[0043] fen = CBcn;
[0044] (化2)根据各通信发射信号时间长度Ten和各通信发射信号的采样频率fen,计算通 信发射信号采样点数Kcn:
[0045] Kcn = TcnfcnD
[0046] 步骤3:计算观测点处各雷达发射机的发射信号幅值Arm和各通信发射机的发射信 号幅值Acji。
[0047] (3a)根据雷达发射机的功率Prm和天线增益Grm,计算观测点处各雷达发射机的发 射信号幅值Arm:
[004引
[0049] 其中,Rrm表示雷达发射机观测点的距离;
[0050] (3b)根据通信发射机的功率Pen和天线增益G。。,计算观测点处各通信发射机的发 射信号幅值Acn:
[0化1 ]
[0052] 其中,Ren表示通信发射机到观测点的距离。
[0053] 步骤4:分别计算观测点处各雷达发射机的发射信号频谱Sr"(f)和各通信发射机的 发射信号频谱Sen(f)。
[0054] (4a)计算观测点处各雷达发射机的发射信号频谱Sr"(f):
[0055] (4al)对各雷达发射机的发射信号Srm(t)w采样频率frm进行采样,并对采样信号 进行傅里叶变换,得到雷达发射机的发射信号采样频谱含w(/);
[0化6] (4a2)根据雷达发射机的发射信号采样频谱心,(/),计算观测点处各雷达发射机 的发射信号频谱Sr"(f):
[0化7]
[005引其中,Trm = RrmZc表不观测点与雷达发射机的时间延迟,Rrm表示雷达发射机到观测 点的距离,C表示光速;
[0059] (4b)计算观测点处各通信发射机的发射信号频谱Stn(f):
[0060] (4bl)对各通信发射机的发射信号Scn(t) W采样频率fen进行采样,并对采样信号 进行傅里叶变换,得到通信发射机的发射信号采样频谱:
[0061] (4b2)根据通信发射机的发射信号采样频谱 ix/),计算各通信发射机的发射信号 频谱 Ssn(f):
[0062]
[0063] 其中,Tcn = RcVc表示观测点与通信发射机的时间延迟,Rcn表示通信发射机到观测 点的距离,C表示光速。
[0064] 步骤5:分别计算观测点处雷达发射机的发射信号叠加频谱Sr(f)和通信发射机的 发射信号叠加频谱Sc(f)。
[0065] (5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f):
[0066] (5al) W雷达发射机的发射信号采样频率范围作为X轴坐标,观测点处雷达发射机 的发射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,该频谱坐标系的横坐标轴单 位长度为
[0067] (5a2)找到各雷达发射机的发射信号频谱Srm(f)中屯、频率点的值在雷达频谱坐标 系横坐标轴上所对应的点yrm;
[0068] (5a3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中屯、 频率点与在雷达频谱坐标系的横坐标轴上所对应的点yrm重合;
[0069] (5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Sr"(f)重叠 部分进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f);
[0070] (5b)计算观测点处通信发射机的发射信号叠加频谱ScXf):
[0071] 巧bl) W通信发射机的发射信号采样频率范围作为X轴坐标,观测点处通信发射机 的发射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,该频谱坐标系的横坐标轴单 位长度另
[0072] (加2)找到各通信发射机的发射信号频谱Scn(f)中屯、频率点的值在通信频谱坐标 系的横坐标轴上所对应的点Zcn ;
[0073] 巧b3)将各通信发射机的发射信号频谱Scn(f)放置在通信频谱坐标系中,使其中屯、 频率点与在通信频谱坐标系的横坐标轴上所对应的点Z。。重合;
[0074] 巧b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Stn(f)重叠部分进行 相加,得到观测点处的通信发射机的发射信号叠加频谱St(f)。
[0075] 步骤6:对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域 的空间电磁信号频谱。
[0076] W上描述仅是本发明的一个具体实例,显然对于本领域的专业人员来说,在了解 了本
【发明内容】
和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的 各种修改和改变,但是运些基于本发明思想的修正和改变仍在本发明的权利要求保护范围 之内。
【主权项】
1. 一种基于变采样率的空间电磁信号频谱获取方法,包括: (1) 在空间放置Μ台雷达发射机和N台通信发射机,设每台雷达发射机的功率prm和天线 增益Grm均大于零,每台雷达发射机的发射信号为s rm(t);设每台通信发射机的功率?。"和天 线增益6。"均大于零,每台通信发射机的发射信号为Sc;n(t),其中,Μ多1,N多1,m表示雷达发 射机序号,m=l,2, · · ·,M,n表示通信发射机序号,η = 1,2, · · ·,N; (2) 根据各雷达发射机的发射信号带宽Brm和各通信发射机的发射信号带宽Βα,分别计 算各雷达发射机的发射信号采样频率f?、采样点数K rm和各通信发射机的发射信号采样频 率f?、采样频率K?; (3) 设置一个观测点,根据各雷达发射机的雷达发射机的功率Prm、天线增益GrjP各通信 发射机的功率P?、天线增益G?,分别计算该观测点处各雷达发射机的发射信号幅值A rm和各 通信发射机的发射信号幅值 (4) 根据各雷达发射机的发射信号采样频率frm和各雷达发射机的发射信号幅值Arm、各 通信发射机的发射信号采样频率hn和各通信发射机的发射信号幅值,分别计算观测点 处各雷达发射机的发射信号频谱S rm(f)和各通信发射机的发射信号频谱S?(f); (5) 根据观测点处各雷达发射机的发射信号频谱Srm(f)和各通信发射机的发射信号频 谱S?(f·),分别计算观测点处雷达发射机的发射信号叠加频谱S r(f)和通信发射机的发射信 号叠加频谱&(f): (5a)计算观测点处各雷达发射机的发射信号叠加频谱Sr(f): (5al)以雷达发射机的发射信号采样频率范围作为X轴坐标,观测点处雷达发射机的发 射信号频谱幅值范围作为y轴坐标,构建雷达频谱坐标系,其横坐标轴单位长度为:(5a2)找到各雷达发射机的发射信号频谱Srm(f)中心频率点的值在雷达频谱坐标系的 横坐标轴上所对应的点yrm; (5a3)将各雷达发射机的发射信号频谱Srm(f)放置在雷达频谱坐标系中,使其中心频率 点与在雷达频谱坐标系的横坐标轴上所对应的点重合; (5a4)在雷达信号频谱直角坐标系中,将各雷达发射机的发射信号频谱Srm(f)重叠部分 进行相加,得到观测点处雷达发射机的发射信号叠加频谱Sr(f); (5b)计算观测点处通信发射机的发射信号叠加频谱&(f): (5bl)以通信发射机的发射信号采样频率范围作为X轴坐标,观测点处通信发射机的发 射信号频谱幅值范围作为y轴坐标,构建通信频谱坐标系,其横坐标轴单位长度为 (5b2)找到各通信发射机的发射信号频谱&η(?·)中心频率点的值在通信频谱坐标系的 横坐标轴上所对应的点Zcn; (5b3)将各通信发射机的发射信号频谱S?(f)放置在通信频谱坐标系中,使其中心频率 点与在通信频谱坐标系的横坐标轴上所对应的点z?重合; (5b4)在通信频谱坐标系中,将各通信发射机的发射信号频谱Sc^f)重叠部分进行相 加,得到观测点处的通信发射机的发射信号叠加频谱&(f); (6) 对整个观测区域的其他观测点重复步骤(3)至步骤(5),即可获得观测区域的空间 电磁信号频谱。2. 根据权利要求1所述的方法,其中步骤(2)中分别计算各雷达发射机的发射信号采样 频率frm、采样点数K r4P各通信发射机的发射信号采样频率f?、采样频率IU,按如下步骤进 行: (2a)计算各雷达发射机的发射信号采样频率frm、采样点数Krm: (2al)根据各雷达发射机的发射信号带宽Brm,计算各雷达发射机的发射信号采样频率 frm : frm- CBrm 其中,ξ表示采样系数,对复信号采样时ξ的取值范围是1.2〈ξ〈1.3,对实信号采样时ξ的 取值范围是2.4〈ξ〈2.6; (2a2)根据各雷达发射机的发射信号时间长度TrdP各雷达发射机的发射信号采样频率 f?,计算雷达发射信号采样点数Krm: Krm - T rmf rm ; (2b)计算各通信发射机的发射信号采样频率fm、采样频率Kcn: (2bl)根据各通信发射机发射信号的带宽Bm,计算各通信发射机发射信号的采样频率 fen: fen- ζΒ〇η ; (2b2)根据各通信发射信号时间长度1^和各通信发射信号的采样频率,计算通信发 射信号采样点数Kcn: Ken - Tenf cn〇3. 根据权利要求1所述的方法,其中步骤(3)中分别计算观测点处各雷达发射机的雷达 信号幅值Arm和各通信发射机的通信信号幅值A?,按如下公式计算:其中,Rrm和R?分别表示雷达发射机和通信发射机到观测点的距离。4. 根据权利要求1所述的方法,其中步骤(4)中分别计算观测点处各雷达发射机的发射 信号频谱Srm(f)和各通信发射机的发射信号频谱s?(f·),按如下步骤进行: (4a)计算观测点处各雷达发射机的发射信号频谱Srm(f): (4al)对各雷达发射机的发射信号srm(t)以采样频率匕》进行采样,并对采样信号进行傅 里叶变换,得到雷达发射机的发射信号采样频谱Ι,Λ/): (4a2)根据雷达发射机的发射信号采样频谱&"(/),计算观测点处各雷达发射机的发射 信号频谱Srm(f):其中,Trm=Rrm/c表示观测点与雷达发射机的时间延迟,R rm表示雷达发射机到观测点的 距离,c表示光速; (4b)计算观测点处各通信发射机的发射信号频谱S?(f): (4bl)对各通信发射机的发射信号S?(t)以采样频率匕"进行采样,并对采样信号进行傅 里叶变换,得到通信发射机的发射信号采样频谱&"(/); (4b2)根据通信发射机的发射信号采样频谱&"(/),计算各通信发射机的发射信号频谱 Scn(f):其中,Tcn = Rcn/C表示观测点与通信发射机的时间延迟,Rcn表示通信发射机到观测点的 距离,C表示光速。
【文档编号】G01R23/16GK105954586SQ201610561491
【公开日】2016年9月21日
【申请日】2016年7月15日
【发明人】刘高高, 蔡晶晶, 苏瑶, 鲍丹, 武斌, 秦国栋, 李鹏
【申请人】西安电子科技大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1