发明名称 一种雷达系统误差估计方法
摘要 一种雷达系统误差估计方法,其步骤如下:⑴网格划分;⑵节点数据集划分;⑶接收目标数据;⑷判定收到的数据种类;⑸存储收到的雷达量测数据;⑹判定高精度信源数据是否存在对应的雷达目标批号;⑺时间配准;⑻计算新的雷达系统误差及其方差;<img file="2011100714008100004DEST_PATH_IMAGE002.GIF" wi="16" he="24" />查找新数据对隶属的节点;<img file="2011100714008100004DEST_PATH_IMAGE004.GIF" wi="16" he="24" />新数据对隶属的节点集合是否为空;<img file="DEST_PATH_IMAGE006.GIF" wi="16" he="24" />计算节点处雷达系统误差;<img file="DEST_PATH_IMAGE008.GIF" wi="16" he="24" />计算网格系统误差函数;<img file="DEST_PATH_IMAGE010.GIF" wi="16" he="24" />计算待估计区域系统误差。其优点在于:由数据驱动的系统误差估计方法克服了雷达系统误差模型复杂,难于建立的困难,采用数值方法得到系统误差的区域分布,提高了估计精度;采用网格节点迭代计算与网格双线性函数插值,计算量小、工程实施方便;本发明可以用于岸基系统的雷达系统误差校准。
申请公布号 CN102221688B 申请公布日期 2013.04.24
申请号 CN201110071400.8 申请日期 2011.03.24
申请人 中国船舶重工集团公司第七○九研究所 发明人 刘颢;肖厚;陈世友;许治华
分类号 G01S7/40(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 武汉金堂专利事务所 42212 代理人 胡清堂
主权项 1.一种雷达系统误差估计方法,其特征在于:其步骤如下:<b>第1步</b>网格划分包括:划定一块扇形区域<img file="2011100714008100001DEST_PATH_IMAGE002.GIF" wi="18" he="18" />覆盖待估计系统误差的区域;<img file="274380DEST_PATH_IMAGE002.GIF" wi="18" he="18" />的距离范围是<img file="2011100714008100001DEST_PATH_IMAGE004.GIF" wi="53" he="25" />,<img file="2011100714008100001DEST_PATH_IMAGE006.GIF" wi="76" he="25" />,其中距离指到雷达的距离,<img file="2011100714008100001DEST_PATH_IMAGE008.GIF" wi="20" he="25" />为最小距离、<img file="2011100714008100001DEST_PATH_IMAGE010.GIF" wi="22" he="25" />为最大距离,<img file="939323DEST_PATH_IMAGE002.GIF" wi="18" he="18" />的方位角范围是<img file="2011100714008100001DEST_PATH_IMAGE012.GIF" wi="54" he="25" />,其中<img file="2011100714008100001DEST_PATH_IMAGE014.GIF" wi="116" he="25" />,表示与正北所成的角度;取<img file="2011100714008100001DEST_PATH_IMAGE016.GIF" wi="202" he="25" />,将距离范围划分为<img file="2011100714008100001DEST_PATH_IMAGE018.GIF" wi="18" he="16" />份,取<img file="DEST_PATH_IMAGE020.GIF" wi="206" he="26" />,将方位角范围划分为<img file="DEST_PATH_IMAGE022.GIF" wi="14" he="16" />份;网格的顶点称为节点,记到雷达距离为<img file="DEST_PATH_IMAGE024.GIF" wi="18" he="25" />,与正北所成方位为<img file="DEST_PATH_IMAGE026.GIF" wi="21" he="26" />的节点为<img file="DEST_PATH_IMAGE028.GIF" wi="25" he="26" />,存储所有节点,网格以其左上角的节点标号作为自己的标号,左上角为<img file="DEST_PATH_IMAGE030.GIF" wi="24" he="26" />的网格记为<img file="DEST_PATH_IMAGE032.GIF" wi="22" he="26" />,存储所有的网格;<b>第2步</b>节点数据集划分包括:设定每个节点的支撑数据集划分方法;数据对是指时间配准后的高精度信源数据与雷达量测数据组成的对,通过划分方法判定数据对是否对<img file="223455DEST_PATH_IMAGE028.GIF" wi="25" he="26" />的系统误差计算有贡献;划分方法为:数据对中的两个数据均满足与<img file="DEST_PATH_IMAGE034.GIF" wi="17" he="25" />的距离差不超过<img file="DEST_PATH_IMAGE036.GIF" wi="28" he="25" />,与<img file="244763DEST_PATH_IMAGE026.GIF" wi="21" he="26" />的方位角差不超过<img file="DEST_PATH_IMAGE038.GIF" wi="30" he="26" />,则属于节点<img file="183769DEST_PATH_IMAGE028.GIF" wi="25" he="26" />的支撑数据集,其中<img file="DEST_PATH_IMAGE040.GIF" wi="54" he="25" />,<img file="DEST_PATH_IMAGE042.GIF" wi="57" he="26" />;<b>第3步</b>接收目标数据包括:接收雷达量测数据、高精度信源数据;其中,雷达量测数据指待估计系统误差的雷达探测到的目标数据,高精度信源数据指目标自身携带的高精度定位设备;<b>第4步</b>判定收到的数据种类:收到的如果是雷达量测数据则执行第5步,如果是高精度信源数据则执行第6步;<b>第5步</b> 存储收到的雷达量测数据;<b>第6步</b>判定高精度信源数据是否存在对应的雷达目标批号:查找目标批号对照表,如果找到第3步收到的高精度信源数据对应的雷达目标批号,则记录此目标批号,执行第7步,否则返回第3步;<b>第7步</b>时间配准包括:数据的时刻是指该数据的产生时刻,记第6步中被判定的高精度信源数据时刻为<img file="DEST_PATH_IMAGE044.GIF" wi="14" he="20" />,找到目标批号<img file="DEST_PATH_IMAGE046.GIF" wi="14" he="20" />索引的雷达量测数据中,时刻与<img file="816351DEST_PATH_IMAGE044.GIF" wi="14" he="20" />相差最小的雷达量测数据,记为<img file="DEST_PATH_IMAGE048.GIF" wi="20" he="25" />,其中<img file="DEST_PATH_IMAGE050.GIF" wi="10" he="17" />为数据的时刻;将第6步中被判定的高精度信源数据对准到<img file="647166DEST_PATH_IMAGE050.GIF" wi="10" he="17" />时刻,利用其速度进行外推或二点线性插值求得;将对准后的高精度信源数据转换到雷达量测坐标系,与<img file="846066DEST_PATH_IMAGE048.GIF" wi="20" he="25" />组成数据对;存储并记雷达目标批号<img file="701895DEST_PATH_IMAGE046.GIF" wi="14" he="20" />在<img file="632942DEST_PATH_IMAGE050.GIF" wi="10" he="17" />时刻的数据对为<img file="DEST_PATH_IMAGE052.GIF" wi="57" he="25" />,高精度信源数据<img file="DEST_PATH_IMAGE054.GIF" wi="103" he="29" />,方差为<img file="DEST_PATH_IMAGE056.GIF" wi="178" he="46" />,雷达量测数据<img file="DEST_PATH_IMAGE058.GIF" wi="97" he="28" />,方差为<img file="DEST_PATH_IMAGE060.GIF" wi="185" he="47" />;其中<img file="DEST_PATH_IMAGE062.GIF" wi="46" he="46" />、<img file="DEST_PATH_IMAGE064.GIF" wi="45" he="45" />取雷达标称精度(RMS),<img file="DEST_PATH_IMAGE066.GIF" wi="44" he="44" />、<img file="DEST_PATH_IMAGE068.GIF" wi="47" he="47" />由高精度信源标称精度(RMS)转至雷达量测坐标系求得;记对准后的高精度信源数据在坐标系<img file="DEST_PATH_IMAGE070.GIF" wi="14" he="16" />中的矢量表示为<img file="DEST_PATH_IMAGE072.GIF" wi="24" he="25" />,方差矩阵为<img file="DEST_PATH_IMAGE074.GIF" wi="20" he="25" />,在雷达量测坐标系<img file="DEST_PATH_IMAGE076.GIF" wi="14" he="20" />中的矢量表示为<img file="DEST_PATH_IMAGE078.GIF" wi="24" he="25" />,方差矩阵为<img file="DEST_PATH_IMAGE080.GIF" wi="20" he="25" />,其中<img file="DEST_PATH_IMAGE082.GIF" wi="20" he="25" />为高精度信源标称精度(RMS)的平方作为对角线元素形成的对角阵;记坐标系<img file="651188DEST_PATH_IMAGE070.GIF" wi="14" he="16" />到坐标系<img file="76615DEST_PATH_IMAGE076.GIF" wi="14" he="20" />的矢量转换函数为<img file="DEST_PATH_IMAGE084.GIF" wi="18" he="18" />,即<img file="DEST_PATH_IMAGE086.GIF" wi="84" he="25" />,那么<img file="DEST_PATH_IMAGE088.GIF" wi="161" he="26" />,其中<img file="DEST_PATH_IMAGE090.GIF" wi="52" he="25" />为函数<img file="662317DEST_PATH_IMAGE084.GIF" wi="18" he="18" />在<img file="783507DEST_PATH_IMAGE072.GIF" wi="24" he="25" />处的雅可比矩阵,<img file="DEST_PATH_IMAGE092.GIF" wi="69" he="26" />为它的转置矩阵,<img file="DEST_PATH_IMAGE094.GIF" wi="47" he="49" />、<img file="DEST_PATH_IMAGE096.GIF" wi="47" he="47" />取<img file="87449DEST_PATH_IMAGE080.GIF" wi="20" he="25" />中对应的对角线元素开方;<b>第8步</b>计算新的雷达系统误差及其方差包括:由第7步得到的数据对计算得:<img file="DEST_PATH_IMAGE098.GIF" wi="170" he="51" />,<img file="DEST_PATH_IMAGE100.GIF" wi="262" he="72" /><img file="DEST_PATH_IMAGE102.GIF" wi="177" he="51" />,<img file="DEST_PATH_IMAGE104.GIF" wi="271" he="72" />其中<img file="DEST_PATH_IMAGE106.GIF" wi="42" he="48" />为新雷达距离系统误差,<img file="DEST_PATH_IMAGE108.GIF" wi="68" he="62" />为<img file="DEST_PATH_IMAGE110.GIF" wi="36" he="40" />的方差,<img file="DEST_PATH_IMAGE112.GIF" wi="42" he="48" />为新雷达方位系统误差,<img file="DEST_PATH_IMAGE114.GIF" wi="64" he="58" />为<img file="DEST_PATH_IMAGE116.GIF" wi="37" he="42" />的方差;<b>第9步</b>查找新数据对隶属的节点包括:遍历第1步存储的节点,按第2步设定的支撑数据集划分方法,判定第7步得到的数据对是否属于其支撑数据集,如果属于则记录此节点,最后存储记录的节点集合;<b>第10步</b>判定新数据对隶属的节点集合是否为空:如果第9步得到的节点集合不为空,则执行第11步,否则返回第3步;<b>第11步</b>计算节点处雷达系统误差包括:对第9步得到的节点集合中每一个节点计算系统误差:假设节点<img file="DEST_PATH_IMAGE118.GIF" wi="51" he="42" />处的距离系统误差为<img file="DEST_PATH_IMAGE120.GIF" wi="53" he="46" />,方差为<img file="DEST_PATH_IMAGE122.GIF" wi="57" he="48" />,距离系统误差为<img file="DEST_PATH_IMAGE124.GIF" wi="55" he="48" />,方差为<img file="DEST_PATH_IMAGE126.GIF" wi="61" he="50" />,如果此节点未计算过系统误差则直接赋值<img file="DEST_PATH_IMAGE128.GIF" wi="138" he="60" />,<img file="DEST_PATH_IMAGE130.GIF" wi="148" he="57" />,<img file="DEST_PATH_IMAGE132.GIF" wi="127" he="55" />,<img file="DEST_PATH_IMAGE134.GIF" wi="148" he="58" />;否则按如下公式迭代计算一次:<img file="DEST_PATH_IMAGE136.GIF" wi="471" he="109" /><img file="DEST_PATH_IMAGE138.GIF" wi="301" he="104" /><img file="DEST_PATH_IMAGE140.GIF" wi="429" he="98" /><img file="DEST_PATH_IMAGE142.GIF" wi="301" he="103" />;其中<img file="DEST_PATH_IMAGE144.GIF" wi="53" he="49" />为距离系统误差新的迭代值,<img file="DEST_PATH_IMAGE146.GIF" wi="58" he="51" />为<img file="DEST_PATH_IMAGE148.GIF" wi="63" he="55" />的方差,<img file="DEST_PATH_IMAGE150.GIF" wi="66" he="58" />为方位系统误差新的迭代值,<img file="DEST_PATH_IMAGE152.GIF" wi="66" he="54" />为<img file="DEST_PATH_IMAGE154.GIF" wi="59" he="52" />的方差;存储迭代计算后的新误差及方差值:<img file="DEST_PATH_IMAGE156.GIF" wi="174" he="64" /><img file="DEST_PATH_IMAGE158.GIF" wi="174" he="58" /><img file="DEST_PATH_IMAGE160.GIF" wi="174" he="62" /><img file="DEST_PATH_IMAGE162.GIF" wi="174" he="59" />;<b>第12步</b>计算网格系统误差函数包括:首先遍历第9步找到的节点集合,找到所有以此集合中节点为顶点的网格集合,然后遍历此网格集合,记其中的网格<img file="DEST_PATH_IMAGE164.GIF" wi="28" he="33" />内的系统误差函数为<img file="DEST_PATH_IMAGE166.GIF" wi="76" he="37" />,<img file="DEST_PATH_IMAGE168.GIF" wi="81" he="38" />,其中<img file="DEST_PATH_IMAGE170.GIF" wi="16" he="20" />为网格<img file="DEST_PATH_IMAGE172.GIF" wi="27" he="32" />内一点到雷达的距离,<img file="DEST_PATH_IMAGE174.GIF" wi="17" he="22" />为此点方位,将网格按照估计出系统误差的顶点数分类,当顶点数为0时,此网格系统误差未知,标记系统误差函数不存在;当仅有一个<img file="DEST_PATH_IMAGE176.GIF" wi="70" he="41" />时,此网格的系统误差<img file="DEST_PATH_IMAGE178.GIF" wi="219" he="99" />;当有两个顶点<img file="DEST_PATH_IMAGE180.GIF" wi="65" he="39" />,<img file="DEST_PATH_IMAGE182.GIF" wi="65" he="37" />时,<img file="DEST_PATH_IMAGE184.GIF" wi="289" he="145" />;当有三个顶点<img file="DEST_PATH_IMAGE186.GIF" wi="63" he="37" />,<img file="DEST_PATH_IMAGE188.GIF" wi="65" he="36" />,<img file="DEST_PATH_IMAGE190.GIF" wi="66" he="38" />时,<img file="DEST_PATH_IMAGE192.GIF" wi="381" he="147" />;当四个顶点均有估计值时,采用拉格朗日双线性插值多项式对四个顶点插值得到网格系统误差函数<img file="DEST_PATH_IMAGE194.GIF" wi="158" he="61" />,其中<img file="DEST_PATH_IMAGE196.GIF" wi="61" he="26" />表示网格<img file="DEST_PATH_IMAGE198.GIF" wi="23" he="26" />内的距离系统误差插值函数,<img file="DEST_PATH_IMAGE200.GIF" wi="62" he="26" />表示网格<img file="730615DEST_PATH_IMAGE198.GIF" wi="23" he="26" />内的方位系统误差插值函数,将后四种情况得到的函数存储于网格<img file="561036DEST_PATH_IMAGE198.GIF" wi="23" he="26" />;<b>第13步</b>计算待估计区域系统误差包括:<b>第13步的第一步</b>输入区域坐标输入区域中待估计点的坐标,记为<img file="DEST_PATH_IMAGE202.GIF" wi="44" he="22" />,其中<img file="522301DEST_PATH_IMAGE170.GIF" wi="16" he="20" />为此点到雷达的距离,<img file="251223DEST_PATH_IMAGE174.GIF" wi="17" he="22" />为此点方位,即雷达到此点的矢量与正北所成角度;<b>第13步的第二步</b>判定系统误差函数是否存在查找此点所属的网格,如果此网格的系统误差函数存在,执行<b>第13步的第三步</b>,反之执行<b>第13步的第四步</b>,所述<b>第13步的第三步</b>是输出系统误差值,包括将第一步输入的<img file="267720DEST_PATH_IMAGE202.GIF" wi="44" he="22" />代入第二步找到的网格的系统误差函数,计算并输出系统误差值;所述<b>第13步的第四步</b>是输出系统误差未知。
地址 430074 湖北省武汉市洪山区珞喻路718号