30MeV/u~(40)Ar+~(115)In反应中的IMF能谱

测量了 3 0MeV/u40 Ar束流轰击115 In反应靶时出射的中等质量碎片(IMF)的能谱 .通过比较在不同出射角度发射的同一种IMF的能谱 ,得到了IMF的发射机制随出射角的演化关系 .在假定运动源速度和出射离子库仑位垒不随出射角变化的条件下 ,通过改变发射源的强度和核温度参量完成了对实验测得能谱的运动源拟合 ,并详细讨论了 3个发射源所占份额和源强随出射角的演化关系 .

阿尔文涨落与地球磁层亚暴

本文讨论了一种地球磁层的亚暴机制。当行星际磁场有大的南向分量时,磁层的位形可由基本闭式转变为开式。磁鞘中的阿尔文波可以携带超过10~(18)尔格/秒的能流传入磁层尾部,并将能量耗散于等离子体片中。等离子体片中的粒子被加热和加速后,注入近地空间,产生环电流和极区亚暴。计算了剪切流场中阿尔文波的传播过程,以及磁层中阿尔文波的耗散。将本文的结算与[4]中的结果合在一起,可以说明当行星际磁场转向南时,容易发生地球磁层亚暴,但这两者并非一一对应的关系,行星际磁场没有南向分量时也可以发生地球磁层亚暴。

用于热核同位旋效应研究的BGO探测器阵列

报道了35MeV/u ~(36)Ar+~(112,124)Sn反应中用于中等质量碎片(IMF)(3≤Z≤5)关联函数测量的ΔE(Si)+E(BGO)望远镜阵列的测试结果。实验表明该探测器阵列的几何安排有利于提高小角关联测量中的动量分辨,可用于IMF发射时标的同位旋依赖性的研究。

35MeV／u~(36)Ar+~(112,124)Sn反应中中等质量碎片关联函数的入射道依赖

测量了35MeV／u36Ar+112,124Sn反应中小角关联出射的中等质量碎片(IMF)约化速度关联函数．结果表明36Ar+124Sn反应系统中的约化速度关联函数在小约化速度处的反关联程度比36Ar+112Sn反应系统中的强,表现出明显的入射道依赖性．考察出射粒子对的单核子总动量时,发现这种差异主要来自于高动量粒子对的贡献．用三体弹道理论模型MENEKA分别计算了两个系统的IMF发射时标,在36Ar+112Sn反应系统中约为150fm／c,而在36Ar+124Sn反应系统中,约为120fm／c．同位旋相关的量子分子动力学计算表明,36Ar+124Sn系统中IMF的发射时间谱比36Ar+112Sn系统略有前移,相应地,其中心密度从最高点随时间的下降亦比36Ar+112Sn系统略快．

中能重离子碰撞中带电粒子多重性的同位旋效应

利用兰州4π带电粒子探测器阵列测量带电粒子多重性,研究了55MeV/u ~(40)Ar+~(58,64)Ni核反应中He和中等质量碎片的产额与反应系统同位旋的关系,以及这种同位旋效应与反应系统的碰撞参数(即碰撞的激烈程度)、系统的激发能的变化关系,对两个反应系统,观察到带电粒子多重性中He的比分随带电粒子多重性的增加而增大,带电粒子多重性中IMF的比分随带电粒子多重性的增加而先增大,后减小的规律.两个反应系统虽然具有相同的核电荷数,但轻粒子He和中等质量碎片在多重性中的比分有明显的同位旋相关性。

中能核反应中的统计发射过程及核温度参量

