Asymptotic Formulae and Divisor Problems

We investigate the asymptotic behaviour of the summatory functions of _z(n, ), _k(n, ) z ⁽ⁿ⁾ and _k(n, ) z ⁽ⁿ⁾.     

Aperture-Angle Optimization Problems in Three Dimensions

Let [a,b] be a line segment with end points a, b and a point at which a viewer is located, all in R ³. The aperture angle of [a,b] from point , denoted by (), is the interior angle at of the triangle (a,b,). Given a convex polyhedron P not intersecting a given segment [a,b] we consider the problem of computing _max() and _min(), the maximum and minimum values of () as varies over all points in P. We obtain two characterizations of _max(). Along the way we solve several interesting special cases of the above problems and establish linear upper and lower bounds on their complexity under several models of computation.     

Parametric stochastic convexity and concavity of stochastic processes

A collection of random variables {X(), } is said to be parametrically stochastically increasing and convex (concave) in if X() is stochastically increasing in , and if for any increasing convex (concave) function , E(X()) is increasing and convex (concave) in whenever these expectations exist. In this paper a notion of directional convexity (concavity) is introduced and its stochastic analog is studied. Using the notion of stochastic directional convexity (concavity), a sufficient condition, on the transition matrix of a discrete time Markov process {X _n(), n=0,1,2,...}, which implies the stochastic monotonicity and convexity of {X _n(), }, for any n, is found. Through uniformization these kinds of results extend to the continuous time case. Some illustrative applications in queueing theory, reliability theory and branching processes are given.Supported by the Air Force Office of Scientific Research, U.S.A.F., under Grant AFOSR-84-0205. Reproduction in whole or in part is permitted for any purpose by the United States Government.     

Minimum Divergence Estimators Based on Grouped Data

The paper considers statistical models with real-valued observations i.i.d. by F(x, ₀) from a family of distribution functions (F(x, ); ), R ^s , s 1. For random quantizations defined by sample quantiles (F _n ^–1 (₁),, F _n ^–1 ( _m–1)) of arbitrary fixed orders 0 < ₁ < _m-1 < 1, there are studied estimators _,n of ₀ which minimize -divergences of the theoretical and empirical probabilities. Under an appropriate regularity, all these estimators are shown to be as efficient (first order, in the sense of Rao) as the MLE in the model quantified nonrandomly by (F ^–1 (₁,₀),, F ^–1 ( _m–1, ₀)). Moreover, the Fisher information matrix I _m (₀, ) of the latter model with the equidistant orders = ( _j = j/m : 1 j m – 1) arbitrarily closely approximates the Fisher information J(₀) of the original model when m is appropriately large. Thus the random binning by a large number of quantiles of equidistant orders leads to appropriate estimates of the above considered type.     

Exact Solution of the Ising Model on the Cayley Tree with Competing Ternary and Binary Interactions

The exact solution is found for the problem of phase transitions in the Ising model with competing ternary and binary interactions. For the pair of parameters =(J) and ₁=₁(J ₁) in the plane (₁,), we find two critical curves such that a phase transition occurs for all pairs (₁,) lying between the curves.     

Prediction Problems for Square-Transformed Stationary Processes

This paper discusses the prediction problems for square-transformed process, Y _t = X _t ², where X _t is a stationary process with spectral density g(). The square-transformation is important in prediction of the volatility of ARCH models. First, we evaluate the mean square prediction error for square-transformed process when the predictor is constructed from the true spectral density g(). However, it is often that the true structure g() is not completely specified. Hence, we consider the problem of misspecified prediction when a conjectured spectral density f (), , is fitted to g(). Then, constructing the best linear predictor based on f (), we can evaluate the prediction error for square-transformed process. Also, we consider a bias adjusted prediction problem for the above two cases. Furthermore, we may suppose that X _t is a non-Gaussian process. Then, we evaluate the mean square prediction errors when the best linear predictor is constructed by the true spectral density g() and the conjectured spectral density f (), respectively. Since is usually unknown we estimate it by a quasi-MLE . The second-order asymptotic approximations of the mean square errors of the predictors based on g() and f () are given. Finally, we provide some numerical examples, which show some unexpected features.     

