CN107315180A - 一种两频点太赫兹全息成像方法及系统 - Google Patents

一种两频点太赫兹全息成像方法及系统 Download PDF

Info

Publication number
CN107315180A
CN107315180A CN201710334334.6A CN201710334334A CN107315180A CN 107315180 A CN107315180 A CN 107315180A CN 201710334334 A CN201710334334 A CN 201710334334A CN 107315180 A CN107315180 A CN 107315180A
Authority
CN
China
Prior art keywords
signal
frequency
way
phase
wave filter
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.)
Pending
Application number
CN201710334334.6A
Other languages
English (en)
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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201710334334.6A priority Critical patent/CN107315180A/zh
Publication of CN107315180A publication Critical patent/CN107315180A/zh
Pending legal-status Critical Current

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
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06—Systems determining position data of a target
    • G01S17/46—Indirect determination of position data
    • G01S17/48—Active triangulation systems, i.e. using the transmission and reflection of electromagnetic waves other than radio waves
    • 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
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/50—Systems of measurement based on relative movement of target

Landscapes

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

Abstract

本发明公开了一种两频点太赫兹全息成像方法及系统。其中，该方法包括如下步骤：构建两频点太赫兹全息成像系统；根据第一发射信号得到第一输出信号和第二输出信号，根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；图像处理单元第一幅度和第一相位和第二幅度和第二相位进行处理得到太赫兹成像全息图。本发明解决了太赫兹波段运动目标难以得到物体准确距离信息的问题，达到的运动目标可以进行三维立体成像，很好的展现了成像的景深效果。

Description

