Fibonacci序列一维光子晶体色散关系解析解和数值解的对比研究

利用转移矩阵方法给出Fibonacci序列一维光子晶体能带，等频率表面的解析公式，同时给出物理量，如等效相折射率、群速度、波矢量的解析式.利用解析公式以及时域有限差分方法对一维光子晶体的色散关系进行对比研究.结果表明：在强调制一维光子晶体色散关系的计算中，数值方法与解析方法得到的结果总体趋势保持一致.但在弱调制一维光子晶体中，等频率表面的数值结果出现了颠覆性的错误.因而，对于弱调制光子晶体色散关系的数值结果，必须通过解析计算或数值模拟计算等方法进行验证.     

Cost Optimized Non-Contacting Experimental Modal Analysis Using a Smartphone

The vibrations behavior analysis is an essential step in the mechanical design process. Several methods such as analytical modelling, numerical analysis and experimental measurements can be used for this purpose. However, the numerical or analytical models should be validated through experimental measurements, usually expensive. This paper introduces an inexpensive smartphone as an accurate, non-intrusive vibrations’ behavior measurement device. An experimental measurement procedure based on the video processing method is presented. This procedure allows the measurement of the natural frequencies and the mode shapes of a vibrating structure, simply by using a smartphone built-in camera. The experimental results are compared to those obtained using an accurate analytical model, where the natural frequencies error is less than 2.7% and the modal assurance criterion is higher than 0.89. In order to highlight the obtained results, a comparison has been done using a high quality and high frame per second (fps) camera-based measurement of material properties. Since the highest recovered natural frequency and its associated mode shape depend on the frame per second rate of the recorded video, this procedure has great potential in low frequencies problems such as for big structures like buildings and bridges. This validated technique re-introduces the personal smartphone as an accurate inexpensive non-contacting vibration measurement tool.     

非均匀平板波导的色散方程

利用转移矩阵理论分析了任意折射率分布平板波导的传输特性，导出了非均匀平板波导的近似解析色散方程，并指出了ＷＫＢ近似的局限性，数值比较的结果表明，文中所得公式的精确度优于传统的ＷＫＢ方法和其它近似方法。     

A finite difference method for a coupled model of wave propagation in poroelastic materials

A computational method for time-domain multi-physics simulation of wave propagation in a poroelastic medium is presented. The medium is composed of an elastic matrix saturated with a Newtonian fluid, and the method operates on a digital representation of the medium where a distinct material phase and properties are specified at each volume cell. The dynamic response to an acoustic excitation is modeled mathematically with a coupled system of equations: elastic wave equation in the solid matrix and linearized Navier-Stokes equation in the fluid. Implementation of the solution is simplified by introducing a common numerical form for both solid and fluid cells and using a rotated-staggered-grid which allows stable solutions without explicitly handling the fluid-solid boundary conditions. A stability analysis is presented which can be used to select gridding and time step size as a function of material properties. The numerical results are shown to agree with the analytical solution for an idealized porous medium of periodically alternating solid and fluid layers.     

Rotating beamforming – motion-compensation in the frequency domain and application of high-resolution beamforming algorithms

A method is presented for the location of rotating sound sources by a microphone array. In contrast to other methods which are formulated in the time domain, this method works completely in the frequency domain and allows the application of advanced, high resolution beamforming techniques. Following the work by Lowis and Joseph for ducted sources, it is shown that a generalized cross-spectral matrix in the rotating frame of reference can be calculated which can serve as a starting point for advanced beamforming techniques. The Green's function of a moving point source under free space conditions is expressed in spherical coordinates, which are, for numerical reasons, more convenient than cylindrical coordinates. The microphones on the array have to be arranged in a ring. As a practical example high resolution beamforming deconvolution methods are applied to simulated data and measured data from rotating sources.     

Love waves in functionally graded piezoelectric materials by stiffness matrix method

A numerical matrix method relative to the propagation of ultrasonic guided waves in functionally graded piezoelectric heterostructure is given in order to make a comparative study with the respective performances of analytical methods proposed in literature. The preliminary obtained results show a good agreement, however numerical approach has the advantage of conceptual simplicity and flexibility brought about by the stiffness matrix method. The propagation behaviour of Love waves in a functionally graded piezoelectric material (FGPM) is investigated in this article. It involves a thin FGPM layer bonded perfectly to an elastic substrate. The inhomogeneous FGPM heterostructure has been stratified along the depth direction, hence each state can be considered as homogeneous and the ordinary differential equation method is applied. The obtained solutions are used to study the effect of an exponential gradient applied to physical properties. Such numerical approach allows applying different gradient variation for mechanical and electrical properties. For this case, the obtained results reveal opposite effects. The dispersive curves and phase velocities of the Love wave propagation in the layered piezoelectric film are obtained for electrical open and short cases on the free surface, respectively. The effect of gradient coefficients on coupled electromechanical factor, on the stress fields, the electrical potential and the mechanical displacement are discussed, respectively. Illustration is achieved on the well known heterostructure PZT-5H/SiO₂, the obtained results are especially useful in the design of high-performance acoustic surface devices and accurately prediction of the Love wave propagation behaviour.     

