An introduction to radial basis functions for system identification: a comparison with other neural network methods

A look is taken at the use of radial basis functions (RBFs), for nonlinear system identification. RBFs are firstly considered in detail themselves and are subsequently compared with a multi-layered perceptron (MLP), in terms of performance and usage.

Necessary and sufficient conditions for realizability of point processes

We give necessary and sufficient conditions for a pair of (generali- zed) functions 1(r1) and 2(r1, r2), ri 2X, to be the density and pair correlations of some point process in a topological space X, for ex- ample, Rd, Zd or a subset of these. This is an infinite-dimensional version of the classical “truncated moment” problem. Standard tech- niques apply in the case in which there can be only a bounded num- ber of points in any compact subset of X. Without this restriction we obtain, for compact X, strengthened conditions which are necessary and sufficient for the existence of a process satisfying a further re- quirement—the existence of a finite third order moment. We general- ize the latter conditions in two distinct ways when X is not compact.

Belowground biomass functions and expansion factors in high elevation Norway spruce

Biomass allocation to above- and belowground compartments in trees is thought to be affected by growth conditions. To assess the strength of such influences, we sampled six Norway spruce forest stands growing at higher altitudes. Within these stands, we randomly selected a total of 77 Norway spruce trees and measured volume and biomass of stem, above- and belowground stump and all roots over 0.5 cm diameter. A comparison of our observations with models parameterised for lower altitudes shows that models developed for specific conditions may be applicable to other locations. Using our observations, we developed biomass functions (BF) and biomass conversion and expansion factors (BCEF) linking belowground biomass to stem parameters. While both BF and BCEF are accurate in belowground biomass predictions, using BCEF appears more promising as such factors can be readily used with existing forest inventory data to obtain estimates of belowground biomass stock. As an example, we show how BF and BCEF developed for individual trees can be used to estimate belowground biomass at the stand level. In combination with existing aboveground models, our observations can be used to quantify total standing biomass of high altitude Norway spruce stands.

Wiener System identification using B-spline functions with De Boor recursion

A simple and effective algorithm is introduced for the system identification of Wiener system based on the observational input/output data. The B-spline neural network is used to approximate the nonlinear static function in the Wiener system. We incorporate the Gauss-Newton algorithm with De Boor algorithm (both curve and the first order derivatives) for the parameter estimation of the Wiener model, together with the use of a parameter initialization scheme. The efficacy of the proposed approach is demonstrated using an illustrative example.

Asymptotic expansions for Taylor coefficients of the composition of two functions