一种两频点太赫兹全息成像方法及系统
技术领域
本发明属于探测成像技术领域，尤其涉及一种两频点太赫兹全息成像方法及系统。
背景技术
太赫兹全息成像可分为点频和宽带两种模式，其中点频只能实现二维全息成像，宽频模式可以实现三维全息成像，并获得比单频成像更高的分辨率，该成像技术不受Fresnel近似条件约束。然而，这两种系统必须预先获知成像目标与扫描孔径之间准确的距离，对运动目标或者远距离目标探测是致命的。Zhaoyang Sun，Chao Li等人在“FastThree-Dimensional Image Reconstruction of Targets Under the Illumination ofTerahertz Gaussian Beams With Enhanced Phase-Shift Migration to ImproveComputation Efficiency”(IEEE transactions on terahertz science andtechnology)，提出了单一频点太赫兹3D增强成像方法，能够显著提升太赫兹成像质量。刘艺青在“一种近距离太赫兹三维全息成像方法及系统”(专利申请号：CN201110038253.4)发射连续波太赫兹信号；在以时间、圆周角及Z轴方向形成三维域中测得回波信号；对最大化的回波信号进行傅里叶变换，再利用相位固定法，实现回波信号的时域向频域的转换；对频域回波信号利用圆柱形傅里叶变换及双线性插值运算进行运动补偿，得到直角坐标系下重构的目标散射强度信号，根据重构的目标散射强度信号进行三维全息成像。吴卫东等“在一种量子级联激光器太赫兹源数字全息成像系统”利用量子级联激光器太赫兹源、光束整形模块、全息光路模块、图像探测与采集模块组成成像系统。以上这些都是对静止物体进行三维全息成像，对运动目标来说，由于目标相对距离的变化，成像将会模糊，全息将难以实现。因而，在全息成像中目标物体距离的探测就成为了一个亟待解决的问题。
发明内容
本发明解决的技术问题是：克服现有技术的不足，提供了一种两频点太赫兹全息成像方法及系统，解决了太赫兹波段运动目标难以得到物体准确距离信息的问题，达到的运动目标可以进行三维立体成像，很好的展现了成像的景深效果。
本发明目的通过以下技术方案予以实现：根据一个方面，本发明提供了一种两频点太赫兹全息成像方法，该方法包括以下步骤：
步骤一：构建两频点太赫兹全息成像系统，其中，两频点太赫兹全息成像系统包括第一信号源、第二信号源、第一倍频器、第二倍频器、开关、本振信号源、第一混频器、第二混频器、第一滤波器、第二滤波器、第三混频器、第四混频器、移相器、第三滤波器、第四滤波器、第一放大器、第二放大器、发射天线和接收天线；
步骤二：第一信号源产生的信号经过第一倍频器和第一放大器后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线后照射目标，经目标反射后得到第一反射信号进入接收天线，第一反射信号与本振信号源产生的本征信号在第一混频器混频后经过第一滤波器形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源产生的本征信号在第二混频器混频后经过第二滤波器形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器混频后经过第三滤波器得到第一输出信号；第一参考中频信号另一路经过移相器与第一接收中频信号另一路在第四混频器混频后经过第四滤波器得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
步骤三：第二信号源产生的信号经过第二倍频器和第二放大器后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线后照射目标，经目标反射后得到第二反射信号进入接收天线，第二反射信号与本振信号源产生的本征信号在第一混频器混频后经过第一滤波器形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源产生的本征信号在第二混频器混频后经过第二滤波器形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器混频后经过第三滤波器得到第三输出信号；第二参考中频信号另一路经过移相器与第二接收中频信号另一路在第四混频器混频后经过第四滤波器得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；其中，频率fRF1与频率fRF2满足关系C表示光速。
步骤四：图像处理单元对步骤二中的第一幅度和第一相位和步骤三中的第二幅度和第二相位进行处理得到太赫兹成像全息图。
上述两频点太赫兹全息成像方法中，在步骤二中，第一发射信号为正弦信号，第一发射信号的公式为：sRF1(t)＝cos(ωRF1t+θ1)，其中，ωRF1＝2πfRF1，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t)为第一发射信号。
上述两频点太赫兹全息成像方法中，在步骤二中，本征信号为如下公式:sLO(t)＝cos(ωLOt+θ2),其中，ωLO为本征信号的角频率，θ2为本征信号的相位，t为时间，sLO(t)为本征信号。
上述两频点太赫兹全息成像方法中，在步骤二中，第一反射信号为如下公式:sRF1(t-Δt)＝cos(ωRF1(t-Δt)+θ1)，其中，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t-Δt)为第一发射信号，Δt为时延。
上述两频点太赫兹全息成像方法中，在步骤二中，第一发射信号另一路与本振信号源产生的本征信号在第二混频器混频后的信号sIF为:sIF＝sRF1(t-Δt)·sLO＝cos(ωRF1(t-Δt)+θ1)cos(ωLOt+θ2)。
上述两频点太赫兹全息成像方法中，第一接收中频信号sIFT的公式如下:sIFT＝cos(ωRF1(t-Δt)+θ1-(ωLOt+θ2))。
上述两频点太赫兹全息成像方法中，第一参考中频信号sIFR的公式如下:sIFR＝cos(ωRF1(t)+θ1-(ωLOt+θ2))。
上述两频点太赫兹全息成像方法中，第一输出信号为sI＝cos(ωRF1Δt)；第二输出信号为sI＝sin(ωRF1Δt)。
上述两频点太赫兹全息成像方法中，在步骤三中，第二发射信号为正弦信号，第一发射信号的公式为：sRF2(t)＝cos(ωRF2t+θ1)，其中，ωRF2＝2πfRF2，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t)为第一发射信号。
上述两频点太赫兹全息成像方法中，在步骤三中，第二反射信号为如下公式:sRF2(t-Δt)＝cos(ωRF2(t-Δt)+θ1)，其中，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t-Δt)为第二发射信号，Δt为时延。
上述两频点太赫兹全息成像方法中，在步骤三中，第二发射信号另一路与本振信号源产生的本征信号在第二混频器混频后的信号sIF'为:sIF'＝sRF2(t-Δt)·sLO＝cos(ωRF2(t-Δt)+θ1)cos(ωLOt+θ2)。
上述两频点太赫兹全息成像方法中，第二接收中频信号sIFT'的公式如下:sIFT'＝cos(ωRF2(t-Δt)+θ1-(ωLOt+θ2))。
上述两频点太赫兹全息成像方法中，第二参考中频信号sIFR'的公式如下:sIFR'＝cos(ωRF2(t)+θ1-(ωLOt+θ2))。
上述两频点太赫兹全息成像方法中，第三输出信号为sI'＝cos(ωRF2Δt)；第四输出信号为sI'＝sin(ωRF2Δt)。
根据另一个方面，本发明还提供了一种两频点太赫兹全息成像系统，其特征在于包括：第一信号源、第二信号源、第一倍频器、第二倍频器、开关、本振信号源、第一混频器、第二混频器、第一滤波器、第二滤波器、第三混频器、第四混频器、移相器、第三滤波器、第四滤波器、第一放大器、第二放大器、发射天线和接收天线；其中，
第一信号源产生的信号经过第一倍频器和第一放大器后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线后照射目标，经目标反射后得到第一反射信号进入接收天线，第一反射信号与本振信号源产生的本征信号在第一混频器混频后经过第一滤波器形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源产生的本征信号在第二混频器混频后经过第二滤波器形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器混频后经过第三滤波器得到第一输出信号；第一参考中频信号另一路经过移相器与第一接收中频信号另一路在第四混频器混频后经过第四滤波器得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
第二信号源产生的信号经过第二倍频器和第二放大器后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线后照射目标，经目标反射后得到第二反射信号进入接收天线，第二反射信号与本振信号源产生的本征信号在第一混频器混频后经过第一滤波器形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源产生的本征信号在第二混频器混频后经过第二滤波器形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器混频后经过第三滤波器得到第三输出信号；第二参考中频信号另一路经过移相器与第二接收中频信号另一路在第四混频器混频后经过第四滤波器得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；其中，频率fRF1与频率fRF2满足关系C表示光速。
本发明与现有技术相比具有如下有益效果：
(1)本发明不用事先预先获知成像目标与T/R天线口径之间关系就可以对目标进行全息成像。
(2)本发明可以得到运动物体距离信息，得到更为清晰的三维图像。
(3)本发明利用太赫兹干涉思想，调整干涉频率，可以进行更为精确的探测。
附图说明
图1是本发明的两频点太赫兹全息成像系统的框图；
图2是本发明的电路仿真图；
图3(a)是本发明的对比图像；
图3(b)是本发明的全息处理的图像。
具体实施方式
下面结合附图对本发明作进一步详细说明：
图1是本发明的两频点太赫兹全息成像系统的框图。结合图1，本发明提供了一种两频点太赫兹全息成像方法，该方法包括以下步骤：
步骤一：构建两频点太赫兹全息成像系统，其中，两频点太赫兹全息成像系统包括第一信号源1、第二信号源2、第一倍频器3、第二倍频器4、开关5、本振信号源6、第一混频器7、第二混频器8、第一滤波器9、第二滤波器10、第三混频器11、第四混频器12、移相器13、第三滤波器14、第四滤波器15、第一放大器16、第二放大器17、发射天线18和接收天线19；
步骤二：第一信号源1产生的信号经过第一倍频器3和第一放大器16后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线18后照射目标，经目标反射后得到第一反射信号进入接收天线19，第一反射信号与本振信号源6产生的本征信号在第一混频器7混频后经过第一滤波器9形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后经过第二滤波器10形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器11混频后经过第三滤波器14得到第一输出信号；第一参考中频信号另一路经过移相器13与第一接收中频信号另一路在第四混频器12混频后经过第四滤波器15得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
步骤三：第二信号源2产生的信号经过第二倍频器4和第二放大器17后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线18后照射目标，经目标反射后得到第二反射信号进入接收天线19，第二反射信号与本振信号源6产生的本征信号在第一混频器7混频后经过第一滤波器9形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后经过第二滤波器10形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器11混频后经过第三滤波器14得到第三输出信号；第二参考中频信号另一路经过移相器13与第二接收中频信号另一路在第四混频器12混频后经过第四滤波器15得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；其中，频率fRF1与频率fRF2满足关系其中，C为光速。
步骤四：图像处理单元对步骤二中的第一幅度和第一相位和步骤三中的第二幅度和第二相位进行处理得到太赫兹成像全息图。具体的，图像处理单元对步骤二中的第一幅度和第一相位和步骤三中的第二幅度和第二相位进行处理得到太赫兹成像全息图。具体包括：由于实际成像目标物体是一个三维的物体，每个物点与成像面的实际距离是一个变化值，而其最大变化量正好是物体可测的最大纵深，传统成像方法必须事先知道成像目标与天线之间的准确关系。对于运动目标来说，成像目标与天线之间的距离是不断变化的，而本发明可以利用两个信号之间的相位差，来得出成像目标与天线之间的距离，根据成像面的实际距离之间的变化，来得到物体的带有“景深”的全息图像。
具体的，在步骤二中，第一发射信号为正弦信号，第一发射信号的公式为：
sRF1(t)＝cos(ωRF1t+θ1)，
其中，ωRF1＝2πfRF1，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t)为第一发射信号。
本征信号为:sLO(t)＝cos(ωLOt+θ2)，其中，ωLO为本征信号的角频率，θ2为本征信号的相位，t为时间，sLO(t)为本征信号。
第一反射信号为:sRF1(t-Δt)＝cos(ωRF1(t-Δt)+θ1)，其中，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t-Δt)为第一发射信号，Δt为时延。
第一发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后的信号sIF为:sIF＝sRF1(t-Δt)·sLO＝cos(ωRF1(t-Δt)+θ1)cos(ωLOt+θ2)。
第一接收中频信号sIFT为:sIFT＝cos(ωRF1(t-Δt)+θ1-(ωLOt+θ2))。
第一参考中频信号sIFR为:sIFR＝cos(ωRF1(t)+θ1-(ωLOt+θ2))。
第一输出信号为sI＝cos(ωRF1Δt)；第二输出信号为sI＝sin(ωRF1Δt)。
具体的，在步骤三中，第二发射信号为正弦信号，第一发射信号的公式为：sRF2(t)＝cos(ωRF2t+θ1)，其中，ωRF2＝2πfRF2，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t)为第一发射信号。
第二反射信号为:sRF2(t-Δt)＝cos(ωRF2(t-Δt)+θ1)，其中，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t-Δt)为第二发射信号，Δt为时延。
第二发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后的信号sIF'为:sIF'＝sRF2(t-Δt)·sLO＝cos(ωRF2(t-Δt)+θ1)cos(ωLOt+θ2)。
第二接收中频信号sIFT'为:sIFT'＝cos(ωRF2(t-Δt)+θ1-(ωLOt+θ2))。
第二参考中频信号sIFR'为:sIFR'＝cos(ωRF2(t)+θ1-(ωLOt+θ2))。
第三输出信号为sI'＝cos(ωRF2Δt)；第四输出信号为sI'＝sin(ωRF2Δt)。
由第一输出信号sI＝cos(ωRF1Δt)得出频点fRF1的相位延迟为由第三输出信号为sI'＝cos(ωRF2Δt)得出频点fRF2的相位延迟为相位延迟还可以表示为2πfRF1Δt＝θ1+2πN1，相位延迟还可以表示为2πfRF2Δt＝θ2+2πN2，其中，N1为整数，N2为整数。
在两个频点的选择上，都处于太赫兹频段，并且不能相差太远。例如一个频点为400GHz，一个为403GHz，两个频点的波长差小于0.1毫米，因此从目标发射回来的电磁波对应的整周期数要么相同要么差1，如果fRF1<fRF2和θ1<θ2,则有N1＝N2＝N或N1＝N＝N2+1，由N1＝N2＝N得到由N1＝N＝N2+1得到
由2πfRF1Δt＝θ1+2πN1和2πfRF2Δt＝θ2+2πN2得到时延Δt，由时延Δt得到最大无模糊探测距离为：R＝c·Δt/2，其中，c为光速。
为了验证设计中的方案，进行建模仿真，具体模型如图2所示。我们所构建的模型中为了验证本方法的正确性，设计了在一定信噪比情况下，进行相位精度测试。因为太赫兹频点的两个频率点相近，我们采用一个频率进行验证。首先倍频产生的射频信号经过高斯白噪声信道，而后经过延时器，进而与本地参考信号混频滤波后形成中频信号，随后，进一步与参考中频信号的1支路和Q支路进行混频滤波，最后获得相位延迟信息，并由示波器显示，并将输出相应的数据进行处理，获得统计不同信噪比下相位均值和方差。
为不失一般性，假定反射信号进入接收机相对从天线口径发射时产生的时延为λ/8，对应相位为45°，设置不同的信噪比，采用图2给出模型，仿真出来的相位均值和方差见表l。从表l中不难看出，随着信噪比的改善，被测延迟相位与设定的相位延迟越接近，同时方差也更小。作为全息成像雷达，如果用于近距离探测，性噪比要求达到10dB，相位探测精度完全由于可以优于λ/10。如果用于目标探测，其作用距离可能较大，同时传统信道会变得复杂，信噪比也变得不可预测，在这种条件下一般信噪比必须大于-5dB，相应的关注不是细节，而是整体，系统依然可以完成成像探测。
表1不同信噪比下相位均值和方差
根据基尔霍夫衍射理论和角谱理论，设定目标表面对某一频率信号的反射系数为U(x，y，0)，传播一段距离z到成像面后为U(x，y，Z)，它们对应的傅里叶变换分别为A(fx,fy,0)和A(fx,fy,Z)，而这两个角谱存在如下关系：
由此可见，如果已有成像面上的空域数据U(x,y,Z)和传播距离Z，对应进行傅里叶变换即可获得角谱A(fx,fy,Z)，根据上式，可以计算出物面角谱A(x,y,0)，然后进行反傅里叶变换即可获得物面反射系数。
U(x,y,0)可用下列公式描述：
由于实际成像目标物体是一个三维的物体，每个物点与成像面的实际距离Z是一个变化值，而其最大变化量正好是物体可测的最大纵深，因此Z值可用取最大值，最小值或者均值，对应到成像后的图像是“景深”会有一个变化，这三个取值对成像质量不会产生明显差异。
由于实际采集的目标数据均是离散数据，实际上采用二维或者三维的快速傅里叶变换算法处理。为了更好说明全息成像算法的有效性，我们采用图像像素大小为512×512，全息成像算法处理前后的图像如图3(a)和图3(b)所示。以400GHz附近频点为例，成像距离100米。可以发现，经过全息成像处理之后，“景深”明显加强，图像更加清晰。
本发明还提供了一种两频点太赫兹全息成像系统，如图1所示，该两频点太赫兹全息成像系统包括第一信号源1、第二信号源2、第一倍频器3、第二倍频器4、开关5、本振信号源6、第一混频器7、第二混频器8、第一滤波器9、第二滤波器10、第三混频器11、第四混频器12、移相器13、第三滤波器14、第四滤波器15、第一放大器16、第二放大器17、发射天线18和接收天线19；其中，
第一信号源1产生的信号经过第一倍频器3和第一放大器16后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线18后照射目标，经目标反射后得到第一反射信号进入接收天线19，第一反射信号与本振信号源6产生的本征信号在第一混频器7混频后经过第一滤波器9形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后经过第二滤波器10形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器11混频后经过第三滤波器14得到第一输出信号；第一参考中频信号另一路经过移相器13与第一接收中频信号另一路在第四混频器12混频后经过第四滤波器15得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
第二信号源2产生的信号经过第二倍频器4和第二放大器17后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线18后照射目标，经目标反射后得到第二反射信号进入接收天线19，第二反射信号与本振信号源6产生的本征信号在第一混频器7混频后经过第一滤波器9形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源6产生的本征信号在第二混频器8混频后经过第二滤波器10形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器11混频后经过第三滤波器14得到第三输出信号；第二参考中频信号另一路经过移相器13与第二接收中频信号另一路在第四混频器12混频后经过第四滤波器15得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元。
以上所述的实施例只是本发明较优选的具体实施方式，本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。

