Derivation of closed-form Green's functions for a generalmicrostrip geometry

The derivation of the closed-form spatial domain Green's functions for the vector and scalar potentials is presented for a microstrip geometry with a substrate and a superstrate, whose thicknesses can be arbitrary. The spatial domain Green's functions for printed circuits are typically expressed as Sommerfeld integrals, which are inverse Hankel transforms of the corresponding spectral domain Green's functions and are time-consuming to evaluate. Closed-form representations of these Green's functions in the spatial domains can only be obtained if the integrands are approximated by a linear combination of functions that are analytically integrable. This is accomplished here by approximating the spectral domain Green's functions in terms of complex exponentials by using the least square Prony's method     

Spherical wave operators and the translation formulas

Translational formulas for both scalar and vector spherical wave solutions of the Helmholtz equation are developed in a straightforward manner using differential operator representations for the modal functions and well-known expressions for the scalar and dyadic free-space Green's functions. The expansion coefficients are given in compact integral or differential operator forms useful for analytic investigation     

Capacitance computations in a multilayered dielectric medium usingclosed-form spatial Green's functions

An efficient method to compute the 2-D and 3-D capacitance matrices of multiconductor interconnects in a multilayered dielectric medium is presented. The method is based on an integral equation approach and assumes the quasi-static condition. It is applicable to conductors of arbitrary polygonal shape embedded in a multilayered dielectric medium with possible ground planes on the top or bottom of the dielectric layers. The computation time required to evaluate the space-domain Green's function for the multilayered medium, which involves an infinite summation, has been greatly reduced by obtaining a closed-form expression, which is derived by approximating the Green's function using a finite number of images in the spectral domain. Then the corresponding space-domain Green's functions are obtained using the proper closed-form integrations. In both 2-D and 3-D cases, the unknown surface charge density is represented by pulse basis functions, and the delta testing function (point matching) is used to solve the integral equation. The elements of the resulting matrix are computed using the closed-form formulation, avoiding any numerical integration. The presented method is compared with other published results and showed good agreement. Finally, the equivalent microstrip crossover capacitance is computed to illustrate the use of a combination of 2-D and 3-D Green's functions     

Multiresolution maximum intensity volume rendering by morphological adjunction pyramids

We describe a multiresolution extension to maximum intensity projection (MIP) volume rendering, allowing progressive refinement and perfect reconstruction. The method makes use of morphological adjunction pyramids. The pyramidal analysis and synthesis operators are composed of morphological 3-D erosion and dilation, combined with dyadic downsampling for analysis and dyadic upsampling for synthesis. In this case the MIP operator can be interchanged with the synthesis operator. This fact is the key to an efficient multiresolution MIP algorithm, because it allows the computation of the maxima along the line of sight on a coarse level, before applying a two-dimensional synthesis operator to perform reconstruction of the projection image to a finer level. For interpolation and resampling of volume data, which is required to deal with arbitrary view directions, morphological sampling is used, an interpolation method well adapted to the nonlinear character of MIP. The structure of the resulting multiresolution rendering algorithm is very similar to wavelet splatting, the main differences being that (i) linear summation of voxel values is replaced by maximum computation, and (ii) linear wavelet filters are replaced by nonlinear morphological filters.     

High-quality image resizing using oblique projection operators

The standard interpolation approach to image resizing is to fit the original picture with a continuous model and resample the function at the desired rate. However, one can obtain more accurate results if one applies a filter prior to sampling, a fact well known from sampling theory. The optimal solution corresponds to an orthogonal projection onto the underlying continuous signal space. Unfortunately, the optimal projection prefilter is difficult to implement when sine or high order spline functions are used. We propose to resize the image using an oblique rather than an orthogonal projection operator in order to make use of faster, simpler, and more general algorithms. We show that we can achieve almost the same result as with the orthogonal projection provided that we use the same approximation space. The main advantage is that it becomes perfectly feasible to use higher order models (e.g. splines of degree n=/>3). We develop the theoretical background and present a simple and practical implementation procedure using B-splines. Our experiments show that the proposed algorithm consistently outperforms the standard interpolation methods and that it provides essentially the same performance as the optimal procedure (least squares solution) with considerably fewer computations. The method works for arbitrary scaling factors and is applicable to both image enlargement and reduction.     

