CN109283498A - 一种线性调频脉冲信号相位误差曲线生成方法 - Google Patents

一种线性调频脉冲信号相位误差曲线生成方法 Download PDF

Info

Publication number
CN109283498A
CN109283498A CN201710598560.5A CN201710598560A CN109283498A CN 109283498 A CN109283498 A CN 109283498A CN 201710598560 A CN201710598560 A CN 201710598560A CN 109283498 A CN109283498 A CN 109283498A
Authority
CN
China
Prior art keywords
phase
phase error
error curve
data sample
pulse 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.)
Granted
Application number
CN201710598560.5A
Other languages
English (en)
Other versions
CN109283498B (zh
Inventor
段云鹏
王平
何靖
王百涛
朱智明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201710598560.5A priority Critical patent/CN109283498B/zh
Publication of CN109283498A publication Critical patent/CN109283498A/zh
Application granted granted Critical
Publication of CN109283498B publication Critical patent/CN109283498B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种线性调频脉冲信号相位误差曲线生成方法，本发明方法首先构建相位误差曲线生成系统，所构建的系统包括：数据样本生成模块、相位信息提取模块和相位误差曲线生成模块，通过数据样本生成模块对线性调频脉冲信号进行采样并生成3个数据样本，相位信息提取模块对这3个数据样本进行离散希尔伯特变换，提取出相位信息，相位误差曲线生成模块利用相位信息生成线性调频脉冲信号的相位误差曲线。本发明方法产生的相位误差曲线平滑，已成功应用在某雷达系统中，满足系统使用要求。

Description

一种线性调频脉冲信号相位误差曲线生成方法
技术领域
本发明涉及一种相位误差曲线生成方法，特别是一种线性调频脉冲信号相位误差曲线生成方法。
背景技术
雷达系统常使用直接波形存储直读技术结合正交调制技术的方法产生线性调频脉冲信号，尽管在设计时对硬件系统作了一些优化处理，但由于器件本身特性，不可避免地在产生线性调频脉冲信号过程中引入相位误差，为了定量分析相位误差，需要产生相位误差曲线。传统的线性调频脉冲信号相位误差曲线生成方法为：对线性调频脉冲信号进行采样，形成单一的数据样本；对数据样本进行离散希尔伯特变换，根据变换结果提取出相位信息；通过相位信息绘制出线性调频脉冲信号的相位误差曲线。此方法比较直观，且计算量小，但由于相位抖动问题，导致生成的相位误差曲线不平滑。
发明内容
本发明目的在于提供一种线性调频脉冲信号相位误差曲线生成方法，解决传统的线性调频脉冲信号相位误差曲线生成方法得到的相位误差曲线不平滑的问题。
一种线性调频脉冲信号相位误差曲线生成方法，其具体步骤为：
第一步搭建相位误差曲线生成系统
相位误差曲线生成系统，包括：数据样本生成模块、相位信息提取模块和相位误差曲线生成模块。所述：
数据样本生成模块的功能为：对线性调频脉冲信号进行采样并生成3个数据样本；
相位信息提取模块的功能为：对这3个数据样本进行离散希尔伯特变换，提取出相位信息；
相位误差曲线生成模块的功能为：利用相位信息生成线性调频脉冲信号的相位误差曲线。
第二步数据样本生成模块对线性调频脉冲信号进行采样并生成数据样本
数据样本生成模块在任意3个脉冲宽度的时间范围内对线性调频信号x(t)进行采样，生成3个数据样本x1(n)、x2(n)和x3(n)，采样频率为fs，采样间隔Ts＝1/fs，t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据样本的长度，N＝fs·τ；线性调频脉冲信号x(t)的脉冲宽度为τ，中心频率为f0，带宽为B，调频斜率μ＝B/τ。
第三步相位信息提取模块对数据样本进行离散希尔伯特变换，提取出相位信息
相位信息提取模块对数据样本x1(n)、x2(n)和x3(n)分别进行离散希尔伯特变换：[x1(n)]h＝hilbert[x1(n)]，[x2(n)]h＝hilbert[x2(n)]，[x3(n)]h＝hilbert[x3(n)]，得到3个数据样本的离散希尔伯特变换结果[x1(n)]h、[x2(n)]h和[x3(n)]h，其中hilbert[·]表示离散希尔伯特变换；
相位信息提取模块提取相位信息：φ1(n)＝angle{[x1(n)]h}，φ2(n)＝angle{[x2(n)]h}，φ3(n)＝angle{[x3(n)]h}，得到数据样本x1(n)的相位φ1(n)、数据样本x2(n)的相位φ2(n)和数据样本x3(n)的相位φ3(n)，其中angle{·}表示复数求相位运算，运算结果是角度值。
第四步相位误差曲线生成模块利用相位信息生成线性调频脉冲信号的相位误差曲线
相位误差曲线生成模块获取相位φ1(n)、φ2(n)和φ3(n)的均值相位 对均值相位进行多项式曲线拟合，得到多项式拟合曲线f1(n)；
相位误差曲线生成模块获取相位误差φ(n)：φ(n)＝f1(n)-[2πf0nTs+μ(nTs)2/2]，对相位误差φ(n)进行多项式曲线拟合，得到线性调频脉冲信号的相位误差曲线f2(n)。
至此，实现了线性调频脉冲信号相位误差曲线的生成。
更优的，其中第四步中对均值相位进行多项式曲线拟合时，多项式的次数为10。
更优的，其中第四步中对相位误差φ(n)进行多项式曲线拟合时，多项式的次数为2。
本发明方法解决了传统的线性调频脉冲信号相位误差曲线生成方法得到的相位误差曲线不平滑的问题，产生的相位误差曲线平滑，此种方法有效、可行，目前本方法已成功应用在某雷达系统中，满足系统使用要求。
具体实施方式
本发明目的在于提供一种线性调频脉冲信号相位误差曲线生成方法，解决传统的线性调频脉冲信号相位误差曲线生成方法得到的相位误差曲线不平滑的问题。
一种线性调频脉冲信号相位误差曲线生成方法，其具体步骤为：
第一步搭建相位误差曲线生成系统
相位误差曲线生成系统，包括：数据样本生成模块、相位信息提取模块和相位误差曲线生成模块。所述：
数据样本生成模块的功能为：对线性调频脉冲信号进行采样并生成3个数据样本；
相位信息提取模块的功能为：对这3个数据样本进行离散希尔伯特变换，提取出相位信息；
相位误差曲线生成模块的功能为：利用相位信息生成线性调频脉冲信号的相位误差曲线。
第二步数据样本生成模块对线性调频脉冲信号进行采样并生成数据样本
数据样本生成模块在任意3个脉冲宽度的时间范围内对线性调频信号x(t)进行采样，生成3个数据样本x1(n)、x2(n)和x3(n)，采样频率为fs，采样间隔Ts＝1/fs，t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据样本的长度，N＝fs·τ；线性调频脉冲信号x(t)的脉冲宽度为τ，中心频率为f0，带宽为B，调频斜率μ＝B/τ。
第三步相位信息提取模块对数据样本进行离散希尔伯特变换，提取出相位信息
相位信息提取模块对数据样本x1(n)、x2(n)和x3(n)分别进行离散希尔伯特变换：[x1(n)]h＝hilbert[x1(n)]，[x2(n)]h＝hilbert[x2(n)]，[x3(n)]h＝hilbert[x3(n)]，得到3个数据样本的离散希尔伯特变换结果[x1(n)]h、[x2(n)]h和[x3(n)]h，其中hilbert[·]表示离散希尔伯特变换；
相位信息提取模块提取相位信息：φ1(n)＝angle{[x1(n)]h}，φ2(n)＝angle{[x2(n)]h}，φ3(n)＝angle{[x3(n)]h}，得到数据样本x1(n)的相位φ1(n)、数据样本x2(n)的相位φ2(n)和数据样本x3(n)的相位φ3(n)，其中angle{·}表示复数求相位运算，运算结果是角度值。
第四步相位误差曲线生成模块利用相位信息生成线性调频脉冲信号的相位误差曲线
相位误差曲线生成模块获取相位φ1(n)、φ2(n)和φ3(n)的均值相位 对均值相位进行多项式曲线拟合，多项式次数设置为10，得到多项式拟合曲线f1(n)；
相位误差曲线生成模块获取相位误差φ(n)：φ(n)＝f1(n)-[2πf0nTs+μ(nTs)2/2]，对相位误差φ(n)进行多项式曲线拟合，多项式次数设置为2，得到线性调频脉冲信号的相位误差曲线f2(n)。
至此，实现了线性调频脉冲信号相位误差曲线的生成。

Claims (3)

1.一种线性调频脉冲信号相位误差曲线生成方法，其特征在于具体步骤为：
第一步 搭建相位误差曲线生成系统
相位误差曲线生成系统，包括：数据样本生成模块、相位信息提取模块和相位误差曲线生成模块；所述：
数据样本生成模块的功能为：对线性调频脉冲信号进行采样并生成3个数据样本；
相位信息提取模块的功能为：对这3个数据样本进行离散希尔伯特变换，提取出相位信息；
相位误差曲线生成模块的功能为：利用相位信息生成线性调频脉冲信号的相位误差曲线；
第二步 数据样本生成模块对线性调频脉冲信号进行采样并生成数据样本
数据样本生成模块在任意3个脉冲宽度的时间范围内对线性调频信号x(t)进行采样，生成3个数据样本x1(n)、x2(n)和x3(n)，采样频率为fs，采样间隔Ts＝1/fs，t为时间参数，n为时域点索引，n＝1,2,…,N，N为数据样本的长度，N＝fs·τ；线性调频脉冲信号x(t)的脉冲宽度为τ，中心频率为f0，带宽为B，调频斜率μ＝B/τ；
第三步 相位信息提取模块对数据样本进行离散希尔伯特变换，提取出相位信息
相位信息提取模块对数据样本x1(n)、x2(n)和x3(n)分别进行离散希尔伯特变换：[x1(n)]h＝hilbert[x1(n)]，[x2(n)]h＝hilbert[x2(n)]，[x3(n)]h＝hilbert[x3(n)]，得到3个数据样本的离散希尔伯特变换结果[x1(n)]h、[x2(n)]h和[x3(n)]h，其中hilbert[·]表示离散希尔伯特变换；
相位信息提取模块提取相位信息：φ1(n)＝angle{[x1(n)]h}，φ2(n)＝angle{[x2(n)]h}，φ3(n)＝angle{[x3(n)]h}，得到数据样本x1(n)的相位φ1(n)、数据样本x2(n)的相位φ2(n)和数据样本x3(n)的相位φ3(n)，其中angle{·}表示复数求相位运算，运算结果是角度值；
第四步 相位误差曲线生成模块利用相位信息生成线性调频脉冲信号的相位误差曲线
相位误差曲线生成模块获取相位φ1(n)、φ2(n)和φ3(n)的均值相位 对均值相位进行多项式曲线拟合，得到多项式拟合曲线f1(n)；
相位误差曲线生成模块获取相位误差 对相位误差进行多项式曲线拟合，得到线性调频脉冲信号的相位误差曲线f2(n)；
至此，实现了线性调频脉冲信号相位误差曲线的生成。
2.如权利要求1所述相位误差曲线生成方法，其特征在于：所述第四步中对均值相位进行多项式曲线拟合时，多项式的次数为10。
3.如权利要求1所述相位误差曲线生成方法，其特征在于：所述第四步中对相位误差进行多项式曲线拟合时，多项式的次数为2。
CN201710598560.5A 2017-07-21 2017-07-21 一种线性调频脉冲信号相位误差曲线生成方法 Active CN109283498B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710598560.5A CN109283498B (zh) 2017-07-21 2017-07-21 一种线性调频脉冲信号相位误差曲线生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710598560.5A CN109283498B (zh) 2017-07-21 2017-07-21 一种线性调频脉冲信号相位误差曲线生成方法

Publications (2)

Publication Number Publication Date
CN109283498A true CN109283498A (zh) 2019-01-29
CN109283498B CN109283498B (zh) 2022-07-15

Family

ID=65185175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710598560.5A Active CN109283498B (zh) 2017-07-21 2017-07-21 一种线性调频脉冲信号相位误差曲线生成方法

Country Status (1)

Country Link
CN (1) CN109283498B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110187310A (zh) * 2019-05-13 2019-08-30 北京遥感设备研究所 一种雷达设备lfm信号瞬时频率曲线拟合方法及系统
CN110221262A (zh) * 2019-07-01 2019-09-10 北京遥感设备研究所 一种雷达设备lfm信号主瓣降低确定平台及方法
CN110554362A (zh) * 2019-08-19 2019-12-10 北京遥感设备研究所 一种雷达设备大时宽带宽积lfm信号相位谱合成系统及方法
CN113884996A (zh) * 2021-09-18 2022-01-04 北京遥感设备研究所 一种雷达专用测试设备多时宽脉冲信号校正方法和系统
CN113933795A (zh) * 2021-10-20 2022-01-14 北京遥感设备研究所 一种雷达设备抗干扰掩护脉冲基带数据产生方法及装置
CN114019460A (zh) * 2021-10-20 2022-02-08 北京遥感设备研究所 一种雷达设备抗干扰掩护脉冲校正序列产生方法及装置
CN115002344A (zh) * 2022-05-11 2022-09-02 昆山丘钛微电子科技股份有限公司 误差曲面拟合方法、装置、设备及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101036068A (zh) * 2004-10-09 2007-09-12 莱卡地球系统公开股份有限公司 用于测距及检测非理想线性调频脉冲曲线的电光方法
US20090257528A1 (en) * 2002-11-27 2009-10-15 Schell Stephan V Quadrature modulation without carrier
CN102506760A (zh) * 2011-11-18 2012-06-20 东南大学 一种物体表面轮廓测量中的相位补偿方法
CN103163511A (zh) * 2013-03-06 2013-06-19 北京理工大学 一种数字阵列雷达频率步进信号相位补偿方法
CN106054159A (zh) * 2016-05-12 2016-10-26 北京航空航天大学 一种多普勒信号的瞬时频率提取方法
CN106170679A (zh) * 2015-09-02 2016-11-30 深圳大学 一种相位误差补偿方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090257528A1 (en) * 2002-11-27 2009-10-15 Schell Stephan V Quadrature modulation without carrier
CN101036068A (zh) * 2004-10-09 2007-09-12 莱卡地球系统公开股份有限公司 用于测距及检测非理想线性调频脉冲曲线的电光方法
CN102506760A (zh) * 2011-11-18 2012-06-20 东南大学 一种物体表面轮廓测量中的相位补偿方法
CN103163511A (zh) * 2013-03-06 2013-06-19 北京理工大学 一种数字阵列雷达频率步进信号相位补偿方法
CN106170679A (zh) * 2015-09-02 2016-11-30 深圳大学 一种相位误差补偿方法及装置
CN106054159A (zh) * 2016-05-12 2016-10-26 北京航空航天大学 一种多普勒信号的瞬时频率提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
彭榕辉等: "一种估计线性调频信号参数的方法", 《遥测遥控》 *
高洪青: "基于希尔伯特解调的LFM信号调频斜率的识别", 《计算机测量与控制》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110187310A (zh) * 2019-05-13 2019-08-30 北京遥感设备研究所 一种雷达设备lfm信号瞬时频率曲线拟合方法及系统
CN110221262A (zh) * 2019-07-01 2019-09-10 北京遥感设备研究所 一种雷达设备lfm信号主瓣降低确定平台及方法
CN110221262B (zh) * 2019-07-01 2021-02-26 北京遥感设备研究所 一种雷达设备lfm信号主瓣降低确定平台及方法
CN110554362A (zh) * 2019-08-19 2019-12-10 北京遥感设备研究所 一种雷达设备大时宽带宽积lfm信号相位谱合成系统及方法
CN113884996A (zh) * 2021-09-18 2022-01-04 北京遥感设备研究所 一种雷达专用测试设备多时宽脉冲信号校正方法和系统
CN113933795A (zh) * 2021-10-20 2022-01-14 北京遥感设备研究所 一种雷达设备抗干扰掩护脉冲基带数据产生方法及装置
CN114019460A (zh) * 2021-10-20 2022-02-08 北京遥感设备研究所 一种雷达设备抗干扰掩护脉冲校正序列产生方法及装置
CN113933795B (zh) * 2021-10-20 2024-07-05 北京遥感设备研究所 一种雷达设备抗干扰掩护脉冲基带数据产生方法及装置
CN115002344A (zh) * 2022-05-11 2022-09-02 昆山丘钛微电子科技股份有限公司 误差曲面拟合方法、装置、设备及介质
CN115002344B (zh) * 2022-05-11 2024-04-19 昆山丘钛微电子科技股份有限公司 误差曲面拟合方法、装置、设备及介质

Also Published As

Publication number Publication date
CN109283498B (zh) 2022-07-15

Similar Documents

Publication Publication Date Title
CN109283498A (zh) 一种线性调频脉冲信号相位误差曲线生成方法
CN107064883B (zh) 一种基于压缩感知的线性调频雷达信号处理方法
Brandwood Fourier transforms in radar and signal processing
CN106560800B (zh) 在雷达和声纳应用中调整定点快速傅立叶变换
WO2018200577A1 (en) Digital communication using lattice division multiplexing
CN107247254B (zh) 一种非线性调频信号去斜处理方法
CN104749432B (zh) 基于聚焦s变换的多分量非平稳信号瞬时频率估计方法
Iglesias et al. Real-time low-complexity automatic modulation classifier for pulsed radar signals
EP2499504A1 (en) A precision measurement of waveforms using deconvolution and windowing
CN104655929B (zh) 一种时域信号的数字时频测量方法及相应的目标识别方法
CN110189743A (zh) 波形拼接中的拼接点平滑方法、装置及存储介质
CN104854473B (zh) 使用多普勒聚焦的亚奈奎斯特雷达处理
JP2020524286A5 (zh)
JP6105286B2 (ja) デジタル信号処理方法、デジタル信号処理装置、およびプログラム
KR20160083276A (ko) 차량 위치 검출 장치 및 방법
CN109188370B (zh) 一种雷达设备lfm脉冲信号包络曲线拟合方法及系统
Batenkov et al. Rethinking super-resolution: the bandwidth selection problem
WO2016197629A1 (en) System and method for frequency estimation
CN110221262B (zh) 一种雷达设备lfm信号主瓣降低确定平台及方法
US7423470B2 (en) Pulse generator
CN105743830B (zh) 一种mfsk信号的解调方法
US8867862B1 (en) Self-optimizing analysis window sizing method
Lin Ridges reconstruction based on inverse wavelet transform
EP3136668B1 (en) Channel equalization and tracking device and method, and receiver
CN110261827B (zh) 一种雷达设备lfm脉冲压缩网络群延时曲线生成平台及方法

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