Thin and thick diffraction gratings: Thin matrix decomposition method

A brief review of the properties of transmission diffraction gratings is presented. Two types of gratings will be analyzed: thin and volume gratings explaining how the efficiency of the different orders that propagate inside the gratings can be calculated in both cases. For thin diffraction gratings the so-called amplitude transmittance method is applied in order to get the amplitude of the different orders, whereas in the case of volume gratings more complex methods are needed, such as Coupled Wave or modal theories. We will comment on the thin matrix decomposition method (TMDM), firstly proposed by Alferness, which gives a very intuitive approach and connects the properties of thin gratings to the properties of volume ones. The thin matrix decomposition method consists in dividing the volume grating in a number of thin gratings and applying the amplitude transmittance method to each thin grating. In this way the output of a grating will be considered as the input of the next and any individual grating can be treated by the amplitude transmittance method. The novelty of this work is that a comparison is made between the analytical expressions obtained by Alferness using the TMDM with the numerical results obtained using the coupled wave (CW) and rigorous coupled wave (RCW) theories for the efficiencies of the zero, first and second order when a plane wave incides onto a sinusoidal diffracion grating at the second on-Bragg replay angular condition.     

啁啾脉冲放大激光系统中展宽器色散的解析算法

时域色散精密控制是超短激光产生及其应用中的关键技术之一,它通过控制各波长的光程产生相对延迟从而改变脉冲宽度.展宽器是啁啾脉冲放大激光系统中对激光脉冲展宽的装置,基于光线追迹法研究光线在展宽器中的传输路径,可计算飞秒脉冲中各波长的光程,进而计算脉冲展宽量并应用于系统设计.由于展宽器的光程表达式复杂,直接对其求导获得色散表达式较困难,目前只能采用数值导数获得近似解,这在计算过程中会引入误差,不利于激光系统精确设计和优化.本文介绍了一种易实现的求解展宽器色散的解析算法,通过归纳展宽器光程表达式特点,引入四个基元函数,将光程表达式分解和反复代换,可得到高阶色散的精确解析值.本文首先对Martinez型展宽器重新光线追迹,获得与Offner型展宽器一致的相位表达式,其次通过解析算法获得了两种展宽器的精确高阶色散值,最后将解析算法与数值算法的结果进行了比较.该解析算法对于啁啾脉冲放大系统的参数设计具有实用价值.     

Solution methods for discrete-state Markovian initial value problems

Solution methods, both numerical and analytical, are considered for solving the Liouville master equation associated with discrete-state Markovian initial value problems. The numerical method, basically a moment (Galerkin) method, is very general and is validated and shown to converge rapidly by comparison with an earlier reported analytical result for the ensemble-averaged transmission of photons through a purely scattering statistical rod. An application of the numerical method to a simple problem in the extended kinetic theory of gases is given. It is also shown that for a certain restricted class of problems, the master equation can be solved analytically using standard Laplace transform techniques. This solution generalizes the analytical solution for the photon transmission problem to a wider class of statistical problems.     

Free vibrations of a mono-coupled periodic system

Vibration problems of periodic systems can be analyzed efficiently by means of the transfer matrix method. The frequency equation for the whole system is shown to be obtained in terms of the eigenvalues, or their natural logarithms, which are often called “propagation constants”, of the transfer matrix for a single periodic subsystem. In case of a mono-coupled system this frequency equation may be solved graphically by using the propagation constant curve, thereby saving a great deal of computational effort. Two types of mono-coupled systems are considered as numerical examples: a spring-mass oscillating system and a continuous Timoshenko beam resting on regularly spaced knife-edge supports. Depending on whether the transfer matrix is derived by an analytical procedure or by the finite element method, the numerical solutions become either exact or approximate.     

数值模拟输运问题的一种破开算子方法

本文通过破开算子方法,把二维输运问题的控制方程破开为对流问题和扩散问题。在任意四边形网格的离散下,用特征线法解对流问题,并采用伽辽金加权余量法,从而有效地减少插值所引起的数值阻尼,提高计算精度。用有限单元法和迭代计算格式解扩散问题。由于采用了辛普生积分公式,在每个时间步长都不需要求逆矩律,节省了计算时间。算例表明,本文数值模拟结果与精确的理论解吻合较好。     

Nonlinear fluctuation,frequency and stability analyses in free vibration of circular sector oscillation systems

