On the forced response of waveguides using the wave and finite element method

The forced response of waveguides subjected to time harmonic loading is treated. The approach starts with the wave and finite element (WFE) method where a segment of the waveguide is modeled using traditional finite element methods. The mass and stiffness matrices of the segment are used to formulate an eigenvalue problem whose solution yields the wave properties of the waveguide. The WFE formulation is used to obtain the response of the waveguide to a convected harmonic pressure (CHP). Since the Fourier transform of the response to a general excitation is a linear combination of the responses to CHPs, the response to a general excitation can be obtained via an inverse Fourier transform process. This is evaluated analytically using contour integration and the residue theorem. Hence, the approach presented herein enables the response of a waveguide to general loading to be found by: (a) modeling a segment of the waveguide using finite element methods and post-processing it to obtain the wave characteristics, (b) using Fourier transform and contour integration to obtain the wave amplitudes and (c) using the wave amplitudes to find the response at any point in the waveguide. Numerical examples are presented.     

Calculating the forced response of two-dimensional homogeneous media using the wave and finite element method

The forced response of two-dimensional, infinite, homogenous media subjected to time harmonic loading is treated. The approach starts with the wave and the finite element (WFE) method where a small segment of a homogeneous medium is modelled using commercial or in-house finite element (FE) packages. The approach is equally applicable to periodic structures with a periodic cell being modelled. This relatively small model is then used, along with periodicity conditions, to formulate an eigenvalue problem whose solution yields the wave characteristics of the whole medium. The eigenvalue problem involves the excitation frequency and the wavenumbers (or propagation constants) in the two directions. The wave characteristics of the medium are then used to obtain the response of the medium to a convected harmonic pressure (CHP). Since the Fourier transform of a general two-dimensional excitation is a linear combination of CHPs, the response to a general excitation is a linear combination of the responses to CHPs. Thus, the response of a two-dimensional medium to a general excitation can be obtained by evaluating an inverse Fourier transform. This is a double integral, one of which is evaluated analytically using contour integration and the residue theorem. The other integral can be evaluated numerically. Hence, the approach presented herein enables the response of an infinite two-dimensional or periodic medium to an arbitrary load to be computed via (a) modelling a small segment of the medium using standard FE methods and post-processing its model to obtain the wave characteristics, (b) formulating the Fourier transform of the response to a general loading, and (c) computing the inverse of the Fourier transform semi-analytically via contour integration and the residue theorem, followed by a numerical integration to find the response at any point in the medium. Numerical examples are presented to illustrate the approach.     

Separation of components from impulses in reassigned spectrograms

Two computational methods for pruning a reassigned spectrogram to show only quasisinusoidal components, or only impulses, or both, are presented mathematically and provided with step-by-step algorithms. Both methods compute the second-order mixed partial derivative of the short-time Fourier transform phase, and rely on the conditions that components and impulses are each well-represented by reassigned spectrographic points possessing particular values of this derivative. This use of the mixed second-order derivative was introduced by Nelson [J. Acoust. Soc. Am. 110, 2575-2592 (2001)] but here our goals are to completely describe the computation of this derivative in a way that highlights the relations to the two most influential methods of computing a reassigned spectrogram, and also to demonstrate the utility of this technique for plotting spectrograms showing line components or impulses while excluding most other points. When applied to speech signals, vocal tract resonances (formants) or glottal pulsations can be effectively isolated in expanded views of the phonation process.     

用改进的传输线算法计算水平层状横向同性地层中海洋可控源电磁响应

利用二维Fourier变换与电磁场分解技术将层状横向同性地层中Maxwell方程转化成两个独立的关于横磁(TM)波和横电(TE)波的传输线方程; 借助传输线理论与叠加原理, 仅利用电流源传输线Green函数得到TM波和TE波的解, 改进传输线算法, 建立横向同性地层中频率-波数域电流源电场和磁场并矢Green函数的新算法与新的解析表达式, 提高海洋可控源电磁响应数值模拟效率. 在此基础上, 利用传输线Green函数的基本解以及边界条件, 推导出广义反射系数与振幅递推公式, 得到各个地层中传输线Green函数的解析解; 然后利用Fourier逆变换与Bessel公式将海洋可控源电磁响应表示为Sommerfeld形式的积分, 借助三次样条插值与Lommel积分公式快速计算其数值解. 通过数值模拟结果考察工作频率以及地层各向异性电阻率变化等对海洋电磁响应的影响. 关键词： 传输线法 横向同性地层 海洋可控源电磁 Sommerfeld积分     

Calculating the forced response of cylinders and cylindrical shells using the wave and finite element method

