参数调制多载波雷达通信共享信号设计

杨慧婷 周宇 谷亚彬 张林让

引用本文:
Citation:

参数调制多载波雷达通信共享信号设计

    作者简介: 杨慧婷(1993–),女,硕士研究生,主要研究方向为雷达通信一体化。E-mail: yanght06@163.com;周 宇(1978–),男,副教授,主要研究方向为雷达通信一体化、STAP、雷达抗干扰技术。E-mail: zhouyu@mail.xidian.edu.cn;谷亚彬(1991–),男,博士研究生,主要研究方向为雷达通信一体化、雷达成像技术。E-mail: yabingu@163.com;张林让(1966–),男,教授,主要研究方向为阵列信号处理、网络化雷达协同抗干扰技术、雷达成像技术和雷达系统建模与仿真技术。E-mail: lrzhang@xidian.edu.cn.
    通讯作者: 周宇, zhouyu@mail.xidian.edu.cn
  • 基金项目:

    国家自然科学基金(61671361, 61731023, 61601343)

  • 中图分类号: TN95

Design of Integrated Radar and Communication Signal Based on Multicarrier Parameter Modulation Signal

    Corresponding author: ZHOU Yu, zhouyu@mail.xidian.edu.cn ;
  • Fund Project: The National Natural Science Foundation of China (61671361, 61731023, 61601343)

    CLC number: TN95

  • 摘要: 雷达通信一体化是减少电子平台体积与电磁干扰的一种有效途径,而共享信号的研究是实现雷达通信一体化的关键技术。该文提出了一种基于Chirp信号的多载波雷达通信共享信号,其主载波采用唯一Chirp信号实现雷达功能,副载波通过改变调频率和初始频率参数组合的Chirp信号调制通信信息。分析了共享信号的模糊函数以及参数设计方法,并对其处理过程及系统性能进行了研究。仿真结果表明,该信号具有较低误码率和高稳健性特性,使用该共享信号可在微量降低雷达性能的前提下实现通信数据的传输。
  • 图 1  64进制数据调制

    Figure 1.  64 system data modulation

    图 2  共享信号参数设计

    Figure 2.  Sharing signal parameters design

    图 3  不同阶次FRFT变换图

    Figure 3.  FRFT graphs of different orders

    图 4  共享信号实现框图

    Figure 4.  Block diagram of realizing sharing signal

    图 5  FRFT幅度峰值随阶数偏移变化

    Figure 5.  Change of FRFT amplitude peak with order deviation

    图 6  共享信号的模糊函数图

    Figure 6.  Ambiguity function of sharing signal

    图 7  主、副载波互相关函数

    Figure 7.  Cross correlation function of main carrier and subcarrier

    图 8  主副载波功率比9:1时FRFT输出

    Figure 8.  FRFT output at main and subcarrier power ratio of 9:1

    图 9  脉压结果

    Figure 9.  Pulse compression result

    图 10  动目标检测结果

    Figure 10.  MTD result

    图 11  检测概率与信噪比关系曲线

    Figure 11.  Curve of relationship between detection probability and signal-to-noise ratio

    图 12  误码率随信噪比变化曲线

    Figure 12.  Change of error rate with SNR

    表 1  主载波与不同主副功率比下共享信号的互相关系数

    Table 1.  Cross-correlation coefficient of shared signal under different power ratios

    主副功率比s13s24s35s46s57
    1:10.73490.74560.73410.73230.7672
    4:10.90820.90530.90220.90150.9115
    9:10.95220.95300.95190.95150.9554
    下载: 导出CSV
  • [1] 张明友. 雷达-电子战-通信一体化概论[M]. 北京: 国防工业出版社, 2010: 87–101.ZHANG Ming-you. The Conspectus of Integrated Radar-EW-Communication[M]. Beijing: National Defend Industry Press, 2010: 87–101.
    [2] TAVIK G C, HILTERBRICK C L, EVINS J B, et al. The advanced multifunction RF concept[J]. IEEE Transactions on Microwave Theory and Techniques, 2005, 53(3): 1009–1020. doi: 10.1109/TMTT.2005.843485
    [3] STURM C and WIESBECK W. Waveform design and signal processing aspects for fusion of wireless communications and radar sensing[J]. Proceedings of the IEEE, 2011, 99(7): 1236–1259. doi: 10.1109/JPROC.2011.2131110
    [4] XU C C and CHEN T Q. Conception of " signal sharing” in integrated radar and jammer system and the integrated signal design[C]. Proceedings of IEEE 2002 International Conference on Communications, Circuits and Systems and West Sino Expositions, Chengdu, China, 2002: 502–505. doi: 10.1109/ICCCAS.2002.1180668.
    [5] 邹广超, 刘以安, 吴少鹏, 等. 雷达-通信一体化系统设计[J]. 计算机仿真, 2011, 28(8): 1–4, 32. doi: 10.3969/j.issn.1006-9348.2011.08.001ZOU Guang-chao, LIU Yi-an, WU Shao-peng, et al. Design of radar-communications integrated system[J]. Computer Simulation, 2011, 28(8): 1–4, 32. doi: 10.3969/j.issn.1006-9348.2011.08.001
    [6] 李晓柏, 杨瑞娟, 程伟. 基于频率调制的多载波Chirp信号雷达通信一体化研究[J]. 电子与信息学报, 2013, 35(2): 406–412. doi: 10.3724/SP.J.1146.2012.00567LI Xiao-bai, YANG Rui-juan, and CHENG Wei. Integrated radar and communication based on multicarrier frequency modulation chirp signal[J]. Journal of Electronics &Information Technology, 2013, 35(2): 406–412. doi: 10.3724/SP.J.1146.2012.00567
    [7] SIT Y L, STURM C, REICHARDT L, et al. The OFDM joint radar-communication system: An overview[C]. Proceedings of the 3rd International Conference on Advances in Satellite and Space Communications, Budapest, 2011: 69–74.
    [8] CHEN X B, WANG X M, XU S F, et al. A novel radar waveform compatible with communication[C]. Proceedings of 2011 International Conference on Computational Problem-Solving, Chengdu, China, 2011: 177–181. doi: 10.1109/ICCPS.2011.6092272.
    [9] 李晓柏, 杨瑞娟, 程伟. 基于Chirp信号的雷达通信一体化研究[J]. 雷达科学与技术, 2012, 10(2): 180–186. doi: 10.3969/j.issn.1672-2337.2012.02.012LI Xiao-bai, YANG Rui-juan, and CHENG Wei. Integrated radar and communication based on chirp[J]. Radar Science and Technology, 2012, 10(2): 180–186. doi: 10.3969/j.issn.1672-2337.2012.02.012
    [10] TAKASE H and SHINRIKI M. A dual-use system for radar and communication with complete complementary codes[C]. Proceedings of the 11th International Radar Symposium, Vilnius, Lithuania, 2010: 1–4.
    [11] ZHAO Z L and JIANG D F. A novel integrated radar and communication waveform based on LFM signal[C]. Proceedings of the 5th International Conference on Electronics Information and Emergency Communication, Beijing, China, 2015: 219–223. doi: 10.1109/ICEIEC.2015.7284525.
    [12] HU F, CUI G L, YE W, et al. Integrated radar and communication system based on stepped frequency continuous waveform[C]. Proceedings of 2015 IEEE Radar Conference, Arlington, VA, USA, 2015: 1084–1087. doi: 10.1109/RADAR.2015.7131155.
    [13] LIU L, TIAN J J, and LIU Z H. A parameter estimation method for moving target in stepped frequency radar[C]. Proceedings of 2010 International Conference on Intelligent Computation Technology and Automation, Changsha, China, 2010: 289–292. doi: 10.1109/ICICTA.2010.170.
    [14] LIU T, CAO N, MAO M H, et al. Waveform design and optimization for OFDM radar signal based on fractional Fourier transform[C]. Proceedings of the 2nd International Conference on Systems and Informatics, Shanghai, China, 2014: 724–729. doi: 10.1109/ICSAI.2014.7009380.
    [15] 陶然, 齐林, 王越. 分数阶Fourier变换的原理与应用[M]. 北京: 清华大学出版社, 2004.TAO Ran, QI Lin, and WANG Yue. Theory and Applications of the Fractional Fourier Transform[M]. Beijing: Tsinghua University Press, 2004.
    [16] 李锋, 朱世华, 荣玫. FRFT-OFDM系统中FRFT域变换角度的估计[J]. 电子与信息学报, 2008, 30(10): 2427–2430.LI Feng, ZHU Shi-hua, and RONG Mei. Estimation of the angle parameter in FRFT domain for FRFT-OFDM systems[J]. Journal of Electronics &Information Technology, 2008, 30(10): 2427–2430.
    [17] 樊昌信, 曹丽娜. 通信原理[M]. 第6版, 北京: 国防工业出版社, 2008: 213–233.FAN Chang-xin and CAO Li-na. Communication Principle[M]. Sixth Edition, Beijing: National Defend Industry Press, 2008: 213–233.
  • [1] 陈芳香易伟周涛孔令讲 . 基于分数阶傅里叶变换的多辐射源被动直接定位算法. 雷达学报, 2018, 7(4): 523-530. doi: 10.12000/JR18027
    [2] 李自琦梅进杰胡登鹏沈序驰李晓柏 . 基于分组格雷编码的OFDM 雷达通信一体化系统峰均功率比抑制. 雷达学报, 2014, 3(5): 548-555. doi: 10.3724/SP.J.1300.2014.14059
    [3] 杨悦叶星炜张方正潘时龙 . 基于微波光子I/Q去斜接收的宽带线性调频雷达成像系统. 雷达学报, 2019, 8(2): 224-231. doi: 10.12000/JR19002
    [4] 黄翀鹏王剑徐保国 . 旁瓣抑制对线性调频脉冲移频干扰的影响. 雷达学报, 2012, 1(3): 270-276. doi: 10.3724/SP.J.1300.2012.20028
    [5] 邵启红万显荣张德磊赵志欣柯亨玉 . 基于OFDM 波形的短波通信与超视距雷达集成实验研究. 雷达学报, 2012, 1(4): 370-379. doi: 10.3724/SP.J.1300.2012.20089
    [6] 朱柯弘王杰梁兴东吴一戎 . 用于SAR与通信一体化系统的滤波器组多载波波形. 雷达学报, 2018, 7(5): 602-612. doi: 10.12000/JR18038
    [7] 杜强宋耀良曹晓健 . 基于Hermite 插值滤波器的直接延时补偿超宽带波束形成技术研究. 雷达学报, 2013, 2(3): 278-283. doi: 10.3724/SP.J.1300.2013.13028
    [8] 易建新万显荣赵志欣程丰柯亨玉 . 单频网CP-OFDM 信号外辐射源雷达的分载波杂波抑制方法 (英文). 雷达学报, 2013, 2(1): 1-13. doi: 10.3724/SP.J.1300.2013.13030
    [9] 路满宋红军罗运华 . 基于调频连续波信号的圆弧式合成孔径雷达成像方法. 雷达学报, 2016, 5(4): 425-433. doi: 10.12000/JR16007
    [10] 陈莹钟菲郭树旭 . 非合作跳频信号参数的盲压缩感知估计. 雷达学报, 2016, 5(5): 531-537. doi: 10.12000/JR15106
    [11] 曹继明李若明杨继尧孙强李王哲 . 基于去调频接收技术的微波光子双波段线性调频连续波雷达. 雷达学报, 2019, 8(2): 189-196. doi: 10.12000/JR18119
    [12] 苑霸姚萍郑天垚 . 基于权重归一化奇异值分解的辐射源信号识别研究. 雷达学报, 2019, 8(1): 44-53. doi: 10.12000/JR18053
    [13] 折小强仇晓兰韩冰雷斌 . 一种基于变换域的滑动聚束SAR 调频率估计方法. 雷达学报, 2014, 3(4): 419-427. doi: 10.3724/SP.J.1300.2014.14008
    [14] 庄旭昇汪玲高瑾池冰清 . 一种基于WiFi信号的运动目标无源雷达成像方法. 雷达学报, 2014, 3(6): 694-701. doi: 10.12000/JR14120
    [15] 霍凯赵晶晶 . 一种基于Bernoulli混沌的四相编码OFDM雷达信号设计方法. 雷达学报, 2016, 5(4): 361-372. doi: 10.12000/JR16050
    [16] 李廉林周小阳崔铁军 . 结构化信号处理理论和方法的研究进展. 雷达学报, 2015, 4(5): 491-502. doi: 10.12000/JR15111
    [17] 熊鹏柳征姜文利 . 复合调制雷达信号时差估计算法. 雷达学报, 2015, 4(4): 460-466. doi: 10.12000/JR15072
    [18] 于彬彬刘畅王岩飞 . 基于chirp-z 变换的斜视FMCW-SAR 非线性尺度变换算法. 雷达学报, 2012, 1(1): 69-75. doi: 10.3724/SP.J.1300.2012.20005
    [19] 邹鲲吴德伟张斌李伟 . 一种针对失配信号的可调自适应检测器. 雷达学报, 2015, 4(4): 411-417. doi: 10.12000/JR14129
    [20] 程丰曾庆平龚子平 . 高频外辐射源雷达一阶海杂波建模及仿真. 雷达学报, 2014, 3(6): 720-726. doi: 10.12000/JR14131
  • 加载中
