一种基于MIMO图像域的BP宽带合成方法
专利汇可以提供一种基于MIMO图像域的BP宽带合成方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于MIMO图像域BP宽带合成方法，其特点是首先利用MIMO体制下多个发射机同时发射中心 频率 步进的LFM 信号 ，对接收到的各个窄带回波信号经BP成像得到各个低 分辨率 的子图像后，再对多个低分辨的子图像进行图像域的相参 叠加 以获得高分辨距离向的SAR图像。这种方法与现有的数据域的宽带合成方法相比，大大简化了宽带合成过程中的通道不一致性误差补偿步骤，在图像域进行宽带合成避免了数据域合成带来的子带重叠误差校正繁琐的问题，在降低运算量的同时实现了SAR距离向的高分辨率成像。本发明可适用于系统 硬件 条件(发射端瞬时发射带宽和接收机 采样 率)受限条件下的SAR距离向高分辨率成像的技术领域。，下面是一种基于MIMO图像域的BP宽带合成方法专利的具体信息内容。
1.一种基于MIMO图像域的BP宽带合成的方法，其特征是它包括以下步骤：
步骤1、用于基于MIMO图像域的BP宽带合成方法相关参数的初始化
初始化的参数均为已知，且初始化的参数如下：所有的坐标都是以笛卡尔直角坐标形式给出；光速为C；雷达发射线性调频信号，脉冲重复频率为PRF，脉冲重复周期为PRI，MIMO阵列天线个数为N，也即每个脉冲重复周期内雷达发射N个中心频率步进的窄带即子带线性调频信号，这N个天线同时也应用于接收各自信号反射的回波；雷达发射的各个子带信号的带宽为Bpn；雷达发射的各个子带信号的时宽为Tpn；发射信号的调频斜率记为γ， 发射脉冲各子带信号的步进频率为Δf；经宽带合成后的宽带信号的中心频率为f0；各子带信号的中心频率记为fc(n)，fc(n)＝f0+nΔf,n＝0,1，…，N-1；发射平台中心Pt0；场景中心Pcenter；雷达参考斜距R0；雷达在距离向上窄带信号的采样频率Fpn；
雷达载机平台速度为V；合成孔径长度为Lsar；合成孔径时间为Tsar， 一个合成孔径时间内距离向采样点数为Nr，便于FFT计算，在此取值为2的整数次幂，方位向采样点数为Na，便于FFT计算，采样点数为取值为2的整数次幂；快时间为一个长度为Nr的行矢量，慢时间ts为一个长度为Na的行矢量，
BP成像场景距离向即斜距平面和方位向点数分别为M和K；BP成像场景距离向即斜距平面和方位向网格间距为dx和dy；目标位置(x0，y0)；
步骤2、获取雷达回波数据
(n)
第n个子带信号回波数据为一个Na×Nr的矩阵，记为S ，其表现形式如下：
其中第n个子带信号的第jj个快时间、第ii个慢时间的雷达回波数据记为
n ＝ 1， …，N；
ii＝1,2，…，Na；jj＝1,2，…，Nr，N、Na、Nr分别为子带数、方位采样点数、距离向采样点数，其中τ(ii)为第ii个方位时刻雷达天线到场景目标点的双程距离延迟时间，ts(ii)为ts的第ii个元素，exp(·)为以自然底数e
为底的指数函数，j为虚数单位，rect(t/Tpn)为矩形窗函数，满足
t、Tpn分别为距离向快时间、子带时宽。
步骤3、对雷达回波数据矩阵每一行即距离向进行脉冲压缩
取出步骤2中所获得的雷达回波数据 构造距离向参考函数Hr，利用传统的脉冲(n) (n) (n)
压缩方法对S 的每一行即距离向进行脉冲压缩，脉冲压缩后的矩阵记为PS ，PS ＝(n) (n) (n) (n)
IFFT(FFT((S )′)·FFT(Hr))，其中(S )′表示S 的转置，PS 的第ii行第jj列记为其中n＝1，…，N；ii＝1,2，…，Na；jj＝1,2，…，Nr。其中Hr为一个Na×Nr的矩阵，其表现形式如下：
2
Hr＝ones(Na，1)×rect(t/Tpn)·exp(jπγt)
其中ones(Na，1)产生Na×1的全1矩阵，FFT(·)、IFFT(·)分别为快速傅里叶变换和快速傅里叶逆变换；
步骤4、对每个距离向的所有回波数据矩阵进行升采样L倍
对步骤3中脉冲压缩后的N个子带、Na个方位向数据矩阵统一做如下处理：
取出步骤3中处理后的第n个子带信号的第ii个方位向的脉压后的数据 n
＝1，…，N；ii＝1,2，…，Na， 表示PS(n)矩阵的第ii行即第ii个方位向数据；
首先将第n个子带信号的第ii个方位向的脉压后的数据 变换到频域得到数据记为 然后将数据 中间插零，得到数据矩阵记为
其中n＝1，…，N；ii＝
1,2，…，Na，其中S0为1行(L-1)·Nr的零矩阵；然后将数据矩阵记为 变换到时域n＝1，…，N；ii＝1,2，…，Na；
步骤5、对每个方位向数据矩阵进行相关叠加进行BP成像
对步骤4中升采样后的N个子带、Na个方位向数据矩阵统一做如下处理：
步骤5.1取出步骤4中第n个子带、第ii个升采样后的方位向数据矩阵 其
中n＝1，…，N；ii＝1，2，…，Na，对第n个子带、第ii个升采样后的方位向数据矩阵做距离向FFT变换到距离-多普勒域，即 其中 是对
的转置操作；然后找到第ii个方位时刻雷达平台的位置Pt＝V·ts与成像场景中的点T(mm，kk)＝[(mm-M/2)·dx，(kk-K/2)·dy]，mm＝1，2，…，M，kk＝1，2，…，K距离历史为记为R(ii)， T(1)表示场景
中第(mm,kk)个点的横坐标，T(2)表示场景中第(mm,kk)个点的纵坐标，其中Rr为斜距平面内的场景参考斜距，并计算第ii个方位向用于BP成像搜索的距离向编号，记为ID，其中round(·)为取四舍五入的函数；最后对
(n) (n)
第ii个方位向、ID对应的距离向数据补偿方位向多普勒相位得到的图像矩阵记为I ，I的第mm行第kk列记为 其中λ(n)＝C/f(n)
为第n个子带信号的波长，
步骤5.2按照步骤5.1对场景每个点(mm,kk)，mm＝1,2，…，M；kk＝1,2，…，K，M和K分别为BP成像场景距离向即斜距平面和方位向点数，补偿Na个方位向的数据ii＝1,2，…，Na；n＝1，…，N的多普勒相位，并做相参叠加就可得到低分辨距离向的子(n) (n)
带图I ，n＝1，…，N，I 为一个M×K的矩阵；
步骤6、对每个距离向低分辨的子图像即图像域作相位补偿
对步骤5中处理后的N个距离向低分辨的子图像矩阵统一做如下处理：
(n)
取出步骤5中第n个距离向低分辨的子图像矩阵I ，n＝1，…，N，对第n个距离向低分辨的子图像矩阵I(n)做相位补偿得到补偿后的子图像记为 其中Φ
为相位补偿函数，其表达式为Φ＝exp(-j4πfc(n)·r/C)，其中r＝C·t/2为斜距平面对应的距离轴；
步骤7、对每个子图像即图像域作频移操作
构造第n个子图像对应的频移参考函数Sref(n)＝exp(j4π(n-(N+1)/2)·Δf·r/C)，其中n＝1，…，N，将该参考函数与子图像矩阵相乘进行频移操作得到频移后的数据记为其中n＝1，…，N；
步骤8、对频移后的子图像作FFT操作变换到波数域
对步骤7中处理后的N个子图像数据矩阵统一做如下处理：
取出步骤7中第n个子数据矩阵 n＝1，…，N，利用快速傅里叶变换方法对 的每一行进行快速傅里叶变换，得到处理后的M行K列的矩阵记为 n
＝1，…，N，其中 是对 的转置矩阵；
步骤9、对波数域各子带图像进行相参叠加
取出步骤8中得到的N个子图像矩阵进行相干叠加得到波数域的等效宽带信号记为Ifull， n＝1，…，N。
步骤10、将合成的波数域的宽带图像变换到空域
对步骤9中合成的波数域的M个距离向行向量统一做如下处理：
取出步骤9中第mm个距离向的行向量Ifull(mm，：)，其中mm＝1,2，…，M，Ifull(mm，：)表示Ifull的第mm行的所有数据，将其进行快速傅里叶逆变换到空域即可得到高分辨距离向的图像记为Image，Image(mm，：)＝IFFT(Ifull(mm，：))，mm＝1，2，…，M，则Image即为最终高分辨距离向图像矩阵。
一种基于MIMO图像域的BP宽带合成方法
技术领域：
雷达载机平台速度为V；合成孔径长度为Lsar；合成孔径时间为Tsar， 一个合成孔径时间内距离向采样点数为Nr(便于FFT计算，在此取值为2的整数次幂)，方位向采样点数为Na(便于FFT计算，在此取值为2的整数次幂)；快时间为一个长度为Nr的行矢量，慢时间ts为一个长度为Na的行矢量，
BP成像场景距离向(斜距平面)和方位向点数分别为M和K；BP成像场景距离向(斜距平面)和方位向网格间距为dx和dy；目标位置(x0，y0)。
为底的指数函数，j为虚数单位(下同)，rect(t/Tpn)为(t、Tpn分别为距离向快时间、子带时宽)矩形窗函数，满足
IFFT(FFT((S )′)·FFT(Hr))，其中(S )′表示S 的转置，PS 的第ii行第jj列记为其中n＝1，…，N；ii＝1，2，…，Na；jj＝1，2，…，Nr。其中Hr为一个Na×Nr的矩阵，其表现形式如下：
2，…，Na，其中S0为1行(L-1)·Nr的零矩阵；
的转置操作；然后找到第ii个方位时刻雷达平台的位置Pt=V·ts与成像场景中的点T(mm,kk)(mm=1,2，…，M；kk=1,2，…，K，其中T(mm,kk)=[(mm-M/2)·dx，(kk-K/2)·dy])的距离历史为记为R(ii)， (T(1)表示场景中
第(mm,kk）个点的横坐标，T(2)表示场景中第(mm，kk)个点的纵坐标)，其中Rr为斜距平面内的场景参考斜距，并计算第ii个方位向用于BP成像搜索的距离向编号，记为ID，其中round(·)为取四舍五入的函数；最后对
(n) (n)
第ii个方位向、ID对应的距离向数据补偿方位向多普勒相位得到的图像矩阵记为I ，I的第mm行第kk列记为 (其中λ(n)=C/f(n)为
第n个子带信号的波长)，
为相位补偿函数，其表达式为Φ=exp(-j4πfc(n)·r/C)，其中r=C·t/2为斜距平面对应的距离轴。
其中 是对 的转置矩阵。
具体实施方式
直角坐标形式给出；光速为C=3×10m/s；雷达发射线性调频信号，脉冲重复频率记为PRF,PRF=500Hz，脉冲重复周期记为PRI， MIMO天线个数为
N=3，也即每个脉冲重复周期内雷达发射3个中心频率步进的窄带(子带)线性调频信号，这3个天线同时也应用于接收各自信号反射的回波；雷达发射的各个子带信号的带宽为Bpn=300MHz；雷达发射的各个子带信号的时宽为Tpn=1μs；发射信号的调频斜率记为γ， 发射脉冲各子带信号的步进频率为
Δf=300MHz；经宽带合成后的宽带信号的中心频率为f0=10GHz；各子带信号的中心频率
9
记为fc(n)，fc(n)=f0+nΔf=(10+n×0.3)×10，n=0，1,2；发射平台中心Pt0=[0,0,4000]；
场景中心Pcener=[3000,0,0]；雷达参考斜距R0=5000；雷达在距离向上窄带信号的的采样频率Fpn=390MHz；雷达载机平台速度为V=[0，100,0]；合成孔径长度为Lsar=400m；合成孔径时间记为 一个合成孔径时间内距离向采样点数为
Nr=512，方位向采样点数为Na=512；快时间为 慢时间为
BP成像场景距离向(斜距平面)和方位向点数分别为M=30和
K=30；BP成像场景距离向(斜距平面)和方位向网格间距为dx=0.5m和dy=0.5m；目标位置(5000，0)。
exp(·)为以自然底数e为底的指数函数，j
为虚数单位(下同)， 为平方根操作(下同)，rect(t/Tpn)(其中t、Tpn分别为距离向快时间、子带时宽)为矩形窗函数，矩形窗函数，满足
IFFT(FFT((S )′)·FFT(Hr))，其中(S )′表示S 的转置，PS 的第ii行第jj列的数据记为 其中n＝1，2，3；ii＝1，2，…，512；jj＝1，2，…，512。其中Hr为一个
512×512的矩阵，其表现形式如下：
的 数 据 n ＝1，2，3；ii＝ 1，2， …，512， 表 示 PS 矩 阵 的 第ii 行(第ii个方位向)数据(下同)。首先将这行数据变换到频域得到的数据记为
然后将该行数据中间插零后得到的数据记为
其中n＝1，2，3；ii
＝1，2，…，512，其中S0为1行3584列的零矩阵；最后将其变换到时域得到的矩阵记为[0089] 步骤5、对每个方位向数据矩阵进行相关叠加进行BP成像
后找到第ii个方位时刻雷达平台的位置 与成
像场景 中的点T(mm，kk)(mm＝ 1，2，…，30；kk＝ 1，2，…，30，其中 T(mm，kk)＝[(mm-15)×0.5，(kk-15)×0.5]) 的距离历史记为R(ii)，并计算第ii个方
位向用于BP成像搜索距离向编号记为ID，
位向多普勒相位后得到的矩阵记为I ，I 的第mm行第kk列记为 则
距离向的子带图I(n)，n＝1，2，3，I(n)为一个M×K的矩阵。
阵相乘进行频移操作得到的数据记为
其中 是对 的转置矩阵。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 等离子体处理装置及其基片直流偏置电压测量方法 | 2021-03-29 | 1 |
| 一种基于决策树模型的车辆强制换道决策方法 | 2022-06-08 | 1 |
| 频率估计的一种多段信号融合方法 | 2022-11-10 | 1 |
| 一种基于克拉克变换的同步相量测量方法 | 2020-08-28 | 0 |
| 一种遥感卫星照射源前视合成孔径雷达成像方法 | 2023-12-30 | 1 |
| 利用频域信息滤除衰减直流的电力信号同步相量测量方法 | 2022-02-11 | 1 |
| 用于补偿用作采样频率的本地时钟的漂移的装置和方法 | 2023-05-29 | 2 |
| 一种基于靶标轨迹跟踪的主轴径向回转误差测量方法 | 2020-06-18 | 2 |
| 基于步态非线性特征的心电RR间隔预测关联方法 | 2020-07-01 | 3 |
| Arrangement for and method of acquiring a monochrome image with a color image capture processor | 2022-12-10 | 2 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
