CN108880648B - 基于微波光子阵列的超宽带信号频率估计方法 - Google Patents
基于微波光子阵列的超宽带信号频率估计方法 Download PDFInfo
- Publication number
- CN108880648B CN108880648B CN201810735983.1A CN201810735983A CN108880648B CN 108880648 B CN108880648 B CN 108880648B CN 201810735983 A CN201810735983 A CN 201810735983A CN 108880648 B CN108880648 B CN 108880648B
- Authority
- CN
- China
- Prior art keywords
- signal
- microwave
- ultra
- frequency
- wideband 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/50—Transmitters
- H04B10/516—Details of coding or modulation
- H04B10/5165—Carrier suppressed; Single sideband; Double sideband or vestigial
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
- H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
- H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
- H04B10/0795—Performance monitoring; Measurement of transmission parameters
- H04B10/07955—Monitoring or measuring power
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/70—Photonic quantum communication
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Optics & Photonics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
Description
技术领域
本发明属于光电通信技术领域,涉及一种超宽带信号频率估计方法,特别涉及一种基于微波光子阵列的超宽带信号频率估计方法,可用于目标检测与无源定位。
背景技术
超宽带信号频率估计是利用空间中的信号接收设备接收多个不同方向的信号源发出的信号,运用现代信号处理方法快速准确的获得信号源频率的技术,在雷达、声纳、无线通信等领域具有重要应用价值。
目前,最新超宽带信号波达方向角估计方法是用微波光子估计超超宽带信号波达方向角,这种方法旨在采用光电器件和光电学方法估计超超宽带信号波达方向角,与传统超超宽带信号波达方向角估计方法相比,具有估计带宽大、损耗小、抗干扰能力强、***小型等优点,能够克服电子瓶颈,适应复杂的电磁环境的优点。例如申请公布号为CN107528638A,名称为“基于微波光子的超宽带接收机装置及实现方法”的专利申请中,公开了一种基于微波光子的超宽带接收机装置及实现方法,该方法将激光源输出的光信号分成两路,一路用于生成频率等间隔的光频梳,采用外调制的方式实现;另一路被待接收射频信号经电光调制器调制,并与光频梳路的光信号合成。该合成信号通过可调谐窄带带通光学滤波器,经该滤波器可将载有待测信号的光频与最临近的光频梳的固定频点滤出,再经过高速探测器变换,实现两个频率拍频合成得到其差频,即得到待测信号下变频的中频信号。本发明将微波光子学中光频梳的技术应用于微波信号接收的设计,结合可调光滤波器,实现了在低频率本振信号(2GHz~5GHz)的条件下,超宽带信号接收和下变频,其中接收带宽十几倍于低频本振信号。但是由于该方法是基于两阵元实现的,因此频率估计范围有限,且估计效率和估计精度较低。
发明内容
本发明的目的在于克服上述现有技术存在的不足,提出了一种基于微波光子阵列的超宽带信号频率估计方法,用于解决现有技术中存在的超宽带信号频率估计效率和估计精度较低的问题。
为实现上述目的,本发明采取的技术方案包括如下步骤:
(1)获取一阶边带信号E0(t):
将天线0接收到的微波信号V0(t)和激光信号源输出的光载波信号E0输入到马赫增德尔调制器MZM-0中进行调制,得到多边带信号,并对多边带信号进行窄带滤波,得到一阶边带信号E0(t);
(2)构建微波光子阵列D:
构建由m个级联的微波时延线组成的微波光子阵列D,相邻微波时延线之间的距离为d,m=1...M,m≥2,第m个微波时延线上级联有马赫增德尔调制器MZM-m、滤波器F-m和光功率计OPM-m,得到微波光子阵元Dm,m个微波光子阵元组成微波光子阵列D;
(3)获取微波光子阵列D输出的超宽带信号Em的光功率Pm:
(3a)将第m个微波时延线接收到的微波时延信号Vm(t)和MZM-0输出的一阶边带信号E0(t)输入到马赫增德尔调制器MZM-m中进行调制,并通过滤波器F-m对调制得到的第m个调制信号进行窄带滤波,得到第m个超宽带信号Em:
其中Vm为马赫增德尔调制器MZM-m的半波电压且V1=V2=...Vm=VM,β1k=β2k=...=βmk...=βMk,τ为相邻微波时延线接收到的微波信号Vm(t)的延迟时间,Ωk表示微波时延线接收到的微波信号的角频率,Ωkm表示m个微波时延线接收到的微波信号的角频率之和,Lm为马赫增德尔调制器MZM-m的***损耗且L1=L2=...Lm=LM,J1(·)表示一阶贝赛尔函数;
(3b)将第m个超宽带信号Em输入到光功率计OPM-m中,计算超宽带信号Em的光功率Pm:
(4)定义光功率比值向量Q:
(4a)计算Pm的相对光功率比值Qm:
(4b)定义M×1维的光功率比值向量:Q=[Q1,Q2,…,Qm…,QM]T,其中[·]T表示矩阵[·]的转置;
(5)定义M×N维矩阵ZD(w)m,k:
(5a)根据微波信号源的空域稀疏特性,采用空间网格划分方法,将可搜索的空间频率域等间隔划分成N份,得到超宽带信号Em的频率的取值空间w为:w=[w1,w2,…,wk…,wN];
(5b)定义与阵列流型向量相关的M×N维的矩阵ZD(w)m,k:
(6)估计微波光子阵列D输出的超宽带信号Em的频率Ωm:
(6b)通过改进的L1-SRACV算法模型,构建基于稀疏重构的约束优化方程式:
其中||·||1和||·||2分别表示l1范数和l2范数,s.t.表示约束关系,ε表示误差的允许值;
本发明与现有技术相比,具有以下优点:
1)本发明在超宽带信号的获取中采用微波光子阵列,无需事先估计入射信号的个数,避免了现有超宽带信号频率估计方法中角度搜索和角度匹配的问题,与现有技术相比,提高了超宽带信号频率的估计效率。
2)本发明在超宽带信号频率的估计中采用稀疏表示的思想将频率度估计问题转化成稀疏重构问题,利用信号源在的空域稀疏特性进行建模,突破了微波信号瑞利限的问题,与现有技术相比,提高了超宽带信号频率的估计精度。
附图说明
图1是实现本发明采用的估计***的结构示意图;
图2是本发明的实现流程图。
具体实施方式
下面结合附图和具体实施例,对本发明作进一步详细描述:
参照图1、实现本发明采用的估计***,包括微波光子阵列和一阶边带信号生成器。在空间等间隔放置m个级联的微波时延线,m=1...M,m≥2,分别定义为微波时延线1,微波时延线2,...,微波时延线m,...,微波时延线M,相邻微波时延线之间的距离为d,m个级联的微波时延线组成微波光子阵列,第m个微波时延线上级联有马赫增德尔调制器MZM-m、滤波器F-m和光功率计OPM-m。一阶边带信号生成器是将激光信号源、天线0和一个滤波器F-0并联到马赫增德尔调制器MZM-0上。一阶边带信号生成器中的滤波器F-0与微波光子阵列D中每一个微波时延线相连接,微波光子阵列D中的每一个光功率计分别级联1个功率比值器,m个光功率比值器并联到稀疏表示模型中。
参照图2、一种基于微波光子阵列的超宽带信号频率估计方法,包括如下步骤:
步骤1)获取一阶边带信号E0(t):
将天线0接收到的微波信号V0(t)和激光信号源输出的光载波信号E0输入到马赫增德尔调制器MZM-0中进行调制成为光电信号,得到多边带光电信号,并对多边带光电信号进行窄带滤波,滤除光载频以外的信号频率分量,得到一阶边带光电信号E0(t):
马赫增德尔调制器是将输入的光载波信号和微波信号调制成为光电信号的光电器件,将光信号和电信号调制成为光电信号,即微波光子。
步骤2)构建微波光子阵列D:
在空间等间隔放置m个微波时延线,m=1...M,m≥2,分别定义为微波时延线1,微波时延线2,...,微波时延线m,微波时延线M,相邻微波时延线之间的距离为d,构建由m个并联的微波时延线组成的微波光子阵列,第m个微波时延线上级联有马赫增德尔调制器MZM-m、滤波器F-m和光功率计OPM-m,得到微波光子阵元Dm,这m个微波光子阵元组成微波光子阵列D,将一阶边带光电信号E0(t)输入到这m个微波光子阵元组成微波光子阵列D中,可同时得到m个超宽带信号,对m个超宽带信号进行频率的估计,这样可以避免传统超宽带信号的频率估计中需要进角度搜索和角度匹配的问题,不用预先知道信号个数,提高了超宽带信号频率的估计效率,具体实施例中,m=50。
步骤3)获取微波光子阵列D输出的超宽带信号Em的光功率Pm:
步骤3a)将第m个微波时延线接收到的微波信号Vm(t)和MZM-0输出的一阶边带光电信号E0(t)输入到马赫增德尔调制器MZM-m中进行调制,并通过滤波器F-m对调制得到的第m个调制信号进行窄带滤波,得到第m个超宽带信号Em:
其中Vm为马赫增德尔调制器MZM-m的半波电压且V1=V2=...Vm=VM,β1k=β2k=...=βmk...=βMk,τ为相邻微波时延线接收到的微波信号Vm(t)的延迟时间,Ωk表示微波时延线接收到的微波信号的角频率,Ωkm表示m个微波时延线接收到的微波信号的角频率之和,Lm为马赫增德尔调制器MZM-m的***损耗且L1=L2=...Lm=LM,J1(·)表示一阶贝赛尔函数;
步骤3b)将第m个超宽带信号Em输入到光功率计OPM-m中,计算超宽带信号Em的光功率Pm:
(4)定义光功率比值向量Q:
(4a)计算Pm的相对光功率比值Qm:
(4b)定义M×1维的光功率比值向量:Q=[Q1,Q2,…,Qm…,QM]T,其中[·]T表示矩阵[·]的转置;
步骤5)定义M×N维矩阵ZD(w)m,k:
步骤5a)根据微波信号源的空域稀疏特性,采用空间网格划分方法,将可搜索的频率域等间隔划分成N份,得到超带微波信号Em的角频率的取值空间w为:w=[w1,w2,…,wk…,wN];
(5b)定义与阵列流型向量相关的M×N维的矩阵ZD(w)m.k:
步骤6)估计微波光子阵列D输出的超宽带信号Em的频率Ωm:
步骤6b)信号源在空域有稀疏特性,应用于超宽带信号的频率可以突破微波信号的瑞利限,进一步提高超宽带信号频率精度。通过改进的L1-SRACV算法模型,可得稀疏重构约束化方程为:
其中||·||1和||·||2分别表示l1范数和l2范数,s.t.表示约束关系,ε表示误差的允许值;
改进的L1-SRACV算法模型是利用阵列接收数据奇异值分解(SVD)得到的信号子空间构造稀疏表示模型,重构稀疏约束化方程。
采用凸优化方法求解稀疏重构约束化方程,凸优化是一种比较特殊的优化,是指目标函数和约束函数均为凸函数的优化问题,凸优化问题有一套非常完备的解决算法,在此采用现有针对凸优化问题的软件包CVX(Grant M,Boyd S.CVX:Matlab software fordisciplined convex programming[J].2008[Online]Available:http://stanfordedu/~boyd/cvx)来求解,通过该方法能够快速地得到稀疏表示系数向量的最优解。
Claims (2)
1.一种基于微波光子阵列的超宽带信号频率估计方法,其特征在于,包括如下步骤:
(1)获取一阶边带信号E0(t):
将天线0接收到的微波信号V0(t)和激光信号源输出的光载波信号E0输入到马赫增德尔调制器MZM-0中进行调制,得到多边带信号,并对多边带信号进行窄带滤波,得到一阶边带信号E0(t);
(2)构建微波光子阵列D:
构建由M个级联的微波时延线组成的微波光子阵列D,相邻微波时延线之间的距离为d,m=1...M,M≥2,第m个微波时延线上级联有马赫增德尔调制器MZM-m、滤波器F-m和光功率计OPM-m,得到微波光子阵元Dm,M个微波光子阵元组成微波光子阵列D;
(3)获取微波光子阵列D输出的超宽带信号Em的光功率Pm:
(3a)将第m个微波时延线接收到的微波时延信号Vm(t)和MZM-0输出的一阶边带信号E0(t)输入到马赫增德尔调制器MZM-m中进行调制,并通过滤波器F-m对调制得到的第m个调制信号进行窄带滤波,得到第m个超宽带信号Em:
其中p表示光载波信号E0的载波功率,K为调制系数,k=K=1,Vm为施加到马赫增德尔调制器MZM-m上的电压,且V1=V2=...Vm=VM,β1k=β2k=…=βmk...=βMk,τ为相邻微波时延线接收到的微波信号Vm(t)的延迟时间,Ωk表示微波时延线接收到的微波信号的角频率,Ωkm表示m个微波时延线接收到的微波信号的角频率之和,Lm为马赫增德尔调制器MZM-m的***损耗且L1=L2=…Lm=LM,J1(·)表示一阶贝赛尔函数,L0为MZM-0的***损耗,ω表示光载波信号E0的载波角频率;
(3b)将第m个超宽带信号Em输入到光功率计OPM-m中,计算超宽带信号Em的光功率Pm:
(4)定义光功率比值向量Q:
(4a)计算Pm的相对光功率比值Qm:
(4b)定义M×1维的光功率比值向量:Q=[Q1,Q2,...,Qm...,QM]T,其中[·]T表示矩阵[·]的转置;
(5)定义M×N维矩阵ZD(w)m,k:
(5a)根据微波信号源的空域稀疏特性,采用空间网格划分方法,将可搜索的空间频率域等间隔划分成N份,N>>M,得到超宽带信号Em的频率的取值空间w为:w=[w1,w2,...,wk...,wN];
(5b)定义与阵列流型向量相关的M×N维的矩阵ZD(w)m,k:
(6)估计微波光子阵列D输出的超宽带信号Em的频率Ωm:
(6b)通过改进的L1-SRACV算法模型,构建基于稀疏重构的约束优化方程式:
其中||·||1和||·||2分别表示l1范数和l2范数,s.t.表示约束关系,ε表示误差的允许值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810735983.1A CN108880648B (zh) | 2018-07-06 | 2018-07-06 | 基于微波光子阵列的超宽带信号频率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810735983.1A CN108880648B (zh) | 2018-07-06 | 2018-07-06 | 基于微波光子阵列的超宽带信号频率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108880648A CN108880648A (zh) | 2018-11-23 |
CN108880648B true CN108880648B (zh) | 2021-04-06 |
Family
ID=64299377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810735983.1A Active CN108880648B (zh) | 2018-07-06 | 2018-07-06 | 基于微波光子阵列的超宽带信号频率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108880648B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693932B (zh) * | 2020-06-08 | 2023-04-21 | 西安电子科技大学 | 基于微波光子的多微波信号波达方向角估计方法 |
CN111698036B (zh) * | 2020-06-08 | 2021-06-08 | 西安电子科技大学 | 基于微波光子的多微波信号频率估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326789A (zh) * | 2013-05-03 | 2013-09-25 | 华中科技大学 | 一种频率可调谐的微波相移***及方法 |
CN107707309A (zh) * | 2017-10-13 | 2018-02-16 | 南京航空航天大学 | 基于级联相位和偏振调制器的微波光子正交混频方法、装置 |
CN108155949A (zh) * | 2017-12-14 | 2018-06-12 | 华东师范大学 | 一种微波光子多频二进制相位编码信号产生装置和方法 |
-
2018
- 2018-07-06 CN CN201810735983.1A patent/CN108880648B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326789A (zh) * | 2013-05-03 | 2013-09-25 | 华中科技大学 | 一种频率可调谐的微波相移***及方法 |
CN107707309A (zh) * | 2017-10-13 | 2018-02-16 | 南京航空航天大学 | 基于级联相位和偏振调制器的微波光子正交混频方法、装置 |
CN108155949A (zh) * | 2017-12-14 | 2018-06-12 | 华东师范大学 | 一种微波光子多频二进制相位编码信号产生装置和方法 |
Non-Patent Citations (2)
Title |
---|
DOA estimation via sparse recovering from the smoothed covariance vector;Jingjing Cai; Dan Bao; Peng Li; Jingjing Cai; Dan Bao; Peng Li;《Journal of Systems Engineering and Electronics》;20160721;第27卷(第3期);第1-7页 * |
基于光子技术的微波频率测量研究进展;邹喜华,卢冰;《数据采集与处理》;20141115;第29卷(第6期);第1-10页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108880648A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109116294B (zh) | 基于微波光子阵列的超宽带信号波达方向角估计方法 | |
CN110061781B (zh) | 使用光频梳的射频光子光谱仪的方法和设备 | |
CN103116170B (zh) | Gnss的基于天线阵列干扰抑制模块的室内测试*** | |
CN106685871A (zh) | 一种iq调制器的控制方法及*** | |
CN109116297B (zh) | 一种被动雷达空间谱估计与合成波束的联合测向方法 | |
CN108880648B (zh) | 基于微波光子阵列的超宽带信号频率估计方法 | |
CN102981048A (zh) | 基于光采样的射频频率测量方法及装置 | |
CN109375201A (zh) | 微波光子雷达探测与测频一体化实现方法及装置 | |
CN111698037B (zh) | 基于微波光子的单微波信号波达方向角估计方法 | |
CN111698036B (zh) | 基于微波光子的多微波信号频率估计方法 | |
CN102007713B (zh) | 用于处理太拉赫波的方法和设备 | |
RU2495449C2 (ru) | Устройство формирования диаграммы направленности активной фазированной антенной решетки | |
Venkatesh et al. | Interferometric direction finding with a metamaterial detector | |
CN103326795B (zh) | 基于时间-频谱卷积原理的宽带射频信号相关检测方法 | |
Yixuan et al. | Parameter estimation of frequency hopping signal based on MWC–MSBL reconstruction | |
CN114696899B (zh) | 基于多频外差原理和光载微波干涉的距离测量方法 | |
Monari et al. | Aperture array for low frequency: The vivaldi solution | |
Sviatkina et al. | On possibility of using of spectral analysis for control the energy distribution of electromagnetic waves in radar channels | |
CN108009373B (zh) | 光控射频波束形成***的时空频多维建模方法 | |
CN111693932B (zh) | 基于微波光子的多微波信号波达方向角估计方法 | |
CN108605029A (zh) | 用于物理和信息信号导引的复合成方法和*** | |
Murugan et al. | BER and eye pattern analysis of 5G optical communication system with filters | |
Song et al. | The progress of terahertz wave source in communication | |
Yang et al. | Photonics-assisted wideband RF source localization method based on synthetic aperture interferometric detection | |
CN114122728B (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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Cai Jingjing Inventor after: Fei Xuming Inventor after: Zong Ru Inventor after: Su Ruolong Inventor after: Che Jinge Inventor before: Cai Jingjing Inventor before: Zong Ru Inventor before: Su Ruolong Inventor before: Che Jinge |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |