CN109471087B - 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 - Google Patents
基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 Download PDFInfo
- Publication number
- CN109471087B CN109471087B CN201811218661.6A CN201811218661A CN109471087B CN 109471087 B CN109471087 B CN 109471087B CN 201811218661 A CN201811218661 A CN 201811218661A CN 109471087 B CN109471087 B CN 109471087B
- Authority
- CN
- China
- Prior art keywords
- array
- vector
- signal
- fourier transform
- elements
- 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.)
- Active
Links
Images
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/42—Diversity systems specially adapted for radar
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,主要解决现有方法计算复杂度较高的问题。其实现步骤是:构建互质MIMO雷达结构;通过雷达接收子阵列接收反射信号并对雷达输出信号建模;构造互质MIMO雷达接收信号对应的二阶差集和集虚拟阵列信号;对二阶差集和集虚拟阵列信号进行二次零填充;对二次零填充后的差集和集虚拟阵列信号进行快速傅里叶变换操作,构建空间谱;根据所得空间谱进行波达方向估计。本发明在降低了波达方向估计计算复杂度的同时,更充分地利用了阵列提供的信息以保证分辨率。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及对雷达信号、声学信号及电磁信号的统计信号处理,具体是一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,可用于有源定位和目标探测。
背景技术
作为阵列信号处理领域的基本问题之一,波达方向(Direction-of-Arrival,DOA)估计利用传感器阵列接收信号,并通过对接收信号进行统计处理,从中提取信号中所包含的波达方向信息,在雷达、声呐、语音、无线通信等领域有着广泛应用。其中,利用多入多出(Multiple-input Multiple-output,MIMO)雷达进行波达方向估计是一个重要分支。
MIMO雷达利用多个传感器阵列,分别发射正交信号并接收反射信号。在传统的MIMO雷达结构中,接收阵列通常被设置为均匀线性阵列,以此满足奈奎斯特采样定理。然而,均匀线性阵列的设置方式限制了阵列孔径,进而限制了空间分辨率等性能。在此背景下,作为互质阵列结构在MIMO雷达上的应用,互质MIMO雷达充分利用了互质阵列的系统化结构,通过构造具有更大孔径的虚拟阵列,实现了DOA估计性能的提升。
现有基于互质MIMO雷达的DOA估计方法大多针对所构建的虚拟阵列对应的二阶等价虚拟信号进行一系列复杂运算,主要包括求逆、特征值分解等复杂的矩阵运算,以及凸优化问题的设计与求解等高复杂度过程,这些运算过程计算复杂度高,在实际系统中的硬件实现较为复杂,且在实时性要求较高的应用场景下面临一定的挑战。此外,由于互质MIMO雷达构造的虚拟阵列中的阵元存在不均匀性,现有的基于互质MIMO雷达的DOA估计方法主要分为两类:一类仅利用虚拟阵列中均匀连续部分的虚拟阵元而舍弃非连续部分,由此导致了一定的性能损失;另一类则通过较为复杂的内插与重建技术对缺失部分进行恢复,由此导致了算法复杂度的上升。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,通过对互质MIMO雷达的二阶差集和集信号进行简单的二次零填充,并进行快速傅里叶变换,利用了虚拟阵列的全部阵元,实现了波达方向估计分辨率性能的提升。此外,本发明所提方法操作简单,计算复杂度低,易于在实际系统的硬件实现。
本发明的目的是通过以下技术方案来实现的:一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,包含以下步骤:
(1)利用M+N个实际阵元构建互质MIMO雷达结构,构建方法如下:选取一组互质的整数M、N,分别用于构造发射子阵列和接收子阵列,其中,发射子阵列包含M个间距为Nd的阵元,其位置为0,Nd,...,(M-1)Nd;接收子阵列包含N个间距为Md的阵元,其位置为0,Md,...,(N-1)Md;d为发射子阵列发射信号的半波长即发射子阵列和接收子阵列部署为单基地雷达的形式;
(2)利用步骤(1)中构建的发射子阵列发射M个正交信号;假设空间远场中有Q个目标,其相对于雷达的方向表示为θ=[θ1,θ2,…,θQ]T,[·]T表示转置操作;该M个信号被空间中Q个目标反射后,其反射信号被接收子阵列的N个阵元接收,经过匹配滤波后,在t时刻雷达输出信号x(t)可建模为:
其中,sq(t)为t时刻第q个目标的反射率,n(t)为与信号无关加性噪声,服从零均值复高斯分布,j为虚数单位;
因此,步骤(2)中的接收信号可以表示为
其中,的理论协方差矩阵, 为(M+N-1)2×Q维虚拟阵列导引矩阵,为包含Q个目标反射功率的向量,为噪声功率,i=vec(IM+N-1),IM+N-1表示(M+N-1)×(M+N-1)维单位矩阵,vec(·)表示向量化操作,即把矩阵中的各列依次堆叠以形成一个新的向量,(·)*表示共轭操作;从向量y中选择对应于中各虚拟阵元位置的等价虚拟信号,得到对应二阶差集和集虚拟阵列信号可表示为
其中,表示与位置上的阵元相对应的元素,max(·)表示一个集合中最大的元素,[·]l表示向量中的第l个元素,l=0,1,...,L-1, 在此基础上,在所得的向量z之后补充一定数量的零,使所得向量中包含P个元素,表示为
其中P为2的整数次幂且P>L,0P-L表示P-L维的零向量;
构建一个空间谱,该空间谱的横轴表示角度θ,其与空间响应第p个元素的关系可表示为:
(6)根据所得空间谱进行波达方向估计:对步骤(5)中空间谱进行谱峰搜索操作,幅度最大的前Q个峰值对应的角度,即为空间中Q个目标的波达方向估计结果。
本发明与现有技术相比具有以下优点:
(1)本发明所提方法充分利用了互质MIMO雷达结构,在阵元数量一定的情况下,获得比传统MIMO雷达更大的阵列孔径,且通过易操作的二次零填充方法,利用了所有虚拟阵元提供的信息,避免了信息损失,从而实现了更高分辨率的波达方向估计。
(2)本发明所提方法在互质MIMO雷达二阶差集和集虚拟阵列信号基础上进行快速傅里叶变换操作,通过所得响应构建空间谱,并通过谱峰搜索方式获得波达方向估计结果,避免了现有方法中常用的矩阵求逆、特征值分解、凸优化问题设计求解等复杂计算过程,具有计算复杂度低、易于在实际系统中实现的特点,且在实时性要求较高的应用场景下具有显著优势。
附图说明
图1是本发明的方法总体流程框图;
图2是本发明中互质MIMO雷达的结构示意图;
图3是用于体现本发明所提方法分辨率性能的空间谱示意图;
图4是用于体现本发明所提方法与传统方法相比在分辨率性能上的优势。
具体实施方式
以下参照附图,对本发明的技术方案和效果作进一步的详细说明。
作为互质阵列在MIMO雷达上的应用,互质MIMO雷达结构充分吸收并体现了互质阵列系统化的结构优势,近年来受到学术界的广泛关注。目前大多基于此的DOA估计方法计算复杂度高,在实时性要求较高的场景下面临挑战,且在实际工程实现时,面临一定的困难。针对以上问题,本发明提供了一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,参照图1,本发明的实现步骤如下:
步骤一:利用M+N个实际阵元构建互质MIMO雷达结构,构建方法如下:选取一组互质的整数M、N,分别用于构造发射子阵列和接收子阵列。其中,发射子阵列包含M个间距为Nd的阵元,其位置为0,Nd,...,(M-1)Nd;接收子阵列包含N个间距为Md的阵元,其位置为0,Md,...,(N-1)Md;d为发射子阵列发射信号的半波长即发射子阵列和接收子阵列部署为单基地雷达的形式。
步骤二:通过雷达接收子阵列接收反射信号并对雷达输出信号建模。利用步骤一中构建的发射子阵列发射M个正交信号。假设空间远场中有Q个目标,其相对于雷达的方向表示为θ=[θ1,θ2,…,θQ]T,[·]T表示转置操作。该M个信号被空间中Q个目标反射后,被接收子阵列的N个阵元接收,经过匹配滤波后,在t时刻雷达输出信号x(t)可建模为:
sq(t)为t时刻第q个目标的反射率,n(t)为与信号无关加性噪声,服从零均值复高斯分布。j为虚数单位。
因此,步骤二中的接收信号可以表示为
其中,为(M+N-1)2×Q维虚拟阵列导引矩阵,为包含Q个目标反射功率的向量,为噪声功率,i=vec(IM+N-1),IM+N-1表示(M+N-1)×(M+N-1)维单位矩阵,vec(·)表示向量化操作,即把矩阵中的各列依次堆叠以形成一个新的向量,(·)*表示共轭操作。从向量y中选择对应于中各虚拟阵元位置的等价虚拟信号,得到对应二阶差集和集虚拟阵列信号可表示为
其中,表示与立置上的阵元相对应的元素,max(·)表示一个集合中最大的元素,[·]l表示向量中的第l个元素,l=0,1,...,L-1, 在此基础上,在所得z向量之后补充一定数量的零,使所得向量中包含P个元素,表示为
其中P为2的整数次幂且P>L,0P-L表示P-L维的零向量。
构建一个空间谱,该空间谱的横轴表示角度θ,其与空间响应第p个元素的关系可表示为:
步骤六:根据所得空间谱进行波达方向估计。对步骤五中空间谱进行谱峰搜索操作,将其峰值按照由高到低排列,则最大的前Q个峰值对应于空间中Q个目标的波达方向估计。
本发明所提波达方向估计方法基于互质MIMO雷达结构,在利用其所接收的反射信号生成的和集差集虚拟信号基础上进行快速傅里叶变换得到其空间响应,并构建空间谱,通过对该空间谱的谱峰搜索过程可得到对空间中多个目标的波达方向估计信息。相比于传统的MIMO雷达结构,本发明所提方法充分利用了互质MIMO雷达结构特性,实现了分辨率性能的提升;相比于现有波达方向估计方法,本发明所提方法计算复杂度仅为在实时性较高的应用场景下具有显著优势,且易于实际系统上的实现。此外,本发明所提方法通过简单的零填充操作利用了全部的虚拟阵元上的信息,避免了信息损失,实现了分辨率性能的进一步提升。
下面结合仿真实例对本发明的效果做进一步的描述。
仿真实例1:利用互质MIMO雷达结构,其参数选取为M=3,N=5,即发射子阵列中包含5个阵元,接收子阵列包含3个阵元。固定采样快拍数T=500,快速傅里叶变换点数P=1024。假设空间中存在3个较为接近的目标,其相对于雷达分别位于-5°,0°,5°,信噪比固定为-10dB。利用本发明所提方法在上述条件下所得空间谱如图3所示。可以看出,本发明所提方法可以准确得到上述条件下三个目标的波达方向估计信息。
仿真实例2:采用与仿真实例1相同的互质MIMO雷达结构,测试本发明所提方法分辨率性能的极限,为体现本发明所提方法的分辨率性能优势,在本仿真中引入具有相同阵元数的传统MIMO结构以及互质MIMO雷达仅利用差集和集虚拟阵列连续部分这两种情况作为对比。假设空间中存在两个目标,相对于雷达分别位于θ1和θ2方向。其中,θ1为满足标准正态分布的随机变量,θ2=θ1+Δθ,Δθ表示两个目标固有的空间角度差。若某次估计中两个目标对应的估计值和分别满足及则认为本次估计成功分辨两个目标。本仿真在每个不同的Δθ下,分别进行500次蒙特卡洛试验,得到在该Δθ下的估计成功率如图4所示。可以看出,相比于其他两种对比情况,本发明所提方法可以在Δθ更小时将两个目标完全分辨,充分体现了本发明所提方法在分辨率性能上的优势。
Claims (3)
1.一种基于互质MIMO雷达差集和集信号快速傅里叶变换的波达方向估计方法,其特征在于,包含以下步骤:
(1)利用M+N个实际阵元构建互质MIMO雷达结构,构建方法如下:选取一组互质的整数M、N,分别用于构造发射子阵列和接收子阵列,其中,发射子阵列包含M个间距为Nd的阵元,其位置为0,Nd,...,(M-1)Nd;接收子阵列包含N个间距为Md的阵元,其位置为0,Md,...,(N-1)Md;d为发射子阵列发射信号的半波长即发射子阵列和接收子阵列部署为单基地雷达的形式;
(2)利用步骤(1)中构建的发射子阵列发射M个正交信号;假设空间远场中有Q个目标,其相对于雷达的方向表示为θ=[θ1,θ2,…,θQ]T,[·]T表示转置操作;该M个信号被空间中Q个目标反射后,其反射信号被接收子阵列的N个阵元接收,经过匹配滤波后,在t时刻雷达输出信号x(t)可建模为:
其中,sq(t)为t时刻第q个目标的反射率,n(t)为与信号无关加性噪声,服从零均值复高斯分布,j为虚数单位;
因此,步骤(2)中的接收信号可以表示为
其中,为x(t)的理论协方差矩阵, 为(M+N-1)2×Q维虚拟阵列导引矩阵,为包含Q个目标反射功率的向量,为噪声功率,i=vec(IM+N-1),IM+N-1表示(M+N-1)×(M+N-1)维单位矩阵,vec(·)表示向量化操作,即把矩阵中的各列依次堆叠以形成一个新的向量,(·)*表示共轭操作;从向量y中选择对应于中各虚拟阵元位置的等价虚拟信号,得到对应二阶差集和集虚拟阵列信号可表示为
其中,表示与位置上的阵元相对应的元素,max(·)表示一个集合中最大的元素,[·]l表示向量中的第l个元素,l=0,1,...,L-1, 在此基础上,在所得的向量z之后补充一定数量的零,使所得向量中包含P个元素,表示为
其中P为2的整数次幂且P>L,0P-L表示P-L维的零向量;
构建一个空间谱,该空间谱的横轴表示角度θ,其与空间响应第p个元素的关系可表示为:
(6)根据所得空间谱进行波达方向估计:对步骤(5)中空间谱进行谱峰搜索操作,幅度最大的前Q个峰值对应的角度,即为空间中Q个目标的波达方向估计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811218661.6A CN109471087B (zh) | 2018-10-18 | 2018-10-18 | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811218661.6A CN109471087B (zh) | 2018-10-18 | 2018-10-18 | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109471087A CN109471087A (zh) | 2019-03-15 |
CN109471087B true CN109471087B (zh) | 2021-01-12 |
Family
ID=65663892
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811218661.6A Active CN109471087B (zh) | 2018-10-18 | 2018-10-18 | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109471087B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115421157B (zh) * | 2022-10-24 | 2023-03-24 | 深圳大学 | 基于无向邻接图构建雷达阵列的方法和装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030058153A1 (en) * | 2001-09-07 | 2003-03-27 | Lockheed Martin Corporation | Adaptive digital beamforming radar method and system for maintaining multiple source angle super-resolution capability in jamming |
CN103323827A (zh) * | 2013-05-27 | 2013-09-25 | 杭州电子科技大学 | 基于快速傅里叶变换的mimo雷达系统角度估计方法 |
CN104849696A (zh) * | 2015-04-27 | 2015-08-19 | 中南大学 | 一种mimo雷达天线阵稀疏布阵方法 |
CN105403856A (zh) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | 基于嵌套式最小冗余阵列的波达方向估计方法 |
CN107576953A (zh) * | 2017-09-12 | 2018-01-12 | 成都理工大学 | 基于互质mimo阵列的相干与非相干混合目标doa估计方法 |
WO2018104468A1 (en) * | 2016-12-08 | 2018-06-14 | Iee International Electronics & Engineering S.A. | Direction of arrival estimation for automotive spread radar systems |
CN108614234A (zh) * | 2018-05-15 | 2018-10-02 | 浙江大学 | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2422210B1 (en) * | 2009-04-23 | 2013-09-11 | Groupe des Ecoles de Telecommunications Ecole National Superieure des Telecommunications | Orientation and localization system |
CN105182293B (zh) * | 2015-08-25 | 2017-07-04 | 西安电子科技大学 | 基于互质阵列mimo雷达doa与dod估计方法 |
CN107104720B (zh) * | 2017-03-01 | 2020-07-10 | 浙江大学 | 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法 |
CN107015190A (zh) * | 2017-03-01 | 2017-08-04 | 浙江大学 | 基于虚拟阵列协方差矩阵稀疏重建的互质阵列波达方向估计方法 |
-
2018
- 2018-10-18 CN CN201811218661.6A patent/CN109471087B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030058153A1 (en) * | 2001-09-07 | 2003-03-27 | Lockheed Martin Corporation | Adaptive digital beamforming radar method and system for maintaining multiple source angle super-resolution capability in jamming |
EP1423726A2 (en) * | 2001-09-07 | 2004-06-02 | Lockheed Martin Corporation | Adaptive digital beamforming radar method and system for maintaining multiple source angle super-resolution capability in jamming |
CN103323827A (zh) * | 2013-05-27 | 2013-09-25 | 杭州电子科技大学 | 基于快速傅里叶变换的mimo雷达系统角度估计方法 |
CN104849696A (zh) * | 2015-04-27 | 2015-08-19 | 中南大学 | 一种mimo雷达天线阵稀疏布阵方法 |
CN105403856A (zh) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | 基于嵌套式最小冗余阵列的波达方向估计方法 |
WO2018104468A1 (en) * | 2016-12-08 | 2018-06-14 | Iee International Electronics & Engineering S.A. | Direction of arrival estimation for automotive spread radar systems |
CN107576953A (zh) * | 2017-09-12 | 2018-01-12 | 成都理工大学 | 基于互质mimo阵列的相干与非相干混合目标doa估计方法 |
CN108614234A (zh) * | 2018-05-15 | 2018-10-02 | 浙江大学 | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 |
Non-Patent Citations (3)
Title |
---|
A Robust and Efficient Algorithm for Coprime Array Adaptive Beamforming;Chengwei Zhou 等;《IEEE Transactions on Vehicular Technology》;20170516;第67卷(第2期);第1099-1112页 * |
基于二次采样的贝叶斯波束成形;顾宇杰 等;《浙江大学学报(工学版)》;20090515;第43卷(第5期);第812-816页 * |
稀疏互质阵列的迭代加权l_1范数约束来波方向估计;刘庆华;《计算机工程》;20170430;第43卷(第4期);第110-115页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109471087A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109932680B (zh) | 一种基于平移互质阵列的非圆信号波达方向估计方法 | |
CN109471086B (zh) | 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法 | |
CN108872929B (zh) | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 | |
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN106972882B (zh) | 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法 | |
CN107104720B (zh) | 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法 | |
CN107092004B (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN110109051B (zh) | 基于频控阵的互耦阵列doa估计方法 | |
CN107290709B (zh) | 基于范德蒙分解的互质阵列波达方向估计方法 | |
WO2022127076A1 (zh) | 基于张量高阶奇异值分解的雷达角度和距离估计方法 | |
CN107561484B (zh) | 基于内插互质阵列协方差矩阵重建的波达方向估计方法 | |
CN106680815B (zh) | 基于张量稀疏表示的mimo雷达成像方法 | |
CN108710102B (zh) | 基于互质阵列二阶等价虚拟信号离散傅里叶逆变换的波达方向估计方法 | |
CN109655799A (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
CN107037392A (zh) | 一种基于压缩感知的自由度增加型互质阵列波达方向估计方法 | |
CN107064926B (zh) | 空域色噪声背景下的双基地mimo雷达角度估算方法 | |
CN107302391A (zh) | 基于互质阵列的自适应波束成形方法 | |
CN108614234B (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
CN113075649B (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
CN104868946B (zh) | 自适应加权的子阵级混合mimo‑相控阵系统的干扰抑制方法 | |
CN109471087B (zh) | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 | |
CN112711000B (zh) | 基于最小化准则的电磁矢量互质面阵张量功率谱估计方法 | |
CN115236589B (zh) | 一种基于协方差矩阵修正的极地冰下doa估计方法 | |
CN108680892B (zh) | 基于角度-空间频率域快速傅里叶变换的互质阵列波达方向估计方法 | |
CN115453450A (zh) | 一种嵌套阵下对抗互藕的稳健doa估计方法 |
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 |