= -0.96043E+03+J-0.29443E+00
0.41849E+03+J-0.29383E+00 ~ 0.44756E+02+J-0.29209E+00 = 0. 13232E+02+J-0 28921£+00 ~ 0.57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.26686E+00 0. 90103E+00+J-0. 25880E+00 m 0.65908E+00+J-0.24988E+00 ~ 0.49402E+00+J-0.24016E+00 0. 3756JE+00+J-0.22974E+00 0.28728E+00+J-0.21868E+00 0. 21923E+00+J-0. 20708E+00 0. 16554E+00+J-0.19502E+00 0. 12239E+00+J-0.18262E+00 e 0.87235E-01+J-0.16995E+00 0. 58352E-01+J-0. 15711E+00 0.345t0E-01+J-0.14419E+00 = 0.14817E-01+J-0.13130E+00 = -0. 13950E-02+J-0. 11852E+00 = 0.41849E+03+J-0.293BJE+00 -0. 96043E+03+J-0. 29443E+00 0. 41849E+03+J-0.29383E+00 = 0.44756E+02+J-0.29209E+00 0. 13232E+02+J-0.28921E+00 0.57993E+01+J-0.28520E+00 0.31J10E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0 26686E+00 = 0 90103E+00+J-0.25880E+00 0. 65908E+00+J-0.24988E+00 = 0.49402E+00+J-0.24016E+00 ~ 0.37563E+00+J-0.22974E+00 0. 28728E+00+J-e.21868E+00 0. 21923E+00+J-0.2070BE+00 0. 16554E+00+J-0. 19502E+00 c 0.12239E+00+J-0.t8262E+00 0.B7235E-01+J-0 16995E+00 o 58352E-01+J-0. 157' 1E+00 o 34510E-01+J-0.14419E+0e 0. 14817E-01+J-0. 13130E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+0J+J-0. 29383E+00 e -0.96043E+03+J-0.29443E+00 ~ 0.41849E+03+J-0 29383E+00 o 44756E+02+J-0 29209E+00 0. 13232E+02+J-0 28921E+00 0.57993E+01+J-0 28520E+00 0.31310E+01+J-0 2R012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.26686E+00 = 0.90103E+00+J-0.25880E+00 0.65908E+09+J-0 24988E+00 0.49402E+00+J-0 24016[+00 o 37563E+00+J-0.22974E+00 0. 2872BE+00+J-0.21868E+00
0. 21923E+00+J-0. 20708E+00 0. 16554E+00+J-0. 19502E+00 0. 12239E+00+J-0.182S2E+00 0. 87235E-01+J-0. 16g95E+00 0.58352E-01+J-0.15711E+00 0.34510E-01+J-0.14419E+00 0. 13232E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 -0. 96043E+03+J-0.29443E+00 0. 41849E+03+J-0.29383E+00 = 0.44756E+02+J-0.29209E+00 - 0. 13232E+02+J-0 28921E+00 0. 57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 o 12790E+01+J-0.26686E+00 0.90103E+00+J-0.25880E+00 - 0.65908E+00+J-0.24988E+00 0. 49402E+00+J-0.24016E+00 0. 37563E+00+J-0. 22974E+00 0. 28728E+00+J-0.21868E+00 0. 21923E+00+J-0. 20708E+00 • 0 16554E+00+J-0.19502E+00 0. 12239E+00+J-0.18262E+00 0.B7235E-01+J-0 16995E+00 0.58352E-01+J-0.15711E+00 0.57993E+01+J-0.28S20E+00 0. 13232E+02+J-0.28921 E+00 0. 44756E+02+J-0.29209E+00 - 0.41849E+03+J-0.29383E+00 ~ -0.96043E+03+J-0.29443E+00 0.4t849E+03+J-0.29383E+00 0. 44756E+02+J-0.29209E+00 0. 132J2E+02+J-0.28921E+00 o 57993E+01+J-0.28520E+00 0. 31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.28686E+00 0.90103E+00+J-0.25B80E+00 = 0.65908E+00+J-0.24988E+00 0. 49402E+00+J-0.240T6E+00 = 0.37563E+00+J-0.22974E+00 0. 28728E+00+J-0.21868E+00 0. 21923E+00+J-0. 2070BE+00 ~ 0.16554E+00+J-0.19502E+00 0. 12239E+00+J-0. 18262E+00 0.87235E-01+J-0.16995E+00 0.31310E+01+J-0.28012E+00 0. 57993E+01+J-0.28520E+00 o 13232E+02+J-0.28921E+00 o 4475SE+02+J-0.29209E+00 0. 41849E+0J+J-0.29383E+00 -0. 96043E+03+J-0.29443E+00 e.41B49E+03+J-0.29383E+00 0.44756E+02+J-0.29209E+00 - 0. 13232E+02+J-0 28921E+00 0. 57993E+01+J-0. 28520E+00 = 0.31310E+01+J-0 28012E+00 3 0 19197E+01+J-0.27398E+00 - 0.12790E+01+J-0.266B6E+00
0. 90103E+00+J-0.25880E+00 0.65908E+00+J-0.249S8E+00 0.49402E+00+J-0 24016£+00 = 0_l7563E+00+J-0.22974E+00 o 28728E+00+J-0.21868E+00 0.21923E+00+J-0.20708E+00 0. 16554E+00+J-0.19502E+00 e.12239E+00+J-0.18262E+00 0. 19197E+01+J-0.27398E+00 o J1310E+01+J-0.28012E+00 0. 57993E+01+J-0.28520E+00 0. 132J2E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29J83E+00 = -0.96043E+03+J-0.29443E+00 0. 41849E+0J+J-0.293B3E+00 0. 44756E+02+J-0.29209E+00 o 13232E+02+J-0.28921E+00 0. 57993E+01+J-0.28520E+00 0.3t310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 = 0.12790E+01+J-0.26686E+00 0. 90103£+00+J-0.25880E+00 0. 65908E+00+J-0.24988E+00 ~ 0.49402E+00+J-0.2401SE+00 = 0 37563E+00+J-0.22974E+00 0. 28728E+00+J-0.2186SE+00 - 0.21923E+00+J-0.20708E+00 0. 16554E+00+J-0.19502E+00 o 12790E+01+J-0 26686E+00 0. 19197E+01+J-0.27398E+00 - 0.31J10E+01+J-0.28012E+00 0.57993E+01+J-0.28520E+00 0. 13232E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29J83E+00 = -e.9604JE+03+J-0.29443E+00 e.41849E+03+J-0.29383E+00 - 0.4475SE+02+J-0.29209E+00 = 0.13232E+02+J-0.28921E+00 0. 57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.2668SE+00 0.90103E+00+J-0.2S880E+00 0. 65908£+00+J-0.2498BE+00 0. 49402E+00+J-0.24016E+00 0. 37S63E+00+J-0.22974E+00 0.2B728E+00+J-0.21868E+00 ~ 0 21923E+00+J-0.20708E+00 o 90103E+00+J-0.25880E+00 0. 12790E+01+J-0.266B6E+00 0. 19197E+01+J-0.27398E+00 0.3t310£+01+J-0.28012E+00 0. 57993E+01+J-0.28520E+00 - 0.13232E+e2+J-0.28921E+00 - 0.44756E+02+J-0.29209E+00 = 0.41849E+03+J-0.29383E+00 - -e.96043E+03+J-0.29443E+00 0.41849E+03+J-0.29383E+00 = 0.44756E+02+J-0.29209E+00
= =
G
~
-
=
e
~ ~
~
z
-
=
= =
~
=
=
m
0. 13232E+02+J-0.28921E+00 0. 57993E+01+J-0. 28520E+00 o 31310E+01+J-0 28012E+00 0. 19197£+01+J-0.27398E+00 0.12790£+01+J-0.26686E+00 0.90103E+00+J-0.25880E+00 0. 65908E+00+J-0.24988£+00 0.49402£+00+J-0.24016£+00 0.37563E+00+J-0.22974E+00 0. 28728E+00+J-0.21868E+00 0.65908E+00+J-0.24988[+00 0.90103E+00+J-0.2S880E+00 0. 12790E+01+J-0.26686E+00 0.19197E+01+J-0.27398E+00 0.31J10E+01+J-0.28012E+00 0. 57993E+01+J-0.28520E+00 0. 13232E+02+J-0.28921E+00 0.44756E+02+J-0.29209E+00 0.41849E+03+J-0.29383E+00 -0.96043E+03+J-0.29443E+00 0.41849E+03+J-0.29383E+00 0.44756E+02+J-0.29209E+00 0. 13232E+02+J-0.28921E+00 0. 57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0.19197E+01+J-0.27398E+00 0. 12790E+01+J-0.26686E+00 o 90103[+00+J-0.25880E+00 0.65908E+00+J-0 24988E+00 0. 49402E+00+J-0. 24016[+00 0.37563E+00+J-0.22974E+00 0. 49402E+00+J-0.24016E+00 0. 65908E+00+J-0.24988E+00 0.90103E+00+J-0.25880E+00 e.12790E+01+J-0.26686E+00 0.19197[+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 0.57993E+01+J-0.28S20E+00 0. 13232[+02+J-0.28921E+00 0.447S6E+02+J-0.29209E+00 0.41849E+03+J-0.29383E+00 -0.96043E+03+J-0.29443E+00 0.41849E+03+J-0.29383E+00 0.44756E+02+J-0.29209E+00 0. 13232E+02+J-0 28921E+00 0.57993E+01+J-0.2B520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.26686E+00 0.90103E+00+J-0 25880[+00 0.65908E+00+J-0.24988E+00 0.49402E+00+J-0.24016E+00 0. 37563E+00+J-0.22974E+00 0.49402E+00+J-0.24016E+00 0.65908E+00+J-0.24988E+00 0.90103[+00+J-0.25880[+00 0. 12790E+01+J-0.26686E+00 e.19197E+01+J-0.2739BE+00 0.31310E+01+J-0.28012E+00 0. 5799JE+01+J-0.2B520E+00 0. 13232E+02+J-0.28921E+00
= 0.44756E+02+J-0.29209E+00 m
5
2
=
s
-
~
c
0.41849E+03+J-0.29383E+00 -0. 96043E+03+J-0.29443E+00 0. 41849E+03+J-0.29383E+00 0. 44756E+02+J-0.29209E+00 o 13232E+02+J-0.28921E+00 o 5799JE+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398ETe0 e.12790E+01+J-0.2668SE+00 0. 90103E+00+J-0.25880E+00 0.S5908E+00+J-0.24988E+00 0. 28728E+00+J-0.21868E+00 0. 375SJE+00+J-0.22974E+00 0. 49402E+00+J-0.2401SE+00 0.S5908E+00+J-0.24988E+00 0. 90103E+00+J-0.25880E+00 0. 12790E+01+J-0.26686E+00 0. 19197E+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 0.S7993E+01+J-0.28520E+00 0.13232E+02+J-0.28921E+00 0. 4475SE+02+J-0.29209E+00 o 41849E+0J+J-0.29383E+00 -0.96043E+03+J-0.29443E+00 0.41849E+03+J-0.29383E+00 0.44756E+02+J...e.29209E+00 0. 13232E+02+J-0.28921E+00 0. 57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12790E+01+J-0.26686E+00 0. 90103E+00+J-0.25880E+00 0. 21923E+00+J-0.20708E+00 e.28728E+00+J-0.21868E+00 0.37563E+00+J-0.22974E+00 0. 49402E+00+J-0.24016E+00 0. 65908E+00+J-0.24988E+00 0. 90103E+00+J-0.25880E+00 0. 12790E+01+J-e.26686E+00 0. 19197E+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 0.57993E+01+J-0.28520E+00 0. 13232E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 -0.9S043E+03+J-0 29443E+00 0. 41849E+03+J-0.293B3E+00 0. 44756E+02+J-0.29209 E+00 0. 13232E+02+J-0.28921E+00 0.57993E+01+J-0.28520E+00 e 31310E+01+J-0.28012E+00 0. 19197E+01+J-0.2739BE+00 0. 12790E+01+J-0.26686E+00 0. 16554E+00+J-0.19502E+00 0. 21923E+00+J-0.20708E+00 e.28728E+00+J-0.21868E+00 0 3756JE+00+J-0.22974E+00 0 49402E+00+J-0.24016E+00 0. 65908E+00+J-0.24988E+00 0. 90103E+00+J-0. 25880E+00
~
-
-
-
= E
0 12790E+01+J-0.26686E+00 0. 19197E+01+J-0.27J98E+00 o 31310E+01+J-0.28012E+00 o 57993E+01+J-0.28520E+00 0. 13232E+02+J-0.28921 E+00 0. 44756E+02+J-0.29209E+00 0.41849E+03+J-0 29383E+00 -0. 9604JE+0J+J-0. 29443£+00 0. 41849E+03+J-0.29383E+00 0. 44756£+02+J-0. 29209E+00 0. 13232E+02+J-0.28921 E+00 0. 57993E+01+J-0.28520E+00 0.31310E+01+J-0.28012E+00 0. 19197E+01+J-0.27398E+00 0. 12239E+00+J-0 18262E+00 0. 16554E+00+J-0.19502E+00 0.21923E+00+J-0.20708E+00 0. 28728E+00+J-0.21868E+00 0.J7563E+00+J-0.22974E+00 0. 49402E+00+J-0. 240'6E+00 0. 6590BE+00+J-e.24988E+00 0.90103E+00+J-0.2S880£+00 0. 12190E+01+J-0.26686E+00 0. 19197E+01+J-0.27398£+00 0.31310E+01+J-0.28012[+00 0. 57993E+01+J-0.28S20E+00 0. 132J2E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 -0.96043E+03+J-0.29443E+00 0.41849E+03+J-0 29383[+00 0. 44756E+02+J-0.29209E+00 0.1 . 28921E+00 0. 28520E+00 0. 31310E+01+J-0.2ae12E+00 0.87235E-01+J-0.16995E+00 0. 12239E+00+J-0.18262E+00 0. 16554E+00+J-0.19502E+00 0. 21923E+00+J-0.20708E+00 0. 2B128E+00+J-e. 21868£+00 0. 37S63E+00+J-0.22974E+0e 0. 49402E+00+J-0.24016E+00 0. 65908£+00+J-0.24988E+00 0. 90103E+00+J-0. 25880E+00 0. 12790E+01+J-0.266B6E+00 0. 19197E+01+J-0.27398E+00 0. 31310E+01+J-0.28012E+00 0. 57993E+01+J-0.28520E+0e 0. 13232E+02+J-0.28921E+00 0. 4475SE+02+J-0. 29209E+00 0. 41849E+03+J-0.29383E+00 -0.96043E+03+J-0.29443E+00 0.41849E+03+J-0.29JB3E+00 0. 44756E+02+J-0.29209E+00 0. 13232E+02+J-0.28921 E+00 0.57993E+01+J-0.28520E+00 0.5B352E-01+J-0.15111E+0e 0.872J5E-01+J-0.16995E+00 0. 12239E+00+J-0. 18262E+00 0. 16554E+00+J-0. 19502E+00 0. 21923E+00+J-0.20708E+00
o 28728E+00+J-0.21868E+00 0. 37563E+00+J-0. 22974E+00 0. 49402E+00+J-0. 24016E+00 0. 65908E+00+J-0.24988E+00 0. 90103E+00+J-0.25880E+00 0. 12790E+01+J-0.26686E+00 0. 19197E+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 ~ 0.57993E+01+J-0.28520E+00 0. 13232E+02+J-0.28921 E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 - -0.96043E+03+J-0.29443E+00 0. 41849E+03+J-0. 29383E+00 0. 44756E+02+J-0.29209E+00 0. 13232E+02+J-0.28921 E+00 0.34510E-01+J-0.14419E+00 0. 58352E-01+J-0. 15711E+00 e.87235E-01+J-0 16995E+00 0. 12239E+00+J-0.18262E+00 0:16554E+00+J-0.19502E+00 0.21923E+00+J-0.20708E+00 0.28728E+00+J-0 21868E+00 0. 37563E+00+J-0. 22974E+00 0. 49402E+00+J-0.24016E+00 0. 65908E+00+J-0.24988E+00 0.90103E+00+J-0.25880E+00 0. 12790E+01+J-0.26686E+00 ~ 0.19197E+01+J-0.27398E+00 ~ 0.31310E+01+J-0.28012E+00 0.57993E+01+J-0 28520E+00 0. 13232E+02+J-0 28921E+00 = 0.44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 = -0.96043E+03+J-0.29443E+00 £ 0.41849E+03+J-0.29383E+00 0. 44756E+02+J-0.29209E+00 e.14817E-01+J-0.13130E+00 0. 34510E-01+J-0. 14419E+00 0. 58352E-01+J-0. 15711E+00 e.8723SE-01+J-0.1699SE+00 = 0.12239E+00+J-0.18262E+00 0. 16554E+00+J-0.19502E+00 = 0.21923E+00+J-0.20708E+00 e 28728E+00+J-0.21868E+00 0. 37563E+00+J-0. 22974E+00 0. 49402E+00+J-e.24016E+00 0. 65908E+00+J-0.24988E+00 0. 90' 03E+00+J-0. 25B80E+00 = 0. 12790E+01+J-0 26686E+00 o 19197E+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 = 0.57993E+01+J-0.28520E+00 o 13232E+02+J-0.28921E+00 0. 44756E+02+J-0.29209E+00 0. 41849E+03+J-0.29383E+00 ~ -0.96043E+03+J-0.29443E+00 o 41849E+03+J-0.29383E+00 -0. 13950E-02+J-0. 1852E+00 0. 14817E-01+J-0.13130E+00 0. 34510E-01+J-0. 14419E+00
z
0. 58352E-01+J-0. 15711E+00
o 87235E-01+J-0.1699SE+00
o 12239E+00+J-0.18262E+00 o 16554E+00+J-0.19502E+00
0.21923E+00+J-0.20708E+00 0.28728E+00+J-0.21868E+00 0.37563E+00+J-0.22974E+00 0.49402E+00+J-0.24016E+00 0. 6590BE+00+J-0.24988E+0e 0. 90103E+00+J-0. 25880E+00 0. 12790E+01+J-0.26686E+00 0. 19197E+01+J-0.27398E+00 0.31310E+01+J-0.28012E+00 0.57993E+01+J-0.28520E+00 ~ 0.13232E+02+J-0.28921E+00 0.44156E+02+J-0.29209E+00 0. 41849E+03+J-0. 29383E+00 = -0.96043E+03+J-0.29443E+00
=
PROGRNtA TORS
c ••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• c • C C C
C
THIS TWO DIMENSIONAL RADIATION AND SCATTERING CODE IS INTENDED • TO GIVE ITS USERS A DEMONSTRATION OF SOLVING THE RADIATION AND • SCATTERING PROBLEMS BY USING THE INTEGRAL EOUATION AND MOMENT • METHOD. TO MAINTAIN SIMPLICITY, EQUAL SEGMENTATION HAS BEEN .. APPLIED TO ALL THE GEOMETRIES AND PIECEWISE PULSE EXPANSION ANO • POINT MATCHING HAVE BEEN SYMMETRIC MATRIX EQUATION • SOLUTION SUBROTUINES HAVE BEEN INCLUDED. IN THEORY THE STRIP • AND THE CIRCULAR CYLINDER PROBLEMS 00 EXHIBIT IMPEDANCE • MATRICES. WHILE THE RECTANGUlAR AND ELLIPTICAL ONES 00 NOT HAVE • SUCH A PROPERTY. HOWEVER, THE DOMINANT CONTRIBUTIONS OF THE .. IMPEDANCE MATRICES ARE SYMMETRIC AND THE SMALLER THE SEGMENT • LENGTH IS. THE BETTER SYMMETRIC THEY HAVE: THEREfORE .. ASSUMING THE S~[TR]C PROPERTY OF THE MATRICES SERVES OUR • PURPOSE WELL. SOME CASES. SUCH AS THE RADIATING SOURCES WHICH .. ARE VERY CLOSE TO THE RECTANGULAR OR ELLIPTICAL BODIES. MAY NOT • GIVE VERY GOOD SOLUTIONS. HOWEVER, THESE CAN BE IMPROVED BY • SPECIFYING A SMALLER SEGMENT LENGTH LIMIT TO ENHANCE THE SYMMETRIC. PROPERTY OF THE IMPEDANCE MATRICES. AND THUS LEAD TO GAIN BETTER • RESULTS. •
C
•
C
C C
C C C C
C C C
C C C C C
c C
C C C
• THE OUTPUT FILE IS ON DEVICE 6, AND THE FOLLOWING PARAMETERS ARE NEEDED TO RUN THE PROGRAM. READ THE FOLLOWING INSTRUCTIONS CAREFULLY TO GIVE THE CORRECT INPUTS.
• • • •
C M THE MAXIMUM SEGMENT NUMBER IS THE MATRIX S • C • C----------------~--------------------------------------------. C
C C
..
THE REMAINING INPUT PARAMETERS SHOULD BE GIVEN IN THE PROGRAM. THOSE ARE AS FOLLOWS
C C C
DL
C
IGEO
C
C C C
C C
C C C
C C C
C C C C C C
C
c c
IPATT
-
THE lARGEST SEGMENT 0.1 IS THE LIMIT. BUT 05 IS RECOMMENDED. THE SMALLER. THE ~E ACCURATE CHOICE OF SCATTERING BODIES 1 STRIP 2 -- CIRCULAR CYLINDER 3 ELLIPTICAL CYLINDER 4 -- RECTANGULAR CYLINDER CHOICES OF THE ELECTROMAGNETIC PROBLEM 1 -- RADIATION PATTERN PROBLEMS 2 -- MONOSTATfC SCATTERING PROBLEM 3 BISTATIC SCATTERING PROBLEM
o
IPOLR --- THE POLARIZATION INDICATOR FOR THE EXCITATION FOR RADIATIOM PATTERN, IT SPECIFIES THE TYPE OF LINE SOURCES 1 -- ELECTRIC LINE SOURCE 2 -- MAGNETIC LINE SOURCE FOR PLANE WAVE SCATTERING, IT SPECIFIES THE POLARIZATION OF THE INCIDENCE PLANE WAVE 1- TRANSVERSE MAGNETIC n ELO TO Z-AXIS 2- TRANSVERSE ELECTRIC FIELD TO Z-AXIS J- ARBITRARY POLARIZATION OF THE TO THE Z-AXIS IrY THE ANGLE IN
• ..
.
...
. ..•... 1III
..
to
•.. ...
...
..... •.. ...
.. .. ...
1III
III
•
C--------------------------------------------------------------.
c
•
C
THE FOLLOWINGS ARE THE INPUTS OF THE GEOMETRY OF THE SCATTERER • THE ORIGIN Of THE COORDINATES IS ALWAYS REFERRED TO THE CENTER OF • THE GEOMETRY. AND THE DIMENSIONS ARE CHOSEN TO BE IN .•
C C C
C
..
C
THE STRIP IS LYING ALONG X-AXIS: AND THE RECTANGULAR AND ELLIPTICAL. CYLINDERS ARE DEFlNED TO HAVE TWO PRINCIPAL AXES LYING ALONG THE .. THE X-AXIS AND Y-AXIS. RESPECTIVELY. •
C C
.. ..
C
C C C C C
C C C
C
C
C C C
C C
THE WIDTH OF THE STRIP fOR STRIP THE RADIUS OF THE CYLINDER THE CYLINDER THE SEMI-AXIAL lENGTH OF ONE OF THE ICAl AXES FOR ELLIPTICAL CYLINDER ALONG THE X-AXIS OR HALF THE WIDTH OF THE RECTANGULAR CYLINDER ALONG THE X-AXIS DI~ENSION. 8 THE SEMI-AXIAL LENGTH Of THE OTHER ELLIPTICAL AXIS fOR ELLIPTICAL CYLINDER ALONG THE V-AXIS OR HALF THE HEIGHT Of THE RECTANGULAR CYLINDER AL0N~ fH~ Y-AXIS D1MENSION. BPHI THE BISTATIC INCIDENT ANGLE- IN PTHETA THE POLARIZATION or THE ELECTRIC FI RESPECT TO Z-AXIS. FOR EXAMPLE PTHETA=0 IMPLIES A TM POlARIZATION AND PTHETA=90 A TE POLARIZATION. THIS PARAMETER ONL V NEEDED WHEN IPOLR IS CHOSe. TO BE J. XS.YS--- THE LOCATION OF THE RADIATING SOURCE WITH RESPECT TO GEOMETRIC CENTER Of THE SCATTERING BODY (IN
W RA A
.. • .. .. • • .. .. III
.. .. .. ..
C
•
C
..
C•••••••••••••••••••••••••••••••••••••• * ••••••••••••••••••••••••••••••••
C ATTENTION!!!
•
C
..
C C C C
FOR A GREAT SAVING OF CPU TI~E AND MEMORIES. THE IMPEDANCE MATRICES fOR STRIP AND CIRCULAR CYLINDER GEOMETRIES ARE DIMENSIONED DIFFERENTly FROM THOSE OF THE RETANGULAR AND ELLIPTICAL CYLENDERS THEREfORE, BE SURE TO SPECIFY ICED IN PAR~EltR STATEMENT
C
• • • • ..
D
EQUIVALENCE DATA RA 1
C .. CRUCIAL PARAMETERS FOR THE PROBLEM DL=0 05
C ..
IPATT=1 THE WIDTH Of THE STRIP IN WAVELENGTH
(FOR STRIP
W=4.
C .. C.. C..
THE RADIUS or THE CIRCULAR CYLINDER IN RA=2. THE LENGTHS OF THE PRINCIPAL AXES OR THE LENGTHS or THE RETANGULAR 80X IN A=.05 8=1.
C. C..
INCIDENT ANGLE FOR 8lSTATIC RCS (IN BPHI=45 POLARIZATION OF THE ELECTRIC FIELD (iN PTHETA=45.
U ...........' ... L . . .
C..
THE LOCATION OF THE RADIATING SOURCE IN WAVELENGTH(S) x~.e
YSo=0.001 C..
SOME CONSTANT
GAM=1.78105 BTA=6 . 2831853 PI-=J.14159265
ETA=120 .• PI 02R-PI/180. J=CMPLX(e.0.1.e) JF(IPATT.EQ.1) THEN WRITE(6.B6} IF(IPOLR.EQ.l) WRITE(6.88) XS,YS
IF(IPOLR.[Q.2) WRITE(6,90) XS,YS ENOIF" IF"(IPATT . EO.2) THEN
WRITE(6,92) IFCIPOLR . EC.lj WRITE(6,94)
If(IPOLR.EC.2 WRITE(6,96) IF(IPOLR.EO.3 WRITE(6.98) PTHETA ENOIF" IF"(IPATT.EQ.3) THEN WRITE(S,100) BPHI IF!IPOLR.EO"j WRTTE!6.94) IF IPOUR.EQ.2 WRITE 6,96) IF IPOLR.EO.3 WRITE 6,9S) PTHETA ENDIF GOTa (1,2,3,4) IGEO
WRITE(6,110) W CALL
STRIP(ZMN,VA,VB,VT,XS,YS,PTHETA.8PHI.XB,W.M,NMA,ET~,IR,WA)
IF(IR.NE.0) WRITE(6,999)
rF(IPATT . NE.2) THEN WRITE(6,112) DOle 1=1 ,NMA
10 WRITE(6,114) X8(1).CA8S(VA(1» ENDIF GOTO S0
2
20
WRITE(6,120) RA CALL CIRCL(ZMN,VA.V8,VT,XS,YS,PTHETA,BPHI,XB,RA,M,NMA,ET~,IR,WA) IF(IR.NE.0) WRITE(6.999) IF"(IPATT.NE.2) THEN WRITE(6.122) DO 20 I=l.~A WRITE(6,124) XB(I).CABS(VA(I» ENDIf GOTO 80
3
WRITE(6,130) A,B CALL ELlIP(ZMN,ZT,VA,VB ,VT,XS,YS,PTHETA.BPHI,XB.YB.M,NMA.ETMM.IR) IF(IR.NE.0) WRITE(6,999) IF(IPATT NE.2) THEN WRlTE( S, 132)
30
WRITE(6.134) XB(I),YB(I),CABS(VA(I» ENDIF
4
WRI1E(S.l40) A.B
DO 30 l=l.NMA
GOTO 80
CALL RECT(ZMN.ZT,VA,V8,VT.XS.YS,PTHETA,BPHI.XB.YB.M,NMA.ET~.IR) If{IR.NE.0) WRITE(6,999) IF(IPATT.NE.2) THEN WR IT E( 6 , '42) DO 40 1-1,Ntw4A
2.65
I)
I»
80
IT SOLVES fOR THE CURRENT THE NORMALIZED RADIATION PATTERN LINE SOURCE LOCATED AT: • , F7 . :5 • • • • • f7 . 3. •
SOURCE lOCATED
'I
FIELD IS • .FB. r"~""DL_iI.-I'Wl
WI TH THE ANGLE Of' I
'j)
150
160
999
•••••••••••••••••••••••••••••••••••••••••••••••••••••••
+5X ' ••• ERROR HAS BEEN FOUND IN THE HANKEL fUNCTION OR ••• +5X.' ••• IN THE INVERSION OF THE IMPEDANCE MATRIX OR ••• +5X.· ••• THE RADIATING SOURCE BEING SHIELDED. • •• +5X.· ••• OUTPUT MAY BE UNRELIABLE. • •• +5X.· •••••••••••••••••••••••••••••••••••••••••••••••••••••• STOP END
c ,VT,XS,YS,PTHETA,BPHI.X.W.NU.NMA.ETMW.
c
STRIP PROfH.El.4 .BTA.ETA.PI.D2R.R2D.DL.J
C
FUNCTIONAL SUBROUTINE FOR ,IPATT
,8
10 1)
1 ) • VT (1) CRT. J • J
SPECIfY M > '.NMA IR-' RETURN
ENOl' 2
ENDfF C.. COMPUTES THE Z MATRIX
c..
NMA
ELE~TS
THE STRIP HAS A TOPLITZ PROPERTY, ONLY
ARE ,........ v ..,v
1
..
•)
*BTA
'"
C .. rJ
IN
MATRIX
00 6 K.... '.NMA
6
C ..
ELSE XO)(=1.
ENDIF
00 B K=1.NMA
8
10
.25
.25
c ..
ODI'\OI:"I:)"'"
TO SAVE TIME
12
lPOLR.[Q.l) RETURN 100
C.. RADIATION PATTERN 20 COOTINUE 00 22
~1.NMA
22
24
RETURN C.. THE BISTATIC CASE 40 CONTINUE C.. GET THE INCIDENT ANGLE fOR BlSTATIC CASE PHI=BPHI.02R 42
DO 42 M-l. NMA J
ELSE XDX=l.
ENDIF
46
)
0.
50
RETURN THEN
60
100 CONTINUE C .. THEN COMPUTES THE TE CASES C .. FIRST THE MATRIX EL~ENTS 10=3
1
102
C.. FILLING IN
~ATRIX
ELEMENTS IN
DO 104 t(l!:1,NMA
104 GOIO 120
GOTO 140 THE MONOSTATJC CASE
c ..
106
108
ELSE
c ..
112
PROPERTY TO SAVE T1ME
RETURN
C. RADIATION PATTERN 120 CONTINUE C,. THE ELECTRIC FIELD DUE TO MAGNETIC LINE SOURCE DO 122 N=1.NW.
122
1J0
RETURN C.. THE BISTATIC CASE 140 CONTINUE C.. GET THE INCIDENT ANGLE FOR BISTATIC CASE PHI=BPHI.02R 00 142 Mz:1 NMA
142
».51
ELSE XDX=1 ,
ENDIF
146
ELSE
150 160
1)
RETURN END
I
c
.VA.VB.VT.XS,YS.PTHETA,BPHI.X.R.NM.NMA.ETMM.
c
SUBROUTINE FOR CIRCULAR CYLINDER PROBLEM .BTA.ETA.PI.D2R,R2D.Dl,J POLR, JPATT B CS.ID 1) l),CRT.J.HANKA DIMENSION EXTERNAL RK-R.BTA •• p I .. '" LT. GT. SIZE. SPECIFY M >' ,NMA ILl • ...n
2
A TOPlITZ PROPERTY
C .. C ..
.»-OC..,STA-ETA6.25
XL.XU.21. .25.CRT IN THE MATRIX C.. fl DO 6 K-l.NMA 4
6
ELECTRIC LINE SOURCE
C ••
.... 20
24
»
30
RETURN C.. THE BISTATIC CASE 40 CONTINUE C.. GET THE INCIDENT ANGLE FOR BISTATIC CASE PHI=BPHI.D2R PTM::=1 •
I
IPOLR. EO. M:=1.
)))
42
50
RETURN THEN 60 100
10=5
C.. THEN COMPUTES THE TE CASES
C.. fIRST THE UATRIX ELEMENTS 1
1 .• -2
102
104
I C ••
THE
LINE SOURCE
10=2 00 120 N==l.NMA
120
124 130
RETURN C.. THE BISTATIC CASE 140 CONTINUE C.. GET THE INCIDENT ANGLE fOR BISTATIC CASE PHI=BPHI.D2R
142
146
ELSE .25 ENDIF I
150
I»
160
RETURN END
c .ZT.VA.VB.VT.XS,YS.PTHETA.8PHI.X,Y.~.NMA.
C
PROBLEM
1)
l).CRT,J.HANKA
c ..
SIZE, SPECIFY M >',NMA
c ..
XAm=A YB=B
XB=0.
A4=.AoA.A-A 84=8·8.8.8
OS-DE.. 1
DO :3 1"'2.NUO
2
3
4
C •.
C..
Z MATRIX OF THE ELLIPSE IT IS ALSO PART Of SO JUST COMPUTE IT ONCE .
. 25)-1.».OE-BTA9ETA-.25
DO 6
DO 6
~1.t.WA N-l.~
IF( ... . EQ.N) THEN ZMN(N. t.I)-=VA( 1) ElSE
DXEA*A*V(N)*(X(M)-X(N~)-B*B.X(N)'(V(M)-Y(N»
6
DY=A*A*(V(N).Y(t.I)-B*B +B*B*X(N)*X(M) R=BTA*SORT~(Y(t.I)-Y(N) •• 2+(X( ... )-X(N» •• 2) CS-DX/SQRT DX.OX+DY.OY) CALL CSINT HANKA.XL.XU.21.CRT) ZWN(N .... )=0.2S.ETA.CRT ENOIF CONTINUE 00 7 No-<1.~
DO 7 ..... 1.NW.
ZT(N.U)-ZMN(N .... ) IF(IPOLR.EQ .2) GOTO 100 IF(IPOLR.EQ.3) THEN PTt.t=COS(PTHETA*02R) PTE-SIN(PTHETA.02R) ENOIF C .. FACTORIZING THE IMPEDANCE MATRIX CALL CROUT(ZMN.VA.0.NMA,MT) IF(IPATT.EQ.l) GOTD 20 IF(IPATT.EO .3) GOTO 40 C .. THE t.«:>NOSTATlC CASE DO 12 1"'1.91 PHI .. (J-l . ) .D2R 7
DO 8
K~l,
Nt.4A
XK-BTA*(X(K).COS(PHI)+Y(K)*SIN(PHI» CRT=CEXP(J.Xl<) VA(K)~RT
8
10
VB(K);DE.CRT CALL CROUT(ZMN,VA,2,NMA.MT) CRT=Ct.4PLX(0 .• e . ) DO 10 t.A=1. Nt.4A CRT~RT+VB(t.I).VA(M)
Ir(IPOLR . EQ.3) THEN VT(I)~RT.CRT.PTM'PTM
ELSE
YA--CABS(CRT'CRT).BTA.ETA.ETA*.25 IF(YA.LE . l . E-9) VA=1.[-9 YA= 10 .• ALOG10(VA) C .. USING THE SYMMETRIC PROPERTY TO SAVE TIME ETMM( r )=VA ElM.4( 1 82- I )=YA E~( 180+I)=YA ETt.4M(362-I )=YA ENOIF 12
CONTINUE
IF(IPOLR.EO . l) RETURN GOTO 100 c .. RADIATION PATTERN 20 CONTINUE I 0=:1
22
DO 22 N=l,NMA XKD=8TA*SORT«Y(N)-YS) •• 2+(X(N)-XS)**2) VA(N)~0.25.8TA.ETA.HANKA(XKD) CROUT(Z~.VA . 2.NMA . MT)
CALL
00 30 K=l,361
PHI=(K-l.)*02R
CRT -c:a.tP l)C ( e. .e _) DO 24 ~1, Nt.tA
X~BTA.(X(~).COS(PHI)+y(~).SIN(PHI»
CRTaCRT+VA(M).CEXP(J'XM) XK %BTA*(XS.COS(PHI)+YS.SIN(PHI» CRT-OE'CRT+CEXP(J.XK) 30 ETMM(K)-cABS(CRT.CRT) CALL D8(ET~,J61,IR) RETURN C .. THE BISTATIC CASE 40 CONTINUE C .. GET THE INCIDENT ANGLE FOR BISTATIC CASE PHlzBPHI*02R 24
42
DO 42 "'1,NMA VA(M)-CEXP(JeBTA'(X(~)*COS(PHI)+Y(M).SIN(PHI»)
CALL CROUT(ZMN. VA, 2. NMA,MT) DO 50 1-1,361 PHI-( 1-1. ).02R CRT-C~LX(0. ,0.) DO 48 M:l,NMA
XM=BTA*(X(M).COS(PHI)+Y(M)-SIN(PHI»
48
CRT~RT+VA(M).OE'CEXP(J'XM)
IF(IPOLR.EO.3) THEN
VT(I)-CRT*CRT.PT~.PTM
ELSE YA~ABS(CRT.CRT)'BTA.ETA'ETA •. 25
IF(YA.LE.1.E-9) YA=1 . E-9
ET~(1)-10 .• ALOG10(YA)
50
60 100
[NDIF CONTINUE If(IPOLR.EO.1) RETURN IF(IPOLR.EO.3) THEN DO 60 1~1,~ V8(1)=VA(I)"2.PTM.PTM ENOIF 10-2
C .. THEN COMPUTES THE TE CASES, THE I~EDANCE MATRIX VA(1)=CMPLX(1 . . -2./PI.(ALOG(BTA'GAM.OE*.25)-1.».OE,BTA,ETA •. 125 ++J*ETA*(0.25*DE-1./(PI*PI*OE)) DO 104 ...... 1 ,NMA DO 104 N=l.NMA IF(M.EQ.N) THEN Z~N(N.M):::::VA(l)
ELSE DX-A4,Y(M).Y(N)+84'X(M)*X(N) DY=(X(M).Y(N)-X{N).Y(M» CT=OX/SQRT(DX.DX+A4.84.0Y.OY) XDc-e.5 eOE*A.A*Y(N)/SQRT(A4.Y(N)-Y(N)+B4.X(N)_X(N» YD-0.S*OE.B.e.X(N)/SORT(A4eY(N),Y(N)+84.X(N).X(N») XA-X(N)-XO YA-Y(N)-YD DX-e.B*X(M)*(Y(M)-YA)-AeA.Y(M)-(X(M)-XA) OY-A.A.(YAeY(M)-8 e B)+B-a.XA.X(M) CT1-0X/SORT(DXeDX+DY.DY) Rl-BTA'SORT«X(M)-XA) •• 2+(Y(M)-YA) •• 2) XA-zX(N)+XD YA--Y(N)+YD DX-e.BeX(M).(Y(~)-YA)-A'A'Y(M).(X(M)-XA) DY-A.A'(YA.Y(~)-B'B)+B.B*XA.X(M) CT2~DX/SQRT(OX.DX+DY.OY)
R2-BTAeSORT«X(M)-XA)*.2+(Y(M)-YA) •• 2)
2.76
ZUN(N.Y)-CT.ZT(N.~)-0 . 25.ETA.(HANKA(R1).CT'-HANKA(R2).CT2)
ENDlf 104
CONT JNUE
C.. NOW FACTORIZE THE IMPEDANCE MATRIX CALL CROUT(ZUN.VA,e,NMA,UT) If(IPATT.EO.1) GOTO 120 IF(IPATT.EO.3) GOTO 140 C.. THEN THIS MUST BE THE MONOSTATIC CASE DO 112 ' .. 1,91
PHI .. ( 1-1 . ) .02R
00 106 K-1,NMA
XK-BTA.(X(K).COS(PHI)+Y(K).SlN(PHI» OX..BeB.X(K).COS(PHI)+AeAeY(K)-STN(PHI) DY-AeAeY(K).COS(PHl)-SeBeX(K)eSIN(PHI)
CRT-CEXP(J.XK)-OX!SORT(DX.DX+OY.DY) VA(K)-CRT 106
VB(K)~E.CRT
CALL
CROUT(ZMN.VA,2.NMA.~T)
CRT-ct.APLX(e. ,e.) DO 108 M=1 .NMA
108
CRT-CRT+VS(~).VA(M)
If(IPOLR.EO.3) THEN YA=CABS(CRT.CRTePTE-PTE+V1(I»eBTAeETAeETAe ,25 ELSE YA=CABS(CRT.CRT)-STAeETA.ETAe . 25 ENDlf If(YA.LE.l.E-9) YA=1.E-9
YA-10 .• ALOG10(YA)
C .. USING THE SYMMETRIC PROPERTY TO SAVE TI~E ET~11 )=YA ETW 182-} )=YA En..-c 180+1 ) ..YA ETt.t.4(J62-} )-YA 112 CONTINUE RETURN C.. RADIATION PATTERN 120 CONTINUE C .. THE ELECTRIC FIELD DUE TO MAGNETIC LINE SOURCE 00 122 N=1,NMA
DX=A.A.Y(N).(X(N)-XS)+8 eee X(N)e(YS-Y(N» DYcAeAe(Y(N)eYS-S.S)+S.S.X(N)_XS RK-BTA.SQRT«Y(N)-YS) •• 2+(X(N)-XS)e.2) VA(N)=-0.25 e J.STA-HANKA(RK).DY/SORT(DX.DX+DY_DY)
122
CONTI NU[ CALL CROUT(ZMN.VA,2.NMA.MT) 00 lJ0 K=I . .36'1
PHI=(K-l.)e02R CRT~L)(e.
,e.)
DO 124 1=1,NMA
X~BTAe(X(I).COS(PHI)+Y(I).SIN(PHI»
DX=S.S.X(I).COS(PHl)+AeA.V(I).SIN(PHI) DY-A.A.Y(I).COS(PHI)-S.S.X(I).SIN(PHI)
124
130
CRT=CRT+VA(I).CEXP(J.X~)-DX/SORT(DX.DX+DY.DY) X~ ~BTA.(XS.COS(PHI)+YS.SIN(PHI» CRT=ETA.CRT-DE+CEXP(J-X~) ET~(K)=CABS(CRT.CRT) CALL DB(ET~,361,IR)
RETuRN
C.. THE BISTATIC CASE 140 CONTINUE C.. GEl THE INCIDENT ANGLE FOR BISTATIC
C~SE
142
148
ELSE
.25
158
1»
168
RETURN £NO
SUB~IN£ RECT(ZMN.ZT.VA.VB.VT.XS.YS.PTHETA.BPHI.X.Y.~.~. +~.IR)
C F"UNCTIONAL SUBROUTINE FOR RECTANGULAR GEOMETRY
COMMON/CST/GAM.BTA,ETA.PI.D2R.R2D.DL.J
C~N/PATH/R.A.B,CS,ID ~N/TYP/JPOLR.IPATT C~L£X Z~N(~T.~T),ZT(~T.MT),VA(1),VB(1),VT(1).CRT.J,HANKA
DIMENSION ETMM(l),X(l).Y(l) EXTERNAL HANKA Nl e 2 .• A/DL+0.5 N2-2 .• B/OL-t-e. 5
If(Nl.LT.l) THEN
Nl .. 1 N2=B/A+.5 ENDIF
If(N2.LT.1) THEN N2-1 N1 a A/B+e.5 ENDfF
NMA-2.(Nl-+N2)
IF(NMA.GT.MT) THEN PRINT *,' MATRIX EXCEEDS PRESET SIZE, SPECIfY M >',NMA IF~·'
RETURN ENDtF
DE=2 .• (A+B)/(Nl+N2)
C.. ASSIGN THE X,Y COORDINATES or THE RECTANGULAR BOX 00 1 l=l,N1 yO )-B .
X(!)S::A-(I-.5je OE Y(NMA-N2+1-1 --Y(I) X(NMA-N2+1-I ~X(I) DO 2 1-1,N2
X(N1+I )-A
Y(N1+1)-S-(I-0.5)eDE X(NMA+l-I)c-X(Nl+I)
Y(NMA+l-I)-Y(Nl+I)
2
NN-W.X0(Nl,N2) NH--Nl+N2 VA(l)-CMPLX(l. ,-2./Pt.(ALOG(BTA.GAMeOE •. 25)-1.».DE*8TA*ETA*.25 ID=l
DO 3 1-2,NN XL=(1-1.5)*OE*STA XU=(I-0.5).OEeBT~
3
CALL CSINT(HANKA,XL,XU,21 ,CRT) VA(I)=0.2S e ETA.CRT 10=6
XL=-e.SeDE.BTA X\J=0.S.0E.STA 00 6 M=l.NMA
DO 6 N=l,NMA IF{M.LE.Nl.AND.N.LE.Nl) THEN ZMN(N,M)=VA(IABS(M-N)+l) GOTO 6
ENOIf
]F{~.CT . Nl.AND.M.LE.NH.AND.N.CT.Nl.ANO.N.LE.NH)
ZMN(N.M)~VA(IABS(M-N)+l)
THEN
GOTO 6 ENDIF If(~.CT.NH.AND.~.lE.Nl+NH.AND.N
Z~(N,M)~VA(IA8S{M-N)+l)
279
GT NH.AND.N.LE.Nl+NH) THEN
OOTO 6
ELSE
D«>IF 6
7
c ..
ZEf«J-OKIJt..H
HANKEL FUNCT ION
100
c ..
c .. )) 8
10
C. .
12
TO SAVE TIME
ENOl CONTINUE
I
IPOlR.EO.1) RETURN 100
C .. RADIATION PATTERN 20
22
CONTINUE 00 22 N:-1.
»)
24 30
RETURN C.. THE BISTATIC CASE 40 CONTINUE:
c .. GET THE INCIDENT ANGLE rOR BISTATIC CASE PHI=BPHI·02R
00 42 M:==1 .NMA
42
»)
48
S0
RETURN THEN
60 100
CONTINUE
C.. THEN COMPUTES THE TE 1
IMPEDANCE MATRIX
.»-DEeBTA.ETAe.125
1. ,-2
.OE.STA .Dr.BTA XL, XU, 21 • 102
DO
)
ENDIF'
281
END)
ENDIF GOTO 104 ENDlf THEN
]
THEN
I
)
ENOlf
ELSE
ENOl ENDff GOTO 104 F
Gl.NH.ANO.M.LE. GT .NH.AND.N.LL
1
ENDlf
THEN
THEN
ELSE
c ..
THE IMPEDANCE MATRIX .VA.0.NMA
GOTO 120 GOTO 140
THE MONOSTATIC CASE
C •.
»
106
108 ELSE
c ..
PROPERTY TO SAVE TIME
CONT RETURN C.. RADIATION PATTERN 120 CONTINUE C .. THE ELECTRIC fIELD DUE TO MAGNETIC LINE SOURCE 112
DO 1
N;;:; 1 • NMA
fIELD AS THE EXITATION .PE CALL DO 130 K=l. . ..02R ,0. )
) 124 130
) 142
148
I».BTA .. ETA .. ETAe.25 ELSE
.STA .. ETA.ETA •. 25
150
I) )
160
RETURN END
c ..
SUBROUTINE
C .• A ROUTINE TO
NUl"UlIIIII"IL
OUTPUT AND TAKE ITS DB
REAL
2
C.. THE WAX IMUM SHOULD BE AROUND 1 TO 2 I
. L T. 0. 1) THEN
••• CHECk THE LOCATION or THE SOURCE
10 1R=1
RETURN (HDlr ..
00 .. 1=1.N I
END
•••• )
C .. COMPLEX FUNCTION HANKEL FUNCTIONS OF
C..
.A.B.CS.ID
GOTO C .• ZERO 1 CALL
HANKA-H RETURN C.. 1ST ORDER HANKEL FUNCTION 2 CAll HANK(XA.H.H1) HANKA-H1 RETURN C . TE STRIP HANKEL fUNCTION J
CALL
.H.Hl)
RETURN C.. TM CIRCULAR HANKEL FUNCTION 4- XA=2. CAll HANKA=H
C
RETURN TE CIRCULAR HANKEL ruNCTION 5 XA=2. CAll HANKA=2.
RETURN
C
TM 6
COMBINATIONS FOR DIFFERENT GEOMETRY
c
c
ARGUMENT ZERO ORDER HANKEL FUNCTION OF SECOND KIND FIRST ORDER HANKEL FUNCTION OF SECOND KIND
c c
8£0.
81-0.
v-e.
Yl-0.
THEN
THEN Xl0==X5.X2 X12=Xle.X2 B=.21E-3eX12-.39444E-2*Xl0+.444479E-1.X8 Y--.24846E-3.X12+.427916E-2eXl0-.4261214E-1.XS 81=. 1109£-4.X12-.31761[-3-X10+.443319E-2.X8 Y1=.27873E-2-X12-.400976E-1·X10+.3123951.XB
ENDIF BsB-.3163866.X6+1.2656208.X4 Y=Y+.25300117.X6-.74350384-X4 81-81-.3954289£-1.)(6+. 21093573.X4 Y1=Y1-1.3164827.X6+2 1682709.)(4 ENDIF 8=8-2.2499997.X2+1. Y=Y+. 69559366.X2+. 36746691+XlN_8 56249985.)(2+. 2212091.X2-. tlSE X1=3.
X2=X1
)(3-)(2.)(1 X4=X3.Xl X5=X4.Xl X6=)(5.)(1 F=.79788456-.77E-S.X1-.55274E-2.X2-.9512E-4.X3+.137237E-2.X4 ~ -.72805E-3.X5+.14476[-3·)(6 T=X-.78539S1S-.4166397E-1.Xl-.3954E-4.X2+.262573[-2·X3 & -.54125£-3.X4-.29333E-3.X5+.13558E-3.X6
& t
Fm. 15GE-56Xl+.1GS9667E-leX2+.17105E-3eX3-.249511E-2eX4 +.11365JE-2.X5-.20033E-3.X6 T=X-2.3561945+.12499612eX1+.5S5E-4.X2-.637879E-2.XJ+ 74348E-3-X4 .29166E-3.XG
ENDIF
RETURN END
c c
C FAST ALGOR I T...... FORM C BY THE AUTHOR
SUBROUTINE CSI
or
THE S I UPSON • S I NT [ORA L ROUT I NE
XL,XU.N.
IMPLICIT
THEN )
ELSE CRT-=CRT+2. ENDIF
28 CONTINUE RETURN [ND
)
c .. SUBROUTINE C .. STRAIGHT FORWARD C .. COMPLEX MATRIX EQUATION
C..
C .. C.. C.. C.
A.,)( - B
JOB =
e
IN SOLVING THE
FACTORING THE MATRIX 1 fACTORING THE MATRIX AND SOLVE THE EQUATION - 2 SOLVE THE EOUATION BASED ON fACTORED MATRIX ~
10
12 20 22
25 30 I==N-L+l 1)",,1+1 I
35
I •
38
END
c c ••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• c
c
c c c c c c c
SUBROUTINE TSLZ
NETLIB
TOEPllTZ PACKAGE. THIS VERSION DATED INPUT: The first row of the T-Matrix fol lowed Its first column Inning with the element. On return is unaltered. The ri hand side vector B. A wor area vector Order of matrix A.
- 1)
c
c
c c
c
c The solution vector. c c PURPOSE: of equations described by a TOEPlITZ motri~. c Solve a c A • X c c SUBROUTINES AND fUNCTIONS: c TOEPlJTZ PACKAGE .. , TSlZl c c ••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• SUBROUTINE
.R
INTEGER
c
COMPLEX
CAll TSlZ1 RETURN END
INE TSlZ1
C SUPPORTING ROUTINE SUBROUTINE TSLZ1 INTEGER t.4 COMPLEX A1, INTEGER 11. .N.
TSLZT .A2.B.X.C1.C2
).
.C1
.N2
COMPLEX R1,R2.R3,RS.R6 R1 -
A1
1
1) 1
GO TO 20
1
11
11:%
1, N2
12 == N - Ii
R5 - R5 + 11).C1 I RS - R6 + Al 11+1) 10 20
CONTINUE
CONTINUE
R2
=
R3 "'" R1 :=: R1
IF
+ RS.RJ
.EO. 2 GO TO 40 = C2(1 .000.0.
11 -
• N1 11
11 11 30 40
= (I1)$R3 + R6 :::; C1(I1) + RS.R2
CONTINUE C2 1
= R3
= (0.000,0.
DO 50 11 :::; 1. Nl 12 "" N 11 R5 R5 + 11)*X(I
50
=
CONTINUE R6 == DO
11) eR6 60 70 80 CONTINUE RETURN END
1
1)
)
)