用户名: 密码: 验证码:
时域有限元的二维导体微波成像及其图形处理器(GPU)加速的探索
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
本文以标量和矢量时域有限元(TDFEM)分别作为电磁散射正问题的求解基础,以全局优化算法——遗传算法和实数微分进化策略来处理电磁逆散射问题(即微波成像问题),在充分应用从时域获得的超宽带信息的基础上实现标量、矢量时域有限元与遗传算法和实数微分进化策略的有机结合,完成对理想导电体(PEC)二维目标几何形状的重构。此外,本文还对使用图形处理器(GPU)对时域有限元运算加速进行了探讨。
     论文首先讨论散射体的描述方法。重点论述文中提出的离散点描述法,并对其优势予以阐述,并最终应用到课题中去。然后,我们将时域MEI(不变量测试方程)方法引入到时域有限元,并以实例说明时域MEI方法应用到时域有限元中是可行的,同时讨论在成像中使用时会产生的问题。
     接下来重点阐述以标量、矢量时域有限元方法(TDFEM)作为电磁散射正问题的求解基础,获得目标散射场的近场和远场信息。并通过对散射信息的频域分析得到目标的大致尺寸,进而为下一步优化提供必要的约束条件,以确定优化算法的搜索范围,使之能够更高效的寻找到全局最优解。在做具体优化之前,将首先对当今较为流行的几种优化算法进行介绍和比较分析,在此基础上,选择优化算法中寻优效果较为稳定的遗传算法以及寻优效率较高的实数微分进化策略作为逆过程的求解方法,并将它们分别和标量、矢量时域有限元结合来进行成像。接下来,通过我们对近场、远场解空间的特性分析,可以看到近场和远场解空间有着很大的差别,远场解空间的情况要远比近场解空间复杂,而这样的情况就必然导致了远场成像过程的复杂性,也必然导致优化算法需要更长的时间来寻找最优解。由文中算例可以看到,我们能够得到较为令人满意的成像效果,只要给逆过程以足够的寻优时间。但是成像过程所耗费的时间过多,这也是不容回避的问题。
     再者,由于时域有限元繁重的运算量,论文对时域有限元运算的GPU加速进行了一个前瞻性的探索。由于GPU具有一定的并行性、高密集的运算能力、高计算精度,通过对GPU这些优良性能的充分利用,使TD-FEM的运算速度得到数倍的提升。
     最后,将对于本课题做一个全面的总结和展望。归纳出本课题的发展前景以及今后需要进一步解决的问题。
In this project,the forward scattering problem is treated via.scalar Time-Domain Finite Element Method(TDFEM) and vector Time-Domain Finite Element Method,while inversion microwave imaging is implemented by global optimization algorithm,genetic algorithm and real differential evolution strategy. And completed reconstruction of the shape information of unknown scatterer is presented by the perfect union of scalar TDFEM,vector TDFEM and genetic algorithm(GA),real differential evolution strategy based on obtaining the ultra-band information through time-domain.Furthermore,we investigate Graphics Processor Unit(GPU) acceleration of TDFEM algorithm.And significant reduction of the computation time of TDFEM is obtained by GPU.
     Firstly,we introduce methods that characterize the shape of scatterer.Our focus is on discrete points method put forward by author and its merits.And then the method is used in this project.After this,we introduce time domain measured equation of invariance(MEI) into TDFEM,and bear out the analysis result by some examples. Finally,discuss some questions as we achieve microwave imaging through it.
     Secondly,focal point is obtaining the scattering information of near-field and far-field by scalar TDFEM and vector TDFEM.We can get the restrictions of optimization algorithm by extrapolating the approximate size of scatterer from the far-field scattering information,and upon that we can obtain the definite search coverage and find the globe optimum solution faster.We should first simply introduce and analysis the popular optimization algorithm in present-day.And then, GA and real differential evolution strategy are choosed to combine with scalar TDFEM and vector TDFEM respectively achieve imaging.And we will analysis the near-field and far-field next.There has big difference between them.As a result,the circumstances of far-firld solution space are more complicated than near-field solution space.For this reason,the far-field imaging is more difficult,and need more time to find optimum solution.In this paper,we will find the shape reconstructed agrees well with the shape of scatterer,so long as time is enough.But the question that can not skirt round is too much time is spent on imaging process.
     In addition,since the complicated computing of the time-domain finite element, GPU accelerated method is applied in the time-domain finite element calculation for prospective study.As the parallelism property of GPU,high-intensive computing capability and high accuracy,through full utilizing the excellent performance of GPU,the TD-FEM computational speed can multiplied.
     At last,we sum up and look forward to vistas of the project.And put forward the question should be solved in the future.