Solution of the nonlinear functional equations representing the roll gap relationships in a cold mill

In the development of a roll force model for cold rolling, techniques were developed for solving the system equations which are of general interest. This paper gives a brief introduction to the physical model but concentrates on the solution of the model equations and the simulation. An unusual feature of the model was that the calculated profiles had to satisfy a number of boundary conditions at different points throughout the roll arc. A new method was developed for calculating these profiles and for determining the gradient functions which satisfied the boundary constraints.Nomenclature p() pressure at roll angle - h() gauge - a() roll radius - y() yield stress - g _i () gradient function on iterationi - e() gauge error - (, ) transition function - H(–) Heaviside unit step function at = - (–) unit impulse function at = - H(, ₁, ₂) defined asH( – ₁) –H( – ₂) - angular position from the roll center line - _T angular limits of roll arc represented - _n angular position of the neutral angle - _i angular position ofith strip elastic-plastic boundary - _pi pressure change at the boundaryi - _i , _i , _i constants defined in Appendix A - k ₁,k ₂ elastic region constants - k total number of strip boundaries (elastic-plastic and entry and exit points) - R undeformed work roll radius - R _s roll separation—distance between roll centers - h ₀₁ unstrained gauge in an elastic region - h _in gauge of the strip at the entry to the roll gap - J gauge error cost function - <x, y> inner product ofx andy - x norm ofx - L ₂[0, _T ] the space of Lebesgue square-integrable functions defined on the interval [0, _T ] - JUVY denotes (Dx)() =dx()/d The author would like to acknowledge the help given by Dr. G. F. Bryant, Director, and Mr. M. A. Fuller, Senior Research Engineer, the Industrial Automation Group, Imperial College of Science and Technology, London. He is also grateful to M. J. G. Henderson of the University of Birmingham for his advice and encouragement during the project. He would like to thank the Directors of GEC Electrical Projects Limited for allowing him to undertake the work and also Mr. J. McTaggart and Mr. C. McKenzie (GEC), Professor H. A. Prime of the University of Birmingham, and Dr. G. F. Bryant for arranging the project.     

Lévy homeomorphic parametrization and exponential families

Summary Let P={P : } be an exponential family of probability distributions with the canonical parameter and consider the one to one mapping : P . It is shown that, under mild regularity assumptions, and ^–1 are continuous with respect to the Lévy metric in P and Euclidean metric in .     

Diffusion Equations and Geometric Inequalities

Let =(₀, ₁) be a fixed vector in R ² with strictly positive components and suppose ₀, ₁ > 0. Set = _{0 0} + _{1 1} and, if x ₀, x ₁ R ⁿ , set x = ₀ x ₀ + ₁ x ₁. Moreover, for any j {0, 1, }, let c _j : R ⁿ R be a continuous, bounded function and denote by p _j , c _j (t, x, y) the fundamental solution of the diffusion equation

On the factor-thickness of regular graphs

For eachr-regular graphG, define a binary sequence(G) = ( ₁, ₂,..., _r-1) by _k = 0 ifG has ak-factor, and _k = 1 otherwise. A binary sequence = ( _i |i = 1, 2,...,r – 1) is said to be realizable if there exists anr-regular graphG such that(G) = . In this paper we characterize all binary sequences which are realizable.     

Inner functions and related spaces of pseudocontinuable functions

Let be an inner function, let C, ¦¦=1. Then the harmonic function [(+)]/(–)] is the Poisson integral of a singular measure D. N. Clark's known theorem enables us to identify in a natural manner the space H² H² with the space L² ( ).Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova Akademii Nauk SSSR, Vol. 170, pp. 7–33, 1989.     

Arithmetic of the Modular Function j1,8

Since the genus of the modular curve X_1 (8) = _1 (8) * is zero, we find a field generator j _1,8(z) = ₃(2z)/₃(4z) (₃(z) := _n ein ^2z ) such that the function field over X ₁(8) is (j _1,8). We apply this modular function j _1,8 to the construction of some class fields over an imaginary quadratic field K, and compute the minimal polynomial of the singular value of the Hauptmodul N(j _1,8) of (j _1,8).     

Radii of convexity and close-to-convexity of certain integral representations

