主权项 |
一种基于集群处理机的SAS频域信号处理的方法,该方法用并行技术、多线程技术和FFTW算法,使得ωk算法能够移植到SAS系统中,能在集群处理机上进行ωk算法,所述的方法包括以下步骤:1)设定SAR系统的参数表;2)距离向脉冲压缩;3)方位向非均匀离散快速傅立叶变换NSFFT;4)距离向和方位向相位补偿及STOLT变换;5)二维傅立叶逆变换;其中,步骤(3)中所述的NSFFT变换方法能够用于多子阵的合成孔径声纳非均匀成像的情况,该变换方法具体步骤如下:设声纳基阵个数为Nc,等效相位中心为间距为dpc,前进速度为v,脉冲重复周期为prt,帧数为P方位向数据点数为M=Nc·P,则有方位向傅立叶变换的离散表示式为: <mrow> <msub> <mi>Ee</mi> <mi>bc</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>u</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>假设等效相位中心做匀速直线运动,则上式可分解为 <mrow> <msub> <mi>Ee</mi> <mi>bc</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>u</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>P</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </msub> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>ω</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>其中 ωp×Nc+m=m×ωpc+p×ωprt则 <mrow> <msub> <mi>Ee</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>P</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </msub> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>ω</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>P</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </msub> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>mω</mi> <mi>pc</mi> </msub> <mo>+</mo> <msub> <mi>pω</mi> <mi>prt</mi> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>mω</mi> <mi>pc</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <munderover> <mi>Σ</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>P</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </msub> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <mi>p</mi> <msub> <mi>ω</mi> <mi>prt</mi> </msub> <mo>)</mo> </mrow> </mrow>令 <mrow> <msub> <mi>Ee</mi> <msub> <mi>g</mi> <mi>m</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>pω</mi> <mi>prt</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mrow> <mi>p</mi> <mo>×</mo> <mi>Nc</mi> <mo>+</mo> <mi>m</mi> </mrow> </msub> </msub> </mrow>则 <mrow> <msub> <mi>Ee</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <mi>m</mi> <msub> <mi>ω</mi> <mi>pc</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <munderover> <mi>Σ</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>P</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>g</mi> <mi>m</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>pω</mi> <mi>prt</mi> </msub> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <msub> <mi>pω</mi> <mi>prt</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>πk</mi> </mrow> <mi>N</mi> </mfrac> <mi>m</mi> <msub> <mi>ω</mi> <mi>pc</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msub> <mi>Ee</mi> <mi>G</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>则 <mrow> <msub> <mi>Ee</mi> <mi>bc</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>Ee</mi> <msub> <mi>b</mi> <mi>m</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <msub> <mi>k</mi> <mi>u</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow> |