Epoch extraction from linear prediction residual for identification of closed glottis interval

In voiced speech analysis epochal information is useful in accurate estimation of pitch periods and the frequency response of the vocal tract system. Ideally, linear prediction (LP) residual should give impulses at epochs. However, there are often ambiguities in the direct use of LP residual since samples of either polarity occur around epochs. Further, since the digital inverse filter does not compensate the phase response of the vocal tract system exactly, there is an uncertainty in the estimated epoch position. In this paper we present an interpretation of LP residual by considering the effect of the following factors: 1) the shape of glottal pulses, 2) inaccurate estimation of formants and bandwidths, 3) phase angles of formants at the instants of excitation, and 4) zeros in the vocal tract system. A method for the unambiguous identification of epochs from LP residual is then presented. The accuracy of the method is tested by comparing the results with the epochs obtained from the estimated glottal pulse shapes for several vowel segments. The method is used to identify the closed glottis interval for the estimation of the true frequency response of the vocal tract system.

An improved iterative finite element solution for pin joints

Accurate, reliable and economical methods of determining stress distributions are important for fastener joints. In the past the contact stress problems in these mechanically fastened joints using interference or push or clearance fit pins were solved using both inverse and iterative techniques. Inverse techniques were found to be most efficient, but at times inadequate in the presence of asymmetries. Iterative techniques based on the finite element method of analysis have wider applications, but they have the major drawbacks of being expensive and time-consuming. In this paper an improved finite element technique for iteration is presented to overcome these drawbacks. The improved iterative technique employs a frontal solver for elimination of variables not requiring iteration, by creation of a dummy element. This automatically results in a large reduction in computer time and in the size of the problem to be handled during iteration. Numerical results are compared with those available in the literature. The method is used to study an eccentrically located pin in a quasi-isotropic laminated plate under uniform tension.

Dynamical Resolution of Redundancy for Robot Manipulators

A large class of work in the robot manipulator literature deals with the kinematical resolution of redundancy based on the pseudo-inverse of the manipulator Jacobian. In this paper an alternative dynamical approach to redundancy resolution is developed which utilizes the mapping between the actuator torques and the acceleration of the end-effector, at a given dynamic state of the manipulator. The potential advantages of the approach are discussed and an example of a planar 3R manipulator following a circular end-effector trajectory is used to illustrate the proposed approach as well as to compare it with the more well-known approach based on the pseudo-inverse.

ECG component delineation by Prony's method

A simple, non-iterative method for component wave delineation from the electrocardiogram (ECG) is derived by modelling its discrete cosine transform (DCT) as a sum of damped cosinusoids. Amplitude, phase, damping factor and frequency parameters of each of the cosinusoids are estimated by the extended Prony method. Different component waves are represented by non-overlapping clusters of model poles in the z plane and thus a component wave is derived by the addition of the inverse transformed (IDCT) impulse responses of the poles in the cluster. Akaike's information criterion (AIC) is used to determine the model order. The method performed satisfactory even in the presence of artifacts. The efficacy of the method is illustrated by analysis of continuous strips of ECG data.

Effect of microstructure and hardness on the grinding abrasive wear resistance of a ball bearing steel

Grinding media wear appears to be non-linear with the time of grinding in a laboratory-scale ball mill. The kinetics of wear can be expressed as a power law of the type w=atb, where the numerical constant a represents wear of a particular microstructure at time t = 1 min and b is the wear exponent which is independent of the particle size prevailing inside a ball mill at any instant of time of grinding. The wear exponent appears to be an indicator of the cutting wear mechanism in dry grinding: a plot of the inverse of the normalised wear exponent (Image ) versusHs (where Hs is the worn surface hardness of the media) yields a curve similar to that of a wear resistance plot obtained in the case of two-body sliding abrasive wear. This method of evaluating the cutting wear resistance of media is demonstrated by employing 15 different microstructures of AISI-SAE 52100 steel balls in dry grinding of quartz in a laboratory-scale ball mill.

1H NMR Study of Molecular Motions and Phase Transitions in Methyl Ammonium Hexabromo Selenate [(CH3NH3)2SeBr6]

The temperature dependence of 1H spin-lattice relaxation time, T1, and that of the second moment, M2, are analysed in the temperature range 390 K to 77 K. A plot of T1 vs inverse temperature shows three phase transitions at 250 K, 167 K and 111 K. At 167 K, T1 displays a large jump while it shows changes in slope at 250 K and 111 K. In the high temperature phase (> 167 K), the correlated motion of CH3 and NH3 groups is found to cause the relaxation while their uncorrelated motion takes over in the low temperature phases (< 167 K). The unusual T1 behaviour in phase II (250 K-167 K) is ascribed to the small angle torsion of the cation. A constant M2 value of ∼ 9.7 G2, throughout the range of temperature studied, indicates the presence of reorientation of CH3 and NH3 groups.

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Synthesis of substituted 8-aminoquinolines and phenanthrolines through a Povarov approach

The synthesis of 8-aminoquinolines and 1,10-phenanthrolines with substituents in alpha of the nitrogen has been performed through an inverse-demanding aza-Diels-Alder (Povarov reaction) in the fluoroalcohols TFE or HFIP. This path involves simple starting materials: 1,2-phenylenediamines, enol ethers and aldehydes.

Interior Photon Absorption Based Adaptive Regularization Improves Diffuse Optical Tomography

An adaptive regularization algorithm that combines elementwise photon absorption and data misfit is proposed to stabilize the non-linear ill-posed inverse problem. The diffuse photon distribution is low near the target compared to the normal region. A Hessian is proposed based on light and tissue interaction, and is estimated using adjoint method by distributing the sources inside the discretized domain. As iteration progresses, the photon absorption near the inhomogeneity becomes high and carries more weightage to the regularization matrix. The domain's interior photon absorption and misfit based adaptive regularization method improves quality of the reconstructed Diffuse Optical Tomographic images.

