利用脑磁图-多重信号分类算法求解真实头模型中磁源定位问题

脑磁图-多重信号分类(MEG-MUSIC)算法与常见的利用时域空域数据定位多偶极子源的全局优化方法相比，具有求解速度快、独立确定各偶极子源的位置、求解难度与偶极子数量无关等优点，在偶极子数量较多的情况下，MEG-MUSIC算法的优势显得更为突出.利用MEG-MUSIC算法求解了真实头模型中脑磁源定位问题. 关键词： 脑磁图-多重信号分类算法 真实头模型 多偶极子     

一种新的求解脑磁逆问题的搜索方法

给出一种新的求解真实头模型下脑磁逆问题的搜索方法.通过不同位置的源的相互关系，由上一个搜索源的计算结果通过简单计算，直接得到下一个搜索源的结果，避免了繁琐耗时的边界元积分方法，简化了求解过程，提高了求解速度. 关键词： 脑磁图 逆问题 搜索方法     

基于三维磁成像的近场磁性体探测

提出一种近场条件下未知磁源的三维磁成像方法.考虑到大多数磁性体不仅受背景磁场磁化，本身也带有较强剩磁，将观测面上的磁场转换为磁场矢量异常模量，并建立目标函数进行最优化求解，以得到符合观测磁场特征的磁性体磁化模型.仿真和实验表明：此方法可有效消除剩磁对反演结果的影响，能够实现对近场多个磁源磁化率分布的成像，验证了所提方法用于探测隐含磁体位置和形状的可行性.     

基于因子分析方法的相位同步脑电源的时-空动力学分析

头皮脑电时间序列的相关性是大脑皮层源的相位同步性的一种体现,因此对相位同步源进行定位,同时找到源对应的时间序列在脑成像研究领域具有重要意义.基于Rössler 模型提出仿真相位同步偶极子源的时间序列的方法,利用时间序列进行同心四层球头模型正演,获得仿真头皮脑电数据.提出了基于最大似然因子分析的相位同步脑电源的时-空动力学分析方法,对仿真和真实头皮脑电数据进行了验证,并与主成分分析法进行对比.仿真实验结果表明:最大似然因子分析法估计的时间序列与仿真源的时间序列具有更高的相关系数,同时估计源与仿真源 关键词： 脑电图 相位同步 因子分析 主成分分析     

非均匀电磁介质中的等效源重构

本文提出了一种用于等效源重构的磁场极值信号法.并通过给定不同的电流偶极子作为信号源,仿真研究了多腔体心脏磁场模型及非均匀介质情况下等效源重构的精度.分析了体电导对心脏磁场的影响,以及磁场极值信号反映的体电导作用.文中将该方法与其他四种源估计方法:磁场梯度极值法、Nelder-Mead单纯形算法、信赖域反射算法和粒子群优化算法作了比较.分析了这些方法的源重构精度和计算所需时间.结果表明,这种针对非均匀介质情况提出的心脏磁场逆问题求解方法有一定的实用价值.     

基于心磁信号的心脏电流偶极子阵列成像及相关性质的研究

超导量子干涉器（SQUID）能探测到微弱的心脏磁场信号. 通过对所得的心磁信号进行分析,可为许多心脏疾病的诊断提供依据. 利用心磁信号,采用极小范数最小二乘法（MNLS）对心脏的电流偶极子阵列进行重建,从而实现了对心脏内部等效电流源的成像. 在使用MNLS进行电流偶极子阵列反演重建的过程中,反演所需的心磁信号,分别由单电流偶极子和电流多极子作为激发源模拟得到,以及由SQUID实际测量得到. 同时,对不同心磁信号反演得到的电流偶极子的分布规律进行了分析. 此外,还给模拟的人体外心磁信号施加了均匀噪声和随机噪声,研究不同信噪比的均匀噪声和随机噪声对电流偶极子阵列重建的影响.     

基于SQUID梯度计的单磁源定位及磁矩反演误差分析

磁源定位及磁矩反演技术在空间探测、无损检测以及目标追踪等领域有着广泛的应用前景。利用磁场梯度张量信息进行磁源定位及磁矩反演的思想是在1975年首次提出的,2006年提出的磁性目标线性方程定位方法以其快速准确求解的特点而备受关注。但是,对于该梯度张量矩阵求逆法进行磁源定位的全方位的误差分析仍有待完善。在简要回顾梯度张量矩阵求逆法进行磁源定位的方法的基础上,提出了一种基于该方法的系统误差和随机误差分析的理论;研究结果表明定位误差是由于将差分近似为微分而引起的系统误差和传感器测量的随机误差构成。以目前具有超高精度、超高灵敏度的超导量子干涉仪作为典型磁场测量传感器,分析了梯度张量矩阵求逆法进行磁源定位及磁矩反演误差与传感器分辨率和相对误差水平之间的关系。     

