基于全变分最小化和交替方向法的康普顿散射成像重建算法

康普顿散射成像技术利用射线与物质作用后的散射光子信息对物质的电子密度进行成像.与传统的透射成像方式相比,康普顿散射成像具有系统结构灵活、成像对比度高、辐射剂量低等优势,在无损检测、医疗诊断、安全检查等领域有着广阔的应用前景.但其重建问题是一个非线性的逆问题,通常是不适定的,其解对噪声和测量误差非常敏感.为解决此问题,本文结合全变分最小化正则化方法和交替方向法提出了一种新的康普顿散射成像重建算法.该算法首先将问题对应的TV模型转化为与之等价的带约束的优化问题,然后利用增广拉格朗日乘子法将优化问题分解为两个具有解析解的子问题,并通过交替求解子问题使增广拉格朗日函数达到最小,进而得到重建的图像.在仿真实验中,通过与主流的ASD-POCS方法进行对比,证明了该算法在重建精度和重建效率方面的优势.     

基于卷积法TV模型的锥束CT散射校正方法

为了解决CT图像中散射伪影问题,提出并使用修正的泊松分布函数作为散射函数,然后利用卷积模型生成多个散射基图像,将基图像与原始重建图像的线性组合作为最终重建图像,最后,通过使重建图像的全变分函数最小求取线性组合系数最优解。该算法具有计算复杂度小、能够灵活应用于不同的成像对象的特点。实验结果表明,算法能够有效抑制散射导致的杯状和黑色带状伪影,提高图像对比度。     

总变差约束的数据分离最小图像重建模型及其Chambolle-Pock求解算法

基于优化的迭代法,可以结合压缩感知和低秩矩阵等稀疏优化技术高精度地重建图像.其中,总变差最小（total variation minimization,TV）模型是一种简单有效的优化模型.传统的约束TV模型,使用数据保真项为约束项,TV正则项为目标函数.本文研究TV约束的、数据分离最小（TV constrained,data divergence minimization,TVcDM）新型TV模型及其求解算法.详细推导了TVcDM模型的Chambolle-Pock（CP）算法,验证了模型及算法的正确性;分析了算法的收敛行为;评估了模型的稀疏重建能力;分析了模型参数的选择对重建的影响及算法参数对收敛速率的影响.研究表明,TVcDM模型有高精度稀疏重建能力;TVcDM-CP算法确保收敛,但迭代过程中有振荡现象;TV限对重建有重要影响,参数值过大会引入噪声而过小会模糊图像细节;算法参数的不同选取会导致不同的收敛速率.     

On the elastic field and image force of dislocations in anisotropic solids and its application to GaN nanostructures

The determination of the elastic field and image force of dislocation in anisotropic media is a nontrivial problem. This work revisits Eshelby’s sextic anisotropic elasticity theory to obtain the stress field of a screw dislocation in an infinite anisotropic solid. The image force of a dislocation in an anisotropic nanowire is then derived by using the concept of ‘image dislocation’. Moreover, this work proposes to study the image force of nanorods by approximating the three-dimensional shape effect as a height-dependent shape function, which could be obtained through curve fitting of the finite element data. The analytical solution is applied to analyse image forces on different dislocations in GaN nanorods oriented along polar (c-axis) and nonpolar (a, m-axis) directions. The result shows the dislocation dissipation could be more effective in a-GaN but less in m-GaN by comparing with the standard growth of c-GaN. The approach developed in this work is applicable to other material systems. Therefore, it could contribute to a wide range of nanostructure design and fabrication for dislocation-free devices.     

Anisotropic localization behavior of graphene in the presence of diagonal and off-diagonal disorders

Anisotropic localization of Dirac fermions in graphene along both the x and y axes was studied using the transfer-matrix method. The two-parameter scaled behavior around the Dirac points was observed along the x axis with off-diagonal disorder. In contrast, the electronic state along the y axis with armchair edges was delocalized, which can be described well by single parameter scaling theory. This implies that the breakdown of the single-parameter scaling is related to the zigzag edge along the x axis. Furthermore, dimerization induced by the substrate suppresses the two-parameter scaling behavior along the x axis and preserves the delocalized state along the y axis. Our results also demonstrate anisotropic localization in graphene with diagonal disorder that can be tuned by dimerization.     

X-ray phase-contrast method and its application to the study of blood vessels with a model object

Studies of the phase-contrast method are reviewed. The principles of phase-contrast image formation are considered. The potentialities of the method are demonstrated using a model object as an example: a capillary tube filled with air and paraffin.     

用动力学夸克并合描述RHIC重离子碰撞中和各相异性流和它们的标度率

We report the results on anisotropic flows and their scaling forφmesons andΩ(Ω~- ■~ )baryons in Au Au collisions at RHIC,obtained from a dynamical quark coalescence model that uses the quark phase- space information from a multi-phase transport(AMPT)model within the string melting scenario and includes the quark structure of hadrons.     