Gain calibration of CCD systems at VBO

The system gain of two CCD systems in regular use at the Vainu Bappu Observatory, Kavalur, is determined at a few gain settings. The procedure used for the determination of system gain and base-level noise is described in detail. The Photometrics CCD system at the 1-m reflector uses a Thomson-CSF TH 7882 CDA chip coated for increased ultraviolet sensitivity. The gain is programme-selected through the parameter 'cgain' varying between 0 and 4095 in steps of 1. The inverse system gain for this system varies almost linearly from 27.7 electrons DN-1 at cgain = 0 to 1.5 electrons DN-1 at cgain = 500. The readout noise is less than or similar 11 electrons at cgain = 66. The Astromed CCD system at 2.3-m Vainu Bappu Telescope uses a GEC P8603 chip which is also coated for enhanced ultraviolet sensitivity. The amplifier gain is selected in discrete steps using switches in the controller. The inverse system gain is 4.15 electrons DN-1 at the gain setting of 9.2, and the readout noise approximately 8 electrons.

Accelerated gradient based diffuse optical tomographic image reconstruction

Purpose: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). Methods: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. Results: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. Conclusions: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]

Concise row-pruning algorithm to invert a matrix

Presented here, in a vector formulation, is an O(mn2) direct concise algorithm that prunes/identifies the linearly dependent (ld) rows of an arbitrary m X n matrix A and computes its reflexive type minimum norm inverse A(mr)-, which will be the true inverse A-1 if A is nonsingular and the Moore-Penrose inverse A+ if A is full row-rank. The algorithm, without any additional computation, produces the projection operator P = (I - A(mr)- A) that provides a means to compute any of the solutions of the consistent linear equation Ax = b since the general solution may be expressed as x = A(mr)+b + Pz, where z is an arbitrary vector. The rank r of A will also be produced in the process. Some of the salient features of this algorithm are that (i) the algorithm is concise, (ii) the minimum norm least squares solution for consistent/inconsistent equations is readily computable when A is full row-rank (else, a minimum norm solution for consistent equations is obtainable), (iii) the algorithm identifies ld rows, if any, and reduces concerned computation and improves accuracy of the result, (iv) error-bounds for the inverse as well as the solution x for Ax = b are readily computable, (v) error-free computation of the inverse, solution vector, rank, and projection operator and its inherent parallel implementation are straightforward, (vi) it is suitable for vector (pipeline) machines, and (vii) the inverse produced by the algorithm can be used to solve under-/overdetermined linear systems.

Cracks emanating from pin-loaded lugs

Pin-loaded lugs were analysed in the presence of cracks emanating from circular holes. The analysis presents a unified treatment of interference, push or clearance fit pins. Both metallic (isotropic) and composite (orthotropic) plates were dealt with. The finite element model used special singular six-noded quadrilateral elements at the crack tip. The non-linear load contact behaviour at the pin-hole interface was dealt with by an inverse technique. A modified crack closure integral (MCCI) technique was used to evaluate the strain energy release rates (SERRs) and stress intensity factors (SIFs) at the crack tips. Numerical results are presented showing the non-linear variation of SIF with applied stress, and the influence of the amount of interference or clearance and the interfacial friction on SIF.

Diffusion of sorbates in zeolites Y and A novel dependence on sorbate size and strength of sor ba t e-Zeoli t e interaction

A systematic investigation of monatomic spherical sorbates in the supercages of zeolites Y and A by molecular dynamics technique is presented. Rates of intercage diffusion, rates of cage visits, and the diffusion coefficients have been calculated as a function of the sorbate-zeolite interaction strength. These properties exhibit markedly different dependences on interaction strength for the two zeolites. The observed behavior is shown to be a consequence of the two principal mechanisms of intercage diffusion and the energetic barrier associated with them. The diffusion coefficient and other properties associated with intercage diffusion are found to be directly proportional to the reciprocal of the square of the sorbate diameter when the sorbate diameter is significantly smaller than the window diameter. As the sorbate diameter increases, a peak is observed in all the transport properties investigated including the diffusion coefficient. We call this surprising effect as the ring or levitation effect and it explains several anomalous results reported in the literature and suggests a breakdown of the geometrical criterion for diffusion of sorbates. It shows that under certain conditions nongeometrical factors play a major role and geometrical factors become secondary in the determination of the molecular sieve property. A generalized parameter has been proposed which suggests conditions under which one can expect the ring or levitation effect in any porous medium. Inverse size selectivity becomes operative under these conditions.

Effect of surface roughness on diffusion-limited charge transfer

A theory is developed for diffusion-limited charge transfer on a non-fractally rough electrode. The perturbation expressions are obtained for concentration, current density and measured diffusion-limited current for arbitrary one- and two-dimensional surface profiles. The random surface model is employed for a rough electrode\electrolyte interface. In this model the gross geometrical property of an electrochemically active rough surface - the surface structure factor-is related to the average electrode current, current density and concentration. Under short and long time regimes, various morphological features of the rough electrodes, i.e. excess area (related to roughness slope), curvature, correlation length, etc. are related to the (average) current transients. A two-point Pade approximant is used to develop an all time average current expression in terms of partial morphological features of the rough surface. The inverse problem of predicting the surface structure factor from the observed transients is also described. Finally, the effect of surface roughness is studied for specific surface statistics, namely a Gaussian correlation function. It is shown how the surface roughness enhances the overall diffusion-limited charge transfer current.