CN110673139B - 一种目标探测雷达系统脉冲压缩系数生成方法及装置 - Google Patents

Abstract

本发明提供一种目标探测雷达系统脉冲压缩系数生成方法及装置，所述方法包括：FPGA主控模块计算出配置参数并输入直接频率合成器；直接频率合成器产生雷达系统要求的线性调频信号，并将其发送至射频前端模块；所述射频前端模块的发射通道将线性调频信号上变频至射频频段，并通过射频线缆环回至射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频并发送至模数转换器；模数转换器将环回信号的数据传输至FPGA主控模块；数据分析模块接收由FPGA主控模块导出的环回DDS数据，运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数。根据本发明的方案，可以提高雷达系统的脉冲压缩性能，获得更高的雷达距离分辨力。

Description

一种目标探测雷达系统脉冲压缩系数生成方法及装置

技术领域

本发明涉及雷达系统探测领域，尤其涉及一种目标探测雷达系统脉冲压缩系数生成方法及装置，特别涉及一种低小慢目标探测雷达系统脉冲压缩系数生成方法及装置。

背景技术

目前，低小慢目标如无人机发展的越来越快，越来越多的低小慢目标投入到军事和民事应用中。对这类低小慢目标的有效探测是当前国防安全和信息安全所急需解决的问题。

雷达系统为了能够获得较高的距离探测精度，需要雷达系统发射宽度足够窄的脉冲，另外雷达系统为了能够探测到远距离的目标，又必须发射能量很大的信号。然而，受限于当前的硬件制造水平，硬件设备极难产生大功率的窄脉冲信号。因此，在传统的低小慢目标探测雷达系统回波信号处理系统中，雷达系统都是发射脉冲宽度较大，达到微秒级的线性调频信号，再回波信号处理时，通过脉冲压缩技术，将大时宽的线性调频信号压缩成脉冲宽度可达纳秒级的窄时宽信号，以使得低小慢目标探测雷达系统系统可以同时获得远探测距离和高测距精度。

脉冲压缩通常由快速傅里叶变换，脉压系数相乘，快速逆傅里叶变换这三个步骤完成，快速傅里叶变换步骤将时域的回波信号变换至频域；脉压系数相乘实现频域回波信号与本地存储的脉冲压缩系数的向量点乘；快速逆傅里叶变换将完成向量点乘的数据从频域变换到时域，完成脉冲压缩，以在时域获得窄时宽的脉冲信号。脉冲压缩的结果受到本地存储的脉冲压缩系数的影响，本地存储的脉冲压缩系数的对雷达系统的测距精度具有至关重要的作用。

现有技术中，脉冲压缩系数一般是根据系统参数设置的脉宽和带宽，生成一组标准的脉冲压缩系数。由于通常由直接频率合成器(DDS)产生由雷达系统发射出的线性调频信号，但DDS产生的信号与标准信号有差距，并且，DDS发出的线性调频信号经过传输、采集和数字下变频等一系列信号处理操作后，与DDS最初产生的信号存在较大的差距，传统的脉冲压缩系数未综合考虑上述因素造成的影响，造成脉冲压缩系数不准确，不能符合实际工作情况。

发明内容

为解决上述技术问题，本发明提出了一种目标探测雷达系统脉冲压缩系数生成方法及装置，所述方法及装置，用以解决现有技术中脉冲压缩系数未综合考虑影响因素，造成脉冲压缩系数不准确，不能符合实际工作情况的技术问题。

根据本发明的第一方面，提供一种目标探测雷达系统脉冲压缩系数生成装置，包括FPGA主控模块、直接频率合成器(DDS)、射频前端模块、模数转换器以及数据分析模块；

所述FPGA主控模块，用于完成直接频率合成器(DDS)的参数配置以及接收由所述模数转换器采集的环回DDS数据，并将其导出；

所述直接频率合成器(DDS)，用于接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块；

所述射频前端模块，用于完成信号上变频、无线发射、无线接收、信号下变频功能；所述射频前端模块的发射通道将DDS输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器；

所述模数转换器，用于完成中频信号的带通采样功能，实现环回信号的数字采集，并将采集到的环回信号的数据传输至所述FPGA主控模块中；

所述数据分析模块，用于接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数。

进一步地，所述直接频率合成器(DDS)为直接频率合成器芯片AD9910。

进一步地，所述FPGA主控模块需配置重频prt、脉宽Pw、起始频率f_start、截止频率f_end、调频率ramp共五个参数。

进一步地，重频prt和脉宽Pw为微秒级。

进一步地，所述直接频率合成器(DDS)的起始频率控制字S，截止频率控制字E，调频率控制字R这三个寄存器控制字受DDS的参考时钟频率Fs影响，其参数计算公式为：

所述FPGA主控模块根据公式(1)-(3)计算出配置参数并输入所述直接频率合成器(DDS)。