Analogue correction method of errors and its application to numerical weather prediction

In this paper, an analogue correction method of errors (ACE) based on a complicated atmospheric model is further developed and applied to numerical weather prediction (NWP). The analysis shows that the ACE can effectively reduce model errors by combining the statistical analogue method with the dynamical model together in order that the information of plenty of historical data is utilized in the current complicated NWP model. Furthermore, in the ACE, the differences of the similarities between different historical analogues and the current initial state are considered as the weights for estimating model errors. The results of daily, decad and monthly prediction experiments on a complicated T63 atmospheric model show that the performance of the ACE by correcting model errors based on the estimation of the errors of 4 historical analogue predictions is not only better than that of the scheme of only introducing the correction of the errors of every single analogue prediction, but is also better than that of the T63 model.     

SOA-fiber ring laser and its application to electric field sensing in frequency domain

We have developed a novel optical fiber ring laser using a semiconductor optical amplifier (SOA) as the gain medium, and taking advantage of polarization anisotropy of its gain. The laser can control the oscillation frequency difference between two counterpropagating lights by incorporating a birefringent medium into the ring resonator. Since the frequency difference is proportional to the birefringence, it can be measured to detect the beat frequency generated by combining two counterpropagating lights and the laser can be applied to an optical fiber sensor whose detecting signal is in frequency domain. Electric field sensor was demonstrated with a lithium niobate crystal as a sensing probe.     

Manifestation and post hoc correction of gradient cross-term artifacts in DTI

Cross-terms between imaging and diffusion gradients, unaccounted for during tensor calculations, can lead to erroneous estimation of diffusivity and fractional anisotropy (FA) in regions of isotropic and anisotropic diffusion. Cross-term of magnitude 136.8±1.6 s/mm(2), artificially introduced in the slice-encode direction, caused an increase in FA in isotropic phantom from 0.0546±0.0001 to 0.0996±0.0001, while the change in chimpanzee brain depended on the orientation of the white matter (WM). Mean diffusivity (MD) remained unchanged in isotropic phantom, but increased by ～20% in the WM due to cross-terms. A bias was observed in the principal eigenvectors in both phantom and chimpanzee brain, resulting in significant increase in midline crossing fibers along the bias than perpendicular to it in tractography in chimpanzee brain. Post hoc correction of these artifacts was achieved by estimating the cross-term factors using calibration scans on an isotropic phantom and modifying the b-matrix before tensor calculation. Upon correction, the FA and MD values closely resembled the values obtained from sequence without cross-terms, and the bias in principal eigenvectors was eliminated. Customized sequences involving large b-values, high-resolution imaging, or long diffusion or echo times should therefore be evaluated and any residual cross-terms corrected before implementation.     

CT实验仪的研制及其在物理教学中的应用

介绍了我们研制的CT实验仪的原理及开发,使用情况,包括CT的物理,数学原理,CT装置的结构,性能、工作原理等,文末给出实验结果,并对装置在物理实验教学中的意义作了讨论。     

平流涡度方程及其在2006年Bilis台风分析中的应用

推导得到气压坐标中的动量叉乘形式的垂直涡度方程,这个动量叉乘形式的涡度方程包含了水平风的平流旋转效应,可称为平流涡度方程.由于水平风场的平流作用可由等压面天气图直观分析得到,因此平流涡度方程可方便用于实际天气分析.对2006年的Bilis台风移动过程中由经典涡度方程和平流涡度方程计算得到的垂直涡度倾向进行对比分析发现,二者计算得到的垂直涡度倾向变化的分布形式接近,但平流涡度方程计算得到的倾向的数值明显大于经典涡度方程的数值,正负涡度倾向区也更集中.对Bilis移动过程中的垂直涡度方程和平流涡度方程中各项的计算分析表明,水平风场的平流旋转作用是Bilis发展移动过程中垂直涡度变化的一个主要因素,是造成垂直涡度增强并发展的主要原因.因此,当水平风场平流旋转效应较强时,平流作用对垂直涡度倾向变化起主导作用,可直接用平流项来近似分析Bilis台风的涡度变化.而平流涡度方程中地转涡度和散度项的变化趋势与Bilis台风的移动路径有较好的一致性,这一项对台风的移动路径预报有更好的指示意义.     

Exact Solution of a(2+1)-Dimensional Anisotropic Star in Finch and Skea Spacetime

We provide a new class of interior solution of a(2+1)-dimensional anisotropic star in Finch and Skea spacetime corresponding to the BTZ black hole. We develop the model by considering the MIT bag model EOS and a particular ansatz for the metric function grrproposed by Finch and Skea [M.R. Finch and J.E.F. Skea, Class. Quantum.Grav. 6(1989) 467]. Our model is free from central singularity and satisfies all the physical requirements for the acceptability of the model.     

