We discuss modelling and simulation of volumetric rainfall in a catchment of the Murray–Darling Basin – an important food production region in Australia that was seriously affected by a recent prolonged drought. Consequently, there has been sustained interest in development of improved water management policies. In order to model accumulated volumetric catchment rainfall over a fixed time period, it is necessary to sum weighted rainfall depths at representative sites within each sub-catchment. Since sub-catchment rainfall may be highly correlated, the use of a Gamma distribution to model rainfall at each site means that catchment rainfall is expressed as a sum of correlated Gamma random variables. We compare four different models and conclude that a joint probability distribution for catchment rainfall constructed by using a copula of maximum entropy is the most effective.

In oscillatory shear experiments, the values of the storage and loss moduli, $G^{\prime }(\unicode[STIX]{x1D714})$ and $G^{\prime \prime }(\unicode[STIX]{x1D714})$, respectively, are only measured and recorded for a number of values of the frequency $\unicode[STIX]{x1D714}$ in some well-defined finite range $[\unicode[STIX]{x1D714}_{\text{min}},\unicode[STIX]{x1D714}_{\text{max}}]$. In many practical situations, when the range $[\unicode[STIX]{x1D714}_{\text{min}},\unicode[STIX]{x1D714}_{\text{max}}]$ is sufficiently large, information about the associated polymer dynamics can be assessed by simply comparing the interrelationship between the frequency dependence of $G^{\prime }(\unicode[STIX]{x1D714})$ and $G^{\prime \prime }(\unicode[STIX]{x1D714})$. For other situations, the required rheological insight can only be obtained once explicit knowledge about the structure of the relaxation time spectrum $H(\unicode[STIX]{x1D70F})$ has been determined through the inversion of the measured storage and loss moduli $G^{\prime }(\unicode[STIX]{x1D714})$ and $G^{\prime \prime }(\unicode[STIX]{x1D714})$. For the recovery of an approximation to $H(\unicode[STIX]{x1D70F})$, in order to cope with the limited range $[\unicode[STIX]{x1D714}_{\text{min}},\unicode[STIX]{x1D714}_{\text{max}}]$ of the measurements, some form of localization algorithm is required. A popular strategy for achieving this is to assume that $H(\unicode[STIX]{x1D70F})$ has a separated discrete point mass (Dirac delta function) structure. However, this expedient overlooks the potential information contained in the structure of a possibly continuous $H(\unicode[STIX]{x1D70F})$. In this paper, simple localization algorithms and, in particular, a joint inversion least squares procedure, are proposed for the rapid recovery of accurate approximations to continuous $H(\unicode[STIX]{x1D70F})$ from limited measurements of $G^{\prime }(\unicode[STIX]{x1D714})$ and $G^{\prime \prime }(\unicode[STIX]{x1D714})$.

In this paper, we show that the cost of an optimal train journey on level track over a fixed distance is a strictly decreasing and strictly convex function of journey time. The precise structure of the cost–time curves for individual trains is an important consideration in the design of energy-efficient timetables on complex rail networks. The development of optimal timetables for busy metropolitan lines can be considered as a two-stage process. The first stage seeks to find optimal transit times for each individual journey segment subject to the usual trip-time, dwell-time, headway and connection constraints in such a way that the total energy consumption over all proposed journeys is minimized. The second stage adjusts the arrival and departure times for each journey while preserving the individual segment times and the overall journey times, in order to best synchronize the collective movement of trains through the network and thereby maximize recovery of energy from regenerative braking. The precise nature of the cost–time curve is a critical component in the first stage of the optimization.

We solve the problem of concept learning using a semi-tensor product method. All possible hypotheses are expressed under the framework of a semi-tensor product. An algorithm is raised to derive the version space. In some cases, the new approach improves the efficiency compared to the previous approach.