在过去的几年里，利用兰州重离子加速器（HIRFL）提供的束流，以及在OUVERTURE合作研究中，利用意大利国家核物理研究院南方实验室（INFN-LNS）超导回旋提供的束流，进行了多次中能区重离子核反应实验研究工作。如，最初的46.7MeV/u ~(12)C+~(58)Ni，~(115)In，~(197)Au的实验及30MeV/u ~(40)Ar+~(58)Ni，~(64)Ni~(115)In和30MeV/u Ni轰击Ni，Au，Al在MULTICS+MEDEA：4π装置上进行的实验工作。此外，本人还从事过一些理论研究工作，包括多粒子散射形式理论和相关数学物理问题研究，量子分子动力学和量子统计模型计算。本文是从事这些核物理研究工作的积累，主要侧重于实验结果的物理内容分析，而不强调实验技术，数据处理的技巧。主要的物理内容有以下几个方面：1.对于利用双同位素产额比提取同位素核温度的方法进行研究推广，使得对于实验中碰到的仅有部分能谱可以实现很好同位素分辩的情况，即使不能得到总的同位素产额，仅仅通过一段能区的同位素产额也可提取核温度。用于具体的实验研究工作中后，对于46.7MeV/u ~(12)C+~(58)Ni，~(115)In，~(197)Au核反应过程，同一体系利用这种方法得到的同位素核温度和利用粒子非稳态布居提取的核温度一致。2.围绕核反应过程中核温度的参量的提取，对于双同位素产额比与核温度的刻度关系进行了分析研究，通过计算考虑中等质量碎片（IMF）内部激发能后的内部配分函数表明，中等质量碎片的内部激发对刻度关系有重要影响。零阶近似下区域密度近似的结果和Gemini模拟计算的结果反映了相同的情况。3.研究核反应机制，多个粒子散射的形式理论的必需的，对于两体散射，其形式理论已经比较成熟，但是对于多个粒子散射问题出现的严重的困难是多体Lippmann-Schwinger方程无唯一收敛的解。作为一种探索性的研究工作，开展了多体散射理论研究工作，发展了一些具有普遍意义的数学物理方法。在本项研究工作中，通过能基础数学中的约当引理的推广，发现一个特例：对非连接图，Lippmann-Schwingwer方程存在收敛的解，因此多体散射形式理论，有可能重新建立。由于核力和多体问题是当今核物理研究的两大难点，世界各国的科学家都在努力以图攻克它们，而且多体问题还是物理学的其它许多领域的难题，因而多体散射还是引起诸多研究学科广泛兴趣的课题。通过发展一些新的数学理论和方法，我们已得到一些有意义的结果。4.将量子分子动力学这种中高能量区域所用的理论分析方法扩展至较低能区，通过对相空间中初始位置和动量抽样增加限制条件。如结合能和实验值要求一致，平均核势，核内Pauli阻塞更强一些，在演化中能量和动量守恒等等。得到一个很稳定的初始基态。均方半径保持不弥散的时间可达1600fm/c，用于研究10.6MeV/u Ne~(20)+Al~(27)的实验分析过程中。另外，量子统计模型（QSM）主要描述中心核-核碰撞，将它和碎裂模型结合，作一些改进后，可以对核-核碰撞进行统一描述。5.在中能核反应研究中发现，核反应过程中有大量的中子，轻带电粒子以及中等质量碎片发射出来，可以将这些粒子发射机制大致分为两大类。其中一类可以归结为动力学发射过程的产物。另一类则可以归于统计发射的产物。在低能核反应中，其发射能谱的斜率的负倒数，可作为复合核的核温度。而在中能重离子核反应中，其发射能谱变得很复杂，不再具有Maxwell分布。通常的三源拟合所给出的温度参数，已不能反映物理实质。提出多阶矩分析方法用于分析中能核反应中统计发射规律及受动力学过程的影响。

35 A×MeV 40Ar + 159Tb/197Au 反应中高激发核性质与反应机制研究

