发明名称 一种基于多频散射建模的昆虫尺寸测量方法
摘要 本发明公开了一种基于多频散射建模的昆虫尺寸测量方法,该方法可以同时测得昆虫的体长和体宽,基于此可以进一步反演昆虫的质量和体长体宽比,该方法所能覆盖的昆虫尺寸范围与传统的昆虫雷达尺寸测量方法相比大大增加,因此也增强了昆虫雷达对昆虫种类识别的能力。
申请公布号 CN105759252A 申请公布日期 2016.07.13
申请号 CN201610171604.1 申请日期 2016.03.24
申请人 北京理工大学 发明人 曾涛;龙腾;刘长江;王锐;胡程
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(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="FDA0000948799000000011.GIF" wi="382" he="151" />其中,<img file="FDA0000948799000000012.GIF" wi="66" he="63" />为水平极化分量,<img file="FDA0000948799000000013.GIF" wi="70" he="63" />为垂直极化分量,<img file="FDA0000948799000000014.GIF" wi="63" he="62" />和<img file="FDA0000948799000000015.GIF" wi="59" 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="FDA0000948799000000016.GIF" wi="1022" he="150" /></maths>得到第n个频点的散射截面积RCS最大值的测量数据σ<sub>max</sub>(dBf<sub>n</sub>);遍历N个频率点,得到各频率点对应的散射截面积RCS最大值的测量数据;步骤二,基于最小二乘意义的搜索分别获得参数dBf<sub>0</sub>和K的估计<img file="FDA0000948799000000017.GIF" wi="91" he="85" />和<img file="FDA0000948799000000018.GIF" wi="106" he="79" />具体为:1)、分别设定参数dBf<sub>0</sub>和K的搜索区间和搜索间隔;其中,dBf<sub>0</sub>表示散射截面积RCS最大值时的对数频率,K是影响散射的绝对强度的变量,表示目标在dBf<sub>0</sub>时的散射截面积RCS;利用电磁仿真软件仿真昆虫在不同波段的散射特性,得到的σ<sub>max</sub>关于对数频率dBf的起伏特性由以下分段函数近似表达:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&sigma;</mi><mi>max</mi></msub><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>d</mi><mi>B</mi><mi>f</mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><msub><mi>dBf</mi><mn>0</mn></msub><mo>,</mo><mi>K</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mn>3.3</mn><mo>&times;</mo><msup><mrow><mo>(</mo><mi>d</mi><mi>B</mi><mi>f</mi><mo>-</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>1.5</mn></msup><mo>+</mo><mi>K</mi></mrow></mtd><mtd><mrow><mi>d</mi><mi>B</mi><mi>f</mi><mo>&le;</mo><msub><mi>dBf</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><msup><mrow><mo>|</mo><mrow><mi>d</mi><mi>B</mi><mi>f</mi><mo>-</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>-</mo><mn>3.3</mn></mrow><mo>|</mo></mrow><mn>1.6</mn></msup><mn>0.9623</mn></mfrac><mo>+</mo><mi>K</mi><mo>-</mo><mn>3.3</mn></mrow></mtd><mtd><mrow><msub><mi>dBf</mi><mn>0</mn></msub><mo>&lt;</mo><mi>d</mi><mi>B</mi><mi>f</mi><mo>&le;</mo><mrow><mo>(</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>+</mo><mn>3.3</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><msup><mrow><mo>|</mo><mrow><mi>d</mi><mi>B</mi><mi>f</mi><mo>-</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>-</mo><mn>3.3</mn></mrow><mo>|</mo></mrow><mn>1.2</mn></msup><mn>1.5513</mn></mfrac><mo>+</mo><mi>K</mi><mo>-</mo><mn>3.3</mn></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>+</mo><mn>3.3</mn><mo>)</mo><mo>&lt;</mo><mi>d</mi><mi>B</mi><mi>f</mi><mo>&le;</mo><mo>(</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>+</mo><mn>4.4</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mo>-</mo><mn>0.35</mn><mo>&times;</mo><mi>d</mi><mi>B</mi><mi>f</mi><mo>+</mo><mi>K</mi><mo>-</mo><mn>4.7607</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mn>0.35</mn><mo>&times;</mo><mrow><mo>(</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>+</mo><mn>4.4</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><msub><mi>dBf</mi><mn>0</mn></msub><mo>+</mo><mn>4.4</mn><mo>)</mo><mo>&lt;</mo><mi>d</mi><mi>B</mi><mi>f</mi></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000948799000000021.GIF" wi="1580" he="722" /></maths>2)、针对每一个频率点对应的对数频率dBf,选取参数dBf<sub>0</sub>和K在各自搜索区间内的一组参数,从公式(5)中找到当前对数频率下σ<sub>max</sub>对应的函数表达式,将当前对数频率、选取的参数dBf<sub>0</sub>和K值代入找到的σ<sub>max</sub>函数表达式中,得到当前对数频率下的σ<sub>max</sub>计算值,将该计算值与步骤一获得的对应频率点的σ<sub>max</sub>测量值相减,得到误差值,再将误差值求平方;遍历所有对数频率点,得到各个频率点对应的误差值的平方,最后,对所有误差值的平方求和;3)、采用2)的方法,按照设定的搜索间隔,在各自的搜索区间内不断更换dBf<sub>0</sub>和K,得到对应的误差值平方的和值,选取使和值最小时对应的一组参数dBf<sub>0</sub>和K的估计<img file="FDA0000948799000000022.GIF" wi="91" he="87" />和<img file="FDA0000948799000000023.GIF" wi="74" he="71" />步骤三,昆虫尺寸计算:将估计结果<img file="FDA0000948799000000024.GIF" wi="97" he="86" />代入下式中:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mn>70.246</mn><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><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><mn>0.482</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000948799000000025.GIF" wi="780" he="78" /></maths>得到昆虫的长度2a;将<img file="FDA0000948799000000026.GIF" wi="48" he="72" />代入到下式中:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><msup><mrow><mo>(</mo><mn>58510</mn><mo>&times;</mo><msup><mn>10</mn><mrow><mover><mi>K</mi><mo>^</mo></mover><mo>/</mo><mn>10</mn></mrow></msup><mo>-</mo><mn>2.89</mn><mo>)</mo></mrow><mrow><mn>4</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000948799000000027.GIF" wi="694" he="103" /></maths>得到昆虫体积V;同时根据椭球的体积公式计算得到昆虫宽度2b:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><msqrt><mfrac><mrow><mn>3</mn><mi>V</mi></mrow><mrow><mn>4</mn><mi>&pi;</mi><mi>a</mi></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000948799000000031.GIF" wi="429" he="137" /></maths>
地址 100081 北京市海淀区中关村南大街5号