基于近似l0范数最小化的NMR波谱稀疏重建算法

在核磁共振（NMR）波谱中,过长的数据采集时间会使很多化学以及分子生物学领域的高分辨率多维谱应用难以实现. 传统的解决办法是使用随机非均匀采样代替奈奎斯特采样,但这样会使谱图质量受损. 压缩传感的出现为此提供了更好的解决办法,合适的压缩传感重建算法可以通过很少的随机非均匀采样将谱图高质量的重建出来. 该文先介绍了一种可用于谱图重建的压缩传感重建算法,名为“平滑l₀范数最小化法”,然后针对该算法对采样噪声鲁棒性较差的缺点进行了改进. 通过将改进后的算法与原算法在一维实数域信号以及NMR波谱信号重建实验中进行对比后表明,改进后的算法对噪声的鲁棒性明显提高,并能获得更好的重建性能.     

The stochastic collocation method for radiation transport in random media

Stochastic spectral expansions are used to represent random input parameters and the random unknown solution to describe radiation transport in random media. The total macroscopic cross section is taken to be a spatially continuous log-normal random process with known covariance function and expressed as a memoryless transformation of a Gaussian random process. The Karhunen-Loève expansion is applied to represent the spatially continuous random cross section in terms of a finite number of discrete Gaussian random variables. The angular flux is then expanded in terms of Hermite polynomials and, using a quadrature-based stochastic collocation method, the expansion coefficients are shown to satisfy uncoupled deterministic transport equations. Sparse grid Gauss quadrature rules are investigated to establish the efficacy of the polynomial chaos-collocation scheme. Numerical results for the mean and standard deviation of the scalar flux as well as probability density functions of the scalar flux and transmission function are obtained for a deterministic incident source, contrasting between absorbing and diffusive media.     

Suppression of sampling artefacts in high-resolution four-dimensional NMR spectra using signal separation algorithm

The development of non-uniform sampling (NUS) strategies permits to obtain high-dimensional spectra with increased resolution in significantly reduced experimental time. We extended a previously proposed signal separation algorithm (SSA) to process sparse four-dimensional NMR data. It is employed for two experiments carried out for a partially unstructured 114-residue construct of chicken Engrailed 2 protein, namely 4D HCCH-TOCSY and 4D C,N-edited NOESY. The SSA allowed us to obtain high-quality spectra using only as little as 0.16% of the available samples, with low sampling artefacts approaching the thermal noise level in most spectral regions. It is demonstrated that NUS 4D HCCH-TOCSY is dominated by sampling noise and requires efficient artefact suppression. On the other hand, 4D C,N-edited NOESY is a particularly attractive experiment for NUS, as the absence of diagonal peaks renders the problem of artefacts less critical. We also present a transverse-relaxation optimized sequence for HMQC that is especially designed for longer evolution periods in the indirectly detected proton dimension in high-dimensional pulse sequences. In conjunction with novel sampling strategies and efficient processing methods, this improvement enabled us to obtain unique structural information about aliphatic-amide contacts.     

Analysis of stochastic bifurcation and chaos in stochastic Duffing--van der Pol system via Chebyshev polynomial approximation

The Chebyshev polynomial approximation is applied to investigate the stochastic period-doubling bifurcation and chaos problems of a stochastic Duffing--van der Pol system with bounded random parameter of exponential probability density function subjected to a harmonic excitation. Firstly the stochastic system is reduced into its equivalent deterministic one, and then the responses of stochastic system can be obtained by numerical methods. Nonlinear dynamical behaviour related to stochastic period-doubling bifurcation and chaos in the stochastic system is explored. Numerical simulations show that similar to its counterpart in deterministic nonlinear system of stochastic period-doubling bifurcation and chaos may occur in the stochastic Duffing--van der Pol system even for weak intensity of random parameter. Simply increasing the intensity of the random parameter may result in the period-doubling bifurcation which is absent from the deterministic system.     

Law of addition in random matrix theory

We discuss the problem of adding random matrices, which enables us to study Hamiltonians consisting of a deterministic term plus a random term. Using a diagrammatic approach and introducing the concept of “gluon connectedness”, we calculate the density of energy levels for a wide class of probability distributions governing the random term, thus generalizing a result obtained recently by Brézin, Hikami and Zee. The method used here may be applied to a broad class of problems involving random matrices.     

Anisotropic hinge model for polarization-mode dispersion in installed fibers