Unified study on characteristic spectrum representation of signals and spectral decomposition for stationary random signals

The unified relationship between the signal characteristic spectrum representation and the spectral decomposition for the stationary random signal was deeply studied. By using the relations among the differential operators, the integral operator and the Green's function of the characteristic differential equation, the inverse relationship between the Hermitian differential operator and the Hermitian integral operator were given, the characteristic differential equation and corresponding characteristic integral equation were demonstrated, and the spectral representations of both Hermitian differential and integral operators and the general spectral representations for both operators were provided. Based on the superposition method of the stochastic simple harmonic vibration and the Hilbert space unitary operator method for the stationary random signal spectral decomposition, the connection and unification on mathematics of the signal characteristic spectral representation and the stationary random signal spectral decomposition are revealed.     

Novel closed-form Green's function in shielded planar layered media

A new method is proposed for the construction of closed-form Green's function in planar, stratified media between two conducting planes. The new approach does not require the a priori extraction of the guided-wave poles and the quasi-static part from the Green function spectrum. The proposed methodology can be easily applied to arbitrary planar media without any restriction on the number of layers and their thickness. Based on the discrete solution of one-dimensional ordinary differential equations for the spectral-domain expressions of the appropriate vector potential components, the proposed method leads to the simultaneous extraction of all Green's function values associated with a given set of source and observation points. Krylov subspace model order reduction is used to express the generated closed-form Green's function representation in terms of a finite sum involving a small number of Hankel functions. The validity of the proposed methodology and the accuracy of the generated closed-form Green's functions are demonstrated through a series of numerical experiments involving both vertical and horizontal dipoles     

A multidimensional Z-transform evaluation of the discrete finite difference time domain Green's function

Analytical expressions for frequency and time domain discrete Green's functions are developed below. Rather than discretizing the continuous Green's functions directly, we derive the discrete Green's functions from first principles, i.e., the difference equations. This is done with the purpose of obtaining discrete integral operators that would replicate results obtained via the finite difference time domain (FDTD) method. The main effort is directed at inversions of multidimensional frequency domain expressions in the Z-transform domain. The Green's functions obtained in this way coincide with time domain developments presented earlier both in the form of direct numerical computations and via combinatorics and the Catalan triangle. The Z-transform methodology, resulting from this development, provides an orderly manner for deriving expressions for multidimensional discrete integral operators that can be hybridized with FDTD-based differential operators in a self consistent manner.     

A new technique for the derivation of closed-form electromagnetic Green's functions for unbounded planar layered media

A closed-form electromagnetic Green's function for unbounded, planar, layered media is derived in terms of a finite sum of Hankel functions. The derivation is based on the direct inverse Hankel transform of a pole-residue representation of the spectral-domain form of the Green's function. Such a pole-residue form is obtained through the solution of the spectral-domain form of the governing Green's function equation numerically, through a finite-difference approximation, rather than analytically. The proposed methodology can handle any number of layers, including the general case where the planar media exhibit arbitrary variation in their electrical properties in the vertical direction. The numerical implementation of the proposed methodology is straightforward and robust, and does not require any preprocessing of the spectrum of the Green's function for the extraction of surface-wave poles or its quasi-static part. The number of terms in the derived closed-form expression is chosen adaptively with the distance between source and observation point as parameter. The development of the closed-form Green's function is presented for both vertical and horizontal dipoles. Its accuracy is verified through a series of numerical examples and comparisons with results from other established methods.     

Chebyshev series for designing RF pulses employing an optimal control approach

Magnetic resonance imaging (MRI) provides bidimensional images with high definition and selectivity. Selective excitations are achieved applying a gradient and a radio frequency (RF) pulse simultaneously. They are modeled by the Bloch differential equation, which has no closed-form solution. Most methods for designing RF pulses are derived from approximation of this equation or are based on iterative optimization methods. The approximation methods are only valid for small tip angles and the optimization-based algorithms yield better results, but they are computationally intensive. To improve the solutions and to reduce processing time, a method for designing RF pulses using a pseudospectral approach is presented. The Bloch equation is expanded in Chebyshev series, which can be solved using a sparse linear algebraic system. The method permits three different formulations derived from the optimal control theory, minimum distance, minimum energy, or minimum time, which are solved as algebraic constrained minimization problems. The results were validated through simulated and real experiments of 90 degrees and 180 degrees RF pulses. They show improvements compared to the corresponding solutions obtained using the Shinnar-Le Roux method. The minimum time formulation produces the best performance for 180 degrees pulses, reducing the excitation length in 4% and the RF pulse energy in 3%.     