基于心磁信号的心脏电流偶极子阵列成像及相关性质的研究

超导量子干涉器（SQUID）能探测到微弱的心脏磁场信号. 通过对所得的心磁信号进行分析,可为许多心脏疾病的诊断提供依据. 利用心磁信号,采用极小范数最小二乘法（MNLS）对心脏的电流偶极子阵列进行重建,从而实现了对心脏内部等效电流源的成像. 在使用MNLS进行电流偶极子阵列反演重建的过程中,反演所需的心磁信号,分别由单电流偶极子和电流多极子作为激发源模拟得到,以及由SQUID实际测量得到. 同时,对不同心磁信号反演得到的电流偶极子的分布规律进行了分析. 此外,还给模拟的人体外心磁信号施加了均匀噪声和随机噪声,研究不同信噪比的均匀噪声和随机噪声对电流偶极子阵列重建的影响. 关键词： 心磁信号 超导量子干涉器 电流偶极子阵列重建 极小范数最小二乘法     

磁刺激穴位复杂脑功能网络构建与分析

本文采用互信息方法对磁刺激内关穴过程中的脑电信 号进行了两两通道间非线性时域关联特性分析, 构建了不同频率刺激前、刺激中、刺激后的脑功能网络, 并基于复杂网络理论对脑功能网络的特征进行了深入研究. 结果表明, 磁刺激频率为3 Hz 时, 大脑功能网络的平均度、平均聚类系数和全局效率与刺激前相比均有显著升高, 平均路径长度显著降低, 并且相应脑功能网络的"小世界"属性有所增强, 信息在大脑各区域间的传递更加高效. 本研究首次开展了磁刺激穴位复杂脑功能网络的构建与分析, 为探索磁刺激穴位对大脑神经调节的作用和机理提供新思路和新方法. 关键词： 复杂网络 磁刺激 脑功能网络 互信息     

Tb0.3Dy0.7Fe2单晶中巨磁致伸缩的逆效应

研究了Terfenol-D材料中巨磁致伸缩的逆效应,即磁机械效应.基于Stoner-Wohlfarth(SW)模型,考虑磁晶各向异性和应力各向异性能,依据自由能极小原理,获得了退磁态下Terfenol-D单晶中磁化强度方向和压应力的关系.采用数值方法求解了平衡条件下的非线性方程组.理论结果表明,Terfenol-D巨磁致伸缩单晶中的磁各向异性取决于磁晶各向异性和应力各向异性之间的竞争.在压应力的作用下,Terfenol-D单晶中的磁各向异性由立方向单轴转变.理论和实验结果的比较表明,存在一个临界压应力,使磁致伸缩效应达到极大值.该理论结果还解释了压应力使得Terfenol-D单晶材料难于磁化和磁致伸缩效应出现极大值的实验事实.理论计算不仅为研究这类问题提供了一个更准确的方法,而且其结果也有助于理解类似材料中的磁化过程.     

Bands of Localized Electromagnetic Waves in 3D Random Media

Anderson localization of electromagnetic waves in three-dimensional disordered dielectric structures is studied using a simple yet realistic theoretical model. An effective approach based on analysis of probability distributions, not averages, is developed. The disordered dielectric medium is modeled by a system of randomly distributed electric dipoles. Spectra of certain random matrices are investigated and the possibility of appearance of the continuous band of localized waves emerging in the limit of an infinite medium is indicated. It is shown that localization could be achieved without tuning the frequency of monochromatic electromagnetic waves to match the internal (Mie-type) resonances of individual scatterers. A possible explanation for the lack of experimental evidence for strong localization in 3D as well as suggestions how to make localization experimentally feasible are also given. Rather peculiar requirements for setting in localization in 3D as compared to 2D are indicated.     

Remote dipolar interactions for objective density calibration and flow control of excitonic fluids

In this Letter we suggest a method to observe remote interactions of spatially separated dipolar quantum fluids, and in particular, of dipolar excitons in GaAs bilayer based devices. The method utilizes the static electric dipole moment of trapped dipolar fluids to induce a local potential change on spatially separated test dipoles. We show that such an interaction can be used for model-independent, objective fluid density measurements, an outstanding problem in this field of research, as well as for interfluid exciton flow control and trapping. For a demonstration of the effects on realistic devices, we use a full two-dimensional hydrodynamical model.     

Simulation of Electromagnetic Wave Logging Response in Deviated Wells Based on Vector Finite Element Method