中能重离子碰撞的反应机制及其形成的高激发核的性质研究是当前中能重离子核物理研究的重要领域。本文通过对35 MeV/u~(40)Ar+~(159)Tb/~(197)Au反应中出射的轻粒子、中等质量碎片(IMF)和复合核裂变碎片的符合测量，研究了中等质量碎片的发射机制、类靶核的形成机制及其裂变时标、复合核的温度与激发能以及两者之间的关系等。利用前角探测器组的测量对前角区出射的IMF的单举和关联数据进行了分析。由弹核碎裂理论中的Gozdhaber关系式拟合了类弹碎片的能谱，由此提取的弹核的动量分布约化宽度为σ_0 ≈ 100MevV/c。通过Gemini模拟与IMF-IMF关联数据的对比揭示出激发的类弹碎片的裂变是关联IMF的一个重要来源。由IMF与裂变碎片以及IMF与后角轻粒子的关联数据显示出，在前角区域，当IMF的能量低于10 MeV/u时，复合核的蒸发对较轻IMF有较大的贡献；当IMF能量在10-24 MeV/u的区域时，IMF能量与复合核平均速度以及激发能之间存在着很强的依赖关系，说明这些IMF主要来源与弹靶之间的耗散反应；当IMF的能量超过24 MeV/u后，IMF能量与复合核平均速度以及激发能之间没有关联或关联很弱，说明在弹核碎裂的同时，靶核也发生了碎裂，它们共同贡献出部分核子形成了"中速源"或"颈部"。通过对复合核裂变时出射的关联裂片以及后角轻带电粒子的测量，提取了复合核的激发能以及在不同激发能能下所对应的核温度。对于中心碰撞，Au靶形成了E~* = 4MeV/u、T_(init) = 7.1 MeV的热核，Tb靶形成了E~* = 4.7MeV/u、T_(init) = 7.7MeV的热核。对其量热曲线的研究显示出，在当前入射能下没有观测到原子核液气相变的特征。由类弹碎片（PLF)与类靶核两个裂片的关联数据，提取了类靶核裂变轴在反应平面内的转动角度，并由此计算出类靶核的、裂变时标。在对称裂变时，类靶核的寿命约5 * 10~(-20) s，在非对称裂变（η =0.6）时降为约1.7 * 10~(-21) S。同时我们认为类靶核的裂变存在两种机制，即动力学裂变和统计裂变，分别对应于裂变轴在反应平面内的非各相同性分别和各相同性分布。

中能重离子碰撞中激发热核性质同位旋效应的实验研究

中能重离子碰撞中的核反应机制及其形成的高激发热核的性质是中能重离子物理研究的重要领域，而高激发热核性质的同位旋效应研究是这一领域的热点之一。选取了具有不同N亿比的反应体系以研究激发热核性质的同位旋效应。本论文涉及的反应系统三对共六个反应体系：55MeV／u~(40)Ar+~(58.64)Ni、30MeV/u~(40)Ar+~(112,124)Sn、35Mev／u~(36)Ar+~(112,124)Sn，这六个反应体系的N／z比分别为1.13，1.26、1.24，1.41、1.18，1.35。分别从带电粒子多重性、相对态布居核温度、关联函数等角度研究了这三对反应体系高激发热核性质的同位旋效应。在55MeV／u 40Ar+58,64Ni核反应中，用兰州4π带电粒子探测器阵列测量带电粒子多重性，研究了He和中等质量碎片的产额与反应系统的同位旋的关系，以及这种同位旋效应与反应系统的碰撞参数（即碰撞的激烈程度）、系统的激发能的变化关系。对两个反应系统，观察到带电粒子多重性中He的比分随带电粒子多重性的增加而增大，带电粒子多重性中IMF的比分随带电粒子多重性的增加而先增大，后减小的规律。两个反应系统虽然具有相同的核电荷数，但轻粒子He和中等质量碎片在多重性中的比分有明显的同位旋相关性。在30Mev／u40Sn、35MeV/u~(40)Ar~(112，124)Sn、35Mev／u 36Ar＋112,124Sn反应中用13单元望远镜探测器阵列测量了小角关联粒子。由价a关联函数提取了30Mev／u 40Ar＋112，12Sn反应系统中激发热核的态布居核温。对于不同同位旋反应系统舜UAr＋112Sn和4VAr+124Sn，提取的相对杰布居核温度分别是4.18+0.28/0.21MEV和4.10士0.22/0.20MeV；考察态布居核温度和粒子能量的关薰时，观察到两个系统的发射温度均随着粒子能量的增加而降低，缺中子系统40Ar＋l12Sn中由低能时的5.13士.30/0.26MEV降低到高能时的3.87士0.37/0.29MeV，丰中子系统40Ar＋124Sn中由低能时的5．39士0.30/0.26MeV降低到高能时的3．32士MeV。讨论了这种布居态核温度的同位旋相关性。在35Mev／u 36Ar＋112，124Sn反应中提取了洲F（3‘25）的约化速度关联函数。相对丰中子36Ar＋124Sn系统的IMF关联函数在小约化速度处反关联程度更强，表明36Ai＋124Sn系统的发射IMF的平均时间更短。用MENEKA程序提取了两个系统IMF的平均发射时间，36Ai＋112sn反应中IMF的发射时间约为150fm/c，而36Ar十124Sn反应IME的发射时间稍短，约为120fm／c。以关联IMF的单核子总能量／动量为窗条件，发现低能IMF关联函数几乎没有差别，而高能IMF关联函数在小约化速度处的差别更为明显，表明两个系统IMF关联函数的同位旋效应可能来自于IMF的早期发射。为了得到进一步的信息，我们提取了高动量窗条件下的IMF发射时间，它们比平均发射时间短，36Ar＋112Sn反应中高能IMF的发射时间约为100蒯c，而36Ai＋124Sn反应中IMF的发射时间则更短，约为50fm/c。