Extended-window interpolation applied to digital filter design

This paper develops a procedure for the design of frequency-selective interpolation operators that can be computed and saved once and for all. These operators are used to design real-time digital operators: interpolators, FIR differentiators, IIR filters, and composed interpolation and filtering operators. Each real-time operator is a matrix relating sets of data points to sets of interpolated values. Since these matrices are characterized by low norms, they permit reduced-word implementations, and are suitable for real-time processing with array processors and massively parallel machines. The design of the interpolation operators uses windows that, unlike traditional approaches, extend beyond the data interval up to the length permitted by the dimensionality theorem. A new form of the dimensionality theorem is used to minimize the minimax interpolation error within a predetermined frequency range, which may be either the passband of the antialiasing filter or the passband of an analog prototype filter. The main application presented in the paper is the design of combined digital filters and interpolators, which will be referred to as interpolating filters. The frequency responses of such filters, as well as the interpolated time responses, almost coincide with those of the corresponding analog prototypes     

Efficient electromagnetic modeling of three-dimensional multilayermicrostrip antennas and circuits

An efficient method-of-moments (MoM) solution is presented for analysis of multilayer microstrip antennas and circuits. The required multilayer Green's functions are evaluated by the discrete complex image method (DCIM), with the guided-mode contribution extracted recursively using a multilevel contour integral in the complex _ρ-plane. An interpolation scheme is employed to further reduce the computer time for calculating the Green's functions in the three-dimensional (3-D) space. Higher order interpolatory basis functions defined on curvilinear triangular patches are used to provide necessary flexibility and accuracy for the discretization of arbitrary shapes and to offer a better convergence than lower order basis functions. The combination of the improved DCIM and the higher order basis functions results in an efficient and accurate MoM analysis for 3-D multilayer microstrip structures     

Concentration maximization and local basis expansions (LBEX) for linear inverse problems

Linear inverse problems arise in biomedicine electroencephalography and magnetoencephalography (EEG and MEG) and geophysics. The kernels relating sensors to the unknown sources are Green's functions of some partial differential equation. This knowledge is obscured when treating the discretized kernels simply as matrices. Consequently, physical understanding of the fundamental resolution limits has been lacking. We relate the inverse problem to spatial Fourier analysis, and the resolution limits to uncertainty principles, providing conceptual links to underlying physics. Motivated by the spectral concentration problem and multitaper spectral analysis, our approach constructs local basis sets using maximally concentrated linear combinations of the measurement kernels.     

Consistent Sampling and Signal Recovery

An attractive formulation of the sampling problem is based on the principle of a consistent signal reconstruction. The requirement is that the reconstructed signal is indistinguishable from the input in the sense that it yields the exact same measurements. Such a system can be interpreted as an oblique projection onto a given reconstruction space. The standard formulation requires a one-to-one relationship between the input measurements and the reconstructed model. Unfortunately, this condition fails when the cross-correlation matrix between the analysis and reconstruction basis functions is not invertible; in particular, when there are less measurements than the number of reconstruction functions. In this paper, we propose an extension of consistent sampling that is applicable to those singular cases as well, and that yields a unique and well-defined solution. This solution also makes use of projection operators and has a geometric interpretation. The key idea is to exclude the null space of the sampling operator from the reconstruction space and to enforce consistency on its complement. We specify a class of consistent reconstruction algorithms corresponding to different choices of complementary reconstruction spaces. The formulation includes the Moore-Penrose generalized inverse, as well as other potentially more interesting reconstructions that preserve certain preferential signals. In particular, we display solutions that preserve polynomials or sinusoids, and therefore perform well in practical applications.     

