亚星游戏官网-www.yaxin868.com

专利名称：用于估计到达信号数的设备和方法
技术领域：
本发明涉及一种雷达设备，具体地说，涉及这样的雷达设备，该雷达设备 使用多个信道从物体接收包含反射雷达波的入射波并且基于指示接收信道间的 接收信号的相互关系的相关矩阵来估计到达信号数(即估计不同反射波的数 目)。
背景技术：
通常，下面的雷达设备是已知的。该雷达设备使用由多个，部件配置的 阵列天线，以及估计同时到达该阵列天线的多个无线电波的到达方向(DOA)。
作为一种估计无线电波的DOA的方法，已知的是多信号分类(MUSIC) 方法、经由旋转恒定技术方法对信号参数的估计(ESPRIT)等等。在这些方法 中，基于指示由每个天线部件(还被称为信道)所接收的接收信号间的相关关 系的相关矩阵，产生方向谱。ilil被扫描的方向谱，来实施高5^,率估计。
下面将描m MUSIC方法的和M。 M3i所谓的线性阵列来形成阵列天线， 在所谓的线性阵列中，N ^h^部件(N为2或大于2的整数)被等间距;NM 齐。
首先，由等式(1)表达的接收向量X(k)被配置用于^样时间KAT经由 阵列織采集的采样数据X!(k)、 X2(k)至XN(k)的片段(AT为采样间隔，以及k 为自然数)。接下来，依照等式(2)使用接收向量X(k)来确定N行和N列的自
相关矩阵Rx^
这里，T是向量转置。H是复共轭转置。
X(k) = {Xl(k)， x2(k)，.. .,XN(k)}T 等式(l)Rxx=X(k)XH(k) 等式(2)
接下来，确定自相关矩阵Rxx中的特征值^至A^ (其中、A^……AN)。 从大于预先设定的噪声阈值TH的特征值个数来估计到达信号的数量L (<N)。 此外，计算对应于特征值^至A^的特征向量A至eN。
然后，由对应于(N-L)个特征值的特征向量组成的噪声特征值向量EN0 由等式G)限定，所述(N-L)等于或小于噪声阈值TH。使用关于方向e表示 的阵列天线的复响应a(e)来确定由等式(4)表示的性能函数PMu(e)。
Eno= {eL+i,eL+2,……,eN〉 等式(3)
膨,
当e匹配入射雷达波的doa时，从性能函数PMu(e)获得的方向谱(称为
"MUSIC谱")分叉并形戯峰。因此，通过搜索MUSIC谱中的峰(也就是， 零点)，可以确定DOA的估计值et至9L。
如上所述，在MUSIC方法中(同样在ESPRIT方法中)，在计算DOA的 过程期间，需要精确估计到达信号的数目L。因此，重要的是，适当地设置噪 声阈值TH。
噪声可以由各种因素引起。例如，在日本专利特开公开物No.200647282 中提出了下面设备。该设备关注于噪声强度随频率改变。该设备基于节拍信号 的频率设定不同噪声阈值TH (也就是，距离物体等的距离)。
当可以确保多个快照(snapshot)时，基于最大似然方法的Akaike信息标 准(AIC)、最小描述长度(MDL)等等是已知的。
特征值^至U与接收长度具有相关性。因此，通过强无线电波等的接收， 与信号元素和噪声元素无关地，可以增强整体接收强度。在这种情况中，如图7 中所示，基于噪声元素的特征值变得大于噪声阈值TH。因此，错误地将特征值 确定为信号元素的那些，导致对到达信号数的估计值增加。
特别地，己知的是，譬如这种状态发生在车辆一安装的雷达设备等上，该 设备被用在包括来自道路等的大量混乱干扰的环境中(也就是，i!31反射产生 的不必要的无线电波)。在不能确保足够数量的快照时及在必须实时使用AIC或 MDL实舰方向的检测时，对到达信号数量的估计精度明显斷氐。
尤其是，如同频率调制的连续波(FMCW)雷达的情况，在单个测量瞬间
一次仅仅可以获得单个物体的一个快照(在节拍信号的频率谱中仅仅有一个频 率峰)时，则在短时期内不能确保大量快照。

发明内容
为了解决战问题，实现了本发明。本发明的目的在于提供一种即使在多 个决照小的时候也能高精度地估计到达信号数的方法，以及一种使用到达信号 数估计方法的雷达设备。
在为了实现上述目的而作出的本发明的方法中，基于相关矩阵确定对应于 指示多个信道间的相互关系的自相关矩阵的特征值，所述多个信道从反射雷达 波的物体接收入射雷达波(即，反射波)。在确定的特征值中选择最大特征值作 为参考特征值。将具有大于预定阈值的参考特征值的比率的(在许多特征值中 的)这些特征值识别为信号空间的特征值，以及比率等于或小于阈值的特征值 被识别为噪声空间中的特征值。被确定为{言号空间的特征值的数量为到达{言号 的数量。
在本发明的用于估计到达信号数的方法中，{顿对参考特征值的比值(即 相对大小)作出识别。因此，即使决照数小时或整体接收强度增加时，也可以
精确地识别特征值，进而可以精确地f古计到达〗言号数L。
接下来，在本发明的雷达设备中，发射及接收装置具有发射雷达波和从反 射雷达波的物体接收入射的雷达波的多个信道。基于从每个信道获得的接收信 号，矩阵产生装置产生指示信道间相互关系的相关矩阵。特征值计算装置计算 对应于由矩阵产生装置产生的相关矩阵的特征值。
然后，在确定的特征值中选择最大特征值作为参考特征值，用于估计到达 信号数的装置在对参考特征值的比值大于预定阈值的特征值中识别特征值作为 信号空间中的特征值，以及识别其比值等于或小于阈值的特征值作为噪声空间 中的特征值。被识别作为信号空间中的特征值的特征微为至哒的信号数。
基于到达信号数和由用于估计到达信号数的装置所识别的噪声空间中的特
征值，方向估计装置估计^h引入雷达波的到达方向。
换言之，本发明的雷达设备使用用于估计到达信号数的方法，并且肖,实 现与该方法一样的效果。再者，该雷达设备可以增强入射雷达波的到达方向的
估计精度。
已知的是，在j細特征值估计弓l入雷达波的至哒方向时，可检观啲到达信 号的最大数比信道数小一。因此，使用下面条件，用于估计引入雷达波的装置 优选限制到达信号数至引入雷达波的最大数。在对参考特征值的比值大于阈值 的特征值的数超过所设定的到达信号数时，该装置限制该数量至小于信道数的 数值。
在本发明的雷达设备中，在M^I寸和接收装置的发射和接收的雷达波为 频率调制连续波时，矩阵产生，ttit地为每个信道从接收信号中确定节拍信 号，进而产生用于每个频率的相关矩阵，在该每个频率处节拍信号的频率谱形 戯峰。
在这种情况中，节拍信号的频率谱形成尖峰处的频率il^存在一个明显的
物体，该物体具有由该频率所确定的距离和相对速度。因此，不管在噪声空间 中的参考特征值(换言之，不管到达信号数为零)，使用这个频率可以阻止到达 信号数被确定。换言之，参考特征值显然为信号空间的特征值。结果，可以提 高检测结果的可靠性。


图1是雷达设备的整体结构的框图2是测量处理细节的流程图3是方向估计处理细节的流程图4是到达信号数估计处理细节的流程图5至6为显示将该雷达设备与常规设备相比较的试验结果的图表；以及
图7是在常规设备中出现的问题的说明图。
具体实施例方式
参考附图将描述本发明的实施例。
图1是应用本发明的雷达设备2的整体结构的框图。 雷达设备2被安装在车辆上。该雷达设备2被配置为用于车辆的物体识别 设备中的一部分，其识别存在于 前方的物体。
如图1中所示，雷达设备2包括数模(D/A)转换器10、压控振荡器(VCO)
14、分配器16和，天线18。根据调制指令C， D/A转换器10产生三角波调 制信号M。经由缓冲器12，将由D/A转换器10产生的调制信号M施加给VCO 14。根据调制信号M， VC014改变振荡频率。分配器16执行功率分配，将来 自VCO 14的输出分为，信号Ss和本地信号LO。根据发射信号Ss，发射天 线18^M雷达波。
雷达设备2还包括接收一端天线部分20、接收开关22、混合器24、放大 器26、模数(A/D)转换器28和信号处理部分30。接收一端天线部分20为由 接收雷达波的N (N为2或大于2的整数)^线配置成的阵列天线。接收开 关22交替选择任一个天线作为接收一端天线部分20。接收开关22将来自所选 择端的信号提供给后续阶段作为接收信号Sr。混合器24将由接收开关22提供 的接收信号Sr与本地信号LO混合，进而产生节拍信号B。放大器26放大由混 合器24产生的节拍信号B。 A/D转换器28对由放大器26放大的节拍信号B进 行采样，以及将被采样的节拍信号B转换为数字数据D。信号处理模块30输出 调制指令C给D/A转换器10，以^tM5i A/D转换器28获得的数字数据D执 行信号处理。
此后，分别榭言道chl至chN分配给配置接收一端天线部分20的*天 线。Sn表/询个信道chi (i = l， 2，……，N )的接收信号。Bi表示基于接收 信号Sn产生的节拍信号。Di表示从被采样的节拍信号Bi转换的数字数据。
在如上所述的配置的雷达设备2中，根据调帝K言号M，分配器16执行由 VCO14所产生的高频信号(即，频率调制连续波[FMCW])的功率分配。结果 是，产生发射信号Ss和本地信号LO。经由，天线18，将，信号SsMt作 为雷达波。
从鄉天线18鄉的且由物体反射之后返回的雷达波(S卩，弓l入的雷达波) 由配置接收一端天线部分20的每个天线(信道chl至cliN)所接收。然而，只 有由接收开关22所选择的信道chi (i=l到N)的接收信号Sri被提供给混合器 24。然后，混合器24耦合来自分配器16的接收信号Sri与本地信号LO，进而 产生节拍信号Bi。 A/D转换器28对由放大器26放大的节拍信号Bi进行釆样。 被采样的节拍信号Bi被加载进信号处理部分30作为数字数据D。
信号处理部分30主要由已知的M计^m构成，该微型计算机包括中央处 理单元(CPU)、只读存储器(ROM)和随机存储器(RAM)。此外，信号
部分30包括被用于对经由A/D转换器28获得的M执行快速傅立叶变换(FFT) 的计算处理设备(诸如数字信号处理器[DSP])等。
基于经由A/D转换器28获得的数字数据D，信号处理部分30中的CPU 执行测量处理。在测量处理中，信号处理部分30计算反射雷达波的物体的距离 和相对速度，进而估计物体所存在的方向。
<测量处理>
参考图2中的流程图，^SMil信号处理模块30中的CPU所重复执行
在开始处理时，首先，在SllO， CPU发送调制指令C给D/A转换器10， 从A/D转换器28获取数字数据D，以及进行到S120。在S120， CPU对在SllO 获取的用于*信道的数据D执行FFT处理，从而计算用于每个信道的节拍信 号的功率谱。
在随后的S130，基于在S120计算的功率谱，使用应用至FMCW雷达的已 知方法(省略了对其细节的解释)，CPU确定与反射雷达波的物体间的距离和相 对皿。然后，CPU进行到S140。
在S140，基于在S120所确定的用于^信道的功率谱，CPU执行方向估 计处理，用于估计存在于S130处提取的物体的方向。然后，CPU结束该处理。
<方向估计处理〉
接下来，参考图3中的流程图，)l繊述在上述S140处实施的方向估计处理 (MUSIC处理)的细节。
首先，在S210处，因为存在基于来自物体的弓l入的雷达波的信号元素，从 在向上扫描或向下扫描的频率调制期间产生的任一分离功率谱中，CPU从在上 述S130中提取的频率(槽(bin))中选择一个频率，该频率仍没有313i方向估 计处理被处理。然后，CPU进行至S220。
在S220， CPU产生接收谱X(i)(参考等式(7))。接收谱X(i)为从全部信 道chl至chN的功率谱中提取的所选择的频率的信号元素阵列(FFT处理结果 数据)。在随后的S230，基于接收向量X(O，根据等式(8)， CPU产生自相关 矩阵Rxx，基于接收向量X(i)，进而进行至S240。
X(i) = {Xl(i),x2(i),... ,XN(i)}T 等式(7 )
R欢二x(i)x"(i) 等式(8)
在步骤S240， CPU计算在S230产生的自相关矩阵Rxx的特征值人!至入n。 特征值、至U被按照最大mJ顿序依次排列。
在随后的S250, CPU在计算的特征值^至^内识另腊号空间中的特征值 和噪声空间中的特征值。CPU还执行到达信号数的估计处理，以估计信号空间 中的特征f微为到达的信号数L。然后，CPU进行至S260。
在S260，基于到达信号数的估计处理的估计结果，CPU计算MUSIC谱， 以及进行至S270。
具体地，基于对应于噪声空间的(N-L)个特征值、w至入M的特征向量 eL+说+2，……,ew，通过等式(9)定义噪声特征向量ENo。然后，j顿a(e沐定义 由等式(10)表达的性能函数PMU(e)，所述a(e)表示接收一端天线部分20相
关于方向e的复响应。根据性能函数pmu(e)确定的方向谱为music谱。
ENo={eL+i，eL+2,...,e^ 等式(9)
尸(0 = (。,_ 等式(10)
"(举w。Aw。，)
在S270， CPU对在S260确定的MUSIC谱执行空扫描。结果是，CPU确 定由^天线所接收的到达信号数L的入射角e,至OL，所述旨天线配置接收 端天线部分20。换言之，CPU确定存在反射入射雷达波的物体的方向。在随后 的S280, CPU判断对于抽取的所有频率(槽)，该处理是否完成，因为存在基 于来自物体的入射雷达波的信号元素。
在存在没有处理的频率(槽)时，CPU返回至S210，并舰没有处理的频 率(槽)重复上述处理(S210至S270)。在所有频率(槽)已经被处理时，CPU 结束该处理。
<到达信号数的估计处理>
这里，参考图4中的流程图，^| 述在上述s250所执行的到达信号数的估 计处理细节。
在开始该处理时，首先，在S310， 4OT等式(11)， CPU计算^^寺征值 入2至A^对参考特征值、的比值R^至R^ (特征值比值)。参考特征值为最大特
征值X"然后，CPU进行至S320。
^^101ogJ夺) 等式(11)
在S320, CPU初始化用于标识特征值、至;w的参数i为1 ，以及初始化到
达信号数L为1 (提前计数参考特征值^)。然后，CPU进行至S330。
在S330， CPU增加参数i (即，i=i+l)。在随后的S340， CPU判断特征值 比值是否大于噪声阈值TH。使用通过试验等预先设定的值作为噪声阈值 TH。
在特征值比值R^大于噪声阈值TH时，将特征值、考虑为信号空间中的特 征值。CPU进行至S350,以及增加到达信号数L (即，L=L+1 )。在随后的S360， CPU判断参数i是否小于天线数N。在参数i小于天线数N时，CPU返回至S330, 并且重复比较特征值比值RXi和噪声阈值TH。
同时，在S340处CPU判iff寺征值比值RAi是否等于或小于噪声阈值TH， 或者在S360处判断参数i是否等于或大于天线数N， CPU进行至S370。 CPU 判断到达信号数L是否大于到达信号最大可识另擞L。基于天线数N，设定 到达信号的最大可识别数Ln^为11^-1。
当到达信号数L等于或小于到达信号的最大可识另擞时，CPU结束该 处理。在到达信号数L大于到达信号的最大可识别数L时，CPU进行至S380。 CPU把到达信号的最大可识别数L^作为到达信号数L (即，IM^)以完成 处理。
依据实施例，VC014、分配器16、发射天线18、接收一端天线部分20和 接收开关22等^^发射和接收^g。 S220至S230等^^矩阵产生装置。S240 等^^特征值计算，。S250 (S310至S380)等j介于引入波估计装置。S260 至S270等价于方向估if^置。
如上所述，雷达设备2产生指示多个信道chl至chN之间的相互关系的相 关矩阵，所述信道接收来自反射雷达波的物体的引入雷达波。然后，雷达设备2 确定相关矩阵的特征值。在确定的特征值之中的最大特征值为参考特征值^。 在其对参考特征值^的比值RA^ (=101OgK)(V^))大于噪声阈值TH中的特征值 入2至Xn中，特征值被识别为信号空间的特征值。等于或小于噪声阈值TH的特 征值被识别为噪声空间的特征值。被识别为信号空间的那些的特征值数为到达
信号数L。
这样，因为雷达设备2使用对参考特征值^的比值R^，也就是相对大小，
识别特征值，所以甚至在快照数小时或在整体接收强度被增加时，雷达设备2 也能够精确识别特征值和估计到达信号数L。结果，可以精确检测到入射雷达 波的DOA (换言之，反射雷达波的物体的位置)。
这里，图5至6为在具有相对速度0km/h的单，##在于安装有雷达设 备的 前大约4米处、在大约4。角度的方向上时，测量物体的相对速度Vr、 垂直位置dist和水平位置x时的结果的图表。
图5显示使用常规设备时的结果，该设备比较特征值自身与阈值且识别特 征值。图6显示4顿依据本实施例的雷达设备2的结果，该设备比较特征值比 值与阈值且识别特征值。
基于Mil FMCW雷达的已知方法获得的距离物体的距离和ilil MUSIC方 法获得的物体的方向，垂直位置dist和水平位置x fflil笛卡尔坐标系统指示物 体的位置。水平位置x为车辆宽度方向上的位置。垂直位置dist为垂直于 宽度方向的方向(沿道路表面的方向)上的位置。
如图5中所示，在常规雷达中，计算的弓l入雷达波数在1和3之间变化。 然而，在雷达设备2中，引入雷达波数L被精确估计为1。清楚的是，可以稳 定地确定垂直位置dist和水平位置x。
<其他实施例>
在上面已经描述了本发明的实施例。然而，本发明不限于上述实施例，在 不脱离本发明的精祌的情况下，可以作出各种改变。
例如，依据实施例，提供一个发射天线和多个接收天线。然而代替之，可 以提供多个发射天线和一个或多个接收天线。在这些情况中，发射天线和接收 天线之间的每一组合可以作为信道。
依据实施例，描述了将本发明应用至FMCW雷达的例子。然而，可以将 本发明应用至具有多个接收入射雷达波的信道的任何设备，以及从指示每个信 道间相互关系的相关函数的特征il^估计入射雷达波的DOA。
权利要求
1、一种用于估计到达信号数的方法，其中该方法产生指示多个信道间的相互关系的相关矩阵以及基于该相关矩阵估计到达信号数，所述多个信道接收来自反射雷达波的物体的入射雷达波，该方法包括计算对应于相关矩阵的特征值；从计算的特征值中选择参考特征值；在计算的特征值中，识别其对参考特征值的比值大于预定阈值的特征值；和设定识别的特征值数作为到达信号数。
2、 依据权利要求l的方法，其中参考特征值为特征值中的最大特征值。
3、 一种雷达设备，包括发射和接收装置，其具有发射雷达波和从反射雷达波的物体接iB(A射雷达波的多个信道；矩阵产生装置，其用于基于从每个信道获得的接收信号产生指示信道间相 互关系的相关矩阵；特征值计算装置，其用于计算对应于由矩阵产生装置所产生的相关矩阵的 特征值；入射波估计装置，其用于在由特征值计算，所计算的特征值内识别作为 第一和第二组特征值的特征值，以及设定所识别的作为第一组特征值的特征值 数为到达信号数；禾口方向估计装置，其基于第二组特征值和到达信号数估计每个入射雷达波的 到达方向，其中入射波估计，从特征值中选择参考特征值，以及在其对参考特征值 的比值大于阈值的特征值中识别特征值作为第一组特征值，以及识别其对参考 特征值的比值等于或小于阈值的特征值作为第二组特征值。
4、 依据权利要求3的雷达设备，其中参考特征值为特征值中的最大特征值。
5、 依据权禾腰求4的雷达设备，其中在其对参考特征值的比值大于阈值的特征值数超过被设定为小于信道数的 数值的最大到达信号数时，入射波估计装置将到达信号数限制为最大到达信号数。
6、 依据权利要求4的雷达设备，其中Mmt和接收装置发射和接收的雷达波为频率调制连续波；禾口根据旨信道的接收信号，所述矩阵产生装置确定节拍信号，以及产生每 个频率的相关矩阵，在所述每个频率处节拍信号的频率谱达到峰值。
7、 依据权利要求5的雷达设备，其中ffl31Mt和接收装置发射和接收的雷达波为频率调制连续波；禾口 根据^M言道的接收信号，所述矩阵产生装置确定节拍信号，以及产生每 个频率的相关矩阵，在所述每个频率处节拍信号的频率谱达到峰值。
8、 一种用于估计到达信号数的设备，该设备产生指示多个信道间相互关系 的相关矩阵以及基于该相关矩阵估计到达信号数，所述多个信道接收来自反射 雷达波的物体的入射雷达波，其中计算对应于相关矩阵的特征值，和从计算的特征值中选择参考特征值； 顿参考特征值的比值大于预定阈值的特征值中，识别特征值；禾口设定识别的特征ttm作为入射波数。
9、 依据权利要求8的到达信号波数估计设备，其中参考特征值为特征值中 的最大特征值。
全文摘要
本发明涉及用于估计到达信号数的设备和方法。在对应于指示多个信道间相互关系的相关矩阵的特征值中，确定最大特征值作为参考特征值λ₁，所述多个信道接收来自反射雷达波的物体的入射雷达波。计算每个特征值λ₂至λ_N对参考特征值λ₁的比值Rλ_i(＝10log₁₀(λ_i/λ₁))。在参考特征值λ₁和其特征值比值Rλ_i大于噪声阈值TH的特征值λ₂至λ_N中，识别特征值作为信号空间中的特征值。其特征值比值Rλ_i等于或小于噪声阈值TH的特征值被识别作为噪声空间中的特征值。将被识别作为信号空间中的特征值的特征值数计数为到达信号数。
文档编号G01S13/00GK101349752SQ20081021470
公开日2009年1月21日 申请日期2008年7月18日 优先权日2007年7月20日
发明者清水耕司, 阿部好浩 申请人:株式会社电装