A numerical method for diffusion-convection equation using high-order difference schemes

In this paper, we propose a simple general form of high-order approximation of O(c²+ch²+h⁴) to solve the two-dimensional parabolic equation αuxx+βuyy=F(x,y,t,u,ux,uy,ut), where α and β are positive constants. We apply the compact form for solving diffusion-convection equation. The results of numerical experiments are presented and compared with analytical solutions to confirm the higher accuracy of the presented scheme.     

Control of singular problem via differentiation of a Min-Max

For Ω a smooth domain in Rⁿ with boundary Λ = Λ₀Λ₁, we are concerned with the wave equation y″ − Δy = S in Q_T =]0, T[ × Ω with = ∂/∂t, at source term satisfying S, S′ ε L¹(0, T L² (Ω)). A Dirichlet condition is imposed on Λ₀ and we consider an absorbing condition ∂y/∂n + uy′ = 0 in [0, T] × gL₁ where u is the control.parameter. We introduce the cost function. and using the Min-Max formulation of J we by-pasas the sensitivity analysis of u → y and obtain the gradient of J with a usual adjoint problem. We first present an abstract frame for this kind of problems. using the differentiability results of a Min-Max [1,2], which we very shortly deduce here, we show that the well posedness of the adjoint equation implies differentiability of the cost function governed by a linear well posed problem.     

Parabolic equations with double variable nonlinearities

The paper is devoted to the study of the homogeneous Dirichlet problem for the doubly nonlinear parabolic equation with nonstandard growth conditions:

ut=div(a(x,t,u)|u^|α(x,t)|∇u^|p(x,t)−2∇u)+f(x,t)     

Two-sided bounds for the asymptotic behaviour of free nonlinear vibration systems with application of the differential calculus of norms

This paper deals with the free nonlinear dynamical system [xdot](t)=A x(t)+h(t, x(t)), t≥t ₀, x(t ₀)=x ₀, where A is an n×n-matrix and h a nonlinear vector function with h(t, u)=o(∥u∥). As the first novel point, a lower bound for the asymptotic behaviour on the solution x(t) is derived. Two methods are applied to determine the optimal two-sided bounds, where one of the methods is the differential calculus of norms. In this context, the second novel point enters; it consists of a new strategy to significantly reduce the computation time for the determination of the optimal constants in the two-sided bounds. The obtained results are especially of interest in engineering and cannot be obtained by the methods used so far.     

On the Autoreducibility of Functions

This paper studies the notions of self-reducibility and autoreducibility. Our main result regarding length-decreasing self-reducibility is that any complexity class C\mathcal{C} that has a (logspace) complete language and is closed under polynomial-time (logspace) padding has the property that if all C\mathcal{C} -complete languages are length-decreasing (logspace) self-reducible then C í P\mathcal{C}\subseteq \mathrm {P} (C í L\mathcal {C}\subseteq \mathrm {L} ). In particular, this result applies to NL, NP and PSPACE. We also prove an equivalent of this theorem for function classes (for example, for #P). We also show that for several hard function classes, in particular for #P, it is the case that all their complete functions are deterministically autoreducible. In particular, we show the following result. Let f be a #P parsimonious function with two preimages of 0. We show that there are two FP functions h and t such that for all inputs x we have f(x)=t(x)+f(h(x)), h(x)≠x, and t(x)∈{0,1}. Our results regarding single-query autoreducibility of #P functions can be contrasted with random self-reducibility for which it is known that if a #P complete function were random self-reducible with one query then the polynomial hierarchy would collapse.     

Improved upper bounds on the L(2,1) -labeling of the skew and converse skew product graphs

An L(2,1)-labeling of a graph G is a function f from the vertex set V(G) to the set of all nonnegative integers such that |f(x)−f(y)|≥2 if d(x,y)=1 and |f(x)−f(y)|≥1 if d(x,y)=2, where d(x,y) denotes the distance between x and y in G. The L(2,1)-labeling number λ(G) of G is the smallest number k such that G has an L(2,1)-labeling with max{f(v):vV(G)}=k. Griggs and Yeh conjecture that λ(G)≤Δ² for any simple graph with maximum degree Δ≥2. This paper considers the graph formed by the skew product and the converse skew product of two graphs with a new approach on the analysis of adjacency matrices of the graphs as in [W.C. Shiu, Z. Shao, K.K. Poon, D. Zhang, A new approach to the L(2,1)-labeling of some products of graphs, IEEE Trans. Circuits Syst. II: Express Briefs (to appear)] and improves the previous upper bounds significantly.     

Oscillation criteria of certain nonlinear partial difference equations

This paper is concerned with the nonlinear partial difference equation with continuous variables

,where a, σ_i, τ_i are positive numbers, h_i(x, y, u) ε C(R⁺ × R⁺ × R, R), uh_i(x, y, u) > 0 for u ≠ 0, h_i is nondecreasing in u, i = 1, …, m. Some oscillation criteria of this equation are obtained.     

Functional expansions for nonlinear discrete-time systems

Given an input-output map associated with a nonlinear discrete-time state equationx(t + 1) =f(x(t);u(t)) and a nonlinear outputy(t) =h(x(t)), we present a method for obtaining a “discrete Volterra series” representation of the outputy(t) in terms of the controlsu(0), ...,u(t − 1). The proof is based on Taylor-type expansions of the iterated composition of analytic functions. It allows us to make an explicit construction of each kernel, that is, each coefficient of the series expansion ofy(t) in powers of the controls. This is achieved by making use of successive directional derivatives associated with a family of vector fields which are deduced from the discrete state equations. We discuss the use of these vector fields for the analysis and control of nonlinear discrete-time systems. This work was carried out while D. Normand-Cyrot was working at the I.A.S.I. (from March to October 1984) and with the financial support of the Italian C.N.R. (Consiglio Nazionale delle Ricerche).     

New explicit filters and smoothers for diffusions with nonlinear drift and measurements

The optimal least-squares filtering of a diffusion x(t) from its noisy measurements {y(τ); 0 τ t} is given by the conditional mean E[x(t)|y(τ); 0 τ t]. When x(t) satisfies the stochastic diffusion equation dx(t) = f(x(t)) dt + dw(t) and y(t) = ∫₀^tx(s) ds + b(t), where f(·) is a global solution of the Riccati equation /xf(x) + f(x)² = f(x)² = αx² + βx + γ, for some , and w(·), b(·) are independent Brownian motions, Benes gave an explicit formula for computing the conditional mean. This paper extends Benes results to measurements y(t) = ∫₀^tx(s) ds + ∫₀^t dx(s) + b(t) (and its multidimensional version) without imposing additional conditions on f(·). Analogous results are also derived for the optimal least-squares smoothed estimate E[x(s)|y(τ); 0 τ t], s < t. The methodology relies on Girsanov's measure transformations, gauge transformations, function space integrations, Lie algebras, and the Duncan-Mortensen-Zakai equation.     

Nonlinear eigenvalue problems for quasilinear systems

The paper deals with the existence of positive solutions for the quasilinear system (Φ(u'))' + λh(t)f(u) = 0,0 < t < 1 with the boundary condition u(0) = u(1) = 0. The vector-valued function Φ is defined by Φ(u) = (q(t)(p(t)u₁), …, q(t)(p(t)u_n)), where u = (u₁, …, u_n), andcovers the two important cases (u) = u and (u) = up > 1, h(t) = diag[h₁(t), …, h_n(t)] and f(u) = (f¹(u), …, fⁿ (u)). Assume that fⁱ and h_i are nonnegative continuous. For u = (u₁, …, u_n), let

, f₀ = maxf¹₀, …, fⁿ₀ and f_∞ = maxf¹_∞, …, fⁿ_∞. We prove that the boundary value problem has a positive solution, for certain finite intervals of λ, if one of f₀ and f_∞ is large enough and the other one is small enough. Our methods employ fixed-point theorem in a cone.     

Filtering of some nonlinear diffusions satisfying the general Bene condition

Under some regularity assumptions and the following generalization of the well-known Bene condition [1]:

Full-size image

, where F(t,z) = g⁻²(t)∫f(t,z)dz, F_t, F_z, F_zz, are partial derivatives of F, we obtain explicit formulas for the unnormalized conditional density q_t(z, x) α Px_t ε dz| y_s, 0 st, where diffusion x_t on R¹ solves x₀ = x, dx_t = [β(t) + α(t)x_t + f(t, x_t] dt + g(t) dw₁, and observation y_t = ∫_o^th(s)x_s ds + ∫_o^t(s) dw_2t, with w = (w₁, w₂) a two-dimensional Wiener process.     

Global state and feedback equivalence of nonlinear systems

The problem of finding global state space transformations and global feedback of the form u(t)= α(x) + ν(t) to transform a given nonlinear system to a controllable linear system on Rⁿ or on an open subset of Rⁿ, is considered here. We give a complete set of differential geometric conditions which are equivalent to the existence of a solution to the above problem.     

Exact finite-dimensional filters via systems realization for a class of discrete-time nonlinear systems

We obtain necessary and sufficient conditions for the existence of a finite-dimensional filter for the discrete-time nonlinear system (ε x_k+1 =φ(x_k), y_k = h(x_k)+η(x_k)ν_k, K=0, 1,…. This system is distinguished by the absence of noise in the dynamic and by the correlation between the state and the intensity of noise in the observations.The necessary and sufficient condition provides an explicit formula for the minimal filter and various system-theoretic properties of (ε) and of the minimal filter.     

On the non-existence of stationary diffusions which satisfy the Bene condition

The one-dimensional diffusion x_t satisfying dx_t = f(x_t)dt + dw_t, where w_t is a standard Brownian motion and f(x) satisfies the Bene condition f′(x) + f²(x) = ax² + bx + c for all real x, is considered. It is shown that this diffusion does not admit a stationary probability measure except for the linear case f(x) = αx + β, α < 0.     

The existence of positive solutions of singular fractional boundary value problems

We discuss the existence of positive solutions for the singular fractional boundary value problem D^αu+f(t,u,u^′,D^μu)=0, u(0)=0, u^′(0)=u^′(1)=0, where 2<α<3, 0<μ<1. Here D^α is the standard Riemann-Liouville fractional derivative of order α, f is a Carathéodory function and f(t,x,y,z) is singular at the value 0 of its arguments x,y,z.     

On the joint nonlinear filtering-smoothing of diffusion processes

The smoothing of diffusions dx_t = f(x_t) dt + σ(x_t) dw_t, measured by a noisy sensor dy_t = h(x_t) dt + dv_t, where w_t and v_t are independent Wiener processes, is considered in this paper. By focussing our attention on the joint p.d.f. of (x_τ x_t), 0 ≤ τ < t, conditioned on the observation path {y_s, 0 ≤ s ≤ t}, the smoothing problem is represented as a solution of an appropriate joint filtering problem of the process, together with its random initial conditions. The filtering problem thus obtained possesses a solution represented by a Zakai-type forward equation. This solution of the smoothing problem differs from the common approach where, by concentrating on the conditional p.d.f. of x_τ alone, a set of ‘forward and reverse’ equations needs to be solved.     

A new unconditionally stable ADI compact scheme for the two-space-dimensional linear hyperbolic equation

We formulate a new alternating direction implicit compact scheme of O(τ²+h ⁴) for the linear hyperbolic equation u _tt +2α u _t +β² u=u _xx +u _yy +f(x, y, t), 0<x, y<1, 0<t≤T, subject to appropriate initial and Dirichlet boundary conditions, where α>0 and β≥0 are real numbers. In this article, we show the method is unconditionally stable by the Von Neumann method. At last, numerical demonstrations are given to illustrate our result.     

Universal stabilization of a class of nonlinear systems with homogeneous vector fields

Under relative-degree-one and minimum-phase assumptions, it is well known that the class of finite-dimensional, linear, single-input (u), single-output (y) systems (A,b,c) is universally stabilized by the feedback strategy u = Λ(λ)y, λ = y², where Λ is a function of Nussbaum type (the terminology “universal stabilization” being used in the sense of rendering /s(0/s) a global attractor for each member of the underlying class whilst assuring boundedness of the function λ(·)). A natural generalization of this result to a class ^k of nonlinear control systems (a,b,c), with positively homogeneous (of degree k 1) drift vector field a, is described. Specifically, under the relative-degree-one (cb ≠ 0) and minimum-phase hypotheses (the latter being interpreted as that of asymptotic stability of the equilibrium of the “zero dynamics”), it is shown that the strategy u = Λ(λ)/vby/vb^k−1y, assures ^k-universal stabilization. More generally, the strategy u = Λ(λ)exp(/vby/vb)y, assures -universal stabilization, where = _{k 1} ^k.     

Stochastic stabilisation of functional differential equations

In this paper we investigate the problem of stochastic stabilisation for a general nonlinear functional differential equation. Given an unstable functional differential equation dx(t)/dt=f(t,x_t), we stochastically perturb it into a stochastic functional differential equation , where Σ is a matrix and B(t) a Brownian motion while X_t={X(t+θ):-τθ0}. Under the condition that f satisfies the local Lipschitz condition and obeys the one-side linear bound, we show that if the time lag τ is sufficiently small, there are many matrices Σ for which the stochastic functional differential equation is almost surely exponentially stable while the corresponding functional differential equation dx(t)/dt=f(t,x_t) may be unstable.     

The asymptotic self-similar behavior for the quasilinear heat equation with nonlinear boundary condition

In this paper the quasilinear heat equation with the nonlinear boundary condition is studied. The blow-up rate and existence of a self-similar solution are obtained. It is proved that the rescaled function

v(y,t)=(T−t)^{1/(2p+α−2)}u((T−t)^{(p−1)/(2p+α−2)}y,t),