主权项 |
高分辨率人体目标运动特征检测方法,其特征在于按照以下步骤进行:步骤1:建立包含有微多普勒特征的雷达解调目标回波信号模型,所述解调目标回波模型具体为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>LO</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>·</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msub><mi>S</mi><mi>Ri</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mn>2</mn><mi>πjαt</mi></mrow></msup><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msup><mi>e</mi><mrow><mn>2</mn><mi>πjα</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>2</mn><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>/</mo><mi>c</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msup><mi>e</mi><mrow><mo>-</mo><mn>2</mn><mi>πj</mi><mrow><mo>(</mo><msup><mrow><mn>2</mn><mi>λ</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></msup></mrow>]]></math><img file="FDA0000582413360000011.GIF" wi="1782" he="183" /></maths>其中,回波信号S(t),S<sub>LO</sub>(t)为雷达本振信号,S<sub>Ri</sub>(t)(i∈[0,N])为人体躯干及各运动关节回波信号,i=0时,S<sub>R0</sub>(t)为人体躯干的回波信号,N为产生回波的人体运动关节总个数,j为虚数单位,α为雷达载波频率,c和λ分别是光速和雷达信号波长,R<sub>0</sub>(t)为人体躯干模型:R<sub>0</sub>(t)=R<sub>0</sub>+vt,R<sub>i</sub>(t)为人体各关节理想条件下的运动简易模型:R<sub>i</sub>(t)=R<sub>0</sub>(t)+r<sub>i</sub>sin(2πf<sub>p</sub>t+θ<sub>i</sub>),其中,R<sub>0</sub>为人体与雷达的距离,v为人体相对雷达的径向速度,r<sub>i</sub>为各关节的幅度,θ<sub>i</sub>为各关节的相位,f<sub>p</sub>为人体走动时各关节的频率;步骤2:对回波信号S(t)过采样,得到离散的信号S[n],对S[n]做时频分析,得到离散化的时频能量分布W<sub>s</sub>(n,f),所述时频能量分布W<sub>s</sub>(n,f)是对<img file="FDA0000582413360000012.GIF" wi="82" he="96" />的近似估计,其中,<img file="FDA0000582413360000013.GIF" wi="286" he="132" />表示人体各个目标相对于雷达的径向速度;步骤3:对时频能量分布矩阵W<sub>s</sub>(n,f)进行处理优化,得到躯干及各关节对应时频谱的曲线点矩阵;步骤4,对曲线点矩阵进行Hough直线检测,得到人体躯干对应直线的斜率k和截距d,根据公式v=f<sub>D</sub>λ/2可求得主体躯干的速度信息,其中f<sub>D</sub>为躯干对应曲线的点对点即时频率,λ为雷达信号波长,再用得到的斜率和截距对曲线点矩阵H<sub>X×Y</sub>进行修正,将各点按照人体躯干的直线在矩阵中心位置进行补偿,即根据求出的截距d将矩阵每列都向下移动<img file="FDA0000582413360000021.GIF" wi="364" he="162" />个位置,使得矩阵第一列表示躯干位置的点在矩阵中间的行,再根据斜率k求出每列应该下移的行数<img file="FDA0000582413360000022.GIF" wi="263" he="94" />得到修正后的曲线点矩阵;步骤5,对修正后的矩阵按列求加权和,并作FFT变换,求得各时频谱曲线的周期T,也就是人体各关节的运动周期,再建立幅度‑相位参数空间(r,θ),其中r为曲线幅度参数,θ为相位参数,建立累加器Σg(r,θ),将参数空间中的每个点,即每一组(r,θ)参数值带入之前的运动简易模型R<sub>i</sub>(t)=r<sub>i</sub>sin(2πft+θ<sub>i</sub>)中对修正后的矩阵H<sub>X×Y</sub>中各曲线点进行计算,将符合该参数下模型的曲线点进行累加,找出累加器g中局部最大点,各局部最大点在累加器中对应的幅度值和相位值即为各条曲线的幅度和相位的具体数值,由于修正后的矩阵H<sub>X×Y</sub>的行表示的频率方向点数,列表示的时间方向点数,在矩阵图像中即为图像的纵坐标和横坐标,根据时频谱曲线上各点在纵坐标方向的位置可以求得该点在频率方向对应的频率,根据v=f<sub>D</sub>λ/2即可得出各个关节的即时速度信息,再根据运动简易模型R<sub>i</sub>(t)=r<sub>i</sub>sin(2πft+θ<sub>i</sub>)可知各关节相对躯干最大速度和幅度关系为v<sub>i_max</sub>=2πfr<sub>i</sub>,可以求得各个关节的幅度信息。 |