A semi-Markov process with an inverse Gaussian distribution as sojourn time

Let X_n denote the state of a device after n repairs. We assume that the time between two repairs is the time τ taken by a Wiener process {W(t), t ? 0}, starting from w₀ and with drift μ < 0, to reach c ∈ [0, w₀). After the nth repair, the process takes on either the value X_n?1 + 1 or X_n?1 + 2. The probability that X_n = X_n?1 + j, for j = 1, 2, depends on whether τ ? t₀ (a fixed constant) or τ > t₀. The device is considered to be worn out when X_n ? k, where k ∈ {1, 2, …}. This model is based on the ones proposed by Rishel (1991) [1] and Tseng and Peng (2007) [2]. We obtain an explicit expression for the mean lifetime of the device. Numerical methods are used to illustrate the analytical findings.     

On Certain C-Test Words for Free Groups

Let F_m be a free group of a finite rank m ≥ 2 and let X_i, Y_j be elements in F_m. A non-empty word w(x₁,…,x_n) is called a C-test word in n letters for F_m if, whenever (X₁,…,X_n) = w(Y₁,…,Y_n) ≠ 1, the two n-typles (X₁,…,X_n) and (Y₁,…,Y_n) are conjugate in F_m. In this paper we construct, for each n ≥ 2, a C-test word v_n(x₁,…,x_n) with the additional property that v_n(X₁,…,X_n) = 1 if and only if the subgroup of F_m generated by X₁,…,X_n is cyclic. Making use of such words v_m(x₁,…,x_m) and v_m + 1(x₁,…,x_m + 1), we provide a positive solution to the following problem raised by Shpilrain: There exist two elements u₁, u₂ ∈ F_m such that every endomorphism ψ of F_m with non-cyclic image is completely determined by ψ(u₁), ψ(u₂).     

Interpolation by Integer-Valued Polynomials

LetRbe a Krull ring with quotient fieldKanda₁,…,a_ninR. If and only if thea_iare pairwise incongruent mod every height 1 prime ideal of infinite index inRdoes there exist for all valuesb₁,…,b_ninRan interpolating integer-valued polynomial, i.e., anf ∈ K[x] withf(a_i) = b_iandf(R) ⊆ R.IfSis an infinite subring of a discrete valuation ringR_vwith quotient fieldKanda₁,…,a_ninSare pairwise incongruent mod allM^k_v ∩ Sof infinite index inS, we also determine the minimald(depending on the distribution of thea_iamong residue classes of the idealsM^k_v ∩ S) such that for allb₁,…,b_n ∈ R_vthere exists a polynomialf ∈ K[x] of degree at mostdwithf(a_i) = b_iandf(S) ⊆ R_v.     

Discrepancy principle for DSM II

Let Ay = f, A is a linear operator in a Hilbert space H, y ⊥ N(A) ≔ {u : Au = 0}, R(A) ≔ {h : h = Au, u ∈ D(A)} is not closed, ∥f_δ − f∥ ⩽ δ. Given f_δ, one wants to construct u_δ such that lim_δ→0∥u_δ − y∥ = 0. Two versions of discrepancy principles for the DSM (dynamical systems method) for finding the stopping time and calculating the stable solution u_δ to the original equation Ay = f are formulated and mathematically justified.     

New exact solitary-wave special solutions for the nonlinear dispersive K(m, n) equations

In this paper, we study the nonlinear dispersive K(m, n) equations: u_t + (u^m)_x − (uⁿ)_xxx = 0 which exhibit solutions with solitary patterns. New exact solitary solutions are found. The two special cases, K(2, 2) and K(3, 3), are chosen to illustrate the concrete features of the decomposition method in K(m, n) equations. The nonlinear equations K(m, n) are studied for two different cases, namely when m = n being odd and even integers. General formulas for the solutions of K(m, n) equations are established.     

Computing best-possible bounds for the distribution of a sum of several variables is NP-hard

In many real-life situations, we know the probability distribution of two random variables x₁ and x₂, but we have no information about the correlation between x₁ and x₂; what are the possible probability distributions for the sum x₁ + x₂? This question was originally raised by A.N. Kolmogorov. Algorithms exist that provide best-possible bounds for the distribution of x₁ + x₂; these algorithms have been implemented as a part of the efficient software for handling probabilistic uncertainty. A natural question is: what if we have several (n > 2) variables with known distribution, we have no information about their correlation, and we are interested in possible probability distribution for the sum y = x₁ + ⋯ + x_n? Known formulas for the case n = 2 can be (and have been) extended to this case. However, as we prove in this paper, not only are these formulas not best-possible anymore, but in general, computing the best-possible bounds for arbitrary n is an NP-hard (computationally intractable) problem.     

Matrix characterization of linear codes with arbitrary Hamming weight hierarchy

