[求助]平面波入射到层状球体介质 求球表面阻抗？

烦请高手指点小弟一下   PuYAoKG  
我想计算平面波入射到层状球体介质模型时  球表面的波阻抗Z   WzDL(~m+Z  
如图所示：

图片:未命名.jpg

  

图片:未命名.jpg

=c8xg/  
我想计算在球面的波阻抗 A]c'`Nf  
可是在计算过程中出现了问题， nc4KeEl  
因为我的球体模型假设为地球，在迭代过程中因为要计算贝塞尔函数 bwG$\Oe6  
用到复宗量 "z=kr"，但是k*r  k是波数 r=每层半径 最外层为地球半径 6371km *6IytWOX5  
这样复宗量z值会很大  而复宗量贝塞尔函数却是发散的  z=kr中当r变化不大 但k太大了就计算结果出现错误了 z.P<)[LUc  
< .. /m"O.17N  
14Y_ oH9  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

我参考的文章是： @oKW$\  
多层球体对电磁场散射响应的求解方法 p!/!
ZIo  
----徐建华  聂在平

这是一个研究问题。见 和里面的参考。 &})Zqc3Lqk  
  author = {O. Pe{\~{n}a} and U. Pal}, ;w;+<Rd  
  title = {Scattering of electromagnetic radiation by a multilayered sphere}, $}EI3a  
  journal = {Computer Physics Communications}, V]Kk=  
  year = {2009}, gI/
SA  
  volume = {180}, e7yn"kd  
  pages = {2348 - 2354}, ,f~)CXNT?  
  number = {11}, bE!z[j]  

球Bessel 函数的计算，有现成的计算机程序，1996年金建铭和张善杰出了一本这个方面的专著。上门有这个方面的程序，不过程序是 Fortran. 应该可以算，大综量的计算可能要用渐近表达式。 Jq:Wt+a  
qFp]jbU  
  表面阻抗是有确定意义的。