30MeV/u~(40)Ar+~(nat)Ag反应中中等质量碎片的发射机制及其发射时间测量

中等质量碎片（IMF）的发射是中能重离子碰撞的一大特点，且随着轰击能量的升高其产额迅速增加。通过测量30MeV/u 40Ar+natAg反应中出射的中等质量碎片，研究了高激发热核的发射机制、发射时标和时空演化规律。 对反应中发射碎片能谱的运动源拟合表明，前角区（11°—22°）中等质量碎片（IMF）来自于三种成分：类靶源、类弹源和中速源成分，其中类弹和中速成分占主导，关联测量的IMF能谱拟合得到，一个来自于类弹源而另一个来自于中速源成分的事件占有相当大的比例。 利用小相对角度内的两碎片关联测量，研究高激发核衰变中中等质量碎片的发射时标和寿命。IMF发射时间随能量的变化很大，从低能碎片的250 fm/c到高能粒子的100 fm/c，表明在此能量下，反应中出射的IMF主要来自于相继两体衰变。通过与其它实验的比较可知，随着束流能量的升高，IMF发射时间由相继两体衰变向多重碎裂过渡。 IMF时空演化研究表明，发射空间的大小对IMF关联函数的影响主要来自于发射源的核物质密度而几乎不依赖于发射源的质量数。

30MeV/u~(40)Ar引起反应中的碎片发射机制

本论文中重点研究了30MeV/u 40Ar+58Ni，64Ni和115In反应中中等质量碎片（IMF）的发射机制。实验中，测量了实验室系5°～140°角度范围内出射碎片的能谱和角分布。对前角区出射的IMF（3≤Z≤13）实现了同位素鉴别，对中后角区出射的碎片在低探测阈（<2MeV/u ）的前提下实现了直到Z～30的元素鉴别。 用运动源模型对不同角度下出射的碎片能谱进行了分析和讨论，并结合角分布特征定性地研究了碎片的三个发射源。通过对各源的贡献随角度以及出射碎片电荷数Z的演化，观察到：类弹源主要发射的是那些前角区出射的、接近束流速度的高能碎片；中等速度源的发射是中角度区出射碎片和前角区低能碎片的主要来源；后角区出射的碎片则主要来自于类熔合源的发射。并观察到相对于类熔合源非平衡源更容易发射较轻的碎片。 通过对前角区出射IMF（3≤Z≤13）的能谱和同位素分布的分析，确定了那些基本保持束流速度的碎片主要来自于弹核碎裂过程。用各种模型对实验同位素分布进行了拟合，发现Sümmerer等人给出的经验公式和abrasion-ablation模型均能比较满意地拟合实验同位素分布的宽度和峰位。同时也观察到abrasion-ablation模型计算对奇Z元素的同位素分布能给出较好的拟合，但对偶Z元素的同位素分布，计算结果与实验值相比出现向丰中子方向的系统性偏移（～lamu）。另外，还着重研究了这些产物的靶核相关性问题。通过系统性分析以及同位旋相关的量子分子动力学（IQMD）模型计算，得出了弹核碎裂产物的靶核相关性是源于靶核表面中子与质子分布的不同和平均场及核子核子相互作用的同位旋效应相关。并且，通过计算还指出了靶核中子皮的厚度对于用弹核碎裂方法产生丰中子同位素的重要性。 通过用统计模型拟合后角区出射碎片的电荷分布，指出了这些碎片主要来自于非完全熔合过程中形成的复合系统的统计发射。 实验中观察到30MeV/u 40Ar轰击Ni和In靶在中角区出射的碎片的电荷分布有不同的特征。相对来说，前者更服从power-law，后者则倾向于服从指数规律。结合核态方程和QMD模型计算分析得出，在30MeV/u 的40Ar引起的反应中，对于弾靶质量接近对称的碰撞对所形成的系统，在一定的碰撞条件下可能已进入spinodal区，而非对称碰撞对的碰撞中压缩能还难以使得系统在膨胀时进入力学不稳定区。从而对观察到的实验现象进行了说明，并认为30MeV/u 40Ar轰击与其质量接近的靶核的反应中出射的碎片中已有相当数量的动力学发射成分的贡献。 利用实验中在前角区出射的接近束流速度的碎片的同位素产额提取了类弹源的温度参数。观察到直接由实验同位素产额比得到的表观核温度与所选择的同位素组合有较强烈的依赖关系，而与靶核和实验室探测角基本无关。利用M.B.Tsang等人给出的修正方法，提取了经边馈效应修正后的发射源温度。得到该温度值与反应靶、探测角以及用于提出温度的同位素组合没有明显的依赖关系，均为～4MeV。并用abrasion-ablation模型计算进行了讨论，得到在假定能级密度参数的倒数k=10MeV时，该温度值与abrasion-ablation模型计算是一致的。