We derive a higher order nonlinear evolution equation for a broader bandwidth three-dimensional capillary–gravity wave packet, in the presence of a surface current produced by an internal wave. Instead of a set of coupled equations, a single nonlinear evolution equation is obtained by eliminating the velocity potential for the wave-induced slow motion. Finally, the equation is expressed in an integro-differential equation form, similar to Zakharov’s integral equation. Using the evolution equation derived here, we show that the two sidebands of a surface capillary–gravity wave get excited as a result of resonance with an internal wave, all propagating in the same direction. It is also shown that surface waves can grow exponentially with time at the expense of the energy of the internal wave.

Selecting important variables and estimating coordinate covariation have received considerable attention in the current big data deluge. Previous work shows that the gradient of the regression function, the objective function in regression and classification problems, can provide both types of information. In this paper, an algorithm to learn this gradient function is proposed for nonidentical data. Under some mild assumptions on data distribution and the model parameters, a result on its learning rate is established which provides a theoretical guarantee for using this method in dynamical gene selection and in network security for recognition of malicious online attacks.

We consider the optimal proportional reinsurance problem for an insurer with two dependent classes of insurance business, where the two claim number processes are correlated through a common shock component. Using the technique of stochastic linear–quadratic control theory and the Hamilton–Jacobi–Bellman equation, we derive the explicit expressions for the optimal reinsurance strategies and value function, and present the verification theorem within the framework of the viscosity solution. Furthermore, we extend the results in the linear–quadratic setting to the mean–variance problem, and obtain an efficient strategy and frontier. Some numerical examples are given to show the impact of model parameters on the efficient frontier.

We would like to present a method to compute the incompatibility operator in any system of curvilinear coordinates (components). The procedure is independent of the metric in the sense that the expression can be obtained by means of the basis vectors only, which are first defined as normal or tangential to the domain boundary, and then extended to the whole domain. It is an intrinsic method, to some extent, since the chosen curvilinear system depends solely on the geometry of the domain boundary. As an application, the in-extenso expression of incompatibility in a spherical system is given.

This paper deals with the flux identification problem for scalar conservation laws. The problem is formulated as an optimization problem, where the objective function compares the solution of the direct problem with observed profiles at a fixed time. A finite volume scheme solves the direct problem and a continuous genetic algorithm solves the inverse problem. The numerical method is tested with synthetic experimental data. Simulation parameters are recovered approximately. The tested heuristic optimization technique turns out to be more robust than classical optimization techniques.

The error of a distributed algorithm for big data classification with a support vector machine (SVM) is analysed in this paper. First, the given big data sets are divided into small subsets, on which the classical SVM with Gaussian kernels is used. Then, the classification error of the SVM for each subset is analysed based on the Tsybakov exponent, geometric noise, and width of the Gaussian kernels. Finally, the whole error of the distributed algorithm is estimated in terms of the error of each subset.

A new approach to jump diffusion is introduced, where the jump is treated as a vertical shift of the price (or volatility) function. This method is simpler than the previous methods and it is applied to the portfolio model with a stochastic volatility. Moreover, closed-form solutions for the optimal portfolio are obtained. The optimal closed-form solutions are derived when the value function is not smooth, without relying on the method of viscosity solutions.

As a new business form, the buy-online and pick-up-in-store (BOPS) mode allows consumers to pay for goods online and pick them up in a physical store. In this paper, an equilibrium model is constructed to formulate an optimal decision-making problem for online and offline retailers under the BOPS mode, where the online retailer determines the retail price of the goods and the consignment quantity in a physical store, while the offline retailer chooses the revenue share of profit by a consignment contract. Different to the existing models, the cost of overstocking and loss of understocking are incorporated into the profit function of the online retailer due to the randomness of demand. For the objective function of the offline retailer, the cross-sale quantity generated by the BOPS mode is taken into account. Then the game between the online and offline retailers is expressed as a stochastic Nash equilibrium model. Based on the analytic properties of the model, necessary conditions for the equilibrium solution are obtained. A case study and sensitivity analysis are employed to reveal the managerial implications of the model, which can provide a number of valuable suggestions on optimizing the strategies for the online and offline retailers under the BOPS mode.

