发明名称 一种基于多波段散射测量的昆虫体长-质量联合测量方法
摘要 本发明公开了一种基于多波段散射测量的昆虫体长‑质量联合测量方法,该方法可以同时测得昆虫的体长和质量,该方法所能覆盖的昆虫尺寸范围与传统的昆虫雷达体型参数反演方法相比大大增加,因此也增强了昆虫雷达对昆虫种类识别的能力。
申请公布号 CN105866770A 申请公布日期 2016.08.17
申请号 CN201610171893.5 申请日期 2016.03.24
申请人 北京理工大学 发明人 胡程;王锐;刘长江;曾涛;龙腾;李卫东
分类号 G01S13/88(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/88(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 李微微;仇蕾安
主权项 一种基于多波段散射测量的昆虫体长‑质量联合测量方法,其特征在于,包括如下步骤:步骤一,多波段RCS数据获取:采用多波段雷达系统探测待测昆虫,获得N个频点的散射信号,其中,每个波段的频点数不小于10个;各频点对应的对数频率为dBf<sub>n</sub>=10log<sub>10</sub>(f<sub>n</sub>),f<sub>n</sub>为第n个频点的频率值,n=1~N;针对第n个频点,获得待测昆虫对应的全极化散射矩阵:<img file="FDA0000948799100000011.GIF" wi="381" he="151" />其中,<img file="FDA0000948799100000012.GIF" wi="64" he="63" />为水平极化分量,<img file="FDA0000948799100000013.GIF" wi="62" he="61" />为垂直极化分量,<img file="FDA0000948799100000014.GIF" wi="63" he="63" />和<img file="FDA0000948799100000015.GIF" wi="66" he="63" />为交叉极化分量;将该雷达全极化散射矩阵S<sub>n</sub>,代入下式中:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mrow><mo>(</mo><msub><mi>dBf</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mo>|</mo><msubsup><mi>S</mi><mrow><mi>h</mi><mi>h</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mrow><mi>v</mi><mi>v</mi></mrow><mi>n</mi></msubsup><mo>|</mo><mo>+</mo><msqrt><mrow><msup><mrow><mo>|</mo><mrow><msubsup><mi>S</mi><mrow><mi>h</mi><mi>h</mi></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>S</mi><mrow><mi>v</mi><mi>v</mi></mrow><mi>n</mi></msubsup></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><mo>|</mo><msubsup><mi>S</mi><mrow><mi>v</mi><mi>h</mi></mrow><mi>n</mi></msubsup><mo>|</mo></mrow></msqrt><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000948799100000016.GIF" wi="1021" he="147" /></maths>得到第n个频点的RCS最大值的测量数据σ<sub>max</sub>(dBf<sub>n</sub>);遍历N个频率点,得到各频率点对应的RCS最大值的测量数据,将各频率点对应的对数频率作为横坐标,将RCS最大值的测量数据作为纵坐标,形成测量数据曲线图;再基于该测量数据曲线图得到对数频率dBf与RCS最大值的测量数据的幂函数关系:σ<sub>max</sub>=α(dBf‑dBf<sub>0</sub>)<sup>β</sup>+σ<sub>0</sub>;其中,dBf<sub>0</sub>表示曲线图中峰值拐点的横坐标,σ<sub>0</sub>为峰值拐点的纵坐标,即为RCS最大值的测量数据;α和β为决定幂函数曲线陡峭程度的参数;步骤二,基于最小二乘意义的搜索获得对峰值拐点参数dBf<sub>0</sub>和σ<sub>0</sub>的估计值<img file="FDA0000948799100000017.GIF" wi="93" he="86" />和<img file="FDA0000948799100000018.GIF" wi="83" he="102" />具体为:1)、选取所述步骤一的测量数据曲线图中峰值拐点的窗口内的RCS测量数据σ<sub>max</sub>(dBf<sub>m</sub>),m=N<sub>s</sub>~N<sub>e</sub>,N<sub>s</sub>≥1,N<sub>e</sub>≤N;2)、分别设定α、β、dBf<sub>0</sub>和σ<sub>0</sub>的遍历搜索区间和遍历搜索间隔;3)、选取α、β、dBf<sub>0</sub>和σ<sub>0</sub>在各自搜索区间范围内的一组参数,代入所述幂函数关系式中,得到其中一个频率点上的计算值,将该计算值与步骤一获得的对应频率点的测量值相减,得到误差值;遍历所述窗口内所有频率点,得到各个频率点对应的误差值,最后,对所有误差值求和;4)、按照设定的搜索间隔,在各自的搜索区间内不断更换参数α、β、dBf<sub>0</sub>和σ<sub>0</sub>,得到对应的误差值的和值,选取使误差值的和值最小时对应的一组参数作为第一峰值拐点参数的估计结果:步骤三,体长和质量反演:将对第一拐点处对数频率的估计结果中的dBf<sub>0</sub>的估计值<img file="FDA0000948799100000021.GIF" wi="98" he="83" />转换为对波长的估计<img file="FDA0000948799100000022.GIF" wi="86" he="85" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>&lambda;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><mi>c</mi><mo>/</mo><mrow><mo>(</mo><msup><mn>10</mn><mrow><mi>d</mi><mover><mi>B</mi><mo>^</mo></mover><msub><mi>f</mi><mn>0</mn></msub><mo>/</mo><mn>10</mn></mrow></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000948799100000023.GIF" wi="358" he="135" /></maths>将估计结果<img file="FDA0000948799100000024.GIF" wi="50" he="85" />直接代入下式中:l=0.521λ<sub>0</sub>+4.7492;得到对昆虫体长的估计<img file="FDA0000948799100000025.GIF" wi="61" he="70" />同时,将估计结果中的σ<sub>0</sub>的估计值<img file="FDA0000948799100000026.GIF" wi="56" he="83" />代入下式中:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>m</mi><mo>^</mo></mover><mo>=</mo><msup><mn>10</mn><mrow><mn>0.125</mn><mover><msub><mi>&sigma;</mi><mn>0</mn></msub><mo>^</mo></mover><mo>+</mo><mn>6.69</mn></mrow></msup></mrow>]]></math><img file="FDA0000948799100000027.GIF" wi="309" he="71" /></maths>获得对昆虫质量的估计<img file="FDA0000948799100000028.GIF" wi="78" he="54" />
地址 100081 北京市海淀区中关村南大街5号
您可能感兴趣的专利