引文
[1]Joachim H GEnder,I Walterscheid,Andreas R Brenner.New aspects of bistatie SAR:processing and experiments[A].Proceedings of Geoscience and Remote Sensing Symposium[C].Alaska:IEEE,Sept 2004.1758-1762.
    [2]汤子跃,张守融.双站合成孔径雷达系统原理[M].北京:科学出版社,2003.
    [3]I G Cumming,J R Bennett.Digital processing of Seasat SAR data[A].Proceedings of IEEE ICASSP 1979[C].Washinton DC:IEEE,April 1979.710-718.
    [4]TRINTINAL IA L C,H Ling.Joint Time-Frequency ISAR using Adaptive Processing[J].IEEE Trans.Antenna Propagat.1997,2(45):221-227.
    [5]Wu H,et al.Transnational Motion Compensation in ISAR Image Processing [G].IEEE Trans.Image Processing,1995,14(11):1561-1571.
    [6]张直中.逆合成孔径雷达(ISAR)成像[J].中国电子科学研究院学报.2006,10(5):391-404.
    [7]莫锦军,刘少斌,袁乃昌.非均匀等离子体覆盖目标隐身研究[J].电波科学学报,2002,17(1):69-73.
    [8]莫锦军,刘少斌,袁乃昌.等离子体覆盖导体柱宽带散射特性分析[J].微波学报,2003,19(1):20-24.
    [9]石建.导体与涂层目标电磁散射特性研究[D].西安:西安电子科技大学硕士论文,2001.
    [10]李民权,李超.波纹表面金属目标散射特性分析[II].合肥学院学报.2007,17(2):41-43.
    [11]WONG D,CARIN L.Analysis and Processing of Ultra-wide-band SAR Imagery for Buried Landmine Detection[J].IEEE Trans.on AP,1998,46(10):1747-1748.
    [12]何建国.冲激雷达反隐身机理及关键技术的理论与实验研究[D].长沙:国防科学技术大学,1994.
    [13]谢洁.基于冲激响应的隐身目标检测和识别方法研究[D].南京:南京理工大学,2003.
    [14]HUANG Chun-lin,SU Yi.A New GPR Calibration Method for High Accuracy Thickness and Permittivity Measurement of Multi-layered Pavement[C]//2005IEEE International Geoscience and Remote Sensing Symposium.2005:1729-1733.
    [15]Devaney A J,"A filtered backpropagation algorithm for diffraction tomography",Ultrosonic Imaging,No.4,1982:pp336-350.
    [16]Bube K P and Burridge R,"The one-dimensional problem of reflection seismology",SIAM Review,Vol.25,1983:pp497-559.
    [17]Lesselier D,"determination of index profiles by time domain reflectometry",J.Optics,Vol.9,1978:pp349-358.
    [18]Devaney A J,"A computer simulation study of diffraction tomography",IEEE Trans.Biomedical Engineering,Vol.30,No.7,1983:pp377-386.
    [19]Born W B,etal."Introduction to the special issue on inverse method in electromagnetics",IEEE Trans.Antennas and Propagat,Vol.29,No.2,1981:pp185-189.
    [20]彭仲秋.二维电磁逆散射变形波恩迭代法的改进.应用科学学报,1992,11(4):297-302.
    [21]Habashy T M and Oristaglio M L,"Simultaneous nonlinear reconstruction of 2-dimensional permittivity and conductivity",Radio Sci.,Vol.29,1994:pp1101-1118.
    [22]Boag A,"A fast multilevel domain decomposition algorithm for radar imaging,"IEEE Trans.Antennas and Propagat,Vol.49,2001:pp666-671.
    [23]Fortuny G J and Lopez S M,"'Extension of the 3-D range migration algorithm to cylindrical and spherical scanning geometries",IEEE Trans.Antennas and Propagat,Vol.49,2001:pp1434-1444.
    [24]Franchois A and C.Pichot,"Microwave imaging-complex permittivity reconstruction with a Levenberg-Marquardt method" IEEE.Trans.,Antennas Propag.,vol.45,1997:203-215.
    [25]程崇虎,张业荣.圆柱面分层介质结构的微波成像[J].电波科学学报.2003,18(1):38-42.
    [26]Weile D S,Michielssen E.."Genetic algorithms optimization applied to electromagnetics:a review".IEEE.Trans.,Antennas Propag.,vol.45,1997:pp343-353.
    [27]Michielssen E,Ranjithan S."Optimal multiplayer filter design using real-coded genetic algorithms".IEE Proc.J,Optoelectron,1992:pp.413-420.
    [28]Chiu C C,Liu P Y."Image reconstruction of a perfectly conduction cylinder by the genetic algorithm".IEE.Proc.,Microw.Antennas Propag,1996:pp249-259.
    [29]丁振宇,钱祖平,洪伟,“快速多极子和遗传算法在电磁成像中的应用”,电波科学学报,Vol.15,No.3,2000:pp347-251.
    [30]Caorsi S and Pastorino M,"Two-dimensional microwave imaging approach based on a genetic algorithm",IEEE Trans.Antennas and Propag.,Vol.48,No.3,2000:pp 370-373.
    [31]Qing A,C.K Lee and L.Jen,"Electromagnetic inverse scattering of two-dimensional perfectly conducting objects by real-coded genetic algorithm",IEEE trans,on Geosci.and remote sensing,Vol.39,No.3,2001:pp.665-676.
    [32]卿安永,李敬,任朗,李政光,钟顺时.多导体柱微波成像的连续编码遗传算法.地球物理学报.1999,42(6):841-848.
    [33]K.A.Michalski,"Electromagnetic imaging of circular-cylindrical conductorsand tunnels using a differential evolution algorithm," Microwave Opt.Tech.Lett.,vol.27,no.5,pp.330-334,Dec.2000.
    [34]A.Qing,"Elecrromagnetic inverse scattering of multiple two-dimensional perfectly conducting objccts by the differential evolution strategy",IEEE Tran.Antennas Propagat.,vol.51,no.6,pp.1251-1262.
    [35]Lee J F.Whitney elements time domain(WETD) methods[J].IEEE Trans Magnetics,1995,31(3):1325-1329.
    [36]Pernet S,Ferrieres X,Cohen G.High spatial order finite element method to solve Maxwell's equations in time domain[J].IEEE Trans on Antenna and Propagate,2005,53(9):2889-2899.
    [37]Jiao Dan,Jin Jianming.Three-dimensional orthogonal vector basis functions for time-domain finite element solution of vector wave equations[J].IEEE Trans on Antenna and Propagate,2003,51(1):59-66.
    [38]Rickard P L E,Jin Jianming.A two-dimensional time-domain finite element formulation for periodic structures[J].IEEE Trans on Antenna and Propagate,2005,53(4):1480-1488.
    [39]金亚秋.轴对称有耗介质电磁问题的时域有限元方法[J].地球物理学报, 1999,42(2):268-276.
    [40]Shao Zhenhai,Hong Wei.Generalized Z-domain absorbing boundary conditions for time-domain finite-element method[C].Antennas and Propagation society International Symposium.Fukuoka,Japan:IEEE,2000,2:765-768.
    [41]李思敏 时域有限元法与模式匹配法的混合法[D].成都:电子科技大学,2003.
    [42]何小祥,刘梅林.SSOR预处理技术在二维电磁特性TDFEM分忻中的应用[J].南京航空航天大学学报,2006,38(6):670-673.
    [43]Krger,J.,Westermann,R.:Linear algebra operators for gpu implementation of numerical algorithms.ACM Transactions on Graphics,22(2003) 908-916
    [44]Bolz,J.,Farmer,I.,Grinspun,E.,Schroder,P.,Schrder,P.:Sparse matrix solvers on the gpu:Conjugate gradients and multigrid.ACM Transactions on Graphics,22(2003) 917-924
    [45]G.Lolan,W.Cliff,etal,"A multi-grid solver for boundary value problems using programmable graphics hardware",Proceedings of Graphics Hardware,Sandiego,2003:pp102-111.
    [46]B.Tie and D.Aubry,"Reliable and Efficient Numerical Modeling of Elastic Wave Propagation in Engineering Structures by Discontinuous Galerkin FEM',The 6~(th) World Congress on Computational Mechanics in Conjunction with the 2nd Asian-Pacific Congress on Computational Mechanics,2004.
    [47]Harris,M.J.:Real-Time Cloud Simulation and Rendering.Dissertaion,University of North Carolina at Chapel Hill,2003.
    [48]Y.Liu,X.Liu,E.Wu,"Real-time 3D fluid simulation on GPU with complex obstacles",Proceedings of 12th Pacific Conference on Computer Graphics and Applications,pp247-256,Oct.2004.
    [49]F.Zhe,Q.Feng,A.Kaufman,etal,"GPU Cluster for High Performance Computing",Proceedings of the ACM/IEEE Supercomputing Conference,pp47-47,2004
    [50]M.J.Inman,A.Elsherbeni and C.Smith."FDTD calculations using graphical processing units",IEEE/ACES International Conference on Wireless Communications and Applied Computational Electromagnetics,pp728-731,April 2005.
    [51]Dominik G"oddeke,"Accelerating Double Precision FEM Simulations with GPUs",Proceedings of ASIM 2005-18th Symposium on Simulation Technique,Sept.2005.
    [52]S.Krakiwsky,L.Turner and M.Okoniewski,"Graphics processor unit(GPU)acceleration of finite-difference time-domain(FDTD) algorithm",Proceedings of the 2004 International Symposium on Circuits and Systems,Vol.5,pp265-268,May 2004.
    [53]http://developer.nvidia.com
    [54]周永祖,“非均匀介质中的场与波”中文版:361-408.
    [55]卫涛,导体目标的超宽带微波成像,西南交通大学博士论文,2008。
    [56]钟选明,廖成,杨丹,杨周炳,孟凡宝.基于实数遗传算法的二维导体目标微波成像[J].西南交通大学学报.2002,12,37(6):689-692.
    [57]Tao Wei,Cheng Liao,YunLin Liu and Wei Chert."Monostatic Microwave local imaging of PEC cylinder based on genetic algorithm".Asia-Pacifie Microwave Conference Proceeding 2005,Vol.4,pp 1606-1607.Dec.2005.
    [58]C.E.Reuter,E.T.Thiele,and A.Taflove,"Validation and extension tO three-dimensions of the Berenger PML absorbing boundary condition for FDTD meshes," IEEE Microwave Guided Wave Lett.Vol.4,pp.268-270,Aug.1994.
    [59]J.Y.Wu,D.M.Kingsland,J.F.Lee,and R.Lee,"A comparison of anisotropie PML to Berenger's PML and its application to the finite-element method for EM scattering," IEEE Trans,Antennas Propagat.,Vol.45,No.1,pp.40-50,Jan 1997.
    [60]D.Jiao,M.Lu,E.Michielssen,and J.M.Jin,"A fast time-domain finite element-boundary integral method for electromagnetic analysis," IEEE Trans.Antennas Propagat.Vol.49,Oct.2001.
    [61]D.Jiao,A.Ergin,B.Shanker,E.Michielssen,and J.M.Jin,"A fast time-domain higher-order finite element-boundary intergral method for three-dimensional electromagnetic scattering analysis," IEEE Trans.Antennas Propagat.,Vol.50,July 2002.
    [62]廖成,卫涛,陈伟.整数微分进化策略及其在微波成像中的应用[J].西南交通大学学报.2007,12,42(6):647-652.
    [63]刘昆,廖成,卫涛,张双文.基于时域有限元的近场单站导体柱微波成像[J]电波科学学报.(己录用)
    [64]K.K.Mei,R Pous,Z Chert,et al."The measured equation of invariance-A new concept in field computation".IEEE Trans.Antennas Propagat,1994,42(3):320-3281.
    [65]J.O.Jevtic and R.Lee."A theoretical and numerical analysis of the measured equation of invariance".IEEE Trans.Antennas Propagat,1994,42(8):1097-1105.
    [66]廖成.时域有限差分法中的截断边界条件研究,西南交通大学博士后研究工作报告,1997。
    [67]K.K.Mei and Y.W.Liu,"On time domain measured equation of invariance",URSI meeting Seattle,Washington State,June 1994.
    [68]Y.W.Liu,K.K.Mei,K.N.Yang,"Measured equation of invariance in time domain",IEEE Asia pacific Microwave conference,1997.Volume 1,Page(s):477-479.
    [69]Reza Pazoki and Ahmad Cheldavi,"A New Technique to Find the MEI Coefficients in the Time Domain:'Self-Metron Technique'",Progress In Electromagnetics Research,PIER 56,53-66,2006.
    [70]廖成,刘昆.运用自测试子法对一阶MEI吸收边界条件的扩展.西南交通大学学报.(已录用)
    [71]Cheng Liao,Liu Meng,Dan Yang and Xu-Ming Zhong,"Analysis of the numerical error for time domain MEI absorbing boundary condition",Proceedings of APMC,Taipei,Taiwan,2001.Page(s):839-842 vol.2.
    [72]K.liu,Cheng Liao,"Application of the Self-Metron Method to first order MEI-ABC",International Symposium on Antennas,Propagation and EM Theory Proceedings,pp.803-806,Guilin,China,2006.
    [73]K.K.Mei,R.Pous,Z.Chen,Y.W.Liu,and M.D.Prouty,Measured equation of invariance:A new concept in field computation[J].IEEE Trans.Antennas Propagat.,1994,42(3):320-327.
    [74]廖成.时域MEI方法初探[J].电波科学学报,2000,15(3):323-327.
    [75]Yaowu Liu,Kang Lan,Cheng Liao and K.K.Mei,Time-domain MEI method for radiation of line source[J].Electronics Letters,1999,35(4):284-285.
    [76]廖成,杨丹.基于MEI吸收边界条件的导体柱散射问题模拟[J].微波学报,2002,18(4):68-70.
    [77]Shuangwen Zhang,Yunlin Liu and Rugui Yang,A new absorbing boundary condition for FEMTD[J].IEEE ISAPE,Guilin China 2006,926-928.
    [78]J.M.Jin.The Finite Element Method in Electromagnetics(the 2nd Edition),New York:John Wiley and Sons,Inc.,2001.
    [79]J.F.Lee,R.Lee,and A.C.Cangellaris,"Time-domain finite element methods," IEEE Trans.Antennas Propagat.,Vol.45,pp.430-442,March 1997.
    [80]J.Berenger,"A perfectly matched layer for the absorption of electromagnetic waves," J.Comput.Phys.,Vol.144,No.2,pp.185-200,Oct.1994.
    [81]C.E.Reuter,R.M.Joseph,E.T.Thiele,D.S.Katz,and A.Taflove,"Ultrawide-band absorbing boundary condition for termination of waveguiding structures in FDTD simulations," IEEE Microwave Guided Wave Lett.,Vol.4,pp.344-346,Oct.1994.
    [82]W.C.Chew and W.Weedon,"A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates," Microwave Op.Tech.Lett.,Vol.7,No.13,pp.599-604,Spet.1994.
    [83]W.C.Chew and J.M.Jin,"Perfectly matched layers in the discretized space An analysis and optimization," Electromagn,Vol.16,pp.325-340,July-Aug 1996.
    [84]C.M.Rappaport,"Interpreting and improving the PML absorbing boundary condition using anisotropic lossy mapping of space," IEEE Trans.Magn.,Vol.32,No.3,pp.968-974,May 1996.
    [85]W.C.Chew,J.M.Jin,and E.Michielssen,"Complex coordinate stretching as a generalized absorbing boundary condition," Microwave Opt.Tech.Lett.,Vol.15,No.6,pp.363-369,Aug.1997.
    [86]Z.S.Sacks,D.M.Kingsland,R.Lee,and J.F.Lee,"A perfectly matched anisotropic absorber for use as an absorbing boundary condition," IEEE Trans.Antennas Propagat.,Vol.43,No.12,pp.1460-1463,Dec.1995.
    [87]S.D.Gedney,"An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices," IEEE Trans.Antennas Propagat.,Vol.44,No.12,pp.1630-1639,Dec.1996.
    [88]M.Kuzuoglu and R.Mittra,"Investigation of nonplanar perfectly matched absorbers for finite-element mesh truncation," IEEE Trans.Antennas Propagat.,Vol.45,No.3,pp.474-486,March.1997.
    [89]L.Zhao and A.C.Cangellaris,"GT-PML:Generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids," IEEE Trans.Microwave Theory Tech.,Vol.44,No.12,pp.2555-2563,Dec.1996.
    [90]T.V.Yioultsis,T.D.Tsiboukis,and E.E.Kriezis,"A generalized nondiagonally anisotropic perfectly matched layer for wide-angle absorption in finite element electromagnetic scattering analysis," IEEE Trans.Magn.,Vol.34,pp.2732-2735,1998.
    [91]D.Jiao and J.M.Jin,"An effective algorithm for implementing perfectly matched layers in time-domain finite element simulation of open-region EM problems," IEEE Trans.Antennas Propagat.,Vol.50,pp.1615-1632,Nov.2002.
    [92]D.Jiao,J.M.Jin,E.Michielssen,and D.Riley,"Time domain finite element simulation of three-dimensional scattering and radiation problems using perfectly matched layers," IEEE Trans.Antennas Propagat.,Vol.51,pp.296-305,Feb.2003.
    [93]Holland J H.Adaptation in Nature and Artificial Systems.MIT Press 1992.
    [94]De Jong K A.An Analysis of the Behavior of a Class of Genetic Adaptive Systems.Ph.D Dissertation,University of Michigan,No.76-9381,1975.
    [95]Goldberg D E.Genetic Algorithms in Search,Optimization and Machine Learning.Adison-Wesley,1989.
    [96]周明,孙树栋.遗传算法原理及应用.国防工业出版社.1999.
    [97]卫涛.优化算法报告.西南交通大学电磁所.2007.
    [98]H.Whitney,Geometric Integration Theory.Princeton,NJ:Princeton Universiry Press,1957.
    [99]J.C.Nedelec,Mixed finite elements in R3,Numer.Meth.,vol.35,pp.315-341,1980.
    [100]A.Bossavit and J.C.Verite,A mixed FFM-BIEM method to solve 3-D eddy current problems,IEEE Trans.Magnetics,vol.MAG-18,pp.431-435,Mar.1982.
    [101]M.Hano,Finite-element analysis of dielectric-loaded waveguides,IEEE Trans.Microwave Theory Tech.,vol.MTT-32,pp.1275-1279,Oct.1984.
    [102]G.Mur and A.T.de Hoop,A finite-element method for computing three-dimensional electromagnetic fields in inhomogeneous media,IEEE Trans.Magnetics,vol.MAG-21,pp.2188-2191,Nov.1985.
    [103]J.S.van Welij,Calculation of eddy current in terms of H on hexhedra,IEEI Trans.Magnetics,vol.MAG-21,pp.2239-2241,Nov.1985.
    [104]A.Kameari,Three dimensional eddy current calculation using edge elernent for magnetic vector potential,in Applied Electromagnetics in Materials.K Miya,Ed.Oxford:Pergamon Press,pp.225-236,1986.
    [105]M.L.Barton and Z.J.Cendes,New vector finite elements for three-dimensional magnetic field computation,J.Appl.Phys.,vol.61,no.8,pp.3919-3921,Apr.1987.
    [106]C.W.Crowley Mixed order covariant projection finite elements for vector fields,Ph.D dissertation,McGill University,Montreal,1988.
    [107]A.Bossavit,A rationale for 'edge-elements' in 3-D fields computations,IEEE Trans.Magnetics,vol.MAG-24,pp.74-79,Jan.1988.
    [108]A.Kameari,Calculation of transient 3D eddy current using edge-elements,IEEE Trans.Magnetics,vol.MAG-26,pp.466-469,MAR.1990.
    [109]A.Chatterjee,J.M.Jin,and J.L.Volakis,Computation of cavity resonances using edge-based finite elements,IEEE Trans.Microwave Theory Tech.,vol.MTT-40,pp.2106-2108,Nov.1992.
    [110]D.Jiao and J.M.Jin,Time-domain finite-element simulation of cavity-backede microstrip patch antennas,Microwave Opt.Tech.Lett.,2002.
    [111]P.P.Silvester and R.L.Ferrari,Finite Elements for Electrical Engineers(2nd edition).Cambridge:Cambridge University Press,1990.
    [112]N.M.Newmark,A method of computation for structural dynamics,J.Eng.Mechanics Div.,Proc.Am.Soc.Civil Eng.,vol.85,no.EM 3,pp.67-94,July 1959.
    [113]葛德彪,闫玉波.电磁波时域有限差分方法.西安电子科技大学出版社,2002.
    [114]Michael L,van Blaricum and Raj Mittra.A Technique for Extracting the Poles and Response.IEEE Trans.Antennas and Propagation,vol.23,no.6,November 1975.
    [115]Kergall J.,Pouliguen P.,Legoff M.,Chevalier Y.,Complex natural resonances extracted from radar signatures measured in time domain with UWB laboratory system.Antennas and Propagation Society International Symposium.vol.2,pp.1346-1349.July 1999.
    [116]Shen J P著,张承义等译.现代处理器设计—超标量处理器基础[M].北京:电子工业出版社,2004.1-2.
    [117]Macedonia M.The GPU Enters Computing's Mainstream[J].IEEE Computer,2003,36(10):106-108.
    [118]Comba J L D.,Dietrich C A,and Christian A P.Computation on GPUs:From a Programmable Pipeline to an Efficient Stream Processor[J].Revista de Informtic Tericae Aplicada,2003,X(2):41-70.
    [119]Geer D.Taking the Graphics Processor beyond Graphics[J].IEEE Computer,2005,38(9):14-16.
    [120]Cg user's manual[OL].http://developer.nvidia.com/Cg.
    [121]张杨,图形处理器并行计算应用研究[D].成都:西南交通大学.2006.
    [122]刘昆,王晓斌,廖成.图形处理器(GPU)加速时域有限元的二维辐射计算.电波科学学报.2008,23(1):111-114.
    [123]Fernando.R.,The Cg Tutorial[M],Addison-Wesley,April,2003.
    [124]K.liu,Xiao-bing Wang,Yang Zhang,Cheng Liao,"Acceleration of Time-Domain Finite Element Method(TD-FEM) Using Graphics Processor Units(GPU)",International Symposium on Antennas,Propagation and EM Theory Proceedings,pp.807-810,Guilin,China,2006.
    [125]吴恩华:图形处理器用于通用计算的技术、现状及其挑战。软件学报,2004,15(10)。
    [126]Wu EH,Liu YQ.General purpose computation on GPU.Chinese Journal of Computer Aided Design & Computer Graphics,2004,16(5):601-612.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700