In the present paper, the modified homotopy perturbation method (MHPM) is employed to investigate about both nonlinear swinging oscillation and the stability of circular sector oscillation systems. The sensitivity study performed for frequency analysis of the mentioned oscillatory circular sector body shows that frequency of nonlinear oscillation depends on some specific parameters and can be optimized. Furthermore onset of the instability is dependent to angle α and initial amplitude.Comparisons made among the results of the present closed-form analytical solution and the traditional numerical iterative time integration solution confirms the accuracy and efficiency of the presented analytical solution.In contrast to the available numerical methods, the present analytical method is free from the numerical damping and the time integration accumulated errors. Moreover, in comparison with the traditional multistep numerical iterative time integration methods, a much less computational time is required for the present analytical method. Responses of the dynamical systems to some extent are affected by the natural frequencies. Results reveal that for nonlinear systems, the natural frequency is remarkably affected by the initial conditions.     

水平变化波导中的全局矩阵耦合简正波方法

提出了一种新的水平变化波导中声场的耦合简正波求解方法,该方法能够处理二维点源和线源问题,提供声场的双向解。该方法利用全局矩阵(DGM)一次性求解耦合模式的系数,消除了传播矩阵递推求解中存在的误差累积问题;此外,改善了现有模型中对距离函数的归一化方法,从而避免了泄露模式指数增长导致的数值溢出问题。本文还给出了绝对软海底理想波导中耦合矩阵的闭合表达式,并分析了单个阶梯下简正波耦合现象。此外,本文还计算了理想楔形波导中的声传播问题(ASA标准问题),并与解析解及COUPLE07计算结果进行了比较,结果表明该方法是一种稳定、精确的水平变化波导中的声场计算方法。      

Energy eigenvalues from an analytical transfer matrix method

A detailed procedure based on an analytical transfer matrix method is presented to solve bound-state problems. The derivation is strict and complete. The energy eigenvalues for an arbitrary one-dimensional potential can be obtained by the method. The anharmonic oscillator potential and the rational potential are two important examples. Checked by numerical techniques, the results for the two potentials by the present method are proven to be exact and reliable.     

一种基于奇异谱的语音激活检测方法

为了提高语音激活检测在低信噪比环境中的检测性能,提出了一种基于奇异谱的语音激活检测方法。首先用多窗口方法计算每一帧语音信号的相关矩阵;然后对相关矩阵进行奇异值分解;利用奇异值可以反映有用信号和噪声分布情况的特性,将每一帧语音信号经过加权处理后的最大奇异值与自适应阈值进行比较进行语音激活检测。该方法原理简单,易于硬件实现,通过实验仿真表明,在低信噪比环境下,和基于对数能量方法相比,本文方法也能够很好的区分语音段和非语音段,有良好的检测性能。     

Analytical methods for AMC and EBG characterizations

AMC and EBG behaviors are respectively characterized by reflection phase method and dispersion diagram. A fast and simple analytical model based on transmission line theory is proposed for each analysis. The validity of these analytical methods is verified by a comparison between circuit representations, measurements and numerical results.     

Exponential decay properties of Wannier functions and related quantities

The spatial decay properties of Wannier functions and related quantities have been investigated using analytical and numerical methods. We find that the form of the decay is a power law times an exponential, with a particular power-law exponent that is universal for each kind of quantity. In one dimension we find an exponent of -3/4 for Wannier functions, -1/2 for the density matrix and for energy matrix elements, and -1/2 or -3/2 for different constructions of nonorthonormal Wannier-like functions.     

Analysis and engineering of coupled cavity waveguides based on coupled-mode theory

The analytical expression for the transmission spectra of coupled cavity waveguides （CCWs） in photonic crystals （PCs） is derived based on the coupled-mode theory （CMT）. Parameters in the analytical expression can be extracted by simple numerical simulations. We reveal that it is the phase shift between the two adjacent PC defects that uniquely determines the flatness of the impurity bands of CCWs. In addition, it is found that the phase shift also greatly affects the bandwidth of CCWs. Thus, the engineering of the impurity bands of CCWs can be realized through the adjustment of the phase shift. Based on the theoretical results, an interesting phenomenon in which a CCW acts as a single PC defect and its impurity band possesses a Lorentz lineshape is predicted. Very good agreement between the analytical results and the numerical simulations based on transfer matrix method has been achieved.     

Tomographic Approach in Reference-Frame-Independent Measurement-Device-Independent Quantum Key Distribution

Recently, a novel reference-frame-independent measurement-device-independent quantum key distribution protocol was proposed, which can remove all detector side channels as well as tolerate unknown and slow variance of reference frame without active alignment. In this paper, we propose a new tomographic method to estimate the key rate in that protocol. We estimate the key rate using conventional method and tomographic method respectively and compare the two methods by numerical simulation. The numerical simulation results show that tomographic approach is equivalent to the conventional approach, which can be used as an alternative method.