A time-series approach to measuring node similarity in networks and its application to community detection

A central concept in network analysis is that of similarity between nodes. In this paper, we introduce a dynamic time-series approach to quantifying the similarity between nodes in networks. The problem of measuring node similarity is exquisitely embedded into the framework of time series for state evolution of nodes. We develop a deterministic parameter-free diffusion model to drive the dynamic evolution of node states, and produce a unique time series for each source node. Then we introduce a measure quantifying how far all the other nodes are located from each source one. Following this measure, a quantity called dissimilarity index is proposed to signify the extent of similarity between nodes. Thereof, our dissimilarity index gives a deep and natural integration between the local and global perspectives of topological structure of networks. Furthermore, we apply our dissimilarity index to unveil community structure in networks, which verifies the proposed dissimilarity index.     

导电平板上任意缝隙填充各向异性介质时TM波散射和传输特性分析

简单介绍了用于分析任意口径问题的通用频域方法——广义网络原理,并利用边界积分法和广义网络原理分析导电平板上任意缝隙填充各向异性介质时TM波的散射及传输特性.由于缝隙填充各向异性介质的情形尚未见公开文献报道,作为验证,将本方法退化计算各向同性介质填充时缝隙的散射和传输特性,并与文献结果进行比较.最后,给出了缝隙填充各向异性介质时的算例.     

超高分辨光学相干层析成像技术与材料检测应用

本文报道了一种超高分辨率谱域光学相干层析成像(SD-OCT)系统. 该系统基于超连续谱激光光源并截取部分光谱作为宽带光源, 其中心波长为665 nm, 光谱半高全宽(FWHM) 230 nm. 系统轴向分辨率0.9 μm, 轴向扫描速率28600行/秒, 横向分辨率3.9 μm, 横向视场1 mm, 最大成像深度0.6 mm(空气中). 利用研制的超高分辨率SD-OCT系统, 对不同型号的工业砂纸精细结构进行了成像, 并与普通SD-OCT的成像结果进行对比, 充分展示了研制系统在材料无损检测中优势.     

拉曼光谱的荧光背景扣除及其用于药物聚类分析

拉曼光谱分析中,由于有机分子或样品中污染物的荧光影响,常会使拉曼光谱产生高背景信号,以致其拉曼光谱吸收信号被淹没。利用自行开发的软件包baselineWavelet,本文对醋酸泼尼松片和格列本脲片的拉曼光谱进行了荧光背景扣除研究,采用主成分分析和随机森林算法对它们进行聚类分析,得到了较好的结果。通过这2种药物的拉曼光谱聚类分析结果,检验了该背景扣除算法的有效性和准确性,并讨论了荧光背景对拉曼光谱聚类分析的影响。结果说明,荧光背景对拉曼光谱聚类分析影响很大,在分析前必须预先扣除。     

Pulsed TV holography and tomography for the study of transient waves in air

Transient waves in air are recorded and reconstructed using pulsed TV holography and computerized tomography (CT). Experiments are performed with an electrical discharge between two electrodes as the acoustic wave source. The free space wave-fronts and pressure fields are reconstructed. Waves reflected and diffracted by different obstacles are also recorded and reconstructed in three dimensions. Speckle averaging and image processing techniques are used to get the high quality projection fields needed for CT reconstruction.     

Exact Solution of a (2+1)-Dimensional Anisotropic Star in Finch and Skea Spacetime

We provide a new class of interior solution of a (2+1)-dimensional anisotropic star in Finch and Skea spacetime corresponding to the BTZ black hole. We develop the model by considering the MIT bag model EOS and a particular ansatz for the metric function g_rr proposed by Finch and Skea [M.R. Finch and J.E.F. Skea, Class. Quantum. Grav. 6 (1989) 467]. Our model is free from central singularity and satisfies all the physical requirements for the acceptability of the model.     

A self-consistent approach to measure preferential attachment in networks and its application to an inherent structure network

Preferential attachment is one possible way to obtain a scale-free network. We develop a self-consistent method to determine whether preferential attachment occurs during the growth of a network, and to extract the preferential attachment rule using time-dependent data. Model networks are grown with known preferential attachment rules to test the method, which is seen to be robust. The method is then applied to a scale-free inherent structure (IS) network, which represents the connections between minima via transition states on a potential energy landscape. Even though this network is static, we can examine the growth of the network as a function of a threshold energy (rather than time), where only those transition states with energies lower than the threshold energy contribute to the network. For these networks we are able to detect the presence of preferential attachment, and this helps to explain the ubiquity of funnels on potential energy landscapes. However, the scale-free degree distribution shows some differences from that of a model network grown using the obtained preferential attachment rules, implying that other factors are also important in the growth process.