基于经验模态分解的探地雷达信号去噪处理
杨建军刘鸿福
(太原理工大学太原 030024
【摘要】探地雷达作为一种先进的地球物理探测方法,具有探测效率高、操作简单、采样迅速、无损伤探测、探测分辨率高等优点。探地雷达的信号的去噪问题已成为一个公认的技术难题。本文用经验模态分解的方法对探地雷达信号进行信号去噪处理,并取得了良好的效果。
【关键词】探地雷达;经验模态分解;信号去噪
1引言
探地雷达又称地质雷达 ,是近几年迅速发展起来的一种高分辨高效率的无损探测技术。探地雷达通过天线向地下发射高频电磁脉冲波 ,电磁波在地下介质传播过程中 ,当遇到存在电性差异的地下目标体,如空洞和分界面时,电磁波便会发生反射,返回到地面时由接收天线所接收。在对接收到的雷达波信号处理和分析的基础上,根据信号的波形、振幅和双程走时等参数便可推断地下目标体的空间位置、结构、电性及几何形态,从而达到对地下隐蔽目标体的探测目的。
信号处理是探地雷达技术中的研究重点之一, 其目的是以高的分辨率在探地雷达显示设备上显示反射波图像,提取反射波的振幅、相位和频率等各种有用的参数,帮助解释地质结构信息。
2固有模态函数
由于大多数信号或数据不是固有模态函数, 在任意时刻数据可能包含多个振荡模式, 这也解释了为什么简单的 Hilbert 变换不能给出一个普通信号的频率内容的完整描述。所以必须把数据分解成固有模态函数,从物理上定义一个有意义的瞬时频率的必要条件是:函数对称于局部零均值,且有相同的极值和过零点。据此,Huang 提出了固有模态函数的定义。一个固有模态函数是满足如下两个条件的函数:
(1在整个数据序列中,极值点的数量与过零点的数量必须相等,或最多相差不能多于一个。
(2在任一时间点上,信号的局部极大值和局部极小值定义的包络平均值为零。
第一个限定条件是非常明显的;它近似于传统的平稳高斯过程关于窄带的定义。第二个条件是一个新的想法;它把传统的全局限定变为局部限定。这种限定是必须的,它可去除由于波形不对称而造成的瞬时频率的波动。采用固有模态函数(以下简称 IMF这个名称是因为它代表了信号数据中的振荡模式。IMF 在按过零点定义的每一个周期中,只包括一个本征模态的振荡,没有复杂的叠加波存在。如此定义,一个基本的 IMF 并不限定为窄带信号,也可以是幅度调制和频率调制的。事实上,它可以是非平稳的。图 1是一个典型的 IMF 。固有模态函数(IMF概念的提出使得用 Hilbert 变换定义的瞬时频率具有实际的物理意义, 而提出 IMF 分量的 EMD 分解方法的出现则使瞬时频率可用于复杂的非平稳信号的分析。图 1所示为一典型的固有模态函数,具有相同数目的过零点和极值点,上下包络关于零值对称。
图 1一个典型的固有模态函数(Huang
3经验模态分解
Huang 认为只有对 IMF 分量求出的瞬时频率才有实际的物理意义, 但是大多数信号不是 IMF 分量, 任何时刻信号中可能包含不只一个 IMF 分量。因此, 必须把信号分解为 IMF 分量。为此,Huang 提出了把信号分解为 IMF 分量的算法—EMD,其具体步骤如下:
设时间序列信号为 X(t,它的上、下包络线分别为和 ,则上、下包络的平均 ( u t ( v t 曲线为:
( m t (3-1
1( ( (]2
m t u t v t =+用 x(r减去 m(t后剩余部分 h 1(t,即:(3-2
1( ( ( h t X t m t =?根据上面的定义,在理论上, 满足:(1极值点(极大值或极小值数目与跨零点数 1( h t 目相等或最多相差一个,(2由局部极大值构成的上包络和由局部极小值构成的下包络的平均值为零;即应该是 IMF。实际上,由干包络线样条逼近的过冲和俯冲作用,会产生新 1( h t 的极值影响原来极值的位置与大小;因此,分解得到的并不完全满足 IMF 条件。
1( h t 用代替 ,与相应的上、下包络线为和 ,重复过程,即:
1( h t ( X t 1( h t 1( u t 1( v t (3-3
1111( ( (]2m t u t v t =+(3-4
211( ( ( h t h t m t =?(3-5
1111( [( (]2k k k m t u t v t ???=+(3-611( ( (
k k k h t h t m t ??=?直到所得的满足 IMF 条件:(1极值点数目与过零点数目相等或最多相差一个, (2( k h t 由局部极大值构成的上包络和由局部极小值构成的下包络的平均值趋近于零。这样就分解得
第一个 IMF, 和信号的剩余部分为 ,即:
1( C t 1( r t (3-7111( (
( ( (
k C t h t r t X t C t ==?对信号的剩余部分 r 1(t继续进行 EMD 分解,直到所得的剩余部分为一单调信号或其值
小于预先给定的值时,分解完毕。最终分解得到所有的 IMF 及余量:(3-8
12( ( ( ..... ( ( n n X t C t C t C t R t =++++如前所述,EMD 分离的本质是筛选。满足 IMF 的第一个条件,可以消除附加波的影响; 而满足第二个条件常常是难以做到的,需要确定一个标准使得这一分离过程能够停下来。 Huang 等提出通过限制标准差 S 的大小来确定,即:(3-9
221(1 11(1 1
[( ( /(]n k k k k S h t h t h t ??==?∑其中:S值定在 0.2到 0.3之间。而本文是用分离结果的上包络和下包络的均值是否小于给定的小数值,来确定是否终止 EMD 分离过程。
事实上,EMD 分离终止标准取的不同,分离出的 IMF 的个数和振幅也各异。
图所示为利用 EMD 方法对某一探地雷达 A-scan 的分解
图 2原始 A-scan 信号050100150200250 300350400450500
采样点幅度 /v 50100150200250300350400450500
IMF 0 50100150200250300350400450500
IMF 0
50100150200250300350400450500
IMF 050
100150200250300350400450500
图 3EMD 分解 A-scan 所得的七个 IMF 分量和残余分量 r7(t
4基于经验模态分解的自适应去噪算法
在受干扰背景下有效地检测信号, 不仅与信号的形式和干扰的性质有关, 也与信号处理的方法有关,对不同类型的信号寻找最佳的处理技术一直是信号处理及检测的主要问题之一。
小波变换通过小波基的伸缩和平移, 实现了信号的时频分析局部化它能够同时保留信号的时域特征和频域特征。由于其多辨特性,在合适的尺度下,非平稳信号中的有效成分会呈现出同噪声截然不同的特性, 利用信号和噪声在多尺度空间中不同的传递特性可以获得干扰背景下信号的有效检测, 这种处理信号的方法, 在获得信噪比增益的同时能够保持对突变信息的良好分辨,在非平稳信号的处理中有自身的优越性。
Hilbert-Huang 变换是最新发展起来的处理非线性非平稳信号的时频分析方法。 Hilbert-Huang 变换吸取了小波变换多分辨的优势, 同时又克服了在小波变换中需要选择小
1-3IMF
IMF
c2
IMF c3
501001502002503003504004505000
IMF c4
波基的困难,因此该方法同样可以用来对非平稳信号进行滤波和去噪 [29]。由于从信号木身的尺度特征出发对信号进行分解,该方法具有良好的局部适应性,加上瞬时频率的引入, 使得可以从时频两方而同时对信号进行分析,增加了处理信号的灵活性和有效性。
基于 Hilbert-Huang 变换的去噪算法如下:首先用经验模态分解方法 EMD (Empiricalmode composition method 获得有限数目的固有模态函数
IMF(Intrinsicmode function, 然后对其中的高频分量进行阈值处理, 处理后的分量叠加便得到去噪后的信号。
实验数据为比利时皇家军事学院(Royal Military Academy所采集,原始数据的
A-scan 有 512个采样点。根据上述思路对实测的某一 A-scan 进行处理并和小波阈值法的处理结果进行了比较。图为两种不同的算法处理后的对比结果:
图4经验模态分解和小波阈值去噪比较
表1是原始信号在小波阈值法和经验模态分解法去噪后得到的均方根误差与信噪比。表 1
两种去噪方法的信噪比(SNR和均方误差(RMSE比较通过以上定量比较可以得出结论, 基于经验模态分解的去噪方法无论在信噪比或均方误差方面均优于小波阈值法(采用sym6小波去噪。
5结论
可以看出无论信噪比还是均方误差经验模态分解法均优于小波阈值法。通过以上例子可估计器
小波阈值法经验模态分解法 SNR
22.668526.7431RMSE 2.1723e-0081.3589e-008
050100150200250300
350400450500-0.02
0.0
2
采样点幅度 /v050100150200250300
350400450500-0.02
0.0
2
采样点幅度 /v050100150200250300
350400450500
-0.020.0
2
采样点幅度 /v
以得出结论基于经验模态分解的去噪方法完全可用于探地雷达的信号去噪处理，而其效果要好于传统的小波阈值方法。参考文献 [1] 孔令讲，浅地层探地雷达信号处理算法的研究[D]，四川成都,电子科技大学，2003. [2] D. J. Daniels, D. J. Gunton, H. F. Scott. Introduction to subsurface radar. IEEE Proceeding, 1998, 133(4, 278~326. [3] 粟毅黄春琳雷文太著探地雷达理论与应用[M] 科学出版社 2006 [4] D.J.Daniels.Surface Penetrating Radar. Published by the Institution of Electrical Engineers, London, United Kingdom, 1996 [5] Vsevolod Kovalenko, Advanced GPR Data Processing Algorithms for Detection of Antipersonnel Landmines [D], Technology
University of Delft, Netherlands, 2006 [6] Harkin Brunzell.Signal Processing Techniques for Detection of Buried Landmine Using Ground Penetrating Radar [D], Sweden, 1998 [7] K. S. Yee. Numerical solution of initial boundary value problems involving Maxwell’s equations initial boundary value problems involving Maxwell equations in isotropic media
现代雷达信号处理技术及发展趋势 摘要：自二战以来，雷达就广泛应用于地对空、空中搜索、空中拦截、敌我识别等领域，后又发展了脉冲多普勒信号处理、结合计算机的自动火控系统、多目标探测与跟踪等新的雷达体制。随着科技的不断进步，雷达技术也在不断发展，现代雷达已经具备了多种功能，如反隐身、反干扰、反辐射、反低空突防等能力，尤其是在复杂的工作环境中提取目标信息的能力不断得到加强。例如，利用雷达系统中的信号处理技术对接收数据进行处理不仅可以实现高精度的目标定位与跟踪, 还能够在目标识别和目标成像、电子对抗、制导等功能方面进行拓展, 实现综合业务的一体化。 一、雷达的起源及应用 雷达，是英文Radar的音译，源于radio detection and ranging的缩写，意思为"无线电探测和测距"，即用无线电的方法发现目标并测定它们的空间位置。因此，雷达也被称为“无线电定位”。雷达是利用电磁波探测目标的电子设备。雷达发射电磁波对目标进行照射并接收其回波，由此获得目标至电磁波发射点的距离、距离变化率（径向速度）、方位、高度等信息。雷达最为一种重要的电磁传感器，在国防和国民经济中应用广泛，最大特点是全天时、全天候工作。雷达由天线、发射机、接收机、信号处理机、终端显示等部分组成。 雷达的出现，是由于二战期间当时英国和德国交战时，英国急需一种能探测空中金属物体的雷达（技术）能在反空袭战中帮助搜寻德国飞机。二战期间，雷达就已经出现了地对空、空对地（搜索）轰炸、空对空（截击）火控、敌我识别功能的雷达技术。二战以后，雷达发展了单脉冲角度跟踪、脉冲多普勒信号处理、合成孔径和脉冲压缩的高分辨率、结合敌我识别的组合系统、结合计算机的自动火控系统、地形回避和地形跟随、无源或有源的相位阵列、频率捷变、多目标探测与跟踪等新的雷达体制。后来随着微电子等各个领域科学进步，雷达技术的不断发展，其内涵和研究内容都在不断地拓展。雷达的探测手段已经由从前的只有雷达一种探测器发展到了红外光、紫外光、激光以及其他光学探测手段融合协作。
脉冲压缩雷达的仿真脉冲压缩雷达与匹配滤波的MATLAB仿真 姓名：-------- 学号：---------- 2014-10-28 西安电子科技大学
一、 雷达工作原理 雷达，是英文Radar 的音译，源于radio detection and ranging 的缩写，原意为"无线电探测和测距"，即用无线电的方法发现目标并测定它们的空间位置。因此，雷达也被称为“无线电定位”。利用电磁波探测目标的电子设备。发射电磁波对目标进行照射并接收其回波，由此获得目标至电磁波发射点的距离、距离变化率（径向速度）、方位、高度等信息。 雷达发射机的任务是产生符合要求的雷达波形（Radar Waveform ），然后经馈线和收发开关由发射天线辐射出去，遇到目标后，电磁波一部分反射，经接收天线和收发开关由接收机接收，对雷达回波信号做适当的处理就可以获知目标的相关信息。 但是因为普通脉冲在雷达作用距离与距离分辨率上存在自我矛盾，为了解决这个矛盾，我们采用脉冲压缩技术，即使用线性调频信号。 二、 线性调频（LFM ）信号 脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率，保证足够大的作用距离；而接受时采用相应的脉冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分辨率之间的矛盾。 脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation ）信号,接收时采用匹配滤波器（Matched Filter ）压缩脉冲。 LFM 信号的数学表达式： （2.1） 其中c f 为载波频率，()t rect T 为矩形信号： （2.2）
地球物理与空间信息学院应用地球物理系 毕业实习报告 题目：探地雷达实习报告 姓名：胡浩 班级：061071-22 学号：20071002609 指导教师：邓世坤 二○一一年四月二十二日
前言 探地雷达是利用超高频脉冲电磁波探测地下介质分布的一种地球物理勘探方法。实践证明，它可以分辨地下1m尺度的介质分布，因此探地雷达方法以其特有的高分辨率在浅层于超浅层地质调查中有着极其广阔的应用前景。 探地雷达利用一个天线发射高频宽带电磁波，另一个天线接收来自地下介质界面的反射波。电磁波在介质中传播的时，其路径、电磁场强度于波形将随所通过介质的电性质及几何形态而变化。因此，根据接收到的波的旅行时间、幅度、与波形资料，可推断介质的结构。 第一章探地雷达的探测原理 探地雷达探测是一种快速、连续、非接触电磁波探测技术，具有采集速度快、分辨率高的特点。探地雷达向地下发送脉冲形式的高频宽带电磁波，电磁波在地下介质传播的过程中，当遇到存在电性目标体时，如空洞、分界面时，电磁波便会发生反射，返回到地面时由接收天线所接收；对接收到的电磁波进行信号处理与分析，根据信号波形、强度、双程走时等参数来推断地下目标体的空间位置、结构、电性及几何形态，从而达到对地下隐蔽目标体的探测。 如图A所示，由发射天线向地下介质中发射一定中心频率的电磁脉冲波，电磁波在地下介质中传播时，遇到介质中的电磁性（电阻率、介电率及磁导率）差异分界面发生反射和透射等现象；被反射的电磁波传回地表，由接收天线接收；通过电脑进行操作和控制；接收天线所接收的地下反射回波信号经由光纤传输到仪器控制台，转换成时间序列信号；这种时间序列即构成每一测点上的雷达波形记录道，它包含该测点处所接收到的雷达波的幅度、相位及旅行时间等信息。由电脑收集并存储每一测点上雷达波形序列，形成一个由若干记录道构成的雷达剖面（见图B）。通过对地质雷达剖面进行处理与推断解释，便可获得探测剖面线下方有关的地质特征与信息（或地下目标体的内部结构特征信息）。
雷达数字信号处理解决方案 1.背景 数字信号处理是现代通信、雷达和电子对抗设备的重要组成部分。在实际应用中，利用数字信号处理技术对接收数据进行处理，不仅可以实现高精准的目标定位和目标跟踪，还能够将目标识别、目标成像、精确制导、电子对抗等功能进行拓展，实现多种业务的一体化集成。 在现代雷达系统中，随着有源相控阵和数字波束形成（DBF）技术的广泛应用，接收前端存在大量的数据需要并行处理，并需要保证高性能和低延迟的特点。雷达日益复杂的应用环境，让雷达系统具备自适应于探测目标和环境的能力，数字信号处理部分也需要使用多种更加复杂的算法，并且可以做到算法模块化，以及通过软件配置功能模块的参数，实现软件定义的功能。更大的数据处理带宽能够使雷达获得更高的分辨率，更高的工作频率使得雷达可以小型化，能够在更小的平台上安装，这样对于硬件平台实现也有低功耗的要求。 在电子对抗设备中，可以在最短的时间内对多个威胁目标进行快速分析和响应，同样需要数字信号处理的相关算法具备高实时，高动态范围和自适应的特点。如何在宽频噪声的环境中寻找到目标的特征数据，如何在宽带范围内制造虚假目标实现全覆盖，数字信号的处理性能是至关重要的设计因素。 加速云的SC-OPS和SC-VPX产品，针对5G通信和雷达的数字信号处理的要求，结合Intel最新14nm 工艺的Stratix10 FPGA系列，提供了一套完整的硬件和软件相结合的解决方案。SC-OPS产品作为单独的硬件加速卡，通过PCIe插卡的方式实现与主机的通信功能，还可以通过多卡级联的方式实现数字信号的分布式处理方案。SC-VPX产品是由FPGA业务单板、主控板和机箱组成的VPX系统。借助于FPGA可编程的特性，加速云提供了高性能数学加速库FBLAS和FFT的RTL级IP，具有高性能和算法参数可配置的特点实现了多重信号分类（MUSIC）和自适应数字波束形成（ADBF）的核心算法，提高了5G通信和雷达在对抗干扰方面的性能。为了方便客户使用高层语言开发，加速云提供基于FPGA完整的OpenCL异构开发环境，快速实现用户自定义的信号处理加速方案。 图1. 加速云SC-OPS和SC-VPX产品
基本雷达信号处理流程 一、脉冲压缩 窄带（或某些中等带宽）的匹配滤波： 相关处理，用FFT 数字化执行，即快速卷积处理，可以在基带实现（脉冲压缩） 快速卷积，频域的匹配滤波 脉宽越小，带宽越宽，距离分辨率越高 ； 脉宽越大，带宽越窄，雷达能量越小，探测距离越近； D=BT （时宽带宽积）； 脉压流程： 频域：回波谱和参考函数共轭相乘 时域：相关 即输入信号的FFT 乘上参考信号FFT 的共轭再逆FFT ； Sc=ifft(fft(Sb).*conj(fft(S))); FFT 输入信号 共轭相乘逆FFT 参考信号的FFT 匹配滤波器 输出 Task1 f0=10e9;%载频tp=10e-6;%脉冲宽度B=10e6;%信号带宽fs=100e6;%采样率 R0=3000;%目标初始距离N=4096;c=3e8;tau=2*R0/c;beita=B/tp;t=(0:N-1)/fs; Sb=rectpuls(t-tp/2-tau,tp).*exp(j*pi*beita*(t-tp/2-tau).^2).*exp(-2j*pi*f0*tau);%回波信号 1000 2000 3000 4000 5000 6000 7000 -1-0.8-0.6-0.4-0.200.20.40.60.81 1000 2000 3000 4000 5000 6000 7000 -1-0.8-0.6-0.4-0.200.20.40.60.81 012345678910 x 10 7 20 40 60 80 100 120
S=rectpuls(t-tp/2,tp).*exp(i*pi*beita*(t-tp/2).^2);%发射信号（参考信号） 0.5 1 1.5 2 2.5 3 3.5 4 4.5x 10 -5 -1-0.8-0.6-0.4-0.200.20.40.60.81 0.5 1 1.5 2 2.5 3 3.5 4 4.5x 10 -5 -1-0.8-0.6-0.4-0.200.20.40.60.81 012345678910x 10 7 20 40 60 80 100 120 So=ifft(fft(Sb).*conj(fft(S)));%脉压 figure(7); plot(t*c/2,db(abs(So)/max(So)))%归一化dB grid on 01000200030004000500060007000 -400 -350-300-250-200-150-100-500
雷达信号处理技术与系统设计 第一章绪论 1.1 论文的背景及其意义 近年来，随着电子器件技术与计算机技术的迅速发展，各种雷达信号处理技术的理论与应用研究成为一大热门领域。 雷达信号的动目标检测(MAD)是利用动目标、地杂波、箔条和气象干扰在频谱上的差别，抑制来自建筑物、山、树、海和雨之类的固定或低速杂波信号。区分运动目标和杂波的基础是它们在运动速度上的差别，运动速度不同会引起回波信号频率产生的多普勒频移不相等，这就可以从频率上区分不同速度目标的回波。固定杂波的中心频率位于零频，很容易设计滤波器将其消除。但对于运动杂波，由于其多普勒频移未知，不能像消除固定杂波那样很容易地设计滤波器，其抑制就变得困难了从本质上来讲，雷达信号的检测问题就是对某一坐标位置上目标信号“有”或“无”的判断问题。最初，这一任务由雷达操作员根据雷达屏幕上的目标回波信号进行人工判断来完成。后来，出现了自动检测技术，一开始为固定或半固定门限检测，这种体制下当干扰和杂波功率水平增加几分贝，虚警概率将急剧增加，以至于显示器画面饱和或数据处理过载，这时即使信噪比很大，也不能作出正确的判断。为克服这些问题进而发展了自适应恒虚警(Constant FalseAlarm Rate，CFAR)检测。CFAR 检测使得雷达在多变的背景信号中能够维持虚警概率的相对稳定，这种虚警概率的稳定性对于大多数的雷达，如搜索警戒雷达、跟踪雷达、火控雷达等。
第二章 雷达信号数字脉冲压缩技术 2.1 引言 雷达脉冲压缩器的设计实际上就是匹配滤波器的设计。根据脉冲压缩系统实 现时的器件不同，通常脉冲压缩的实现方法分为两类，一类是用模拟器件实现的 模拟方式，另一类是数字方式实现的，主要采用数字器件实现。 脉冲压缩处理时必须解决降低距离旁瓣的问题，否则强信号脉冲压缩的旁瓣 会掩盖或干扰附近的弱信号的反射回波。这种情况在实际工作中是不允许的。采 用加权的方法可以降低旁瓣，理论设计旁瓣可以达到小于-40dB 的量级。但用模拟技术实现时实际结果与理论值相差很大，而用数字技术实现时实际输出的距离旁瓣与理论值非常接近。数字脉压以其许多独特的优点正在或已经替代模拟器件进行脉冲压缩处理。 2.2 数字脉压实现方法 用数字技术实现脉冲压缩可采用时域方法或频域方法。至于采用哪种方法。 要根据具体情况而定，一般而言，对于大时宽带宽积信号，用频域脉压较好；对 于小时宽带宽积信号，用时域脉压较好。 2.2.1 时域卷积法实现数字脉压 时域脉冲压缩的过程是通过对接收信号)(t s 与匹配滤波器脉冲响应)(t h 求卷积的方法实现的。根据匹配滤波理论，)()(0*t t s t h -=，即匹配滤波器是输入信号的共轭镜像，并有响应的时移0t 。 用数字方法实现时，输入信号为)(n s ，起匹配滤波器为)(n h ，即匹配滤波器的输出为输入离散信号)(n s 与其匹配滤波器)(n h 的卷积
一种雷达通用信号处理系统的实现与应用 一种雷达通用信号处理系统的实现与应用 FPGA是一种现场可编程器件，设计灵活方便可以反复修改内部逻辑，适用于算法结构比较简单、处理速度较高的情况。DSP是一种基于指令集的处理器，适于大信息、复杂算法的信息处理场合。鉴于两种处理器件自身优势，FPGA+DSP信号处理架构，已成为信号处理系统的常用结构。但当前FPGA+DSP的信号处理平台或者是基于某些固定目的，实现某些固定功能，系统的移植性、通用性较差。或者仅仅简要介绍了平台的结构没有给出一些具体的实现。本文提出的基于FPGA+DSP通用信号处理平台具有两种处理器的优点，兼颐速度和灵活性，而且可以应用在不同雷达信号处理系统中，具有很强的通用性。本文举例说明该系统在连续波雷达和脉冲雷达中的典型应用。1系统资源概述1．1处理器介绍本系统FPGA选择Altera公司的EP2S60F1020。Stratix II FPGA采用TSMC的90nm 低k绝缘工艺技术。Stratix II FPGA支持高达1Gb·s-1的高速差分I／O信号，满足新兴接口包括LVDS，LNPECL和HyperTransport标准的高性能需求，支持各种单端I／O接口标准。EP2S60系列内部有48352个ALUT；具有2544192bit的RAM 块，其中M512RAM(512bit)329个，M4K RAM(4kbit)255个，M-RAM(512kbit)2个。具有嵌入式DSP块36个，等效18bit×18bit乘法器144个；具有加强型锁相环EPLL4个，
快速锁相环FPLL8个。这些锁相环具有高端功能包括时钟切换，PLL 重新配置，扩频时钟，频率综合，可编程相位偏移，可编程延迟偏移，外部反馈和可编程带宽等。本系统DSP选择ADI公司的ADSP TS201。它有高达600MHz的运行速度，1．6ns的指令周期；有24MB的片内DRAM；双运算模块，每个计算块包含1个ALU，一个乘法器，1个移位器，1个寄存器组和1个通信逻辑单元(CLU)；双整数ALU，提供数据寻址和指针操作功能；集成I／O接口，包括14通道的DMA控制器，外部端口，4个链路口，SDRAM控制器，可编程标识引脚，2个定时器和定时器输出引脚等用于系统连接；IEEE1149．1兼容的JTAG端口用于在线仿真；通过共享总线可以无缝连接多达8个TigerSHARC DSP。1．2FPGA+DSP结构由于FPGA和DSP各自的自身优势，FPGA+DSP信号处理架构已成为信号处理系统的常用结构。一般情况下FPGA+DSP的拓扑结构会根据需要进行不同的连接，这就导致这种结构的专用性，缺乏灵活性。对于一个通用处理平台要考虑到各种不同的信号通路，因此大部分通用FPGA+DSP平台都采取各个处理器间均有通路的方式。这种拓扑结构灵活方便，可以满足各种不同的通路需求，这种结构的缺点就是硬件设计的复杂以及可能会有资源浪费。对于这种通用FPGA+DSP 结构，FPGA与各个DSP之间均有连接，不同之处便是DSP之间的拓扑结构。一般分两种，一是高速外部总线口耦合结构组成多DSP 系统，这种结构可以实现多DSP共享系统内的资源，系统内的个处理器可以共享RAM，SDRAM和主机等资源，还可共享其他处理器核
ReflexW读取GSSI劳雷雷达数据和处理方法 Dr.Zhang/ 2020.2.21 1雷达数据 1.1雷达数据文件 探地雷达的采集原理类似于地震，老一点的数据格式一般跟地震一样都是segy。后来 GSSI 公司的 TerraSIRch SIR3000地质雷达系统(简称SIR-3000) 在国内普及起来后，格式dzt就成了主流格式（美国SIR系列探地雷达的数据格式）。现在最新的情况是国内生产单位还是美国最多，其次是加拿大和意大利的。 国内其他常用的GPR数据还包括：DT格式（意大利RIS系列探地雷达的数据格式），DT1格式（加拿大Pulse-Ekko 系列探地雷达的数据格式）以及RD3格式(瑞典MALA系列探地雷达的数据格式)。 1.2 雷达文件dzt数据 目前主要是劳雷GSSI SIR 3000和SIR 4000的数据格式。 2REFLEXW软件简介 2.1 REFLEXW软件简介 本文基于MATLAB开发的GUI界面，主要实现了REFLEXW软件的相应功能。 REFLEXW是地质雷达数据（类地震）数据处理及解释软件，应用于地质雷达的数据处理以及资料解释。Reflexw软件兼容了世界上大多数雷达的数据格式，在欧美地区，Reflexw已经成为了地质雷达数据处理的标准软件。随着地质雷达行业的发展，在国内也越来越多的人开始使用Reflexw软件。 软件特点：功能强大，可做多种滤波处理可对个雷达数据进行批量处理导入GPS数据，可绘制测线轨迹、修正地形可显示测线中的标记对不需要的雷达数据可进行删除可做2D剖
面处理和3D时间切片处理。 2.2 REFLEXW软件功能分析 通过REFLEXW软件，可以实现对探地雷达数据的读入，一维滤波，校正，二维滤波，波形图观察等功能，REFLEXW提供了较为全面的滤波手段，可以将探地雷达图谱处理的更加容易观察 在REFLEXW软件中，在显示数据方面，也同样提供了大量的处理方式。其中，最主要的处理方式为Plot Options。在Plot Options选项中，包括了绘图模式，点模式比例，能量衰减和振幅比例。其中，点模式比例又分为了XY比例绘图，每样点像素，每道像素，能量衰减，振幅比例。以上的这些选项，在用户导入探地雷达图谱之后，可以对图谱进行一个基本的预处理，或者通过不同的方式来观察图谱。 3 读取.dzt文件 3.1 建立预处理的工程文件 雷达信号分析时会产生一系列的文件，因此最好为每个雷达文件建立一个工程，并指定其存放位置。ReflexW会自动生成相应的文件。 图3.1 建立工程(project -选择文件-confirm)
连续波雷达及信号处理技术初探 摘要：连续波雷达，主要就是连续发生电磁波的雷达，可以根据不同发射信号的形式，将其划分成为非调制单频与调频两种类型。在连续波雷达系统实际应用的过程中，应当科学使用信号处理技术开展相关处理工作，在实际观测的过程中，解决收发开关中存在的问题，保证雷达信号接收与发射工作效果。关键词：连续波雷达；信号处理技术；应用措施在使用信号处理技术对连续波雷达进行控制的过程中，应当建立多元化的管理机制，明确各方面工作要求，创新信号处理工作形式，保证能够提升信号处理技术的应用水平，创建专门的管理机制。一、连续波雷达定义与特征分析对于连续波雷达而言，主要是针对电磁波进行连续的发射，根据发射信号形式将其划分成为非调制单频与调频两种类型。在1924年的时候，英国就开始通过连续波课调频测距相关分析，对电离层开展观测工作。且在第二次世界大战的过程中，已经使用连续波雷达开展飞机观测与地面观测工作。然而，在实际使用的过程中，经常会出现收发隔离的现象，难以保证工作效果，因此，使用收发开关对此类问题进行了解决。当前，在使用连续波雷达的过程中，已经能够通过同一天线开展信号接收与发射工作，产生良好的工作效果。在使用连续波雷达发射机设备的过程中，不需要高压的支持，也不会出现打火的现象，能够利用多元化的方式开展信号调制工作，有利于提升信号的发射效率，增强雷达处理效果，因此，在相同体积、重量的雷达设备中，连续波雷达受到广泛关注与重视，应用于世界的各个国家。同时，连续波雷达的体积很小，重量很轻，馈线的损耗最低，使用流程简单，与其他雷达相较可以得知，连续波雷达在接收机方面，所使用的宽带脉冲较窄，有利于抵抗杂波问题，提升电磁干扰的抵抗能力。在应用连续波雷达对距离与速度进行测量的过程中，其测量准确性较高，不会受到其他因素的干扰。对于连续波雷达而言，其特点主要表现为以下几点： （一）发射机的运行功率较低连续波雷达的发射机运行功率很低，有利于应用在侦查工作中。一般情况下，在使用侦查接收机的过程中，可以利用连续波雷达对其进行处理，提升工作效率，加快侦查速度，保证瞬时频率符合相关规定。同时，在使用连续波雷达的过程中，还要使用伪随机码调相方式对其进行处理，减少外界带来的干扰，做好反侦察工作，保证可以符合实际发展需求。（二）接收机的宽带很窄连续波雷达在实际运行的过程中，接收机的宽带很窄，在杂波环境中，能够实行检测工作，提升自身抗干扰能力。且在电磁干扰的环境中，可以提升自身的抗干扰性能，满足实际处理需求[1]。（三）对小目标进行检测连续波雷达设备的使用，可以提升发射机的功率，增加收发天线的收益，且可以减少噪音问题，在一定程度上，能够减少微波损耗问题，更好的对隐身目标进行检测，合理开展雷达探测等工作，提升相关信号的处理效果，满足实际发展需求。二、连续波雷达的相关工作园林分析连续波雷达的运行，需要明确实际工作原理，通常情况下，雷达发射线性三角调频的相关连续性信号，那么，雷达设备的载频就在f0的数值之上，在此过程中，可以将调频宽带设置成为A，将调频间隔设置成为C。在对信号频率与时间进行计算的过程中，应当明确相关原理，创新管理工作形式，对具有代表性的内容进行合理分析，保证可以提升自身分析工作效果。在信号处理工作中，应当重点关注发射信号与目标回波信号，通过合理的计算方式，创建多
雷达系统中的信号处理技术 摘要本文介绍了雷达系统及雷达系统信号处理的主要内容，着重介绍与分析了雷达系统信号处理的正交采样、脉冲压缩、MTD和恒虚警检测几种现代雷达技术，雷达系统通过脉冲压缩解决解决雷达作用距离和距离分辨力之间的矛盾，通过MTD来探测动目标，通过恒虚警（CFAR）来实现整个系统对目标的检测。 关键词雷达系统正交采样脉冲压缩MTD 恒虚警检测 1雷达系统概述 雷达是Radar（Radio Detection And Ranging）的音译词，意为“无线电检测和测距”，即利用无线电波来检测目标并测定目标的位置，这也是雷达设备在最初阶段的功能。雷达的任务就是测量目标的距离、方位和仰角，还包括目标的速度，以及从目标回波中获取更多有关目标的信息。典型的雷达系统如图1，它主要由雷达发射机、天线、雷达接收机、收发转换开关、信号处理机、数据处理机、终端显示等设备组成。 图1雷达系统框图
随着现代电子技术的不断发展，特别是数字信号处理技术、超大规模集成数字电路技术、计算机技术和通信技术的告诉发展，现代雷达信号处理技术正在向着算法更先进、更快速、处理容量更大和算法硬件化方向飞速发展，可以对目标回波与各种干扰、噪声的混叠信号进行有效的加工处理，最大程度低剔除无用信号，而且在一定的条件下，保证以最大发现概率发现目标和提取目标的有用信息。 雷达发射机产生符合要求的雷达波形，然后经馈线和收发开关由发射天线辐射出去，遇到目标后，电磁波一部分反射，经接收天线和收发开关由雷达接收机接收，然后对雷达回波信号依次进行信号处理、数据处理，就可以获知目标的相关信息。 雷达信号处理的流程如下： 图 2 雷达信号处理流程 2雷达信号处理的主要内容 雷达信号处理是雷达系统的主要组成部分。信号处理消除不需要的杂波，通过所需要的目标信号，并提取目标信息。内容包括雷达信号处理的几个主要部分：正交采样、脉冲压缩、MTD和恒虚警检测。 正交采样是信号处理的第一步，担负着为后续处理提供高质量数据的任务。采样的速率和精度是需要考虑的首要问题，采样系统引起的失真应当被限定在后续信号处理任务所要求的误差范围内，直接中频数字正交采样是当代雷达的主要技术之一。脉冲压缩技术在现代雷达系统中得到了广泛的应用。脉冲压缩雷达既能保持窄脉冲雷达的高距离分辨力，又能获得脉冲雷达的高检测力，并且抗干扰能力强。现在，脉冲压缩雷达使用的波形正在从单一的线性调频发展到时间、频率、编码混合调制，在尽可能不增加整机复杂度的条件下实现雷达性能的提升。杂波抑制是雷达需要具备的重要功能之一。动目标指示与检测是通过回波多普勒频移的不同来区分动目标和固定目标，通过设计合理的滤波器（组），就可以把目标号和杂波分开。
浅谈探地雷达的原理与特点 摘要：地下管线系统的建立是城市现代化建设的重要因素，但由于地下管线中的非金属管线的大量存在以及城市建设快速安全的需要，探地雷达探测技术的独特优势就显现出来，本文通过对探底雷达和地下管线的分析，为应用探底雷达在城市地下管线建设提供参考。 关键字：探地雷达；地下管线；探测技术 0 引言 随着城市现代化的发展，地下管线的密集程度也在不断地扩大。地下管线作为城市的重要基础设施之一，它一方面关系着城市居民生活及城市工业的发展，担负着巨大的社会责任，另一方面又由于它深埋于地下，具有不透明性，纵横交错、结构复杂。近年来，在许多大城市出现施工时挖断通信、电力电缆导致通讯中断、区域性停电、停产事故，这些事故给该地区经济和人们的生产生活带来了巨大的损失。因此，地下工程在施工时如何避免破坏这些地下管线就变得越来越重要，建立完整的城市地下管线系统成为现代城市快速建设的关键因素。 探地雷达(Ground Penetrating Radar，简称GPR)是一种对地下或物体内不可见的目标或界面进行定位的电磁法，并以其探测的高分辨率和高效率而成为地球物理勘探的一种有力工具。随着信号处理技术和电子技术的不断发展及工程实践的增多和经验的不断积累，探地雷达技术也得到极大提高，仪器也不断更新，探地雷达检测技术具有分辨率高、采集速度快、后期数据处理简便等特点。因此在铁路、公路、建筑、市政、考古等领域得到广泛的应用,并受到广大现场技术人员的认可和喜爱。 1 探地雷达的发展 国外探底雷达技术最早可追溯到二十世纪初，西方国家以专利形式提出将雷达原理用于探地，正式提出了探地雷达的概念。但是直到50年代后期探地雷达技术才被慢慢重视起来。探地雷达在矿井、冰层厚度、地下粘土属性、地下水位等方面的得到了应用。1967年，一个与Stern最初用于冰川探测的仪器类似的系统被设计研制出来，1972年，Rex Morcy和Art Drake开创了GSSI公司，主要从事商业探底雷达的销售。随着电子技术的发展，电子存储设备的问世，加之现代数据处理技术的应用，特别是拟反射地震处理的应用，探底雷达的应用领域迅速扩大，其中有：石灰岩地区采石场的探测、淡水和沙漠地区的探测、工程地质探测、煤矿井探测、泥灰调查、放射性废弃物处理调查、埋设物探测、水文地质调查、地基和道路下空洞及裂缝调查、水坝的缺陷检测、隧道及堤坝探测等。 自70年代以来，许多商业化的通用数字探地雷达系统先后问世，其中有代表性的有：美国GSSI公司的SIR和MK系列，加拿大Sensor&Software公司的Pulse Ekko系列。这些雷达的基本原理大同小异主要功能有多通道采集、多维显示、实时处理、变频天线、多次叠加、多波形处理等，另外还有井中雷达系统，
目录 目录.......................................................................................................................... I 表目录............................................................................................................... III 图目录................................................................................................................ I V 摘要 (i) ABSTRACT ................................................................................................................ i i 第一章 绪论 .. (1) 1.1 研究背景及意义 (1) 1.2 国内外研究现状 (3) 1.3 本文主要工作及内容安排 (5) 第二章 探地雷达的基本原理 (7) 2.1 麦克斯韦方程组和波动方程 (7) 2.2 探地雷达的探测方式 (9) 2.3 探地雷达的数据形式 (10) 2.4 探地雷达的分辨率 (11) 2.5 探地雷达的信号模型 (13) 2.5.1 探地雷达A-scan信号模型 (13) 2.5.2 探地雷达B-scan信号模型 (15) 2.6 小目标回波的双曲线特征 (16) 2.6.1 小目标回波双曲线特征分析 (16) 2.6.2 基于双曲线对称性的小目标定位与检测 (17) 2.7 本章小结 (22) 第三章 探地雷达杂波抑制算法研究 (23) 3.1 引言 (23) 3.2 三种常用的杂波抑制算法分析 (24) 3.2.1 均值法 (24) 3.2.2 SaS(Scale and Shift)法 (25) 3.2.3 自适应对消法 (26) 3.3 基于奇异值分解的探地雷达杂波抑制算法研究 (30) 3.3.1 基于矩阵内积的B-scan数据分析 (30) 3.3.2 奇异值分解带通滤波算法理论分析 (33)
雷达信号处理和数据处理技术 定价: ￥89.00元金桥价: ￥84.55元节省: ￥4.45元 内容简介 雷达信号处理和数据处理技术是雷达的神经中枢。信号处理通过对雷达回波信号的处理来发现目标和测定目标的坐标和速度等，形成目标点迹，数据处理通过对目标点迹的处理形成目标的航迹供指挥决策使用。 本书的主要内容包括雷达信号的形式、雷达杂波抑制、雷达脉冲压缩、雷达信号检测、雷达抗干扰、雷达目标识别、雷达点迹处理和雷达航迹处理等。 全书共14章，第1章为概论，第2章到第10章为雷达信号处理技术，第11章到第14章为雷达数据处理技术。全部内容既包含处理理论，也包含设计技术。 本书可以帮助雷达工程技术人员和雷达使用人员掌握有关雷达信号处理和数据处理技术，解决有关应用问题；同时还可以作为高等学校电子工程相关专业高年级本科生和研究生的参考用书。 雷达信号处理基础 定价: ￥55.00元金桥价: ￥52.25元节省: ￥2.75元
内容简介 本书译自国际著名雷达信号处理专家Mark A. Richards教授编写的教科书。该书介绍了雷达系统与信号处理的基本理论和方法，主要内容包括：雷达系统导论、雷达信号模型、脉冲雷达信号的采样和量化、雷达波形、多普勒处理、检测基础原理、恒虚警率检测、合成孔径雷达成像技术、波束形成和空-时二维自适应处理导论。书中包含了大量反映雷达信号处理最新研究成果和当前研究热点的补充内容，提供了大量有助于读者深入的示例。该书对基础理论和方法进行了详尽的介绍与深入严谨的论述，是一本雷达信号处理领域中高水平的教科书。 本书适合于从事雷达成像、检测、数据处理及相关信号处理的研究生作为教材使用，也是相关专业研究人员不可多得的一本参考书。Mark A.Richards。博士，佐治亚理工学院（Georgia Institute of Technology）的首席研发工程师和兼职教授。他具有20余年在学术界、工业界及政府部门从事雷达信号处理和嵌入式计算方面研究的经历。他曾被聘为美国国防高级研究计划署项目经理、IEEE 2001年雷达会议的总主席，以及IEEE图像处理和IEEE信号处理期刊的副编辑。Eichards博士长期从事关于雷达信号处理、雷达图像处理及相关学科的研究生教育和职业教育。这本严谨的著作源自于一位该领域令人尊敬的领导者，它提供了其他文献中所没有的关于雷达DSP基础及其应用的详细内容。对于那些不只想从普通雷达系统的书籍中粗略学习信号处理，还想学到更多关于信号模型、波形、干扰抑制、探测，以及诸如SAR和SFAP等高级雷达信号处理主题的人而言，本书是非常合适的。经过多年研究生和职业教育的完善与检验，这本深入介绍雷达DSP技术的书籍，以现有的先进雷达技术为基础，全面讨论了以下几方面的问题，并提供了详尽的例子：多域信号获取和采样、目标和干扰模型、常见雷达波形、干扰抑制技术、检测算法和工具、合成孔径成像和自适应阵列处理基础。 信息传输与正交函数 定价: ￥28.00元金桥价: ￥26.60元节省: ￥1.40元 内容简介 本书叙述了非正弦正交函数理论和以之为基础的信息传输系统，主要内容包括正交函数系、信息传输的基本思想和方法，移动通信与正交函数之间的关系，沃尔什函数的复制生成理论，一般复制生成理论及桥函数的概念，沃尔什函数及桥函数的相关函数的定义及其特性，序率分割制多路传输系统，信息传输系统的统一模型等。 本书可供从事通信、遥控、遥测和雷达工作的技术人员、科研人员以及高等院校师生参考。 DSP开发应用技术
探地雷达的发展与现状 探地雷达的历史最早可追溯到20世纪初。1904年，德国人Hülsmeyer首次将电磁波信号应用于地下金属体的探测。1910年，Leimback和L?wy以专利形式提出将雷达原理用于探地，他们用埋设在一组钻孔中的偶极天线探测地下相对高导电性质的区域，正式提出了探地雷达的概念。1926年Hülsenbeck第一个提出应用脉冲技术确定地下结构的思路，他指出介电常数不同的介质交界面会产生电磁波反射。由于地下介质具有比空气强得多的电磁衰减特性，加之地下介质情况的多样性，电磁波在地下的传播比空气中复杂的多，之后二三十年尽管在美国出现过一些相关的专利，这项技术很少被运用到其它领域,直到50年代后期，探地雷达技术才慢慢重新被人们所重视。探地雷达在矿井(1960,J.C.Cook)、冰层厚度(1963,S.Evans)、地下粘土属性(1965,Barringer)、地下水位(1966,Lundien)的探测方面得到了应用。1967年，一个与stern最初用于冰川探测的仪器类似的系统被设计研制出来，1972年Procello将其于探测月球表面结构。同样在1972年，Rex Morcy和Art Drake开创了GSSI(Geophysical Survey Systems Inc.)公司，主要从事商业探地雷达的销售。随着电子技术的发展，数字磁带记录问世，加之现代数据处理技术的应用，特别是拟反射地震处理的应用，探地雷达的实际应用范围在70年代以后迅速扩大，其中有：石灰岩地区采石场的探测(1971,Takazi;1973,kithara;)、淡水和沙漠地区的探测(1974,R.M.Morey;1976,P.K.Kadaba)、工程地质探测(1976,A.P.Annan和J.L.Davis;1978,G.R.Olhoeft,L.T.Dolphin)、煤矿井探测(1975,J.C.Cook)、泥炭调查(1982,C.P.F.Ulriken)、放射性废弃物处理调查(1982,D.L.Wright;1985,O.Olsson)、以及地面和井中雷达用于地质构造填图(1997,M.Serzu )、水文地质调查(1996,A.Chanzy ;1997,Chieh-Hou Yang )、地基和道路下空洞及裂缝调查、埋设物探测、水坝的缺陷检测、隧道及堤岸探测等。 自70年代以来、许多商业化的通用数字探地雷达系统先后问世，其中有代表性的有：美国Geophysical Survey System Inc公司的SIR系统、Microwave Associates 的MK系列，加拿大Sensor & Software的Pulse Ekko系列，瑞典地质公司(SGAB)的RAMAC/GPR系列，日本应用地质株式会社OYO公司的GEORADAR系列及一些国内产品(电子工业部LTD系列，北京爱迪尔公司CR-20、CBS-900等)。这些雷达仪器的基本原理大同小异，主要功能有多通道采集、多维显示、实时处理、变频天线、多次叠加、多波形处理等，另外还有井中雷达系统，多态雷达系统，层析成像雷达系统等。 国内探地雷达的研究始于70年代初。当时，地矿部物探所、煤炭部煤科院，以及一些高校和其他研究部门均做过探地雷达设备研制和野外试验工作，但由于种种原因，这些研究未能正式用于实际。90年代以来，由于大量国外仪器的引进，探地雷达得到了广泛的应用与研究。1990-1993年，中国地质大学(武汉)在国家自然科学基金资助下，开展了大量的理论研究和工程实践，取得了不少成果。探地雷达主要应用领域有隧道(1998，隋景峰；2001，刘敦文等)、水利工程设施(1997，赵竹占等)、混凝土基桩(2000，李梁等)、煤矿(1998，刘传孝等)、公路(1996牛一雄等；1997，沈飚等)；岩溶(1994，王传雷，祁明松；1995，李玮，梁晓园)；工程地质(1994，胡晓光；1999，刘红军，贾永刚)；钻孔雷达(1999，宋雷，黄家会)等。
ReflexW读取意大利IDS雷达数据DT和处理方法 Dr.Zhang/ 2020.2.23 1雷达数据 1.1雷达数据文件 .DT格式（意大利RIS系列探地雷达的数据格式） 国内其他常用的GPR数据还包括：，DT1格式（加拿大Pulse-Ekko 系列探地雷达的数据格式）以及RD3格式(瑞典MALA系列探地雷达的数据格式)。GSSI 公司的TerraSIRch SIR3000地质雷达系统(简称SIR-3000) ，格式dzt（美国SIR系列探地雷达的数据格式）。 1.2 雷达文件.DT数据 目IDS的探地雷达（GPR）数据格式。
2REFLEXW软件简介 2.1 REFLEXW软件简介 本文基于MATLAB开发的GUI界面，主要实现了REFLEXW软件的相应功能。 REFLEXW是地质雷达数据（类地震）数据处理及解释软件，应用于地质雷达的数据处理以及资料解释。Reflexw软件兼容了世界上大多数雷达的数据格式，在欧美地区，Reflexw已经成为了地质雷达数据处理的标准软件。随着地质雷达行业的发展，在国内也越来越多的人开始使用Reflexw软件。 软件特点：功能强大，可做多种滤波处理可对个雷达数据进行批量处理导入GPS数据，可绘制测线轨迹、修正地形可显示测线中的标记对不需要的雷达数据可进行删除可做2D剖面处理和3D时间切片处理。 2.2 REFLEXW软件功能分析 通过REFLEXW软件，可以实现对探地雷达数据的读入，一维滤波，校正，二维滤波，波形图观察等功能，REFLEXW提供了较为全面的滤波手段，可以将探地雷达图谱处理的更加容易观察 在REFLEXW软件中，在显示数据方面，也同样提供了大量的处理方式。其中，最主要的处理方式为Plot Options。在Plot Options选项中，包括了绘图模式，点模式比例，能量衰减和振幅比例。其中，点模式比例又分为了XY比例绘图，每样点像素，每道像素，能量衰减，振幅比例。以上的这些选项，在用户导入探地雷达图谱之后，可以对图谱进行一个基本的预处理，或者通过不同的方式来观察图谱。 3 读取.DT文件 3.1 建立预处理的工程文件 雷达信号分析时会产生一系列的文件，因此最好为每个雷达文件建立一个工程，并指定其存放位置。ReflexW会自动生成相应的文件。
附录一——MATLAB信号处理程序 %% 1、准备工作 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 开始 clc; clear; close all; clear vars; %% 雷达波形参数定义及说明 f1=1e3; % 最低频率 f2=11e3; % 最高频率 B=f2-f1; % 信号带宽 T=1e-2; % 信号扫频时宽(10ms) c=3e8; % 电磁波空间传播速度 f0=(f1+f2)/2; % 雷达工作频率（中心频率）(3kHz) fs=1e5; % 采样率(100kHz) N_signal_T=round(fs*T); % 单周期信号的数据点数 number_of_signal_period=400; % 脉冲信号的周期个数 duty_ratio=0.5; % 信号占空比 T_signal=T/duty_ratio; % 脉冲信号周期 %% 导入AD数据时频分析 [FileName,PathName] = uigetfile('C:\Users\XYB\Desktop\课程设计之雷达信号分析处理\AD数据\USB (3).dat','Select the USB.dat file'); f = fullfile(PathName,filesep,FileName); fid = fopen(f,'r'); data = fscanf(fid,'%x'); fclose(fid); data = data(1:2:end)*256 + data(2:2:end); %将16进制转换为10进制 datsgn = data./1000; %单位换算（mV->V） %转化为有符号数（去直流） datsgn=datsgn-mean(datsgn); %时域波形 figure; plot([0:1/fs:(length(datsgn)-1)/fs],datsgn); xlabel('时间/s') ylabel('振幅/V') title('LFMCW时域波形') %频谱图 N=1024; datfft = (2/N)*fftshift(fft(datsgn(1:N))); nordat = abs(datfft)/max(abs(datfft)); %对信号做FFT并归一化 figure; plot([-length(datfft)/2:(length(datfft)/2- 1)].*(fs/N),20*log10(abs(nordat)));