中能区~(40)Ar + ~(197)Au/~(159)Tb/~(nat)Ag反应中高激发热核性质研究

通过对中能区Ar+Au/Tb/Ag反应中高激发核的发射时标、发射次序、发射机制、核温度、集体运动等衰变特性的研究，提取了轻粒子和中等质量碎片（IMF）的发射时间，IMF发射时间随束流能的升高而变短，发射机制逐渐由相继衰变过渡至多重碎裂。研究了轻粒子和碎片间的发射次序，对高能粒子和碎片，轻粒子先于碎片发射，而低能时，则为碎片先于轻粒子发射。IMF发射成分与角度和碰撞参数有关，前角区来自于弹核碎裂，后角区来自于类靶热核的蒸发。在平面和出平面研究表明，中速粒子和碎片为在平面发射占主导，即存在类转动效应；对轻粒子，转动效应随粒子质量增加而增加；对中速产物均观测到该效应随碰撞参数饿增大而增大。采用几种不同的方法提取了热核的核温度，研究了不同方法之间的区别。

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

HHT在电磁信号处理中的应用

In modem signal Processing，non－linear，non－Gaussian and non－stable signals are usually the analyzed and Processed objects，especially non－stable signals. The convention always to analyze and Process non－stable signals are： short time Fourier transform，Wigner－Ville distribution，wavelet Transform and so on. But the above three algorithms are all based on Fourier Transform，so they all have the shortcoming of Fourier Analysis and cannot get rid of the localization of it. Hilbert－Huang Transform is a new non－stable signal processing technology，proposed by N. E. Huang in 1998. It is composed of Empirical Mode Decomposition (referred to as EMD) and Hilbert Spectral Analysis (referred to as HSA). After EMD Processing，any non－stable signal will be decomposed to a series of data sequences with different scales. Each sequence is called an Intrinsic Mode Function (referred to as IMF). And then the energy distribution plots of the original non－stable signal can be found by summing all the Hilbert spectrums of each IMF. In essence，this algorithm makes the non－stable signals become stable and decomposes the fluctuations and tendencies of different scales by degrees and at last describes the frequency components with instantaneous frequency and energy instead of the total frequency and energy in Fourier Spectral Analysis. In this case，the shortcoming of using many fake harmonic waves to describe non－linear and non－stable signals in Fourier Transform can be avoided. This Paper researches in the following parts： Firstly，This paper introduce the history and development of HHT，subsequently the characters and main issues of HHT. This paper briefly introduced the basic realization principles and algorithms of Hilbert－Huang transformation and confirms its validity by simulations. Secondly, This paper discuss on some shortcoming of HHT. By using FFT interpolation, we solve the problem of IMF instability and instantaneous frequency undulate which are caused by the insufficiency of sampling rate. As to the bound effect caused by the limitation of envelop algorithm of HHT, we use the wave characteristic matching method, and have good result. Thirdly, This paper do some deeply research on the application of HHT in electromagnetism signals processing. Based on the analysis of actual data examples, we discussed its application in electromagnetism signals processing and noise suppression. Using empirical mode decomposition method and multi-scale filter characteristics can effectively analyze the noise distribution of electromagnetism signal and suppress interference processing and information interpretability. It has been founded that selecting electromagnetism signal sessions using Hilbert time-frequency energy spectrum is helpful to improve signal quality and enhance the quality of data.