Claims (10)

1.一种两频点太赫兹全息成像方法，其特征在于，所述方法包括以下步骤：
步骤一：构建两频点太赫兹全息成像系统，其中，两频点太赫兹全息成像系统包括第一信号源(1)、第二信号源(2)、第一倍频器(3)、第二倍频器(4)、开关(5)、本振信号源(6)、第一混频器(7)、第二混频器(8)、第一滤波器(9)、第二滤波器(10)、第三混频器(11)、第四混频器(12)、移相器(13)、第三滤波器(14)、第四滤波器(15)、第一放大器(16)、第二放大器(17)、发射天线(18)和接收天线(19)；
步骤二：第一信号源(1)产生的信号经过第一倍频器(3)和第一放大器(16)后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线(18)后照射目标，经目标反射后得到第一反射信号进入接收天线(19)，第一反射信号与本振信号源(6)产生的本征信号在第一混频器(7)混频后经过第一滤波器(9)形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后经过第二滤波器(10)形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器(11)混频后经过第三滤波器(14)得到第一输出信号；第一参考中频信号另一路经过移相器(13)与第一接收中频信号另一路在第四混频器(12)混频后经过第四滤波器(15)得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
步骤三：第二信号源(2)产生的信号经过第二倍频器(4)和第二放大器(17)后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线(18)后照射目标，经目标反射后得到第二反射信号进入接收天线(19)，第二反射信号与本振信号源(6)产生的本征信号在第一混频器(7)混频后经过第一滤波器(9)形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后经过第二滤波器(10)形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器(11)混频后经过第三滤波器(14)得到第三输出信号；第二参考中频信号另一路经过移相器(13)与第二接收中频信号另一路在第四混频器(12)混频后经过第四滤波器(15)得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；其中，频率fRF1与频率fRF2满足关系C表示光速；
步骤四：图像处理单元对步骤二中的第一幅度和第一相位和步骤三中的第二幅度和第二相位进行处理得到太赫兹成像全息图。
2.根据权利要求1所述的两频点太赫兹全息成像方法，其特征在于：在步骤二中，第一发射信号为正弦信号，第一发射信号的公式为：
sRF1(t)＝cos(ωRF1t+θ1)，
其中，ωRF1＝2πfRF1，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t)为第一发射信号；
本征信号为如下公式:
sLO(t)＝cos(ωLOt+θ2)，
其中，ωLO为本征信号的角频率，θ2为本征信号的相位，t为时间，sLO(t)为本征信号；
第一反射信号为如下公式:
sRF1(t-Δt)＝cos(ωRF1(t-Δt)+θ1)，
其中，ωRF1为第一发射信号的角频率，θ1为初始相位，t为时间，sRF1(t-Δt)为第一发射信号，Δt为时延。
3.根据权利要求2所述的两频点太赫兹全息成像方法，其特征在于：在步骤二中，第一发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后的信号sIF为:
sIF＝sRF1(t-Δt)·sLO＝cos(ωRF1(t-Δt)+θ1)cos(ωLOt+θ2)。
4.根据权利要求3所述的两频点太赫兹全息成像方法，其特征在于：在步骤二中，第一接收中频信号sIFT的公式如下:
sIFT＝cos(ωRF1(t-Δt)+θ1-(ωLOt+θ2))。
5.根据权利要求4所述的两频点太赫兹全息成像方法，其特征在于：在步骤二中，第一参考中频信号sIFR的公式如下:
sIFR＝cos(ωRF1(t)+θ1-(ωLOt+θ2))。
6.根据权利要求5所述的两频点太赫兹全息成像方法，其特征在于：在步骤二中，第一输出信号为sI＝cos(ωRF1Δt)；第二输出信号为sI＝sin(ωRF1Δt)。
7.根据权利要求6所述的两频点太赫兹全息成像方法，其特征在于：在步骤三中，第二发射信号为正弦信号，第二发射信号的公式为：
sRF2(t)＝cos(ωRF2t+θ1)，
其中，ωRF2＝2πfRF2，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t)为第一发射信号；
第二反射信号为如下公式:
sRF2(t-Δt)＝cos(ωRF2(t-Δt)+θ1)，
其中，ωRF2为第二发射信号的角频率，θ1为初始相位，t为时间，sRF2(t-Δt)为第二发射信号，Δt为时延；
第二发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后的信号sIF'为:
sIF'＝sRF2(t-Δt)·sLO＝cos(ωRF2(t-Δt)+θ1)cos(ωLOt+θ2)；
第二接收中频信号sIFT'的公式如下:
sIFT'＝cos(ωRF2(t-Δt)+θ1-(ωLOt+θ2)) 。
8.根据权利要求7所述的两频点太赫兹全息成像方法，其特征在于：第二参考中频信号sIFR'的公式如下:
sIFR'＝cos(ωRF2(t)+θ1-(ωLOt+θ2))。
9.根据权利要求8所述的两频点太赫兹全息成像方法，其特征在于：第三输出信号为sI'＝cos(ωRF2Δt)；第四输出信号为sI'＝sin(ωRF2Δt)。
10.一种两频点太赫兹全息成像系统，其特征在于包括：第一信号源(1)、第二信号源(2)、第一倍频器(3)、第二倍频器(4)、开关(5)、本振信号源(6)、第一混频器(7)、第二混频器(8)、第一滤波器(9)、第二滤波器(10)、第三混频器(11)、第四混频器(12)、移相器(13)、第三滤波器(14)、第四滤波器(15)、第一放大器(16)、第二放大器(17)、发射天线(18)和接收天线(19)；其中，
第一信号源(1)产生的信号经过第一倍频器(3)和第一放大器(16)后转换为频率为fRF1的第一发射信号，第一发射信号分为两路，第一发射信号一路传输至发射天线(18)后照射目标，经目标反射后得到第一反射信号进入接收天线(19)，第一反射信号与本振信号源(6)产生的本征信号在第一混频器(7)混频后经过第一滤波器(9)形成第一接收中频信号，第一接收中频信号分成两路；第一发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后经过第二滤波器(10)形成第一参考中频信号，第一参考中频信号分为两路；第一参考中频信号一路与第一接收中频信号一路在第三混频器(11)混频后经过第三滤波器(14)得到第一输出信号；第一参考中频信号另一路经过移相器(13)与第一接收中频信号另一路在第四混频器(12)混频后经过第四滤波器(15)得到第二输出信号；根据第一输出信号和第二输出信号得到第一幅度和第一相位，然后将第一幅度和第一相位传输至图像信息处理单元；
第二信号源(2)产生的信号经过第二倍频器(4)和第二放大器(17)后转换为频率为fRF2的第二发射信号，第二发射信号分为两路，第二发射信号一路传输至发射天线(18)后照射目标，经目标反射后得到第二反射信号进入接收天线(19)，第二反射信号与本振信号源(6)产生的本征信号在第一混频器(7)混频后经过第一滤波器(9)形成第二接收中频信号，第二接收中频信号分成两路；第二发射信号另一路与本振信号源(6)产生的本征信号在第二混频器(8)混频后经过第二滤波器(10)形成第二参考中频信号，第二参考中频信号分为两路；第二参考中频信号一路与第二接收中频信号一路在第三混频器(11)混频后经过第三滤波器(14)得到第三输出信号；第二参考中频信号另一路经过移相器(13)与第二接收中频信号另一路在第四混频器(12)混频后经过第四滤波器(15)得到第四输出信号；根据第三输出信号和第四输出信号得到第二幅度和第二相位，然后将第二幅度和第二相位传输至图像信息处理单元；其中，频率fRF1与频率fRF2满足关系C表示光速。
CN201710334334.6A 2017-05-12 2017-05-12 一种两频点太赫兹全息成像方法及系统 Pending CN107315180A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710334334.6A CN107315180A (zh) 2017-05-12 2017-05-12 一种两频点太赫兹全息成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710334334.6A CN107315180A (zh) 2017-05-12 2017-05-12 一种两频点太赫兹全息成像方法及系统