The dynamic response of circular cylinders can be obtained analytically in very few (and simple) cases. For complicated (thick or anisotropic) circular cylinders, researchers often resort to the finite element (FE) method. This can lead to large models, especially at higher frequencies, which translates into high computational costs and memory requirements. In this paper, the response of axially homogenous circular cylinders (that can be arbitrarily complex through the thickness) is obtained using the wave and finite element (WFE) method. Here, the homogeneity of the cylinder around the circumference and along the axis are exploited to post-process the FE model of a small rectangular segment of the cylinder using periodic structure theory and obtain the wave characteristics of the cylinder. The full power of FE methods can be utilised to obtain the FE model of the small segment. Then, the forced response of the cylinder is posed as an inverse Fourier transform. However, since there are an integer number of wavelengths around the circumference of a closed circular cylinder, one of the integrals in the inverse Fourier transform becomes a simple summation, whereas the other can be resolved analytically using contour integration and the residue theorem. The result is a computationally efficient technique for obtaining the response to time harmonic, arbitrarily distributed loads of axially homogenous, circular cylinders with arbitrary complexity across the thickness.     

A subfringe integration method for multiple-beam Fizeau fringe analysis

Some aspects concerning the subfringe integration method in interferogram analysis have been investigated and modified. The modified algorithm, introduced in this paper, is capable of reconstructing the phase in the presence of noise or errors in carrier frequency. The subfringe integration method was applied to analyze two computer simulated patterns of equispaced Fizeau fringes using N bucket integration. Also, it is used to analyze the multiple-beam Fizeau fringe. The refractive index profile of polyethylene fiber is obtained by using two methods, subfringe integration method, and Fourier transform method. A comparison between the obtained results using the maintained methods is presented.     

The measurement of the free field impulse response of microphones in a laboratory environment

A method is presented of measuring the free field frequency and impulse response of microphones by using a pulse technique in an ordinary laboratory environment. The pressure transient generated by exciting a loudspeaker with a narrow pulse is detected at some point in the loudspeaker's far field by a “reference” microphone whose response is assumed flat over the frequency range of interest. The microphone to be tested is then substituted in exactly the same position and its response to the transient measured. The outputs of the two microphones are accurately sampled and deconvolved by using a discrete Fourier transform technique to give the magnitude and phase parts of the “test” microphone's frequency response, and hence, via the discrete Fourier transform, its response to a delta function of pressure propagating in the free field. The computed impulse responses are presented and the accuracy of the method discussed.To illustrate the use of the method, the free field frequency response and free field correction curves of a one inch instrumentation microphone are measured. The method is then used to measure the pressure which occurs at the centre of the flat end-face of a long a cylinder when excited by an impulse of acoustic pressure propagating in the free field from various angles of incidence.     

Fourier analysis for rotating-element ellipsometers

We introduce a Fourier analysis of the waveform of periodic light-irradiance variation to capture Fourier coefficients for multichannel rotating-element ellipsometers. In this analysis, the Fourier coefficients for a sample are obtained using a discrete Fourier transform on the exposures. The analysis gives a generic function that encompasses the discrete Fourier transform or the Hadamard transform, depending on the specific conditions. Unlike the Hadamard transform, a well-known data acquisition method that is used only for conventional multichannel rotating-element ellipsometers with line arrays with specific readout-mode timing, this Fourier analysis is applicable to various line arrays with either nonoverlap or overlap readout-mode timing. To assess the effects of the novel Fourier analysis, the Fourier coefficients for a sample were measured with a custom-built rotating-polarizer ellipsometer, using this Fourier analysis with various numbers of scans, integration times, and rotational speeds of the polarizer.     

双重分数傅里叶变换计算全息

提出了双重分数傅里叶变换计算全息，在这种方法中，将两个图像的信息分别经不同阶的分数傅里叶变换后，记录在同一张分数傅里叶变换计算全息图上，它需要两个特定的分数傅里叶变换系统才能再现出所记录的图像信息，利用其再现方式的特殊性，可制成一种新的安全认证系统。     

Propagation characteristics of Bessel-like beams through ABCD optical system

Based on the Collins diffraction integral formula, this work focus to the study of the propagation properties of a Bessel-like beam through an ABCD optical system. As it is not easy to develop an analytical expression describing the transformation of the Bessel-like beams by an ABCD optical system, numerical integration method was indispensable to conduct this research. The characteristics of the transformation of Bessel-like beams by three optical systems: free space, Fourier transform, Fractional Fourier transform system and quadratic graded-index system are numerically examined and treated as examples of optical systems in this investigation. Some numerical results are carried-out to show how the superposition of multiple Airy beams propagates and transforms by an ABCD optical system.     

X-ray interferometric Fourier holography

An X-ray interferometric Fourier holography was proposed and theoretically investigated. X-ray Interferometric Young fringes and the reconstruction of an object image were investigated by the Fourier transform method. It was shown that on the output surface of the analyzer crystal (the third plate of the interferometer) the interference pattern of two slits gives X-ray interferometric Young fringes. An expression for the period of X-ray interferometric Young fringes was obtained. The subsequent reconstruction of the slit image as an object is conducted by means of the Fourier transform method of intensity distribution on the hologram. Three methods for reconstruction of the function of complex transmission of the object are presented: an analytical one–the approximate method, the iteration method and the step-by-step approach. As examples a recording of X-ray interferometric Fourier hologram and the reconstruction of the function of complex amplitude transmission of a beryllium circular cylinder are given.     