Strict upper bounds are determined for ¦s(z)¦, ¦Re s(z)¦, and ¦Im s(z) ¦ in the class of functions s(z)=a _nzⁿ+a_n+1zⁿ⁺¹+... (n¹) regular in ¦z¦<1 and satisfying the condition ¦u (₁) –u (₂) ¦K¦ ₁- ₂¦, where U()=Re s (eⁱ ), K>0, and ₁ and ₂ are arbitrary real numbers. These bounds are used in the determination of radii of convexity and close-to-convexity of certain integral representations.Translated from Matematicheskie Zametki, Vol. 7, No. 5, pp. 581–592, May, 1970.The author wishes to thank L. A. Aksent'ev for his guidance in this work.     

On the Embedding of an Affine Space into a Projective Space

Let k, K be fields, and assume that |k| 4 and n, m 2, or |k| = 3 and n 3, m 2. Then, for any embedding of AG(n, k) into PG(m, K), there exists an isomorphism from k into K and an (n+1) × (m+1) matrix B with entries in K such that can be expressed as (x₁,x₂,...,x_n) = [(1,x₁ ,x₂ ,...,x_n )B], where the right-hand side is the equivalence class of (1,x₁ ,x₂ ,...,x_n )B. Moreover, in this expression, is uniquely determined, and B is uniquely determined up to a multiplication of element of K*. Let l 1, and suppose that there exists an embedding of AG(m+l, k) into PG(m, K) which has the above expression. If we put r = dim _k K, then we have r 3 and m > 2 l-1)/(r-2). Conversely, there exists an embedding of AG(l+m, k) into PG(m, K) with the above expression if K is a cyclic extension of k with dim _k K=r 3, and if m 2l/(r-2) with m even or if m 2l/(r-2) +1 with m odd.     

A Comparison of Methods for Estimating the Extremal Index

The extremal index, (01), is the key parameter when extending discussions of the limiting behavior of the extreme values from independent and identically distributed sequences to stationary sequences. As measures the limiting dependence of exceedances over a threshold u, as u tends to the upper endpoint of the distribution, it may not always be informative about the extremal dependence at levels of practical interest. Therefore we also consider a threshold-based extremal index, (u). We compare the performance of a range of different estimators for and (u) covering processes with < 1 and = 1. We find that the established methods for estimating actually estimate (u), so perform well only when (u) . For Markov processes, we introduce an estimator which is as good as the established methods when (u) but provides an improvement when (u) < = 1. We illustrate our methods using simulated data and daily rainfall measurements.     

Das Zeitverhalten von einfachen offenen exponentiellen Warteschlangensystemen mit unendlich vielen Warteplätzen

Zusammenfassung Die zeitabhängige (instationäre) Lösung für die Zustandswahrscheinlichkeiten und für einige Kenngrößen von Warteschlangensystemen mit einer Bedienungsstation, unendlich vielen Warteplätzen, exponentiellem Zu- und Abgang und beliebigem Anfangszustand wird bestimmt. Die ZustandswahrscheinlichkeitenP _v (), d. h. die Wahrscheinlichkeiten für Einheiten im System zur Zeit, ergeben sich als Integrale, in denen modifizierteSessel-Funktionen 1. Art auftreten. Der ErwartungswertL () und die VarianzV() der Zahl von Einheiten im System lassen sich als Integrale darstellen, in denen nur die ZustandswahrscheinlichkeitP ₀() auftritt.Für<1 und erreichen die Systeme einen stationären Zustand (für den die Lösung bekannt ist); für1 und giltP _v ()0 für alle, L(),V().Ist>1, dann wachsenL() undV() für große linear mit; ihre Asymptoten werden berechnet. Ist=1, dann wachsenL() und die Standardabweichung() für große mit ; einfache Näherungsformeln werden gefunden.

---

Summary The time dependent solution is determined for the state probabilities and for some characteristic values of queuing systems with a single server, an infinite number of waiting places, exponentially distributed inter-arrival and service times, and any initial state. The state probabilitiesP _v (), i.e. the probabilities for units in the system at time, are given in the form of integrals in which modifiedBessel functions of the first kind occur. Integrating the state probalityP ₀() over leads to the meanL() and the varianceV() of the number of units in the system.For<1 and the systems tend to a steady state (for which the solution is known); for1 and we haveP _v ()0 for all, L(),V().If>1 asymptotic expansions for large are found givingL() andV() proportional to. If=1 simple approximate formulas for large are obtained givingL() and the standard deviation() proportional to .