We discuss a generalized waveplate hinge model to characterize anisotropic effects associated with the hinge model of polarization-mode dispersion in installed systems. In this model, the action of the hinges is a random time-dependent rotation about a fixed axis. We obtain the probability density function of the differential group delay and the outage probability of an individual wavelength band using a combination of importance sampling and the cross-entropy method, and we then compute the noncompliant capacity ratio by averaging over wavelength bands. The results show that there are significant differences between the outage statistics predicted by isotropic and anisotropic hinge models.     

Stationary response of nonlinear magneto-piezoelectric energy harvester systems under stochastic excitation

Recent years have shown increasing interest of researchers in energy harvesting systems designed to generate electrical energy from ambient energy sources, such as mechanical excitations. In a lot of cases excitation patterns of such systems exhibit random rather than deterministic behaviour with broad-band frequency spectra. In this paper, we study the efficiency of vibration energy harvesting systems with stochastic ambient excitations by solving corresponding Fokker-Planck equations. In the system under consideration, mechanical energy is transformed by a piezoelectric transducer in the presence of mechanical potential functions which are governed by magnetic fields applied to the device. Depending on the magnet positions and orientations the vibrating piezo beam system is subject to characteristic potential functions, including single and double well shapes. Considering random excitation, the probability density function (pdf) of the state variables can be calculated by solving the corresponding Fokker-Planck equation. For this purpose, the pdf is expanded into orthogonal polynomials specially adapted to the problem and the residual is minimized by a Galerkin procedure. The power output has been estimated as a function of basic potential function parameters determining the characteristic pdf shape.     

基于耗损失效模型的金属化膜脉冲电容器可靠性评估

 高储能密度金属化膜脉冲电容器是惯性约束聚变装置的关键元器件，由于其“自愈”特性，在短时间内很难得到它的失效数据。通过分析电容器的失效机理，给出了金属化膜脉冲电容器的一个耗损失效模型，推导了该模型的失效概率密度函数和分布函数，并利用电容器的性能衰退数据对其进行了可靠性分析。所选的某型金属化膜脉冲电容器未知参数估计值为0.000 119 4和0.006 7，将该值代入失效分布函数和概率密度函数中，从而确定电容器的失效模型，由此模型求得该型电容器充放电10 000次的可靠度为0.988 5，预计寿命为23 461次充放电。在工程实践中使用该模型对该型电容器进行可靠性分析可节约大量的试验成本。     

Precise measurement of the top-quark mass in the lepton+jets topology at CDF II

