发明名称 一种无线网络系统中基于信号到达时间的非合作定位方法
摘要 本发明公开了一种无线网络系统中基于信号到达时间的定位方法,特点是根据布置好的发射机以及每个接收机的位置、每个基站的测量噪声功率以及信号从发射机到每个接收机的传播时间,构造一个加权最小二乘关系式,然后构造拉格朗日函数,并利用二分法得到最优拉格朗日乘子,最后利用最优拉格朗日乘子得到目标源位置的最终定位值,优点在于可以确保得到全局最优解而不受局部收敛的影响,因而定位精度高且非常稳健,并且只需要求解一个参数的非线性方程,因此其计算复杂度低。
申请公布号 CN103969620B 申请公布日期 2016.04.27
申请号 CN201410154272.7 申请日期 2014.04.17
申请人 宁波大学 发明人 何杉;王刚;金明;李有明
分类号 G01S5/00(2006.01)I 主分类号 G01S5/00(2006.01)I
代理机构 宁波奥圣专利代理事务所(普通合伙) 33226 代理人 程晓明
主权项 一种无线网络系统中基于信号到达时间的定位方法,其特征在于包括以下步骤:(1)首先在无线网络环境中建立一个平面或空间坐标系作为参考坐标系,设一个目标源、一个发射机和N个接收机存在于该无线网络中,发射机在参考坐标系的坐标记为s<sub>0</sub>,接收机在参考坐标系的坐标记为s<sub>1</sub>,...,s<sub>N</sub>,目标源在参考坐标系的坐标记为x,发射机向目标源发射信号,经目标源反射后由接收机接收,测量信号从发射机发出到接收机接收所经历的时间t<sub>i</sub>,进而得到信号传输的距离d<sub>i</sub>=ct<sub>i</sub>,其中c为光速,i=1,...,N,N≥3;(2)根据接收机和发射机的坐标以及信号传输的距离d<sub>i</sub>,采用非线性加权最小二乘方法估计目标源的位置,表示为以下关系式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>y</mi><mo>=</mo><msup><mrow><mo>[</mo><msup><mi>x</mi><mi>T</mi></msup><msub><mrow><mo>,</mo><mi>r</mi></mrow><mn>0</mn></msub><mo>]</mo></mrow><mi>T</mi></msup></mrow></munder><msup><mrow><mo>(</mo><mi>Ay</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow><mi>T</mi></msup><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>Ay</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000492252360000011.GIF" wi="606" he="103" /></maths>s.t.||x‑s<sub>0</sub>||<sup>2</sup>=r<sub>0</sub><sup>2</sup>,其中y=[x<sup>T</sup> r<sub>0</sub>]<sup>T</sup>为优化变量,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>2</mn><msup><mrow><mo>(</mo><msub><mi>s</mi><mn>1</mn></msub><mo>-</mo><msub><mi>s</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><mo>-</mo><mn>2</mn><msub><mi>d</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msup><mrow><mn>2</mn><mrow><mo>(</mo><msub><mi>s</mi><mi>N</mi></msub><mo>-</mo><msub><mi>s</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mi>T</mi></msup></mtd><mtd><mo>-</mo><mn>2</mn><msub><mi>d</mi><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000492252360000012.GIF" wi="549" he="246" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>s</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>s</mi><mn>0</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>d</mi><mn>1</mn><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>s</mi><mi>N</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>s</mi><mn>0</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>d</mi><mi>N</mi><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000492252360000013.GIF" wi="495" he="263" /></maths><img file="FDA0000492252360000014.GIF" wi="422" he="246" /><img file="FDA0000492252360000015.GIF" wi="191" he="79" />为接收机测量的噪声功率,“||·||”是欧几里德范数,“T”为矩阵转置,“s.t.”表示“受约束为”,“min”表示“使最小化”;(3)将步骤(2)中的非线性加权最小二乘关系式等价写为如下关系式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munder><mi>min</mi><mi>y</mi></munder><msup><mrow><mo>(</mo><mi>Ay</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow><mi>T</mi></msup><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>Ay</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000492252360000016.GIF" wi="550" he="103" /></maths>s.t.y<sup>T</sup>Dy+2f<sup>T</sup>y+||s<sub>0</sub>||<sup>2</sup>=0,其中<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mi>n</mi></msub></mtd><mtd><msub><mn>0</mn><mrow><mi>n</mi><mo>&times;</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mn>0</mn><mrow><mi>n</mi><mo>&times;</mo><mn>1</mn></mrow><mi>T</mi></msubsup></mtd><mtd><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000492252360000017.GIF" wi="358" he="159" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msub><mi>s</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000492252360000018.GIF" wi="219" he="157" /></maths>n为参考坐标系的坐标维数,n=2或3,“I<sub>n</sub>”为n×n单位矩阵,“0<sub>n×1</sub>”为n×1零向量;然后构造上述加权最小二乘关系式的拉格朗日函数L(λ)=(Ay‑b)<sup>T</sup>Q<sup>‑1</sup>(Ay‑b)+λ(y<sup>T</sup>Dy+2f<sup>T</sup>y+||s<sub>0</sub>||<sup>2</sup>),其中λ为拉格朗日乘子;最后,利用二分法得到最优拉格朗日乘子,记为λ<sup>*</sup>;(4)将步骤(3)得到的最优拉格朗日乘子λ<sup>*</sup>代入关系式y<sup>*</sup>=(A<sup>T</sup>Q<sup>‑1</sup>A+λ<sup>*</sup>D)<sup>‑1</sup>(A<sup>T</sup>Q<sup>‑1</sup>b‑λ<sup>*</sup>f)中,得到上述加权最小二乘关系式的全局最优解y<sup>*</sup>,将y<sup>*</sup>代入公式y=[x<sup>T</sup> r<sub>0</sub>]<sup>T</sup>,根据y<sup>*</sup>=[x<sup>*T</sup> r<sub>0</sub><sup>*</sup>]<sup>T</sup>,得到目标源在参考坐标系的坐标x<sup>*</sup>,即得到目标源位置的最终定位值。
地址 315211 浙江省宁波市江北区风华路818号