Publications (1)

Publication Number Publication Date
CN107315180A true CN107315180A (zh) 2017-11-03

Family

ID=60181427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710334334.6A Pending CN107315180A (zh) 2017-05-12 2017-05-12 一种两频点太赫兹全息成像方法及系统

Country Status (1)

Country Link
CN (1) CN107315180A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113608175A (zh) * 2021-08-03 2021-11-05 上海无线电设备研究所 一种基于量子级联的rcs测量收发系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901498A (zh) * 2012-12-26 2014-07-02 中国电子科技集团公司第五十研究所 增强被动太赫兹成像效果的系统
CN105676218A (zh) * 2016-01-08 2016-06-15 中国人民解放军国防科学技术大学 一种基于双频干涉的太赫兹雷达目标三维成像方法
CN105974491A (zh) * 2016-05-06 2016-09-28 中国工程物理研究院电子工程研究所 一种基于矢量网络分析仪的太赫兹安检成像系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901498A (zh) * 2012-12-26 2014-07-02 中国电子科技集团公司第五十研究所 增强被动太赫兹成像效果的系统
CN105676218A (zh) * 2016-01-08 2016-06-15 中国人民解放军国防科学技术大学 一种基于双频干涉的太赫兹雷达目标三维成像方法
CN105974491A (zh) * 2016-05-06 2016-09-28 中国工程物理研究院电子工程研究所 一种基于矢量网络分析仪的太赫兹安检成像系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
汤燕 等: "一种毫米波全息成像系统方案及其仿真", 《现代雷达》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113608175A (zh) * 2021-08-03 2021-11-05 上海无线电设备研究所 一种基于量子级联的rcs测量收发系统
CN113608175B (zh) * 2021-08-03 2023-09-19 上海无线电设备研究所 一种基于量子级联的rcs测量收发系统