---

---

Vorgel. v.:J. Nitsche.     

Similarity solutions of radially symmetric two-phase flow

When both the diffusivityD and fractional flow functionf have a power law dependence on the water content , i.e.D=D _o andf=⁺¹, the nonlinear transport equation for radially symmetric two phase flow can, in certain circumstances, be reduced to a weakly coupled system of two first order nonlinear ordinary differential equations. Numerical solutions of these equations for a constant flux boundary conditionV _wo and comparison with experimental data are given. In particular, when the fluxV _wo and a are related byV _wo( + 1)/D _o=2, a new fully explicit analytical solution is found as (r, t)=(1 – r ²/4D _ot)^1/ forr ² < 4D _ot/ and (r, t)=0 forr ² 4D _ot/ We show that the existence of this exact soution is due to the presence of a Lagrangian symmetry.     

A graded scale of parametric families of distributions,and parameter estimates based on the sample mean

Summary Denote by _k a class of familiesP={P} of distributions on the line R¹ depending on a general scalar parameter , being an interval of R¹, and such that the moments µ₁()=xdP ,...,µ_2k ()=x ^2k dP are finite, ₁ (), ..., _k (), _k+1 () ..., _k () exist and are continuous, with ₁ () 0, and _j +1 ()= ₁ () _j () +[₂() -₁()²] _j ()/ ₁ (), J=2, ..., k. Let ₁=¯x=x ₁ + ... +x _n/n, ₂=x ₁ ² + ... +x _n ²/n, ..., _k =(x ₁ ^k + ... +x _n ^k/n denote the sample moments constructed for a sample x₁, ..., x_n from a population with distribution Pg. We prove that the estimator of the parameter by the method of moments determined from the equation ₁= ₁() and depending on the observations x₁, ..., x_n only via the sample mean ¯x is asymptotically admissible (and optimal) in the class _k of the estimators determined by the estimator equations of the form ₀ () + ₁ () ₁ + ... + _k () _k =0 if and only ifP _k .The asymptotic admissibility (respectively, optimality) means that the variance of the limit, as n (normal) distribution of an estimator normalized in a standard way is less than the same characteristic for any estimator in the class under consideration for at least one 9 (respectively, for every ).The scales arise of classes _{1 2}... of parametric families and of classes _{1 2} ... of estimators related so that the asymptotic admissibility of an estimator by the method of moments in the class _k is equivalent to the membership of the familyP in the class _k .The intersection consists only of the families of distributions with densities of the form h(x) exp {C₀() + C₁() x } when for the latter the problem of moments is definite, that is, there is no other family with the same moments ₁ (), ₂ (), ...Such scales in the problem of estimating the location parameter were predicted by Linnik about 20 years ago and were constructed by the author in [1] (see also [2, 3]) in exact, not asymptotic, formulation.Translated from Problemy Ustoichivosti Stokhasticheskikh Modelei, pp. 41–47, 1981.     

On accuracy and unconditional stability of linear multistep methods for second order differential equations

Linear multistep methods for solution of the equationy=f(t, y) are studied by means of the test equationy=–² y, with real. It is shown that the order of accuracy cannot exceed 2 for an unconditionally stable method.This work was supported by the NASA-Ames Research Center, Moffett Field, California, under Interchange No. NCA2-OR745-712, while the author was a visitor at the Computer Science Department, Stanford University, Stanford, California.     

On dense sequences of polynomials in several variables

In the present paper a form of generalization of Gelfond's lemma on dense sequences of polynomials is proposed. For a set of complex numbers ₁, ..., _s we define the coefficients_gk( ₁, ..., _s ) (0ks) and give the relations between them and the transcendental degrees or the transcendence types of the field © ( ₁, ..., _s ) or its subfields.This work was completed at the Dept. of Math., Univ. of Southern Mississippi, USA in 1987.