The support of an [n, k] linear code C over a finite field F_q is the set of all coordinate positions such that at least one codeword has a nonzero entry in each of these coordinate position. The rth generalized Hamming weight d_r(C), 1 ⩽ r ⩽ k, of C is defined as the minimum of the cardinalities of the supports of all [n, r] subcodes of C. The sequence (d₁(C), d₂(C), … , d_k(C)) is called the Hamming weight hierarchy (HWH) of C. The HWH, d_r(C) = n − k + r;  r = 1, 2 , …, k, characterizes maximum distance separable (MDS) codes. Therefore the matrix characterization of MDS codes is also the characterization of codes with the HWH d_r(C) = n − k + r; r = 1, 2, … , k. A linear code C with systematic check matrix [I∣P], where I is the (n − k) × (n − k) identity matrix and P is a (n − k) × k matrix, is MDS iff every square submatrix of P is nonsingular. In this paper we extend this characterization to linear codes with arbitrary HWH. Using this result, we characterize Near-MDS codes, Near-Near-MDS (N²-MDS) codes and A^μ-MDS codes. The MDS-rank of C is the smallest integer η such that d_η+1 = n − k + η + 1 and the defect vector of C with MDS-rank η is defined as the ordered set {μ₁(C), μ₂(C), μ₃(C), … , μ_η(C), μ_η+1(C)}, where μ_i(C) = n − k + i − d_i(C). We call C a dually defective code if the defect vector of the code and its dual are the same. We also discuss matrix characterization of dually defective codes. Further, the codes meeting the generalized Greismer bound are characterized in terms of their generator matrix. The HWH of dually defective codes meeting the generalized Greismer bound are also reported.     

Pattern Frequency Sequences and Internal Zeros

Let q be a pattern and let S_n, q(c) be the number of n-permutations having exactly c copies of q. We investigate when the sequence (S_n, q(c))_c ≥ 0 has internal zeros. If q is a monotone pattern it turns out that, except for q = 12 or 21, the nontrivial sequences (those where n is at least the length of q) always have internal zeros. For the pattern q = 1(l + 1)l…2 there are infinitely many sequences which contain internal zeros and when l = 2 there are also infinitely many which do not. In the latter case, the only possible places for internal zeros are the next-to-last or the second-to-last positions. Note that by symmetry this completely determines the existence of internal zeros for all patterns of length at most 3.     

Electrostatics in fractal geometry: Fractional calculus approach

An electric field in a composite dielectric with a fractal charge distribution is obtained in the spherical symmetry case. The method is based on the splitting of a composite volume into a fractal volume V_d ～ r^d with the fractal dimension d and a complementary host volume V_h = V₃ ? V_d. Integrations over these fractal volumes correspond to the convolution integrals that eventually lead to the employment of the fractional integro-differentiation.     

Differential invariants of the one-dimensional quasi-linear second-order evolution equation

We consider evolution equations of the form u_t = f(x, u, u_x)u_xx + g(x, u, u_x) and u_t = u_xx + g(x, u, u_x). In the spirit of the recent work of Ibragimov [Ibragimov NH. Laplace type invariants for parabolic equations. Nonlinear Dynam 2002;28:125–33] who adopted the infinitesimal method for calculating invariants of families of differential equations using the equivalence groups, we apply the method to these equations. We show that the first class admits one differential invariant of order two, while the second class admits three functional independent differential invariants of order three. We use these invariants to determine equations that can be transformed into the linear diffusion equation.     

Statistical properties of nucleotides in human chromosomes 21 and 22

In this paper the statistical properties of nucleotides in human chromosomes 21 and 22 are investigated. The n-tuple Zipf analysis with n = 3, 4, 5, 6, and 7 is used in our investigation. It is found that the most common n-tuples are those which consist only of adenine (A) and thymine (T), and the rarest n-tuples are those in which GC or CG pattern appears twice. With the n-tuples become more and more frequent, the double GC or CG pattern becomes a single GC or CG pattern. The percentage of four nucleotides in the rarest ten and the most common ten n-tuples are also considered in human chromosomes 21 and 22, and different behaviors are found in the percentage of four nucleotides. Frequency of appearance of n-tuple f(r) as a function of rank r is also examined. We find the n-tuple Zipf plot shows a power-law behavior for r < 4ⁿ⁻¹ and a rapid decrease for r > 4ⁿ⁻¹. In order to explore the interior statistical properties of human chromosomes 21 and 22 in detail, we divide the chromosome sequence into some moving windows and we discuss the percentage of ξη (ξ, η = A, C, G, T) pair in those moving windows. In some particular regions, there are some obvious changes in the percentage of ξη pair, and there maybe exist functional differences. The normalized number of repeats N₀(l) can be described by a power law: N₀(l) ∼ l^−μ. The distance distributions P₀(S) between two nucleotides in human chromosomes 21 and 22 are also discussed. A two-order polynomial fit exists in those distance distributions: log P₀(S) = a + bS + cS², and it is quite different from the random sequence.     