Frequency-domain Hadamard spectroscopy

A new technique is proposed for multichannel excitation and detection of NMR signals in the frequency domain, an alternative to the widely used pulse-excited Fourier transform method. An extensive array of N radiofrequency irradiation channels covers the spectrum of interest. A selective radiofrequency pulse sequence is applied to each channel, generating a steady-state NMR response acquired one-point-at-a-time in the intervals between pulses. The excitation pattern is repeated N times, phase-encoded according to a Hadamard matrix, and the corresponding N composite responses are decoded by reference to the same matrix. This multiplex technique offers the same sensitivity advantage as conventional Fourier transform spectroscopy. The irradiation pattern may be tailored to concentrate on interesting spectral regions, to facilitate homonuclear double resonance, or to avoid exciting strong solvent peaks. As no free induction decay is involved, the new method avoids problems of pulse breakthrough or lineshape distortion by premature termination of the time-domain signal.     

Implementing compressive fractional Fourier transformation with iterative kernel steering regression in double random phase encoding

This paper proposes a novel approach in double random phase encryption based on compressive fractional Fourier transform along with the kernel steering regression. The method increases the complexity of the image by using fractional Fourier transform and taking fewer measurements from the image data. Numerical results are given to analyze the validity of this technique. Considering natural images to be sparse in some domain, we apply a compressive sensing (CS) approach by using a TwIST algorithm. The encryption process has kernel steering regression algorithm for denoising and compressive sensing technique for image compression along with the fractional Fourier transform that makes the image in more complex form.     

Optical image encryption by random shifting in fractional Fourier domains

A number of methods have recently been proposed in the literature for the encryption of two-dimensional information by use of optical systems based on the fractional Fourier transform. Typically, these methods require random phase screen keys for decrypting the data, which must be stored at the receiver and must be carefully aligned with the received encrypted data. A new technique based on a random shifting, or jigsaw, algorithm is proposed. This method does not require the use of phase keys. The image is encrypted by juxtaposition of sections of the image in fractional Fourier domains. The new method has been compared with existing methods and shows comparable or superior robustness to blind decryption. Optical implementation is discussed, and the sensitivity of the various encryption keys to blind decryption is examined.     

System identification from multiple input/output data

For stationary random or transient data representing multicorrelated (multicoherent) input/output data occurring in physical systems, iterative computational algorithms are developed to identify the frequency response functions of optimum constant parameter linear systems connecting this data. Results are obtained from Fourier transform methods and optimum least-squares prediction techniques by changing original arbitrary input records into ordered sets of conditioned input records. These procedures provide the basis for efficient digital computer analysis of general multiple input/output problems. The Appendix contains useful error analysis results for the optimum frequency response estimates determined from measured data.     

The performance of methods based on the fractional Fourier transform for detecting marine mammal vocalizations

The analysis of cetacean vocalizations is considered using Fourier-based techniques that employ chirp functions in their decomposition. In particular, the paper considers a short-time methods based on the fractional Fourier transform for detecting frequency modulated narrow-band signals, such as dolphin whistles, and compares this to the classical short-time Fourier methods. The fractional Fourier technique explored computes transforms associated with a range of chirp rates and automatically selects the rate for the final analysis. This avoids the need for prior knowledge of signal's chirp rate. An analysis is presented that details the performance of both methods as signal detectors and allows one to determine their detection thresholds. These thresholds are then used to measure the detectability of synthetic signals. This principle is then extended to measure performance on a set of recordings of narrow-band vocalizations from a range of cetacean species.     

Pulse sequences for realizing the quantum Fourier transform on multilevel systems

Sequences of selective pulses of an RF magnetic field for realizing the quantum Fourier transform by NMR methods on systems with four, six, and eight nonequidistant levels are found using the virtual spin formalism. The results can be applied to other quantum systems when laser pulses are used.     

Study on the simplified phase-shifting digital holographic microscopy

In this paper, a new simplified technique for effectively eliminating the zero order and the conjugate virtual image in digital holographic microcopy, which makes use of two-step phase-shifting method of just recording two holograms and an intensity image of object wave, is proposed. Meanwhile, combined with the principle of making full use of spatial bandwidth of the CCD sensor by in-line lens-less Fourier holographic recording geometry, the theory and experimental methods to increase the resolution of the reconstructed image in digital holography by using phase-shifting technique are detailedly analyzed. At end, the validity and availability of this technique has been demonstrated through the off-axis and in-line Fourier transform recording geometry. The study provides some theoretical and experimental guidance for the design and operation of a digital holographic microscopy system.     

Estimation of generalized frequency response functions for quadratically and cubically nonlinear systems

A new method of estimating the Volterra kernels in the frequency-domain is introduced based on a non-parametric algorithm. Unlike the traditional non-parametric methods using the frequency-domain formulations based on discrete Fourier transform (DFT) data, this new approach uses the time-domain measurements directly to estimate the frequency-domain response functions.