图(12)表(1)
计量
  • 文章访问数:  218
  • HTML浏览量:  110
  • PDF下载量:  99
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-02
  • 录用日期:  2018-07-13
  • 网络出版日期:  2018-08-27
  • 刊出日期:  2019-02-28

参数调制多载波雷达通信共享信号设计

    通讯作者: 周宇, zhouyu@mail.xidian.edu.cn
    作者简介: 杨慧婷(1993–),女,硕士研究生,主要研究方向为雷达通信一体化。E-mail: yanght06@163.com;周 宇(1978–),男,副教授,主要研究方向为雷达通信一体化、STAP、雷达抗干扰技术。E-mail: zhouyu@mail.xidian.edu.cn;谷亚彬(1991–),男,博士研究生,主要研究方向为雷达通信一体化、雷达成像技术。E-mail: yabingu@163.com;张林让(1966–),男,教授,主要研究方向为阵列信号处理、网络化雷达协同抗干扰技术、雷达成像技术和雷达系统建模与仿真技术。E-mail: lrzhang@xidian.edu.cn
  • 西安电子科技大学雷达信号处理国家重点实验室   西安   710071
基金项目:  国家自然科学基金(61671361, 61731023, 61601343)

摘要: 雷达通信一体化是减少电子平台体积与电磁干扰的一种有效途径,而共享信号的研究是实现雷达通信一体化的关键技术。该文提出了一种基于Chirp信号的多载波雷达通信共享信号,其主载波采用唯一Chirp信号实现雷达功能,副载波通过改变调频率和初始频率参数组合的Chirp信号调制通信信息。分析了共享信号的模糊函数以及参数设计方法,并对其处理过程及系统性能进行了研究。仿真结果表明,该信号具有较低误码率和高稳健性特性,使用该共享信号可在微量降低雷达性能的前提下实现通信数据的传输。

