基于矩量法的二維金屬體散射（內(nèi)含matlab程序）

1、基于矩量法的二維金屬體散射計(jì)算基于矩量法的二維金屬體散射計(jì)算1問(wèn)題的描述問(wèn)題的描述本題是用矩量法計(jì)算二維金屬圓柱體的散射場(chǎng)，如圖所示為一圓柱體和一個(gè)橢圓柱的截面，為了計(jì)算簡(jiǎn)單，選入射波為垂直z軸入射的TM或TE平面波yxy2izEizE?x2?2矩量法求解過(guò)程矩量法求解過(guò)程2.1電場(chǎng)積分方程電場(chǎng)積分方程2.1.1問(wèn)題的分析問(wèn)題的分析由麥克斯韋方程組(1)HjE??????(2)JEjH??????可得電場(chǎng)積分方程為(3)20)()(4)

2、(dsKHJKZExz?????????表示在圓柱表面的面電流在遠(yuǎn)處產(chǎn)生的總場(chǎng)。設(shè)入射場(chǎng)為E，散射場(chǎng)為E，由金屬izsz表面的邊界條件=0(4)szizzEEE??得(5)20)()(4)(dlKHJKZECziz???????2.1.2離散化離散化設(shè)入射波為將散射體截面C分為N份△C，用點(diǎn)匹配法對(duì)上述積分)sincos(??yxjkizeE???n式子進(jìn)行離散化，即基函數(shù)可取(6)?上在其它Cnf??10)(?可得下列離散方程：[P]

3、J=b(7)y????22)()()()(tytytxtxR????ytdlnxx(18))(sin)(cos)(tytxtt????其中t表示邊界上的一點(diǎn)，是和X的夾角。)(t?dl根根據(jù)前面的過(guò)程，圓柱邊界分成N分。等效電流密度可以近似為一些脈沖函數(shù)的迭加：(19)???NnnntttpjJ1)()(其中(20)?????其它上在如果01)(nnCttp則得到(21)??????HjZincznnmn?矩陣非對(duì)角元(22))2(1)

4、())()(cos)()((sin4dtkRtyttxtjkmmmCmmmnHRyRxZn????????(23)????22)()(tyytxxRmm????在上認(rèn)為是常量故nC?(24))()cos(sin)2(1mmmnmnmnmnkRynxHRyRxZ??????對(duì)角元(25)21??Zmn其中????22nmnmmnyyxxR????回波寬度的近似公式為：(26)21)sincos()sin(4)(???????Nnyxjkn