进一步地，所述FPGA主控模块计算出配置参数并输入所述直接频率合成器(DDS)之后，以示波器和频谱仪观测所述直接频率合成器(DDS)输出的线性调频信号是否达标，再根据实际情况微调各参数，以使所述直接频率合成器(DDS)输出适当的线性调频波形。

进一步地，模数转换器的时钟和所述直接频率合成器(DDS)的参考时钟时钟同源。

进一步地，所述数据分析模块使用MATLAB分析软件进行数据分析及计算。

根据本发明第二方面，具有如前所述的脉冲压缩系数生成装置，其特征在于，所述方法执行以下步骤：

步骤S201：所述FPGA主控模块计算出配置参数并输入所述直接频率合成器(DDS)；

步骤S202：所述直接频率合成器(DDS)接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块；

步骤S203：所述射频前端模块的发射通道将所述直接频率合成器(DDS)输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器；

步骤S204：所述模数转换器将采集到的环回信号的数据传输至所述FPGA主控模块；

步骤S205:所述数据分析模块接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数。

进一步地，所述先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数，包括：

步骤S2051：运行雷达系统的数字下变频算法，将环回DDS信号转换成数字基带复信号；

步骤S2052：获取采样率Fs，计算ts＝1/Fs，设置矩阵行值K＝400、列值L＝400；设置k＝1；

步骤S2053：判断k是否小于所述矩阵行值K，若是，进入步骤S2054；若否，进入步骤S2058；

步骤S2054：设置l＝1；

步骤S2055：判断l是否小于所述矩阵列值L，若是，进入步骤S2056；若否，进入步骤S2057；

步骤S2056：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；

ref_f＝conj(fft(ref_t))；pc＝ifft(fft(s)*ref_f)；max_pc(K,L)＝max(abs(pc))；l赋值为l+1；进入步骤S2055；

其中，以标准DDS参数为基础，在雷达系统参数上下浮动20％的范围内，以参数的1‰为步长进行搜索，exp为指数运算符，conj为共轭运算符，fft为快速傅里叶变换，ifft为逆快速傅里叶变换，abs为取绝对值运算符，max为求向量最大值函数；Pw为脉宽，是脉冲波形的时域持续时间；所述t＝-Pw/2:ts:Pw/2-ts表示生成从-Pw/2到Pw/2-ts、以ts为间隔的数值；Bw为带宽，取值可以为截止频率-起始频率；t为时间变量，是脉冲波型持续时间内的时间变量；ref_f为回波信号s进行脉冲压缩处理时需要用的脉冲压缩系数；pc为脉冲压缩后的时域波形值；

步骤S2057：k赋值为k+1，进入步骤S2053；

步骤S2058：对于得到的K*L的矩阵，找出矩阵所有元素中的最大值，记录最大值对应的二维坐标，将K、L的值分别赋值为最大值对应的坐标值；

步骤S2059：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；

ref_f＝conj(fft(ref_t))，计算得到脉冲压缩性能最好的脉冲压缩系数ref_f。

根据本发明的上述方案，能够替代现有技术中通常采用的脉冲压缩系数的生成方法，考虑了雷达系统的实际工作情况，环回采集经过实际雷达系统处理的线性调频信号，依据该线性调频信号迭代处理计算出最优的脉冲压缩系数，可以提高雷达系统的脉冲压缩性能，获得更高的雷达距离分辨力。

上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。

附图说明

构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明提供如下附图进行说明。在附图中：

图1为一个实施方式的目标探测雷达系统脉冲压缩系数生成装置的总体架构图；

图2为本发明一个实施方式的目标探测雷达系统脉冲压缩系数生成方法的流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

首先结合图1说明本发明的实现目标探测雷达系统脉冲压缩系数生成装置的总体架构，图1示出了根据本发明的一个实施方式的目标探测雷达系统脉冲压缩系数生成装置的总体架构图。如图1所示：

所述目标探测雷达系统脉冲压缩系数生成装置模拟实际的雷达系统而构建，包括FPGA主控模块、直接频率合成器(DDS)、射频前端模块、模数转换器以及数据分析模块。

所述FPGA主控模块，用于完成直接频率合成器(DDS)的参数配置以及接收由所述模数转换器采集的环回DDS数据，并将其导出。

所述直接频率合成器(DDS)，用于接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块。

所述射频前端模块，用于完成信号上变频、无线发射、无线接收、信号下变频功能。所述射频前端模块的发射通道将DDS输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器。

所述模数转换器，用于完成中频信号的带通采样功能，实现环回信号的数字采集，并将采集到的环回信号的数据传输至所述FPGA主控模块中。

所述数据分析模块，用于接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数。

