发明名称 TOPSAR系统参数工程设计方法
摘要 本发明提供一种TOPSAR系统参数工程设计方法。技术方案是:已知卫星参数、雷达参数、波位参数,在满足整个测绘带的图像连续而均匀的前提下,根据设计指标和待设计系统参数之间的定量关系,设计出适于工程设计的系统参数。TOPSAR待设计系统参数包括方位处理角、方位分辨率、扫描速率、各子测绘带的工作时间和扫描周期。本发明可使设计结果中的各子测绘带有一致的方位分辨率,并且各子测绘带的中心方位模糊度,方位模糊度恶化,NESZ恶化和方位向重叠率均满足指标要求。本发明设计的TOPSAR系统参数,可使TOPSAR工作中整个测绘带的图像连续而均匀。
申请公布号 CN102967861B 申请公布日期 2014.02.12
申请号 CN201210395168.8 申请日期 2012.10.17
申请人 中国人民解放军国防科学技术大学 发明人 黄海风;陈祺;董臻;张永胜;何峰;余安喜;杜湘瑜;孙造宇;金光虎;何志华;李力;马喜乐
分类号 G01S13/90(2006.01)I;G01S7/02(2006.01)I;G06F17/50(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 郭敏
主权项 1.一种方位向电扫描合成孔径雷达系统参数工程设计方法,其特征在于,包括下述步骤:已知卫星参数和雷达参数:卫星速度V<sub>s</sub>,星载合成孔径雷达信号波长λ,波束切换时间T<sub>G</sub>,天线扫描角为<img file="FDA00002269098200011.GIF" wi="28" he="33" />时对应的双程天线方向图函数<img file="FDA00002269098200012.GIF" wi="145" he="60" />天线波束3dB宽度θ<sub>03dB</sub>;已知波位参数:脉冲重复频率<img file="FDA00002269098200013.GIF" wi="91" he="61" />波束脚印速度<img file="FDA00002269098200014.GIF" wi="91" he="61" />中心斜距<img file="FDA00002269098200015.GIF" wi="91" he="55" />上标i表示第i个子测绘带,i=1~N<sub>s</sub>,N<sub>s</sub>为子测绘带数目;系统设计要求:各子测绘带有一致的方位分辨率,各子测绘带的中心方位模糊度小于等于AASR<sub>0</sub>,各子测绘带的扫描角度小于等于<img file="FDA00002269098200016.GIF" wi="104" he="36" />各子测绘带的方位模糊度恶化小于等于AAST<sub>det0</sub>,各子测绘带的等效噪声系数NESZ恶化小于等于NESZ<sub>det0</sub>,方位向重叠率不小于e;第一步:方位处理角设计对第i个子测绘带,进行下述处理:令<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>&theta;</mi><msub><mn>0</mn><mi>j</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>2</mn><mi>arcsin</mi><mrow><mo>(</mo><mfrac><msubsup><mi>&lambda;f</mi><mi>p</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><msub><mrow><mn>4</mn><mi>V</mi></mrow><mi>s</mi></msub></mfrac><mo>)</mo></mrow><mfrac><mi>j</mi><msub><mi>N</mi><mi>j</mi></msub></mfrac><mo>,</mo></mrow></math>]]></maths>j=1,2,3,…,N<sub>j</sub>其中,N<sub>j</sub>为一正整数,与需要的方位处理角精度有关,N<sub>j</sub>越大,方位处理角计算精度越高;对每一个<img file="FDA00002269098200018.GIF" wi="78" he="62" />用下式计算多普勒带宽<img file="FDA00002269098200019.GIF" wi="91" he="67" /><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>B</mi><msub><mi>a</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><msub><mrow><mn>4</mn><mi>V</mi></mrow><mi>s</mi></msub><mi>&lambda;</mi></mfrac><mi>sin</mi><mrow><mo>(</mo><mfrac><msubsup><mi>&theta;</mi><msub><mn>0</mn><mi>j</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>再用下式计算对应的中心方位模糊度<img file="FDA000022690982000111.GIF" wi="177" he="67" /><img file="FDA000022690982000112.GIF" wi="739" he="366" />其中,<img file="FDA00002269098200021.GIF" wi="301" he="62" />表示对<img file="FDA00002269098200022.GIF" wi="301" he="62" />取平方,<img file="FDA00002269098200023.GIF" wi="143" he="54" />为双程天线方向图<img file="FDA00002269098200024.GIF" wi="120" he="61" />当<img file="FDA00002269098200025.GIF" wi="102" he="43" />时的函数;按下式计算方位处理角θ<sub>0amb</sub>;<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>j</mi><mrow><mn>0</mn><mi>amb</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><munder><mi>max</mi><mrow><msubsup><mi>AAST</mi><mrow><mi>c</mi><mo>,</mo><mi>j</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>&le;</mo><msub><mi>AASR</mi><mn>0</mn></msub></mrow></munder><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&theta;</mi><mrow><mn>0</mn><mi>amb</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msubsup><mi>&theta;</mi><msub><mn>0</mn><msubsup><mi>j</mi><mrow><mn>0</mn><mi>amb</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mrow></math>]]></maths>则,最终方位处理角<img file="FDA00002269098200028.GIF" wi="60" he="56" />和多普勒带宽<img file="FDA00002269098200029.GIF" wi="64" he="56" />为:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>&theta;</mi><mn>0</mn><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mi>min</mi><mo>{</mo><msubsup><mi>&theta;</mi><mrow><mn>0</mn><mi>amb</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msub><mi>&theta;</mi><mrow><mn>03</mn><mi>dB</mi></mrow></msub><mo>}</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>B</mi><mi>a</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><msub><mrow><mn>4</mn><mi>V</mi></mrow><mi>s</mi></msub><mi>&lambda;</mi></mfrac><mi>sin</mi><mrow><mo>(</mo><mfrac><msubsup><mi>&theta;</mi><mn>0</mn><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>第二步,方位分辨率设计令<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>&rho;</mi><msub><mi>a</mi><mi>k</mi></msub><mo>&prime;</mo></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>s</mi></msub></munderover><mfrac><mi>&lambda;</mi><msubsup><mrow><mn>2</mn><mi>&theta;</mi></mrow><mn>0</mn><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup></mfrac><mfrac><msubsup><mi>V</mi><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><msub><mi>V</mi><mi>s</mi></msub></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mi>k</mi><msub><mi>N</mi><mi>k</mi></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>k=1,2,3,…,N<sub>k</sub>其中,N<sub>k</sub>为正整数,与需要的方位分辨率精度有关,N<sub>k</sub>越大,方位分辨率精度越高;对每个<img file="FDA000022690982000213.GIF" wi="82" he="57" />进行下述处理:首先,计算:<img file="FDA000022690982000214.GIF" wi="583" he="149" />i=1,2,3,…,N<sub>s</sub>再解线性方程组:<img file="FDA000022690982000215.GIF" wi="881" he="223" />i=1,2,3,…,N<sub>s</sub>可得到各子测绘带工作时间<img file="FDA000022690982000216.GIF" wi="61" he="63" />和扫描周期T<sub>k</sub>;然后,用下式计算需要的最大扫描角度<img file="FDA000022690982000217.GIF" wi="49" he="57" />及其对应的中心频率<img file="FDA000022690982000218.GIF" wi="91" he="63" /><img file="FDA000022690982000219.GIF" wi="392" he="151" /><img file="FDA000022690982000220.GIF" wi="457" he="151" />再根据下式计算各子测绘带的方位模糊度恶化<img file="FDA00002269098200031.GIF" wi="194" he="63" /><img file="FDA00002269098200032.GIF" wi="1023" he="326" />并用下式计算各子测绘带的NESZ恶化<img file="FDA00002269098200033.GIF" wi="199" he="62" /><img file="FDA00002269098200034.GIF" wi="962" he="144" />最后,按下式计算方位分辨率ρ<sub>a</sub>:<img file="FDA00002269098200035.GIF" wi="426" he="181" /><maths num="0008"><![CDATA[<math><mrow><msub><mi>&rho;</mi><mi>a</mi></msub><mo>=</mo><msubsup><mi>&rho;</mi><msub><mi>a</mi><msub><mi>k</mi><msub><mi>&rho;</mi><mi>a</mi></msub></msub></msub><mo>&prime;</mo></msubsup></mrow></math>]]></maths>第三步:时序设计首先根据下式求出扫描速率<img file="FDA00002269098200037.GIF" wi="83" he="62" /><img file="FDA00002269098200038.GIF" wi="554" he="149" />然后解下列方程组:<img file="FDA00002269098200039.GIF" wi="870" he="213" />i=1,2,3,…,N<sub>s</sub>即可求得各子测绘带工作时间<img file="FDA000022690982000310.GIF" wi="62" he="55" />和扫描周期T。
地址 410073 湖南省长沙市开福区德雅路109号