An extended rational interpolation method

To interpolate function, f(x), a ? x ? b, when we have some information about the values of f(x) and their derivatives in separate points on {x₀, x₁, … , x_n} ? [a, b], the Hermit interpolation method is usually used. Here, to solve this kind of problems, extended rational interpolation method is presented and it is shown that the suggested method is more efficient and suitable than the Hermit interpolation method, especially when the function f(x) has singular points in interval [a, b]. Also for implementing the extended rational interpolation method, the direct method and the inverse differences method are presented, and with some examples these arguments are examined numerically.     

A Bijective Approach to the Area of Generalized Motzkin Paths

For fixed positive integer k, let E_n denote the set of lattice paths using the steps (1, 1), (1, − 1), and (k, 0) and running from (0, 0) to (n, 0) while remaining strictly above the x-axis elsewhere. We first prove bijectively that the total area of the regions bounded by the paths of E_n and the x-axis satisfies a four-term recurrence depending only on k. We then give both a bijective and a generating function argument proving that the total area under the paths of E_n equals the total number of lattice points on the x-axis hit by the unrestricted paths running from (0, 0) to (n − 2, 0) and using the same step set as above.     

The stochastic soliton-like solutions of stochastic KdV equations

By means of a generalized method and symbolic computation, we consider a stochastic KdV equation U_t + f(t)U ♢ U_x + g(t)U_xxx = W(t) ♢ R^♢(t, U, U_x, U_xxx). We construct new and more general formal solutions. At the same time, we recover all the solutions found by Xie [Phys. Lett. A 310 (2003) 161]. The solutions obtained include the nontravelling wave and coefficient function’s stochastic soliton-like solutions, singular stochastic soliton-like solutions, stochastic triangular functions solutions.     

Solution of the Burgers equation using an implicit linearizing transformation

The initial-boundary-value problem on the semi-infinite interval and on a finite interval for the Burgers equation u_t = u_xx + 2u_xu is solved using a stream function ? and a linearizing transformation w = e^?. The transformation reduces the equation to a heat equation with appropriate initial and homogeneous time-dependent linear boundary conditions. One advantage of this method is that we never need to find an explicit expression for ? in our computations.     

Complex dynamics in Duffing–Van der Pol equation

Duffing–Van der Pol equation with fifth nonlinear-restoring force and two external forcing terms is investigated. The threshold values of existence of chaotic motion are obtained under the periodic perturbation. By second-order averaging method and Melnikov method, we prove the criterion of existence of chaos in averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 1, 3, 5, and cannot prove the criterion of existence of chaos in second-order averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 2, 4, 6, 7, 8, 9, 10, where σ is not rational to ω₁, but can show the occurrence of chaos in original system by numerical simulation. Numerical simulations including heteroclinic and homoclinic bifurcation surfaces, bifurcation diagrams, Lyapunov exponent, phase portraits and Poincaré map, not only show the consistence with the theoretical analysis but also exhibit the more new complex dynamical behaviors. We show that cascades of interlocking period-doubling and reverse period-doubling bifurcations from period-2 to -4 and -6 orbits, interleaving occurrence of chaotic behaviors and quasi-periodic orbits, transient chaos with a great abundance of period windows, symmetry-breaking of periodic orbits in chaotic regions, onset of chaos which occurs more than one, chaos suddenly disappearing to period orbits, interior crisis, strange non-chaotic attractor, non-attracting chaotic set and nice chaotic attractors. Our results show many dynamical behaviors and some of them are strictly departure from the behaviors of Duffing–Van der Pol equation with a cubic nonlinear-restoring force and one external forcing.     

Complexity Estimates for the Schmüdgen Positivstellensatz

LetKbe a closed basic set inRⁿgiven by the polynomial inequalities φ₁≥ 0, . . . , φ_m≥ 0 and let Σ be the semiring generated by the φ_kand the squares inR[x₁, . . . ,x_n]. Schmüdgen has shown that ifKis compact then any polynomial function strictly positive onKbelongs to Σ. Easy consequences are (1)f≥ 0 onKif and only iff∈R⁺+ Σ (Positivstellensatz) and (2) iff≥ 0 onKbutf∈ Σ then asdtends to 0⁺, in any representation off + das an element of Σ in terms of the φ_k, the squares and semiring operations, the integerN(d) which is the minimum over all representations of the maximum degree of the summands must become arbitrarily large. A one-dimensional example is analyzed to obtain asymptotic lower and upper bounds of the formcd^−1/2≤N(d) ≤Cd^−1/2log (1/d).     

On Superheight Conditions for the Affineness of Open Subsets

In this paper we consider the open complement U of a hypersurface Y = V(a) in an affine scheme X. We study the relations between the affineness of U, the intersection of Y with closed subschemes, the property that every closed surface in U is affine, the property that every analytic closed surface is Stein, and the superheight of a defining ideal a.