本实施例中，所述直接频率合成器(DDS)为直接频率合成器芯片AD9910，所述FPGA主控模块需配置重频prt、脉宽Pw、起始频率f_start、截止频率f_end、调频率ramp共五个参数。其中重频prt和脉宽Pw根据实际的系统参数而定，一般为微秒级，例如重频75微秒，脉宽7.5微秒。

所述直接频率合成器(DDS)的起始频率控制字S，截止频率控制字E，调频率控制字R这三个寄存器控制字受DDS的参考时钟频率Fs影响，其参数计算公式为：

所述FPGA主控模块根据公式(1)-(3)计算出配置参数并输入所述直接频率合成器(DDS)。之后，以示波器和频谱仪观测所述直接频率合成器(DDS)输出的线性调频信号是否达标，再根据实际情况微调各参数，以使所述直接频率合成器(DDS)输出适当的线性调频波形。

为了确保雷达系统的测距可靠性，务必保证模数转换器的时钟和所述直接频率合成器(DDS)的参考时钟时钟同源。本实施例中，所述数据分析模块使用MATLAB分析软件进行数据分析及计算。

下面结合图2说明本发明的实现目标探测雷达系统脉冲压缩系数生成方法，图2示出了根据本发明的一个实施方式的目标探测雷达系统脉冲压缩系数生成方法的流程图。如图2所示：

步骤S201：所述FPGA主控模块计算出配置参数并输入所述直接频率合成器(DDS)；

步骤S202：所述直接频率合成器(DDS)接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块；

步骤S203：所述射频前端模块的发射通道将所述直接频率合成器(DDS)输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器；

步骤S204：所述模数转换器将采集到的环回信号的数据传输至所述FPGA主控模块；

步骤S205:所述数据分析模块接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数。

所述步骤S205:所述数据分析模块接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数，包括：

步骤S2051：运行雷达系统的数字下变频算法，将环回DDS信号转换成数字基带复信号；

步骤S2052：获取采样率Fs，计算ts＝1/Fs，设置矩阵行值K＝400、列值L＝400；设置k＝1；

步骤S2053：判断k是否小于所述矩阵行值K，若是，进入步骤S2054；若否，进入步骤S2058；

步骤S2054：设置l＝1；

步骤S2055：判断l是否小于所述矩阵列值L，若是，进入步骤S2056；若否，进入步骤S2057；

步骤S2056：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；ref_f＝conj(fft(ref_t))；pc＝ifft(fft(s)*ref_f)；max_pc(K,L)＝max(abs(pc))；l赋值为l+1；进入步骤S2055；

其中，以标准DDS参数为基础，在雷达系统参数上下浮动20％的范围内，以参数的1‰为步长进行搜索，exp为指数运算符，conj为共轭运算符，fft为快速傅里叶变换，ifft为逆快速傅里叶变换，abs为取绝对值运算符，max为求向量最大值函数；Pw为脉宽，是脉冲波形的时域持续时间；所述t＝-Pw/2:ts:Pw/2-ts表示生成从-Pw/2到Pw/2-ts、以ts为间隔的数值；Bw为带宽，取值可以为截止频率-起始频率；t为时间变量，是脉冲波型持续时间内的时间变量；ref_f为回波信号s进行脉冲压缩处理时需要用的脉冲压缩系数；pc为脉冲压缩后的时域波形值。

步骤S2057：k赋值为k+1，进入步骤S2053；

步骤S2058：对于得到的K*L的矩阵，找出矩阵所有元素中的最大值，记录最大值对应的二维坐标，将K、L的值分别赋值为最大值对应的坐标值；

步骤S2059：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；ref_f＝conj(fft(ref_t))，计算得到脉冲压缩性能最好的脉冲压缩系数ref_f，其中“j*π*Ka*t∧2”中的j表示这是一个虚数。

以下实施例结合说明然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数的方式。

以采样率Fs为40MHz为例，以标准DDS参数为基础，在雷达系统参数上下浮动20％的范围内，以参数的1‰为步长进行搜索，按照步骤S2056中的公式计算得到400*400的矩阵，从所述矩阵中找出最大值对应的二维坐标(k₀,l₀)，将其代入公式Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；

Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t∧2)，其中“j*π*Ka*t∧2”中的j表示这是一个虚数；

ref_f＝conj(fft(ref_t))，计算得到脉冲压缩性能最好的脉冲压缩系数。

需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。

在本发明所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。

上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机装置(可以是个人计算机，实体机服务器，或者网络云服务器等，需安装Windows或者Windows Server操作系统)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (2)

1.一种低小慢目标探测雷达系统脉冲压缩系数生成装置，其特征在于，所述目标探测雷达系统脉冲压缩系数生成装置模拟实际的雷达系统而构建，包括FPGA主控模块、直接频率合成器(DDS)、射频前端模块、模数转换器以及数据分析模块；

所述FPGA主控模块，用于完成直接频率合成器(DDS)的参数配置以及接收由所述模数转换器采集的环回DDS数据，并将其导出；