The vector finite element method of tetrahedral elements is used to model 3D electromagnetic wave logging response. The tangential component of the vector field at the mesh edges is used as a degree of freedom to overcome the shortcomings of node-based finite element methods. The algorithm can simulate inhomogeneous media with arbitrary distribution of conductivity and magnetic permeability. The electromagnetic response of well logging tools are studied in dipping bed layers with the borehole and invasion included. In order to simulate realistic logging tools, we take the transmitter antennas consisting of circular wire loops instead of magnetic dipoles. We also investigate the apparent resistivity of inhomogeneous formation for different dip angles.     

Tunable Wigner states with dipolar atoms and molecules

We study the few-body physics of trapped atoms or molecules with electric or magnetic dipole moments aligned by an external field. Using exact numerical diagonalization appropriate for the strongly correlated regime, as well as a classical analysis, we show how Wigner localization emerges with increasing coupling strength. The Wigner states exhibit nontrivial geometries due to the anisotropy of the interaction. This leads to transitions between different Wigner states as the tilt angle of the dipoles with the confining plane is changed. Intriguingly, while the individual Wigner states are well described by a classical analysis, the transitions between different Wigner states are strongly affected by quantum statistics. This can be understood by considering the interplay between quantum-mechanical and spatial symmetry properties. Finally, we demonstrate that our results are relevant to experimentally realistic systems.     

From the Anderson Model on a Strip to the DMPK Equation and Random Matrix Theory

We study weakly disordered quantum wires whose width is large compared to the Fermi wavelength. It is conjectured that such wires display universal metallic behavior as long as their length is shorter than the localization length (which increases with the width). The random matrix theory that accounts for this behavior—the DMPK theory—rests on assumptions that are in general not satisfied by realistic microscopic models. Starting from the Anderson model on a strip, we show that a twofold scaling limit nevertheless allows to recover rigorously the fundaments of DMPK theory, thus opening a way to settle some conjectures on universal metallic behavior.     

Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations

A new type of ultra-high resolution atmospheric global circulation model is developed. The new model is designed to perform “cloud resolving simulations” by directly calculating deep convection and meso-scale circulations, which play key roles not only in the tropical circulations but in the global circulations of the atmosphere. Since cores of deep convection have a few km in horizontal size, they have not directly been resolved by existing atmospheric general circulation models (AGCMs). In order to drastically enhance horizontal resolution, a new framework of a global atmospheric model is required; we adopted nonhydrostatic governing equations and icosahedral grids to the new model, and call it Nonhydrostatic ICosahedral Atmospheric Model (NICAM). In this article, we review governing equations and numerical techniques employed, and present the results from the unique 3.5-km mesh global experiments—with O(10⁹) computational nodes—using realistic topography and land/ocean surface thermal forcing. The results show realistic behaviors of multi-scale convective systems in the tropics, which have not been captured by AGCMs. We also argue future perspective of the roles of the new model in the next generation atmospheric sciences.     

Finite difference computation of head-related transfer function for human hearing

Modeling the head-related transfer function (HRTF) is a key to many applications in spatial audio. To understand and predict the effects of head geometry and the surrounding environment on the HRTF, a three-dimensional finite-difference time domain model (3D FDTD) has been developed to simulate acoustic wave interaction with a human head. A perfectly matched layer (PML) is used to absorb outgoing waves at the truncated boundary of an unbounded medium. An external source is utilized to reduce the computational domain size through the scattered-field/total-field formulation. This numerical model has been validated by analytical solutions for a spherical head model. The 3D FDTD code is then used as a computational tool to predict the HRTF for various scenarios. In particular, a simplified spherical head model is compared to a realistic head model up to about 7 kHz. The HRTF is also computed for a realistic head model in the presence of a wall. It is demonstrated that this 3D FDTD model can be a useful tool for spatial audio applications.     

Simulation study of a magnetocardiogram based on a virtual heart model: effect of a cardiac equivalent source and a volume conductor

In this paper,we present a magnetocardiogram(MCG) simulation study using the boundary element method(BEM) and based on the virtual heart model and the realistic human volume conductor model.The different contributions of cardiac equivalent source models and volume conductor models to the MCG are deeply and comprehensively investigated.The single dipole source model,the multiple dipoles source model and the equivalent double layer(EDL) source model are analysed and compared with the cardiac equivalent source models.Meanwhile,the effect of the volume conductor model on the MCG combined with these cardiac equivalent sources is investigated.The simulation results demonstrate that the cardiac electrophysiological information will be partly missed when only the single dipole source is taken,while the EDL source is a good option for MCG simulation and the effect of the volume conductor is smallest for the EDL source.Therefore,the EDL source is suitable for the study of MCG forward and inverse problems,and more attention should be paid to it in future MCG studies.