We present a measurement of the mass of the top quark from proton-antiproton collisions recorded at the CDF experiment in Run II of the Fermilab Tevatron. We analyze events from the single lepton plus jets final state (tt-->W(+)bW(-)b-->lnubqq'b). The top-quark mass is extracted using a direct calculation of the probability density that each event corresponds to the tt final state. The probability is a function of both the mass of the top quark and the energy scale of the calorimeter jets, which is constrained in situ by the hadronic W boson mass. Using 167 events observed in 955 pb(-1) of integrated luminosity, we achieve the single most precise measurement of the top-quark mass, 170.8+/-2.2(stat.)+/-1.4(syst.) GeV/c(2).     

基于低秩矩阵的非均匀采样NMR波谱重建进展

多维核磁共振（Nuclear Magnetic Resonance，NMR）利用多维波谱来分析分子结构，被广泛用于化学、生物学和医学等领域，但信号采样时间随波谱维度和采样点数增加而迅速增长．非均匀采样通过降低间接维采样点数来加速数据采集，并引入合理的重建方法获得完整的NMR波谱．如何快速重建高质量的波谱，是NMR信号处理研究的前沿．本文主要综述近年来基于低秩矩阵的NMR波谱重建方法的发展．首先介绍了低秩矩阵的相关数学基础；然后从一般低秩矩阵和结构化低秩汉克尔矩阵两个角度来论述重建模型，并讨论相关的NMR波谱应用；最后分析了该技术存在的不足，并展望其未来发展的趋势．     

Stochastic period-doubling bifurcation analysis of stochastic Bonhoeffer--van der Pol system

In this paper, the Chebyshev polynomial approximation is applied to the problem of stochastic period-doubling bifurcation of a stochastic Bonhoeffer--van der Pol (BVP for short) system with a bounded random parameter. In the analysis, the stochastic BVP system is transformed by the Chebyshev polynomial approximation into an equivalent deterministic system, whose response can be readily obtained by conventional numerical methods. In this way we have explored plenty of stochastic period-doubling bifurcation phenomena of the stochastic BVP system. The numerical simulations show that the behaviour of the stochastic period-doubling bifurcation in the stochastic BVP system is by and large similar to that in the deterministic mean-parameter BVP system, but there are still some featured differences between them. For example, in the stochastic dynamic system the period-doubling bifurcation point diffuses into a critical interval and the location of the critical interval shifts with the variation of intensity of the random parameter. The obtained results show that Chebyshev polynomial approximation is an effective approach to dynamical problems in some typical nonlinear systems with a bounded random parameter of an arch-like probability density function.     

Measurement of the Bs0-Bs0 oscillation frequency

We present the first precise measurement of the Bs0-Bs0 oscillation frequency Deltams. We use 1 fb-1 of data from pp collisions at sqrts=1.96 TeV collected with the CDF II detector at the Fermilab Tevatron. The sample contains signals of 3600 fully reconstructed hadronic Bs decays and 37,000 partially reconstructed semileptonic Bs decays. We measure the probability as a function of proper decay time that the Bs decays with the same, or opposite, flavor as the flavor at production, and we find a signal consistent with Bs0-Bs0 oscillations. The probability that random fluctuations could produce a comparable signal is 0.2%. Under the hypothesis that the signal is due to Bs0-Bs0 oscillations, we measure Deltams=17.31(-0.18)+0.33(stat)+/-0.07(syst) ps-1 and determine |Vtd/Vts|=0.208(-0.002)+0.001(expt)-0.006(+0.008)(theor).     

A non-adapted sparse approximation of PDEs with stochastic inputs

We propose a method for the approximation of solutions of PDEs with stochastic coefficients based on the direct, i.e., non-adapted, sampling of solutions. This sampling can be done by using any legacy code for the deterministic problem as a black box. The method converges in probability (with probabilistic error bounds) as a consequence of sparsity and a concentration of measure phenomenon on the empirical correlation between samples. We show that the method is well suited for truly high-dimensional problems.     

The Milne problem in a continuous stochastic planar medium with Gaussian statistics

The Milne problem of radiative transfer in a planar medium, with isotropic scattering is considered. The medium is assumed to be continuous stochastic medium, with fluctuations described as Gaussian field. Pomraning-Eddington method is used to obtain an explicit form for the radiation energy density in the deterministic case. It depends on two random variables, namely the optical space variable and the optical thickness of the medium. The Gaussian joint probability density function of these two random variables is defined and used to find the ensemble-averaged energy density and the linear extrapolation distance. It is shown that the statistical nature of the medium leads to two quite different solutions of the Milne problem. Numerical results are implemented for the sake of clarification.     

On stochastic complex beam-beam interaction models with Gaussian colored noise

This paper is to continue our study on complex beam-beam interaction models in particle accelerators with random excitations Y. Xu, W. Xu, G.M. Mahmoud, On a complex beam-beam interaction model with random forcing [Physica A 336 (2004) 347-360]. The random noise is taken as the form of exponentially correlated Gaussian colored noise, and the transition probability density function is obtained in terms of a perturbation expansion of the parameter. Then the method of stochastic averaging based on perturbation technique is used to derive a Fokker-Planck equation for the transition probability density function. The solvability condition and the general transforms using the method of characteristics are proposed to obtain the approximate expressions of probability density function to order ε.Also the exact stationary probability density and the first and second moments of the amplitude are obtained, and one can find when the correlation time equals to zero, the result is identical to that derived from the Stratonovich-Khasminskii theorem for the same model under a broad-band excitation in our previous work.     

Modification of instability processes by multiplicative noises

We study two dynamical systems submitted to white and Gaussian random noise acting multiplicatively. The first system is an imperfect pitchfork bifurcation with a noisy departure from onset. The second system is a pitchfork bifurcation in which the noise acts multiplicatively on the non-linear term of lowest order. In both cases noise suppresses some solutions that exist in the deterministic regime. Besides, for the first system, the imperfectness of the bifurcation reduces the regime of on-off intermittency. For the second system, the unstable mode can achieve a jump of finite amplitude at instability but without hysteresis. We finally identify a generic property that is verified by the stationary probability density function of the dynamical variable when a control parameter is varied.