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

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

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

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

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