所述直接频率合成器(DDS)，用于接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块；

所述射频前端模块，用于完成信号上变频、无线发射、无线接收、信号下变频功能；所述射频前端模块的发射通道将DDS输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器；

所述模数转换器，用于完成中频信号的带通采样功能，实现环回信号的数字采集，并将采集到的环回信号的数据传输至所述FPGA主控模块中；

所述数据分析模块，用于接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数；

所述FPGA主控模块需配置重频prt、脉宽Pw、起始频率f_start、截止频率f_end、调频率ramp共五个参数；

所述直接频率合成器(DDS)的起始频率控制字S，截止频率控制字E，调频率控制字R这三个寄存器控制字受DDS的参考时钟频率Fs影响，其参数计算公式为：

所述FPGA主控模块根据公式(1)-(3)计算出配置参数并输入所述直接频率合成器(DDS)；

所述FPGA主控模块计算出配置参数并输入所述直接频率合成器(DDS)之后，以示波器和频谱仪观测所述直接频率合成器(DDS)输出的线性调频信号是否达标，再根据实际情况微调各参数，以使所述直接频率合成器(DDS)输出适当的线性调频波形；

所述直接频率合成器(DDS)为直接频率合成器芯片AD9910；

重频prt和脉宽Pw为微秒级；

模数转换器的时钟和所述直接频率合成器(DDS)的参考时钟同源；

所述数据分析模块使用MATLAB分析软件进行数据分析及计算。

2.一种低小慢目标探测雷达系统脉冲压缩系数生成方法，具有如权利要求1所述脉冲压缩系数生成装置，其特征在于，所述方法执行以下步骤：

步骤S201：所述FPGA主控模块计算出配置参数并输入所述直接频率合成器(DDS)；

步骤S202：所述直接频率合成器(DDS)接受所述FPGA主控模块的控制，产生雷达系统要求的线性调频信号，并将其发送至所述射频前端模块；

步骤S203：所述射频前端模块的发射通道将所述直接频率合成器(DDS)输入的线性调频信号上变频至射频频段，并通过射频线缆环回至所述射频前端模块的接收通道，接收通道将射频频段的信号下变频至中频后，发送至所述模数转换器；

步骤S204：所述模数转换器将采集到的环回信号的数据传输至所述FPGA主控模块；

步骤S205:所述数据分析模块接收由所述FPGA主控模块导出的所述环回DDS数据，先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数；

所述先运行雷达系统的数字下变频算法，然后运行脉冲压缩系数优化方法，计算出最优的脉冲压缩系数，包括：

步骤S2051：运行雷达系统的数字下变频算法，将环回DDS信号转换成数字基带复信号；

步骤S2052：获取采样率Fs，计算ts＝1/Fs，设置矩阵行值K＝400、列值L＝400；设置k＝1；

步骤S2053：判断k是否小于所述矩阵行值K，若是，进入步骤S2054；若否，进入步骤S2058；

步骤S2054：设置l＝1；

步骤S2055：判断l是否小于所述矩阵列值L，若是，进入步骤S2056；若否，进入步骤S2057；

步骤S2056：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；ref_f＝conj(fft(ref_t))；pc＝ifft(fft(s)*ref_f)；max_pc(K,L)＝max(abs(pc))；l赋值为l+1；进入步骤S2055；

其中，以标准DDS参数为基础，在雷达系统参数上下浮动20％的范围内，以参数的1‰为步长进行搜索，exp为指数运算符，conj为共轭运算符，fft为快速傅里叶变换，ifft为逆快速傅里叶变换，abs为取绝对值运算符，max为求向量最大值函数；Pw为脉宽，是脉冲波形的时域持续时间；所述t＝-Pw/2:ts:Pw/2-ts表示生成从-Pw/2到Pw/2-ts、以ts为间隔的数值；Bw为带宽，取值可以为截止频率-起始频率；t为时间变量，是脉冲波型持续时间内的时间变量；ref_f为回波信号s进行脉冲压缩处理时需要用的脉冲压缩系数；pc为脉冲压缩后的时域波形值；

步骤S2057：k赋值为k+1，进入步骤S2053；

步骤S2058：对于得到的K*L的矩阵，找出矩阵所有元素中的最大值，记录最大值对应的二维坐标，将K、L的值分别赋值为最大值对应的坐标值；

步骤S2059：计算Pw1＝0.8*Pw+K*(Pw/1000)；Bw1＝0.8*Bw+L*(Bw/1000)；Ka＝Bw1/Pw1；t＝-Pw/2:ts:Pw/2-ts；ref_t＝exp(j*π*Ka*t^2)；

ref_f＝conj(fft(ref_t))，计算得到脉冲压缩性能最好的脉冲压缩系数ref_f。

Images