We address the construction and approximation for feed-forward neural networks (FNNs) with zonal functions on the unit sphere. The filtered de la Vallée-Poussin operator and the spherical quadrature formula are used to construct the spherical FNNs. In particular, the upper and lower bounds of approximation errors by the FNNs are estimated, where the best polynomial approximation of a spherical function is used as a measure of approximation error.

We generate an algebra on blood phenotypes with multiplication based on the human ABO-blood group inheritance pattern. We assume that gametes are not chosen randomly during meiosis. We investigate some of the properties of this algebra, namely, the set of idempotents, lattice of ideals and the associative enveloping algebra.

A theoretical study of an unsteady two-layered blood flow through a stenosed artery is presented in this article. The geometry of a rigid stenosed artery is assumed to be $w$-shaped. The flow regime is assumed to be laminar, unsteady and uni-directional. The characteristics of blood are modelled by the generalized Oldroyd-B non-Newtonian fluid model in the core region and a Newtonian fluid model in the periphery region. The governing partial differential equations are derived for each region by using mass and momentum conservation equations. In order to facilitate numerical solutions, the derived differential equations are nondimensionalized. A well-tested explicit finite-difference method (FDM) which is forward in time and central in space is employed for the solution of a nonlinear initial boundary value problem corresponding to each region. Validation of the FDM computations is achieved with a variational finite element method algorithm. The influences of the emerging geometric and rheological parameters on axial velocity, resistance impedance and wall shear stress are displayed graphically. The instantaneous patterns of streamlines are also presented to illustrate the global behaviour of the blood flow. The simulations are relevant to haemodynamics of small blood vessels and capillary transport, wherein rheological effects are dominant.

Image inpainting methods recover true images from partial noisy observations. Natural images usually have two layers consisting of cartoons and textures. Methods using simultaneous cartoon and texture inpainting are popular in the literature by using two combined tight frames: one (often built from wavelets, curvelets or shearlets) provides sparse representations for cartoons and the other (often built from discrete cosine transforms) offers sparse approximation for textures. Inspired by the recent development on directional tensor product complex tight framelets ($\text{TP}\text{-}\mathbb{C}\text{TF}$s) and their impressive performance for the image denoising problem, we propose an iterative thresholding algorithm using tight frames derived from $\text{TP}\text{-}\mathbb{C}\text{TF}$s for the image inpainting problem. The tight frame $\text{TP}\text{-}\mathbb{C}\text{TF}_{6}$ contains two classes of framelets; one is good for cartoons and the other is good for textures. Therefore, it can handle both the cartoons and the textures well. For the image inpainting problem with additive zero-mean independent and identically distributed Gaussian noise, our proposed algorithm does not require us to tune parameters manually for reasonably good performance. Experimental results show that our proposed algorithm performs comparatively better than several well-known frame systems for the image inpainting problem.

The central projection transform can be employed to extract invariant features by combining contour-based and region-based methods. However, the central projection transform only considers the accumulation of the pixels along the radial direction. Consequently, information along the radial direction is inevitably lost. In this paper, we propose the Mellin central projection transform to extract affine invariant features. The radial factor introduced by the Mellin transform, makes up for the loss of information along the radial direction by the central projection transform. The Mellin central projection transform can convert any object into a closed curve as a central projection transform, so the central projection transform is only a special case of the Mellin central projection transform. We prove that closed curves extracted from the original image and the affine transformed image by the Mellin central projection transform satisfy the same affine transform relationship. A method is provided for the extraction of affine invariants by employing the area of closed curves derived by the Mellin central projection transform. Experiments have been conducted on some printed Chinese characters and the results establish the invariance and robustness of the extracted features.