English Abstract

    • 随着科学技术的发展,现代战争中的装备必须趋向于综合化发展,但同时也不能让过多的设备加剧恶化平台周围的电磁环境、增加负荷,例如无人机,就需在小体积平台上集成多种装备功能,并保持平台的机动性和综合性,雷达和通信系统是平台广泛配备的两种电子系统[1],若能实现雷达通信一体化[2,3],将大大提高电子系统的综合利用率。

      雷达通信一体化的理念在20世纪60年代出现以后,对其研究主要分为分时、分波束和同时3种体制。分时体制在通信时不能兼顾雷达探测,即在通信时存在雷达探测盲区,但相对最易实现,故研究的较多;分波束体制将相控阵面划分为不同区域,利用划分的各个阵面实现不同功能;同时体制将雷达信号和通信信号融合在一起,在同一平台同时实现探测与传输功能。其一体化程度最高,是未来雷达通信一体化的发展方向。这种体制的关键技术更多地集中在共享信号设计,而共享信号设计主要需要解决通信数据传输和雷达探测之间的关系。现有的共享信号设计方法基本可分为3类:①雷达与通信信号各自独立产生后叠加[4],②基于通信信号,将其改造成雷达探测波,③基于雷达信号,在其上调制通信数据[5]。文献[6]中研究了利用线性调频信号(Chirp信号)上、下扫频分别作为雷达波形和通信波形,叠加产生共享信号,接收时利用正交性将其分离的一体化系统,但其通信速率受到很大的限制,雷达性能有所降低;文献[7]中研究了利用正交频分复用(Orthogonal Frequency Division Multiplexing, OFDM)信号实现一体化波形,但是OFDM信号不是恒包络,峰均比较高不利于在雷达的C类放大器中放大,且对多普勒频移较敏感,仅适用于短距离通信与探测;文献[8]中研究了将通信最小频移键控(Minimum Shift Keying, MSK)调制到线性调频信号上实现一体化波形,该波形能在实现雷达检测动能的同时完成通信功能;文献[9]研究了通过键控Chirp信号的初始频率来调制通信数据从而实现雷达通信一体化的方法,但是雷达检测处理的匹配滤波器要随着发射信号的改变而改变。

      本文提出一种基于Chirp信号参数调制的多载波雷达通信共享信号,主载波用于雷达检测功能,副载波的调频率与初始频率参数可选,从而携带数据实现通信信息调制。在设计共享信号时,通信数据的随机性常使不同脉冲间的信号相关性减弱,而雷达探测为进行相干积累,需在接收端使用与之对应的匹配滤波器,大大增加雷达系统负担,本文所设计信号利用主载波的确定性提高脉冲的相关性,雷达处理系统不需要增加额外单元,采用同原始雷达相同的处理流程;而不同起始频率、不同调频率的Chirp信号能在带宽利用率及正交性之间提供平衡。在文中对所设计共享波形的模糊函数、主副载波之间的正交性等性能进行了分析;在接收端通过分数傅里叶变换,根据检测点的能量聚集位置进行解调。

    • 共享信号设计中,在雷达探测波形上调制通信信息后,由于通信数据的随机性,使脉冲波形产生差异性,需要增加额外的雷达信号处理单元,造成负担。为减少脉冲差异性,便于雷达目标检测处理,设计主副载波的共享信号形式,主载波作雷达目标检测,副载波调制通信信息[10,11]

      副载波由待传输码元从一组Chirp信号组$\left\{ {{s_{kl}}\left( t \right)} \right\}$中选取,表达式为

      $ \begin{aligned} {s_{kl}}\left( t \right) =&{A_2}\exp\left[ {{\rm{j}}\left( {{{π}}{\mu _k}{t^2} + 2{{π}}{f\!_l}t} \right)} \right], \;\; \\ {\rm {}}&k = 0,1, ·\!·\!· ,{N_1} - 1; \; l = 0,1, ·\!·\!· ,{N_2} - 1 \end{aligned} $

      式中,等间隔调频率${\mu _k} = {\mu _{\rm{0}}} + k{\rm{\Delta }}\mu $;等间隔初始频率${f_l} = {f_{\rm{0}}} + l\Delta f$$n = {n_1} + {n_2}$位二进制数据中,${n_1}$位数据映射${N_1}$个调频率,${n_2}$位数据映射${N_2}$个初始频率,由通信数据键控映射得调频率为${\mu _k}$、初始频率为${f_l}$的Chirp信号,单个Chirp信号可携带n bit的数据[12]。以${n_1} = 3$, ${n_2} = 3$的8 bit调制为例,示意图如图1所示,$\mu $$f$两方向间没有约束关系。

      图  1  64进制数据调制

      Figure 1.  64 system data modulation

      主载波确定为调频率大于副载波调频率选取范围的Chirp信号,即${\mu _{\rm{r}}} > {\mu _k}$,带宽覆盖副载波的可用带宽,使得共享信号的带宽始终保持不变,

      ${s_{\rm{r}}}\left( t \right) = {A_1}\exp \left[ {{\rm{j}}\left( {{{π}}{\mu _{\rm{r}}}{t^2} + 2{{π}}{f\!_{\rm{r}}}t} \right)} \right]$

      利用调频率的多样性,给主、副载波提供良好的准正交性[13]。参数选取范围如图2所示,主载波为确定的Chirp信号,副载波为众多参数组合中选取的某一Chirp信号,共享信号表示为

      图  2  共享信号参数设计

      Figure 2.  Sharing signal parameters design

      $s\left( t \right) = {s_{\rm{r}}}\left( t \right) + {s_{kl}}\left( t \right), \quad - {\tau / 2} \le t \le {\tau / 2}$

    • 共享信号的通信码元信息利用分数阶傅里叶变换(FRactional Fourier Transform, FRFT)解调,Chirp信号的FRFT变换为

      $ \begin{align} {S_\alpha }\left( u \right) =& A\sqrt {1 - {\rm{j}}\cot \alpha } \exp\left( {{\rm{j{{π}} }}{u^2}\cot \alpha } \right) \\ {\rm{}}& \cdot \int_{ - {\tau / 2}}^{{\tau / 2}} {{{\rm{e}}^{{\rm{j2{{π}} }}\left( {\mu + \cot \alpha } \right){t^2} + {\rm{j{{π}} }}\left( {{f_0} - u\csc \alpha } \right)t}}} {\rm{d}}t \end{align} $

      式中,旋转角度$\alpha = {{p{{π}}} / 2}$; $p$为变换阶次。FRFT变换与Chirp信号的调频率、初始频率参数间的关系满足

      $ \left\{\!\!\! \begin{array}{l} \mu = - \cot \alpha \\ f = u\csc \alpha \\ \end{array} \right. \Rightarrow \left\{\!\!\! \begin{array}{l} p = {{ - 2{\rm{arccot}}\mu } / {{π}}} \\ u = - f\sin \left( {{\rm{arccot}}\mu } \right) \\ \end{array} \right. $

      时,出现幅度峰值$\left| {{S_\alpha }\left( u \right)} \right| = {{A\tau } / {\sqrt {\left| {\sin \alpha } \right|} }}$。即Chirp信号的调频率$\mu $决定了FRFT唯一的最优变换阶次${p_{\rm{m}}}$,载频${f_0}$决定了在${p_{\rm{m}}}$阶分数阶傅里叶域上最优的能量聚集位置[14]

      图3所示,只有在最优阶次的FRFT变换时,在分数阶傅里叶域上才有峰值输出。接收端对接收到的共享信号进行FRFT变换,检测峰值,得到峰值所在的变换阶次与分数阶傅里叶域坐标,由FRFT变换与Chirp信号参数之间的关系,解得副载波的调频率与初始频率,从而映射出调制的码元数据[15]

      图  3  不同阶次FRFT变换图

      Figure 3.  FRFT graphs of different orders

    • 首先根据设定的映射规则,将通信数据映射到对应的初始频率、调频率的Chirp信号序列中,在接收端通过分数阶傅里叶变换进行解调。一体化框图如图4所示,在发送端,将通信数据串并转换后分成$n$位一组,前${n_1}$位键控得调频率,后${n_2}$位键控得初始频率,产生一个特定的Chirp信号,与主载波确定Chirp信号叠加组成共享信号,送入高斯白噪声信道。

      图  4  共享信号实现框图

      Figure 4.  Block diagram of realizing sharing signal

      在接收端,雷达处理系统与常规雷达相同,不会增加额外处理单元。通信处理系统,对回波进行分数阶傅里叶变换,依次进行${p_0} \sim {p_{{N_1} - 1}}$阶FRFT变换,将变换后FRFT域${u_0} \sim {u_{{N_2} - 1}}$位置处的采样点设为检测点,对其数据进行门限判定,得到超过阈值的检测点的阶次${p_k}$和位置${u_l}$,并根据映射关系解调出通信码元数据。

    • 为满足通信解调准确率与雷达目标检测分辨率,需对共享信号各参数进行设计。

    • 副载波$\left\{ {{\mu _k}} \right\}$直接影响数据传输速率及解调性能;固定信号脉宽$\tau $, $\mu $越大Chirp信号占用的带宽越大,$\left\{ {{\mu _k}} \right\}$的最大取值被带宽所限制;$\left\{ {{\mu _k}} \right\}$对应${n_1}$位二进制通信信息,初始频率${f_{\rm r}}$相同时相邻的${\mu _k}$${\mu _{k + 1}}$对应的两个最优分数阶傅里叶变换阶次${p_k}$${p_{k + 1}}$的间隔$\Delta p$决定了采用分数阶傅里叶变换进行解调时对相邻符号的区分度,即${\rm{\Delta }}p$决定了相邻符号间的干扰程度。由式(5)推导可得$\Delta p$${\mu _k}$, ${\mu _{k + 1}}$具有以下关系:

      ${\rm{\Delta }}p = \left| {{p_k} - {p_{k + 1}}} \right| = {{2\left| {{\rm{arccot}}{\mu _k} - {\rm{arccot}}{\mu _{k + 1}}} \right|} / {{π}}}$

      对Chirp信号进行非最优阶FRFT时,分数傅里叶域谱不具有聚集性质,而且随着变换阶次偏离$p$的程度$\Delta p$增大,FRFT的峰值明显下降;同时调频率越大,随$\Delta p$增加下降得更快[16]图5仿真了不同$\Delta p$下归一化峰值的变化趋势,归一化调频率分别为0.1, 0.3, 0.5, 0.7, 0.9,脉宽相同,从图中看出,随$\Delta p$增大峰值幅度降低;要使得FRFT解调输出能唯一确定峰值点对应参数,可根据所需要求设定峰值幅度门限,若设定FRFT峰值幅度门限为–10 dB,则$\left\{ {{\mu _k}} \right\}$的设计要使得$\Delta p$至少为0.03。

      图  5  FRFT幅度峰值随阶数偏移变化

      Figure 5.  Change of FRFT amplitude peak with order deviation

    • 为使具有相同调频率${\mu _k}$、不同初始频率${f_l}$的Chirp 信号在分数傅里叶域上可以区分,相邻的${f\!_l}$${f\!_{l + 1}}$对应的分数阶傅里叶变换域上的两个谱峰位置${u_l}$${u_{l + 1}}$的间隔$\Delta u$决定了采用FRFT进行解调时对${n_2}$位数据中相邻符号的区分度,即$\Delta u$决定了相邻符号间的干扰程度。由式(5)推导可得$\Delta u$${f\!_l}$, ${f\!_{l + 1}}$存在如下关系:

      $ \Delta u = \left| {{u_l} - {u_{l + 1}}} \right| = \Bigr| {\left( {{f_l} - {f_{l + 1}}} \right) \cdot \sin \alpha } \Bigr| \\ \;\;\;\;\;\;= \left| {\Delta f \cdot \sin \alpha } \right| $

      对Chirp信号进行最优阶FRFT变换时,需要在$u$轴上能根据峰值幅度区分出${f\!_l}$${f\!_{l + 1}}$,由Chirp信号的分数傅里叶域的幅度谱

      $\Big| {{S_\alpha }\left( u \right)} \Big| = \left| {A\tau \sqrt {1 - {\rm{j}}\cot \alpha } \cdot } \right. {{\rm{Sa}}\left[ {{{π}} \left( {{f\!_0} - u\csc \alpha } \right)\tau } \right]} \Bigr|$

      (8)

      第1零点间距离为$\left| {{{2\sin \alpha } / \tau }} \right|$, $\left\{ {{f\!_l}} \right\}$的设计应使得$\Delta u$满足

      $\Delta u = \left| {\Delta f \cdot \sin \alpha } \right| > \left| {{{2\sin \alpha } / \tau }} \right|$

      $\Delta f > {2 / \tau }$

    • 模糊函数表征了波形的距离与多普勒分辨率等特性。共享信号$s\left( t \right) = {s_{\rm r}}\left( t \right) + {s_{kl}}\left( t \right)$, $ - {{\tau '} / 2} \le t \le$$ {{\tau '} / 2}$的模糊函数为

      $ \begin{align} \chi \left( {\tau ,{f\!_{\rm{d}}}} \right) =& \int_{ - \infty }^\infty {s\left( t \right){s^*}\left( {t - \tau } \right){{\rm e}^{{\rm{j2{{π}} }}{f\!_{\rm{d}}}t}}{\rm{d}}t} \\ =&\int_{{t_1}}^{{t_2}} {\underbrace {\Bigr[ {{s_{\rm{r}}}\left( t \right)s_{\rm{r}}^*\left( {t - \tau } \right) + {s_{kl}}\left( t \right)s_{kl}^*\left( {t - \tau } \right)} }_{{\chi _{_{\rm M}}}}} \\ {\rm {}} & + \underbrace { {{s_{\rm{r}}}\left( t \right)s_{kl}^*\left( {t - \tau } \right) + {s_{kl}}\left( t \right)s_{\rm{r}}^*\left( {t - \tau } \right)} \Bigr]}_{{\chi _{_{\rm I}}}} \\ {\rm {}} & {\cdot {\rm e}^{{\rm{j2{{π}} }}{f_{\rm{d}}}t}}{\rm d}t \end{align} $

      $\left\{ \begin{array}{l} \!\!\!0 \le \tau \le \tau '{{时,}}\ {t_1} = - \tau '/2 + \tau ,{t_2} = \tau '/2\\ \!\!\!- \tau ' \le \tau \le 0{{时,}}\ {t_1} = - \tau '/2,{t_2} = \tau '/2 + \tau \end{array} \right.$

      由表达式可看出,可将共享信号的模糊函数分为主瓣区域${\chi _{\rm_{M}}}$与邻道干扰项${\chi _{\rm_{I}}}$。主瓣区域${\chi _{_\rm{M}}}$为主、副载波的自模糊函数之和,表示为

      $ \begin{align} {\chi _{_{\rm{M}}}} =& A_1^2{{\rm e}^{{\rm{j{{π}} }}\left( {2{f\!_{\rm{r}}} + {f\!_{\rm{d}}}} \right)\tau }}\frac{{\sin \Bigr[ {{{π}}\left( {{\mu _{\rm{r}}}\tau + {f\!_{\rm{d}}}} \right)\bigr( {\tau ' - \left| \tau \right|} \bigr)} \Bigr]}}{{{{π}}\bigr( {{\mu _{\rm{r}}}\tau + {f\!_{\rm{d}}}} \bigr)}} \\ {\rm {}}& + A_2^2{{\rm e}^{{\rm{j{{π}} }}\left( {2{f\!_l} + {f\!_{\rm{d}}}} \right)\tau }}\frac{{\sin \Bigr[ {{{π}}\left( {{\mu _k}\tau + {f\!_{\rm{d}}}} \right)\bigr( {\tau ' - \left| \tau \right|} \bigr)} \Bigr]}}{{{{π}}\left( {{\mu _k}\tau + {f\!_{\rm{d}}}} \right)}} \end{align} $

      而邻道干扰项${\chi _{_{\rm{I}}}}$由主、副载波间的互模糊函数之和表示,是应该尽量抑制的部分。

      $ \begin{align} {\chi _{\rm_{I}}} =& {A_1}{A_2}\sum\limits_{\displaystyle i=1 ,\ s = 1\atop {{s \ne i}}}^2 {\exp \left( {{\rm j}2{{π}} {f\!_i}\tau - {\rm j}{{π}} {\mu _i}{\tau ^2}} \right)} \\ {\rm {}}& \cdot \int_{{t_1}}^{{t_2}} {\exp \left[ { - {\rm j}{{π}} \left( {{\mu _i} - {\mu _s}} \right){t^2}} \right]} \\ {\rm {}}& \cdot \exp \Bigr[ { - {\rm j}2{{π}} \left( {{f\!_i} - {\mu _i}\tau - {f\!_s} - {f\!_{\rm{d}}}} \right)t} \Bigr]{\rm d}t\\ =& {A_1}{A_2}\sum\limits_{\displaystyle i=1 , \ s = 1\atop {{s \ne i}}}^2 {\exp \left( \alpha \right)} \cdot \int_{{t_1}}^{{t_2}} {\exp \left( { - {\beta ^2}} \right){\rm d}t} \\ =& {A_1}{A_2}\sum\limits_{\displaystyle i=1, \ s = 1\atop {{s \ne i}}}^2 {\frac{1}{{2\sqrt {{\rm j}\left( {{\mu _i} - {\mu _s}} \right)} }}\exp \left( \alpha \right)} \\ {\rm {}}& \cdot \Bigr\{ {{\rm{erf}}\left[ {\beta \left( {{t_2}} \right)} \right] - {\rm{erf}}\left[ {\beta \left( {{t_1}} \right)}\right]} \Bigr\} \end{align} $

      式中,${f\!_1},{\mu _1}$为主载波参数,${f_2},{\mu _2}$为副载波参数,

      $\alpha = \displaystyle\frac{{{\rm j}{{π}} {{\left( {{f\!_i} - {\mu _i}\tau - {f\!_s} - {f\!_{\rm{d}}}} \right)}^2}}}{{{\mu _s} - {\mu _i}}} + {\rm j}{\rm{2}}{{π}} {f\!_i}\tau - {\rm j}{{π}} {\mu _i}{\tau ^2}$,

      $\beta = \sqrt {{\rm{{\rm j}{{π}} }}} \left[ {\displaystyle\frac{{\left( {{f\!_i} - {\mu _i}\tau - {f\!_s} - {f\!_{\rm{d}}}} \right)}}{{\sqrt {{\mu _i} - {\mu _s}} }}+ } {\sqrt {{\mu _i} - {\mu _s}} t} \right]$, ${\rm{erf}}\left( x \right) = \displaystyle\frac{2}{{\sqrt {{π}} }}\int_0^x {\exp \left( { - {z^2}} \right)} {\rm d}z$

      由于无法求得模糊函数的具体表达式,故对其模糊函数进行了统计意义上的仿真分析,仿真参数同第5节,主瓣区域为两Chirp信号模糊函数叠加,而邻道干扰性是由两项主、副载波的互模糊函数之和,其幅度相较于信号模糊函数的峰值,幅度较低,由多次仿真得到,邻道干扰项幅度峰值的平均值仅为模糊函数峰值的2.2%,方差为0.000115,且峰值不位于速度-距离平面原点,故认为模糊函数主要由主瓣区域决定。随调制数据的改变,多普勒切片的主瓣宽度变化范围均低于${1 / \tau }$,第1旁瓣峰值随调制数据在–13 dB上下变化;时延切片的主瓣宽度在${1 / B}$上下变化,第1旁瓣峰值在–12 dB上下变化;故当利用发射的共享信号进行匹配滤波时,性能将有所下降。图6仿真了调制某一数据的共享信号的主瓣区域模糊函数图和邻道干扰项的模糊函数图,主瓣位于速度-距离平面原点处的峰值幅度最高。

      图  6  共享信号的模糊函数图

      Figure 6.  Ambiguity function of sharing signal

    • 主副载波间的互相关性决定了在接收端进行匹配滤波时副载波的剩余量。主载波与副载波信号表达式如下:

      ${s_{\rm{r}}}\left( t \right) = {A_1}\exp \left( {{\rm j}{{π}} {\mu _{\rm{r}}}{t^2} + {\rm j}2{{π}} {f\!_{\rm{r}}}t} \right)$

      ${s_{kl}}\left( t \right) = {A_2}\exp \left( {{\rm j}{{π}} {\mu _k}{t^2} + {\rm j}2{{π}} {f\!_l}t} \right)$

      其中,$ - {{\tau '} / 2} \le t \le {{\tau '} / 2}$, $\tau '$为信号脉冲宽度。相关函数为

      $ \begin{aligned} {\rm{}}&{R_{{s_{\rm{r}}}, {s_{kl}}}}\left( \tau \right) \\ {\rm{}}&\quad = \int_{ - \infty }^{ + \infty } {s_{\rm{r}}^*\left( t \right)} {s_{kl}}\left( {t + \tau } \right){\rm d}t \\ {\rm{}} &\quad ={A_1}{A_2}\!\!\int_{{t_1}}^{{t_2}} {\!\!\!\!\!\exp \left\{ {{\rm{j2{{π}} }}{f\!_l}\tau - {\rm{j{{π}} }}\frac{{{{\left( {{f\!_l} + {\mu _k}\tau - {f\!_{\rm{r}}}} \right)}^2}}}{{{\mu _k} - {\mu _{\rm{r}}}}} } \right.} \\ {\rm{}}&\quad \quad +{\rm{j\!\!}}\frac{{{π}}}{{\rm{2}}}\!{{\left( {\frac{{2\left( {{f\!_l} \!+\! {\mu _k}\tau \!-\! {f_{\rm{r}}}} \right)}}{{\sqrt {2\left( {{\mu _k} \!-\! {\mu _{\rm{r}}}} \right)} }} \!+\! \sqrt {2\left( {{\mu _k} \!-\! {\mu _{\rm{r}}}} \right)} \!\cdot \!t\!} \right)}^2} \\ {\rm{}}&\quad \quad +{\rm{j{{π}} }}{\mu _k}{\tau ^2} \Biggr\}{\rm d}t \end{aligned} $

      式中,积分区间$\left[ {{t_1},{t_2}} \right]$取决于$\tau $,当$0 \le \tau \le \tau '$时,${t_1} = \, -\, {{\tau '} / 2, }$ ${t_2} \,=\, {{\tau '} / \,2} \,-\, \tau\,; $$ \;-\; \tau ' \;\le \;\tau \;\le 0$时,${t_1} = - {{\tau '} / 2} - \tau $, ${\;\;t_2} = {{\tau '} / 2\;\;}$。令$\gamma \left( t \right) = \sqrt {2\left( {{\mu _k} - {\mu _{\rm{r}}}} \right)} $ $ \cdot t +\displaystyle\frac{{2\left( {{f_l} + {\mu _k}\tau - {f_{\rm{r}}}} \right)}}{{\sqrt {2\left( {{\mu _k} - {\mu _{\rm{r}}}} \right)} }}$,得${\rm d}\gamma = \sqrt {2\left( {{\mu _k} - {\mu _{\rm{r}}}} \right)} \cdot {\rm d}t$,则式(16)变为

      $ \begin{align} {R_{{s_{\rm{r}}},\; {s_{kl}}}}\left( \tau \right) =& {A_1}{A_2}\frac{{\exp \left( {{\rm{j2{{π}} }}{f\!_l}\;\tau + {\rm{j{{π}} }}{\mu _k}{\tau ^2}} \right)}}{{\sqrt {2\left( {{\mu _k} - {\mu _{\rm{r}}}} \right)} }} \\ {\rm {}}&\cdot \exp \left[ { - {\rm{j{{π}} }}\frac{{{{\left( {{f\!_l} + {\mu _k}\tau - {f_{\rm{r}}}} \right)}^2}}}{{{\mu _k} - {\mu _{\rm{r}}}}}} \right] \\ {\rm {}}& \cdot \int_{\gamma \left( {{t_1}} \right)}^{\gamma \left( {{t_2}} \right)} {\exp \left[ {{\rm j}\frac{{{π}} }{2}{\gamma ^2}\left( t \right)} \right]} {\rm d}\gamma \end{align} $

      其 中,$ \displaystyle\int_{\gamma \left( {{t_1}} \right)}^{\gamma \left( {{t_2}} \right)} {\exp\!\left[ {{\rm{j}}\frac{{{π}}}{{\rm{2}}}{\gamma ^2}\left( t \right)} \right]} {\rm d}\gamma = C\Bigr( {\gamma \left( {{t_2}} \right)} \Bigr) +{\rm{j}}S \Bigr( {\gamma \left( {{t_2}} \right)}\Bigr)$         $ \quad\quad -C\Bigr( {\gamma \left( {{t_1}} \right)} \Bigr) - {\rm{j}}S\Bigr( {\gamma \left( {{t_1}} \right)} \Bigr), $ $\displaystyle C\left( \gamma \right) = \int_0^\gamma {\cos \left( {\displaystyle\frac{{{{π}}{v^2}}}{2}} \right)} {\rm d}v $$ \approx \displaystyle\frac{1}{2} + \frac{1}{{{{π}}\gamma }}\sin \left( {\displaystyle\frac{{{π}}}{2}{\gamma ^2}} \right), $$\displaystyle S\left( \gamma \right) = \int_0^\gamma {\sin \left( {\displaystyle\frac{{{{π}}{v^2}}}{2}} \right)} {\rm d}v$$ \approx \displaystyle\frac{1}{2} - \displaystyle\frac{1}{{{{π}}\gamma }}\sin \left( {\displaystyle\frac{{{π}}}{2}{\gamma ^2}} \right)$

      互相关值取决于主副载波的调频率差值以及频率差值与调频率差值的比值。根据前述参数设计主载波${\mu _{\rm{r}}} = 120 \; {{{\rm{MHz}}} / {{\rm{{{μ}} s}}}}$, ${f_{\rm{r}}} = 0 \; {\rm{ MHz}}$,副载波${\mu _{{k}}} = 99 \; {{{\rm{MHz}}} / {{\rm{{{μ}} s}}}}$, ${f_l} = 21 \; {\rm{ MHz}}$,对载波自相关特性及互相关特性进行了仿真,由图7可看出,主副载波的互相关值比主载波的自相关值低35 dB,表明在接收端进行匹配滤波时,副载波剩余量很小。模糊函数表征利用发射信号进行匹配滤波的输出,性能有明显下降,但本文中,主载波保持不变,仅利用主载波进行匹配滤波时,探测性能的降低量将变小。

      图  7  主、副载波互相关函数

      Figure 7.  Cross correlation function of main carrier and subcarrier

    • 共享信号由主载波与副载波叠加得到,则用于雷达探测的功率会有所下降,但主、副载波功能相互独立,故可调整主载波与副载波的不同功率配比,增加用于雷达探测的功率。

      雷达探测目标由以主载波为参考信号的匹配滤波器进行脉压处理,处理结果基于主载波与各脉冲共享信号的相关性,相关性表示为主载波的自相关函数与主、副载波的互相关函数之和:

      $ \begin{align} {R_{s,\; {s_{\rm{r}}}}}\left( \tau \right) =& \int_{ - \infty }^{ + \infty } {{s^*}\left( t \right)} {s_{\rm{r}}}\left( {t + \tau } \right){\rm d}t \\ =& \int_{ - \infty }^{ + \infty } {\Bigr[ {s_{\rm{r}}^*\left( t \right) + s_{kl}^*\left( t \right)} \Bigr]} \cdot {s_{\rm{r}}}\left( {t + \tau } \right){\rm d}t \\ =& {R_{{s_{\rm{r}}} ,\; {s_{\rm{r}}}}}\left( \tau \right) + {R_{{s_{kl}},\; {s_{\rm{r}}}}}\left( \tau \right) \end{align} $

      由4.2节分析可知,主、副载波的互相关函数相较于主载波的自相关函数幅度很低,故雷达探测结果受副载波分量影响很小。

      主副载波的功率分配决定了用于雷达探测的功率,可在适当范围内提高主载波的功率以用于雷达探测,表1列出了在不同主、副载波功率比时,主载波与几组不同参数下的共享信号之间的互相关系数,从表1中可以看出,主载波所占功率越大,相关系数越接近于1,雷达探测性能越好。

      主副功率比s13s24s35s46s57
      1:10.73490.74560.73410.73230.7672
      4:10.90820.90530.90220.90150.9115
      9:10.95220.95300.95190.95150.9554

      表 1  主载波与不同主副功率比下共享信号的互相关系数

      Table 1.  Cross-correlation coefficient of shared signal under different power ratios

      增加了主载波的功率后,副载波的功率必然会下降,利用FRFT变换的解调性能会有所下降,如图8所示,在主、副载波的功率比为9:1时,FRFT变换的旁瓣在略微升高后,依然能保持在–10 dB左右,能检测到明显峰值,解调出码元数据,但是主副功率比不能无限制增大,主载波功率过高时,在FRFT解调输出谱中会覆盖掉峰值,无法解调出数据,功率比越高,误码率越差,可根据应用条件选择主副载波的功率比。

      图  8  主副载波功率比9:1时FRFT输出

      Figure 8.  FRFT output at main and subcarrier power ratio of 9:1

    • 雷达与通信接收端之间的相对运动会存在多普勒频移${f\!_{\rm{d}}}$,此时接收Chirp信号形式为

      $ \begin{align} s\left( t \right) =& A\exp \Bigr[ {{\rm{j}}\!\left({2{{π}}f't + {{π}}\mu {t^2} + \varphi } \right)} \Bigr],\\ {\rm{}}&\;\;\;t \in \left[ { - {\tau / 2},{\tau / 2}} \right],\ f' = {f\!_{\rm{c}}} + {f\!_{\rm{d}}} \end{align} $

      多普勒频移${f\!_{\rm{d}}}$可看作是初始频率${f\!_{\rm{c}}}$偏移,调频率$\mu $不受影响,对应的是分数阶傅里叶域上的峰值位置,而峰值幅度不变,则解调器输出峰值位置偏移量和检测点幅度平方输出分别为

      $ {\rm{\Delta }}u = {f\!_{\rm{d}}}\sin \Bigr( {{\rm{arccot}}\left( { - \mu } \right)} \Bigr) $

      ${\left| {{S_\alpha }\left( u \right)} \right|^2} = \frac{{{A^2}}}{{\sin \alpha }} \cdot \frac{{{{\sin }^2}\left( {{{π}}{f\!_{\rm{d}}}\tau } \right)}}{{{{\left( {{{π}}{f\!_{\rm{d}}}} \right)}^2}}}$

      多普勒频移${f\!_{\rm{d}}}$带来的幅度平方衰减系数为

      $\gamma = {{\frac{{{A^2}{{\sin }^2}\left( {{{π}}{f\!_{\rm{d}}}\tau } \right)}}{{{{\left( {{{π}}{f\!_{\rm{d}}}} \right)}^2}\sin \alpha }}} \biggr/ {\frac{{{A^2}{\tau ^2}}}{{\sin \alpha }} = }}\frac{{{{\sin }^2}\left( {{{π}}{f\!_{\rm{d}}}\tau } \right)}}{{{{\left( {{{π}}{f\!_{\rm{d}}}\tau } \right)}^2}}}$

      由式(22)可看出衰减系数只与多普勒频移${f\!_{\rm{d}}}$和信号脉宽$\tau $有关,${{π}}{f\!_{\rm{d}}}\tau $达到0.5时幅度衰减尚不到0.1,可知多普勒频移对检测点幅值影响较小,说明本文所设计的共享信号对多普勒是稳健的。

    • 本文调制方式需要考虑调频率与频率的配比,设比特宽度为${T_{\rm{b}}}$,仅对调频率键控时,$M = {2^n}$进制符号宽度为${T_{\rm{s}}} = {T_{\rm{b}}}n$,由于调频率${\mu _k}$是变化的,不同符号键控输出的Chirp信号的扫频带宽${\mu _k}{T_{\rm{s}}}$也是变化的,其最大宽度由最大的${\mu _{k\max }}$决定。则调频率键控方式的带宽效率${\eta _{{μ}}}$满足:

      ${\eta _{\rm{{{μ}} }}} = \frac{{{R_{\rm{b}}}}}{B} = \frac{1}{{{T_{\rm{b}}}B}} = \frac{n}{{{T_{\rm{s}}}B}} \ge \frac{n}{{{\mu _{k\max }}T_{\rm{s}}\;\!^2}}$

      MFSK (Multiple Frequency-Shift Keying)的信道带宽理论值为${{{R_{\rm{b}}}\left( {M + 3} \right)} / {2n}}$ (相干MFSK),则MFSK的带宽效率为

      ${\eta _{\rm{f}}} = \frac{{{R_{\rm{b}}}}}{B} = \frac{{2n}}{{M + 3}}$

      因此,当${\mu _{k\max }}T_{\rm{s}}\,\!^2 < {{\left( {M + 3} \right)} / 2}$即用于调频率键控的Chirp信号的最大时宽带宽积$\tau {B_{\max }} \!<\! {{\left( {M \!+\! 3} \right)}\! / 2}$时,调频率键控有优于MFSK的带宽效率。而MFSK的误码率性能优于调频率键控,故调频率键控与频率键控同时使用时,可以通过调整2种调制方式的配比,在误码率性能与带宽效率间折中选择。

    • 在仿真实验中,设定二进制数据对8调频率与8初始频率的64个Chirp信号进行调制,根据第2部分要求设计仿真参数为:射频${f\!_{\rm{p}}} = 10 \; {\rm{ GHz}}$,时宽$\tau = 1 \; {{{μ}}{\rm s}}$,占空比10%,带宽$B = 120 \; {\rm{ MHz}}$,主载波频率${f\!_{\rm{r}}} = 0 \; {\rm{ MHz}}$,调频率${\mu _{\rm{r}}} = 120 \; {{{\rm{MHz}}} /}$μs;副载波调频率组${\mu _k} = \left[ {{\rm{15,27,39,51,63,75,87,99}}} \right] $ MHz/μs,频率组${f\!_l} = \left[ {{\rm{0, 3, 6, 9, 12, 15, 18, 21}}} \right] \; {\rm{ MHz}}$,主副载波功率比为1:1,目标参数$\left[ {1000 \; {\rm{ m}},200 \; {{\rm{m}} / {\rm{s}}}} \right]$

    • 因主载波Chirp信号的初始频率${f\!_{\rm{r}}} = {f\!_1}$,调频率${\mu _{\rm{r}}} > {\mu _k}$在可用范围内达到最大,则主载波的频谱带宽内包含了副载波用于通信的所有频谱带宽,主载波用于雷达目标检测,接收机利用匹配滤波器对回波进行脉压处理,结果如图9

      图  9  脉压结果

      Figure 9.  Pulse compression result

      分析知主、副载波互相关性较低,不影响雷达的目标检测,第1旁瓣依然在–13 dB左右,经过匹配滤波器后副载波剩余量很小,幅度保持在–20 dB以下,而且增加主载波的功率后,剩余量幅度会更低,经过匹配滤波器后的回波脉冲串有很高的相关性,经过多普勒滤波器组进行动目标检测(Moving Target Detection, MTD)处理即可得出目标相对速度,表明在接收端仅需使用单一滤波器即可完成,测速结果如图10所示。

      图  10  动目标检测结果

      Figure 10.  MTD result

      在虚警概率等于10–4条件下,本文共享信号进行脉压处理和不同脉冲数积累MTD处理后,检测概率与信噪比(SNR)的关系曲线如图11所示,脉压与MTD处理利用相干积累提高了SNR,由于本文共享信号存在通信副载波分量,故与同参数下的单Chirp信号相比,检测概率有所下降;但在进行MTD相干处理后,提高了SNR,从而改善检测概率,而且相干积累的脉冲数越多,检测概率越优,故采用较多脉冲积累来弥补共享信号雷达检测性能降低的不足。

      图  11  检测概率与信噪比关系曲线

      Figure 11.  Curve of relationship between detection probability and signal-to-noise ratio

    • 副载波通过键控Chirp信号的调频率与初始频率来调制通信数据,通信接收端对回波进行$\left\{ {{\mu _k}} \right\}$对应的${2^{{n_1}}}$个固定阶次的FRFT处理,在$\left\{ {{f\!_l}} \right\}$对应的${2^{{n_2}}}$个固定位置处检测幅度值,得到上述${2^{{n_1}}} \times {2^{{n_2}}}$个检测点中幅度值高于设定阈值的检测点对应的阶次$p$和位置$u$,即可解调出调制进制数据。

      调制数据时,将通信数据串并转换并分组后,根据数据组的大小排列方式可在共享信号中叠加多个副载波,若后一组数据大于前一组数据,则可将这两组数据调制到同一共享信号中,则此脉冲就有多个副载波,通信接收时不需要改变解调方式,每个副载波携带的数据均可解调出,只需将数据组按大小排列,以此提高通信速率;若后一组数据不大于前一组数据,则后组数据在下一脉冲调制。因此,此共享信号的通信传输速率在$n{\rm{PRF}} \sim {2^n} \cdot n{\rm{PRF}}$范围内变化,$n = {n_1} + {n_2}$为单个副载波携带的比特位数。若仅采用主载波叠加单个副载波的形式,此仿真参数下的通信速率为$n{\rm{PRF}} = 600 \; {{{\rm{kb}}} / {\rm{s}}}$,通过改变参数增加调制位数可得更高传输速率。

      AWGN信道中,本文共享信号的误码率仿真曲线如图12所示,从上往下第3, 4, 5条曲线为64进制调制的3种不同配比,即16K-4F, 8K-8F和4K-16F, K表示调频率,F表示载频,由图12可见这3种方式的误码率性能逐渐改善,根据4.5节分析得知,这是由于FSK的误码率性能优于调频率调制,通过改变调频率与初始频率的不同配比可以调整本文共享信号的抗干扰性能与带宽效率。图中给出数字调制中的键控法MFSK, MASK, MPSK的理论误码率曲线作为对比参考,随着调制位数M的增大,MASK和MPSK的抗噪声性能下降,频带利用率上升,而MFSK抗噪声性能更好,有更好的误码率性能,但频带利用率较差[17]

      图  12  误码率随信噪比变化曲线

      Figure 12.  Change of error rate with SNR

      $\left\{ {\begin{array}{*{20}{l}} {{P_{{\rm{MASK}}}} = \left( {1 - \frac{1}{M}} \right){\rm{erfc}}\left( {\sqrt {\frac{{3r}}{{{M^2} - 1}}} } \right)}\\ {{P_{{\rm{MFSK}}}} = \frac{{M - 1}}{2}{\rm{erfc}}\left( {\sqrt {\frac{r}{2}} } \right)}\\ {{P_{{\rm{MPSK}}}} \approx {\rm{erfc}}\left( {\sqrt {2r} \sin \frac{{{π}} }{{2M}}} \right)} \end{array}} \right.$

    • 本文设计并研究了一种多载波雷达通信共享信号,通过对副载波Chirp信号的调频率与初始频率键控调制通信数据,利用主载波进行雷达目标检测。对共享信号的模糊函数及主副载波间的正交性进行了分析,对Chirp信号参数间关系进行设计,在通信接收端采用FRFT变换进行解调,并对共享信号的抗多普勒性能进行了分析。共享信号的设计实现了复杂集成电子装备平台中,雷达和通信信号能量和时间的一体化,这将是未来一体化电子战系统的一个重要的发展方向。

参考文献 (17)

目录

    /

    返回文章
    返回