We give an asymptotic expansion for the Taylor coe±cients of L(P(z)) where L(z) is analytic in the open unit disc whose Taylor coe±cients vary `smoothly' and P(z) is a probability generating function. We show how this result applies to a variety of problems, amongst them obtaining the asymptotics of Bernoulli transforms and weighted renewal sequences.

Periodic time series modeling of environmental proxy records with guaranteed positive growth rate estimation

Identifying a periodic time-series model from environmental records, without imposing the positivity of the growth rate, does not necessarily respect the time order of the data observations. Consequently, subsequent observations, sampled in the environmental archive, can be inversed on the time axis, resulting in a non-physical signal model. In this paper an optimization technique with linear constraints on the signal model parameters is proposed that prevents time inversions. The activation conditions for this constrained optimization are based upon the physical constraint of the growth rate, namely, that it cannot take values smaller than zero. The actual constraints are defined for polynomials and first-order splines as basis functions for the nonlinear contribution in the distance-time relationship. The method is compared with an existing method that eliminates the time inversions, and its noise sensitivity is tested by means of Monte Carlo simulations. Finally, the usefulness of the method is demonstrated on the measurements of the vessel density, in a mangrove tree, Rhizophora mucronata, and the measurement of Mg/Ca ratios, in a bivalve, Mytilus trossulus.

Biomass functions and expansion factors in young Norway spruce (Picea abies [L.] Karst) trees

Roots, stems, branches and needles of 160 Norway spruce trees younger than 10 years were sampled in seven forest stands in central Slovakia in order to establish their biomassfunctions (BFs) and biomassexpansionfactors (BEFs). We tested three models for each biomass pool based on the stem base diameter, tree height and the two parameters combined. BEF values decreased for all spruce components with increasing height and diameter, which was most evident in very young trees under 1 m in height. In older trees, the values of BEFs did tend to stabilise at the height of 3–4 m. We subsequently used the BEFs to calculate dry biomass of the stands based on average stem base diameter and tree height. Total stand biomass grew with increasing age of the stands from about 1.0 Mg ha−1 at 1.5 years to 44.3 Mg ha−1 at 9.5 years. The proportion of stem and branch biomass was found to increase with age, while that of needles was fairly constant and the proportion of root biomass did decrease as the stands grew older.

Well-posed two-point initial-boundary value problems with arbitrary boundary conditions

We study initial-boundary value problems for linear evolution equations of arbitrary spatial order, subject to arbitrary linear boundary conditions and posed on a rectangular 1-space, 1-time domain. We give a new characterisation of the boundary conditions that specify well-posed problems using Fokas' transform method. We also give a sufficient condition guaranteeing that the solution can be represented using a series. The relevant condition, the analyticity at infinity of certain meromorphic functions within particular sectors, is significantly more concrete and easier to test than the previous criterion, based on the existence of admissible functions.

Application of cylindrical distribution functions to wide-angle X-ray scattering from oriented polymers

A systematic approach is presented for obtaining cylindrical distribution functions (CDF's) of noncrystalline polymers which have been oriented by extension. The scattering patterns and CDF's are also sharpened by the method proposed by Deas and by Ruland. Data from atactic poly(methyl methacrylate) and polystyrene are analysed by these techniques. The methods could also be usefully applied to liquid crystals.

Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations

In this paper, we extend to the time-harmonic Maxwell equations the p-version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284] for Trefftz-discontinuous Galerkin approximations of the Helmholtz problem. While error estimates in a mesh-skeleton norm are derived parallel to the Helmholtz case, the derivation of estimates in a mesh-independent norm requires new twists in the duality argument. The particular case where the local Trefftz approximation spaces are built of vector-valued plane wave functions is considered, and convergence rates are derived.

System identification of Wiener systems with B-spline functions using De Boor recursion

In this article a simple and effective algorithm is introduced for the system identification of the Wiener system using observational input/output data. The nonlinear static function in the Wiener system is modelled using a B-spline neural network. The Gauss–Newton algorithm is combined with De Boor algorithm (both curve and the first order derivatives) for the parameter estimation of the Wiener model, together with the use of a parameter initialisation scheme. Numerical examples are utilised to demonstrate the efficacy of the proposed approach.

Rhinocetin, a venom-derived integrin-specific antagonist inhibits collagen-induced platelet and endothelial cell functions

Snaclecs are small non-enzymatic proteins present in viper venoms reported to modulate haemostasis of victims through effects on platelets, vascular endothelial and smooth muscle cells. In this study, we have isolated and functionally characterised a snaclec which we named rhinocetin from the venom of West African gaboon viper, Bitis gabonica rhinoceros. Rhinocetin was shown to comprise α and β chains with the molecular masses of 13.5 and 13kDa respectively. Sequence and immunoblot analysis of rhinocetin confirmed this to be a novel snaclec. Rhinocetin inhibited collagen-stimulated activation of human platelets in dose dependent manner, but displayed no inhibitory effects on glycoprotein VI (collagen receptor) selective agonist, CRP-XL-, ADP- or thrombin-induced platelet activation. Rhinocetin antagonised the binding of monoclonal antibodies against the α2 subunit of integrin α2β1 to platelets and coimmunoprecipitation analysis confirmed integrin α2β1 as a target for this venom protein. Rhinocetin inhibited a range of collagen induced platelet functions such as fibrinogen binding, calcium mobilisation, granule secretion, aggregation and thrombus formation. It also inhibited integrin α2β1 dependent functions of human endothelial cells. Together, our data suggest rhinocetin to be a modulator of integrin α2β1 function and thus may provide valuable insights into the role of this integrin in physiological and pathophysiological scenarios including haemostasis, thrombosis and envenomation.

Smoothness of scale functions for spectrally negative Lévy processes

Scale functions play a central role in the fluctuation theory of spectrally negative Lévy processes and often appear in the context of martingale relations. These relations are often require excursion theory rather than Itô calculus. The reason for the latter is that standard Itô calculus is only applicable to functions with a sufficient degree of smoothness and knowledge of the precise degree of smoothness of scale functions is seemingly incomplete. The aim of this article is to offer new results concerning properties of scale functions in relation to the smoothness of the underlying Lévy measure. We place particular emphasis on spectrally negative Lévy processes with a Gaussian component and processes of bounded variation. An additional motivation is the very intimate relation of scale functions to renewal functions of subordinators. The results obtained for scale functions have direct implications offering new results concerning the smoothness of such renewal functions for which there seems to be very little existing literature on this topic.