Similar Documents

Publication Publication Date Title
CN106405541B (zh) 全相参连续波多普勒雷达及其测距测速方法
CN107340269A (zh) 近距离无损探测的超宽带太赫兹三维成像系统及方法
JP4665590B2 (ja) 干渉型レーダ
CN102426360B (zh) 空中微动旋转目标的二维isar成像方法
CN100538394C (zh) 一种采用多发多收频分信号的宽带雷达及其成像方法
CN1932551B (zh) 雷达
CN102288963B (zh) 基于子孔径参数估计的双基地isar图像融合方法
CN107238866B (zh) 基于合成孔径技术的毫米波视频成像系统及方法
JP2006177907A (ja) 干渉型レーダー
CN102707273B (zh) 一种THz雷达的收发前端
CN207163907U (zh) 近距离无损探测的超宽带太赫兹三维成像系统
CN104181531B (zh) 一种基于相控阵雷达的三维关联成像方法
CN102812379B (zh) 超分辨率成像雷达
CN112505719A (zh) 多普勒频率校正二次补偿的激光测风雷达测风方法及系统
CN115508820A (zh) 一种线性调频脉冲雷达的目标探测方法
Chen et al. Wave measurements with multi-frequency HF radar in the East China Sea
CN109001729A (zh) 太赫兹成像中调频连续波线性度实时校准方法及其系统
CN110554377B (zh) 基于多普勒中心偏移的单通道sar二维流场反演方法及系统
CN202693789U (zh) 一种THz雷达的收发前端
CN105866743B (zh) 基于混响室的雷达杂波频谱特性模拟系统及方法
CN107315180A (zh) 一种两频点太赫兹全息成像方法及系统
CN111750972A (zh) 一种适用于双频连续激光雷达的目标振动测量方法
CN106526547A (zh) 基于InSAR技术的直线扫描近场RCS测试杂波抑制方法
CN106405551B (zh) 一种v-fm信号体制isar二维成像方法
CN105842681B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171103