Optimal nonnegative definite approximations of estimated movingaverage covariance sequences

The problem of finding the closest nonnegative definite moving average covariance sequence to a given estimate which may not be nonnegative definite is considered. An algorithm is developed which is based on a set of constrained minimization problems, each parameterized by the zero frequencies of the spectral density function corresponding to the optimal solution. The algorithm entails first solving a simple minimization problem with linear constraints whose closed-form solution is given by a projection onto a subspace. These solutions lie either outside the set of nonnegative definite sequences, or on its boundary; if the solution lies on the boundary, it is the optimal solution. The problem is considered directly in the space of covariance sequence elements. As a result, the nonlinear maximization step is performed on sets of low dimension. By considering the minimization problem in this space, it is possible to characterize some of the geometrical properties of the optimal solution in terms of the locations of its zero frequencies     

微分求积法求解高速大规模集成电路互连线的瞬态响应

本文将微分求积法（ＤＱ方法）应用于高速大规模集成电路互连线的瞬态模拟。ＤＱ方法是一种直接的数值方法，与差分和有限元法相比，它的计算量可以大大降低，具有较高的精度。ＤＱ方法的主要思想是将某坐标方向上的微分算子用该方向上一系列适当的离散点的函数值加仅逼近，将偏微分方程化为常微分方程或代数方程求解。ＤＱ方法用于高速大规模集成电路互连线系统的瞬态模拟非常有效，其适用范围也相当广泛。     

Finite-Element Solution of Unbounded Field Problems

An unbounded region is divided into local picture-frame regions where a partial differential-equation solution is obtained, with the remaining unbounded region represented by an integral equation. (The method permits the use of free-space Green's functions, and thus special problem-dependent Green's functions need not be found.) The integral equation is formulated as a constraint upon the local picture-frame solutions, whence these local solutions are solved directly by a variational method, using finite elements, in a manner such that the problem of the Green's-function singularity is side-stepped. The technique is applicable where sources and media inhomogeneities and anisotropies are local, and can all be placed within one or several picture frames. It is in these cases that the integral-equation approach is at a particular disadvantage, and the use of a partial differential-equation technique is advisable if not necessary. Examples presented include the static and harmonic fields of a parallel-plate capacitor, a microstrip line on a dielectric substratum, and a radiating antenna with dielectric obstacles.     

New improved fractional order integrators using PSO optimisation

This article investigates the optimal results of new improved fractional order integrators (FOIs) of different orders. Mathematical models of FOIs have been first developed by a single-step procedure of direct linear interpolation of fractional integrators based on Al-Alaoui operator in fractional domain itself, instead of using three steps of the well-known conventional method, namely, digital interpolation, series expansion and truncation. Later, these transfer functions (TFs) are optimised for their coefficient values for finding a minimum error function by particle swarm optimisation (PSO) algorithm. Simulation results of magnitude responses, phase responses and relative magnitude errors (dB) for all the proposed half integrators have validated the effectiveness of this new technique of interpolation of fractional order operators, mixed with PSO algorithm. A parallel comparison has been also drawn between the proposed optimised half integrators and those obtained by discretisation of PSO optimised integer order digital integrators (DIs) to properly support the proposed novel combination of interpolation and PSO, both applied together in fractional domain.     

用电磁场算子理论分析TEM室的传输特性

本文用电磁场算子理论通过多个虚拟边界的电场和磁场的耦合求解TEM室(TEM cell)的本征值,在此基础上讨论了TEM室的的传输特性,文中所采用的并矢格林函数没有奇异项,可以化为标量格林函数来计算,并与数值计算的结果进行了比较.     

An Efficient Application of the Discrete Complex Image Method for Quasi-3-D Microwave Circuits in Layered Media

An efficient means of evaluating reactions arising from the mixed-potential integral equation for layered media for quasi-3D microwave circuits is presented. Analytical formulations for the z-integration in the spectral domain are derived, thus avoiding expensive 2D evaluation and interpolation of the layered Green's function. In this paper, closed-form formulations for the resultant Sommerfeld integrals are evaluated via the robust two-step discrete complex image method. The overall computational time can consequently be greatly reduced when analyzing quasi-3D circuits in layered media using the proposed method.