共查询到20条相似文献,搜索用时 31 毫秒
1.
《国际计算机数学杂志》2012,89(7):932-943
In this paper, a split least-squares characteristic mixed finite element method is proposed for solving nonlinear nonstationary convection–diffusion problem. By selecting the least-squares functional property, the resulting least-squares procedure can be split into two independent symmetric positive definite sub-schemes. The first sub-scheme is for the unknown variable u, which is the same as the standard characteristic Galerkin finite element approximation. The second sub-scheme is for the unknown flux σ. Theoretical analysis shows that the method yields the approximate solutions with optimal accuracy in L 2(Ω) norm for the primal unknown and in H(div; Ω) norm for the unknown flux, respectively. Some numerical examples are given to confirm our theory results. 相似文献
2.
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) = ∫0tx(s) ds + b(t), where f(·) is a global solution of the Riccati equation /xf(x) + f(x)2 = f(x)2 = αx2 + β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) = ∫0tx(s) ds + ∫0t 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. 相似文献
3.
Jue Yan 《Journal of scientific computing》2013,54(2-3):663-683
We propose a discontinuous Galerkin finite element method for convection diffusion equations that involves a new methodology handling the diffusion term. Test function derivative numerical flux term is introduced in the scheme formulation to balance the solution derivative numerical flux term. The scheme has a nonsymmetric structure. For general nonlinear diffusion equations, nonlinear stability of the numerical solution is obtained. Optimal kth order error estimate under energy norm is proved for linear diffusion problems with piecewise P k polynomial approximations. Numerical examples under one-dimensional and two-dimensional settings are carried out. Optimal (k+1)th order of accuracy with P k polynomial approximations is obtained on uniform and nonuniform meshes. Compared to the Baumann-Oden method and the NIPG method, the optimal convergence is recovered for even order P k polynomial approximations. 相似文献
4.
We present a study of the local discontinuous Galerkin method for transient convection–diffusion problems in one dimension. We show that p-degree piecewise polynomial discontinuous finite element solutions of convection-dominated problems are O(Δx
p+2) superconvergent at Radau points. For diffusion- dominated problems, the solution’s derivative is O(Δx
p+2) superconvergent at the roots of the derivative of Radau polynomial of degree p+1. Using these results, we construct several asymptotically exact a posteriori finite element error estimates. Computational results reveal that the error estimates are asymptotically exact.This revised version was published online in July 2005 with corrected volume and issue numbers. 相似文献
5.
In this paper, we derive a high-order compact finite difference scheme for solving the reaction–subdiffusion equation with Neumann boundary value condition. The L1 method is used to approximate the temporal Caputo derivative, and the compact difference operator is applied for spatial discretization. We prove that the compact finite difference method is unconditionally stable and convergent with order O(τ2?α+h4) in L2 norm, where τ, α, and h are the temporal step size, the order of time fractional derivative and the spatial step size, respectively. Finally, some numerical experiments are carried out to show the effectiveness of the proposed difference scheme. 相似文献
6.
In this paper, we study the sample complexity of weak learning. That is, we ask how many data must be collected from an unknown distribution in order to extract a small but significant advantage in prediction. We show that it is important to distinguish between those learning algorithms that output deterministic hypotheses and those that output randomized hypotheses. We prove that in the weak learning model, any algorithm using deterministic hypotheses to weakly learn a class of Vapnik-Chervonenkis dimension d(n) requires Ω ([formula]) examples. In contrast, when randomized hypotheses are allowed, we show that Θ (1) examples suffice in some cases. We then show that there exists an efficient algorithm using deterministic hypotheses that weakly learns against any distribution on a set of size d(n) with only O(d(n)2/3) examples. Thus for the class of symmetric Boolean functions over n variables, where the strong learning sample complexity is Θ (n), the sample complexity for weak learning using deterministic hypotheses is Ω ([formula]) and O(n2/3), and the sample complexity for weak learning using randomized hypotheses is Θ (1). Next we prove the existence of classes for which the distribution-free sample size required to obtain a slight advantage in prediction over random guessing is essentially equal to that required to obtain arbitrary accuracy. Finally, for a class of small circuits, namely all parity functions of subsets of n Boolean variables, we prove a weak learning sample complexity of Θ(n). This bound holds even if the weak learning algorithm is allowed to replace random sampling with membership queries, and the target distribution is uniform on {0, 1}n. 相似文献
7.
Let f(xθ) = αθαx−(α+1)I(x>θ) be the pdf of a Pareto distribution with known shape parameter α>0, and unknown scale parameter θ. Let {(Xi, θi)} be a sequence of independent random pairs, where Xi's are independent with pdf f(xαi), and θi are iid according to an unknown distribution G in a class
of distributions whose supports are included in an interval (0, m), where m is a positive finite number. Under some assumption on the class
and squared error loss, at (n + 1)th stage we construct a sequence of empirical Bayes estimators of θn+1 based on the past n independent observations X1,…, Xn and the present observation Xn+1. This empirical Bayes estimator is shown to be asymptotically optimal with rate of convergence O(n−1/2). It is also exhibited that this convergence rate cannot be improved beyond n−1/2 for the priors in class
. 相似文献
8.
《国际计算机数学杂志》2012,89(8):1644-1663
A fully discrete two-grid finite-volume method (FVM) for a nonlinear parabolic problem is studied in this paper. This method involves solving a nonlinear parabolic equation on coarse mesh space and a linearized parabolic equation on fine grid. Both L 2 and H 1 norm error estimates of the standard FVM for the nonlinear parabolic problem are derived. Compared with the standard FVM, the two-level method is of the same order as the one-level method in the H 1-norm as long as the mesh sizes satisfy h=𝒪(H 3/2). However, the two-level method involves much less work than the standard method. Numerical results are provided to demonstrate the effectiveness of our algorithm. 相似文献
9.
Stability of differential inclusions
F(x(t)) is studied by using minorant and majorant mappings F
– and F
+, F
–(x)F(x)F
+(x). Properties of F
–,F
+ are developed in terms of partial orderings, with the condition that F
–, F
+ are either heterotone or pseudoconcave. The main results concern asymptotically stable absorbing sets, including the case of a single equilibrium point, and are illustrated by examples of control systems. 相似文献
10.
Eric Bach 《Information and Computation》1991,90(2)
Pollard's “rho” method for integer factorization iterates a simple polynomial map and produces a nontrivial divisor of n when two such iterates agree modulo this divisor. Experience and heuristic arguments suggest that a prime divisor p should be detected in
steps, but this has never been proved. Indeed, nothing seems to be have been rigorously proved about the probability of success that would improve the obvious lower bound of 1/p. This paper shows that for fixed k, this probability is at least (2k)/p + O(p−3/2) as p → ∞. This leads to an Ω(log2p)/p estimate of the success probability. 相似文献
11.
Baruch Awerbuch Margrit Betke Ronald L. Rivest Mona Singh 《Information and Computation》1999,152(2):155
We study how a mobile robot can learn an unknown environment in a piecemeal manner. The robot's goal is to learn a complete map of its environment, while satisfying the constraint that it must return every so often to its starting position (for refueling, say). The environment is modeled as an arbitrary, undirected graph, which is initially unknown to the robot. We assume that the robot can distinguish vertices and edges that it has already explored. We present a surprisingly efficient algorithm for piecemeal learning an unknown undirected graph G=(V, E) in which the robot explores every vertex and edge in the graph by traversing at most O(E+V1+o(1)) edges. This nearly linear algorithm improves on the best previous algorithm, in which the robot traverses at most O(E+V2) edges. We also give an application of piecemeal learning to the problem of searching a graph for a “treasure.” 相似文献
12.
Beg M. Faisal Miller Michael I. Trouvé Alain Younes Laurent 《International Journal of Computer Vision》2005,61(2):139-157
This paper examine the Euler-Lagrange equations for the solution of the large deformation diffeomorphic metric mapping problem studied in Dupuis et al. (1998) and Trouvé (1995) in which two images I
0, I
1 are given and connected via the diffeomorphic change of coordinates I
0○ϕ−1=I
1 where ϕ=Φ1 is the end point at t= 1 of curve Φ
t
, t∈[0, 1] satisfying .Φ
t
=v
t (Φ
t
), t∈ [0,1] with Φ0=id. The variational problem takes the form
where ‖v
t‖
V
is an appropriate Sobolev norm on the velocity field v
t(·), and the second term enforces matching of the images with ‖·‖L
2 representing the squared-error norm.In this paper we derive the Euler-Lagrange equations characterizing the minimizing vector fields v
t, t∈[0, 1] assuming sufficient smoothness of the norm to guarantee existence of solutions in the space of diffeomorphisms. We describe the implementation of the Euler equations using semi-lagrangian method of computing particle flows and show the solutions for various examples. As well, we compute the metric distance on several anatomical configurations as measured by ∫0
1‖v
t‖
V
dt on the geodesic shortest paths. 相似文献
13.
Plant template generation is the key step in applying quantitative feedback theory (QFT) to design robust control for uncertain systems. In this paper we propose a technique for generating plant templates for a class of linear systems with an uncertain time delay and affine parameter perturbations in coefficients. The main contribution lies in presenting a necessary and sufficient condition for the zero inclusion of the value set f(T,Q)={f(τ,q): τT[τ−,τ+], qQ∏k=0m−1[qk−,qk+]}, where f(τ,q)=g(q)+h(q)e−jτω*, g(q) and h(q) are both complex-valued affine functions of the m-dimensional real vector q, and ω* is a fixed frequency. Based on this condition, an efficient algorithm which involves, in the worst case, evaluation of m algebraic inequalities and solution of m2m−1 one-variable quadratic equations, is developed for testing the zero inclusion of the value set f(T,Q). This zero-inclusion test algorithm allows one to utilize a pivoting procedure to generate the outer boundary of a plant template with a prescribed accuracy or resolution. The proposed template generation technique has a linear computational complexity in resolution and is, therefore, more efficient than the parameter gridding and interval methods. A numerical example illustrating the proposed technique and its computational superiority over the interval method is included. 相似文献
14.
This paper first shows how the Bézier coefficients of a given degree n polynomial are perturbed so that it can be reduced to a degree m (<n) polynomial with the constraint that continuity of a prescribed order is preserved at the two endpoints. The perturbation vector, which consists of the perturbation coefficients, is determined by minimizing a weighted Euclidean norm. The optimal degree n−1 approximation polynomial is explicitly given in Bézier form. Next the paper proves that the problem of finding a best L2-approximation over the interval [0,1] for constrained degree reduction is equivalent to that of finding a minimum perturbation vector in a certain weighted Euclidean norm. The relevant weights are derived. This result is applied to computing the optimal constrained degree reduction of parametric Bézier curves in the L2-norm. 相似文献
15.
For Ω a smooth domain in Rn with boundary Λ = Λ0Λ1, we are concerned with the wave equation y″ − Δy = S in QT =]0, T[ × Ω with = ∂/∂t, at source term satisfying S, S′ ε L1(0, T L2 (Ω)). A Dirichlet condition is imposed on Λ0 and we consider an absorbing condition ∂y/∂n + uy′ = 0 in [0, T] × gL1 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. 相似文献
16.
17.
Shunsheng Guo Hongbiao Jiang Qiulan Qi 《Computers & Mathematics with Applications》2007,54(11-12):1387-1394
In this paper, we give direct, inverse and equivalence approximation theorems for the Bézier type of Meyer–König and Zeller operator with unified Ditzian–Totik modulus ωφλ(f,t) (0≤λ≤1). 相似文献
18.
《国际计算机数学杂志》2012,89(4):445-463
An algorithm for approximating a non negative solution of inconsistent systems of linear equations is presented. We define a best approximate solution of a system Ax = b x≥0 to be the vector x≥0 which minimizes the norm of the residual r(x) = b ? Ax, for a smooth and strictly convex norm. The algorithm is shown to be feasible and globally convergent. The special case of the ? p norm is included. In particular, the method converges for 1 < p < 2. A generalization of this algorithm is also given. Numerical results are included. 相似文献
19.
Existence of triple positive solutions of two-point right focal boundary value problems on time scales 总被引:2,自引:0,他引:2
We consider the following boundary value problem, (−1)n−1yΔn(t)=(−1)p+1F(t,y(σn−1(t))),t[a,b]∩T, yΔn(a)=0,0≤i≤p−1, yΔn(σ(b))=0,p≤i≤n−1,where n ≥ 2, 1 ≤ p ≤ n - 1 is fixed and T is a time scale. By applying fixed-point theorems for operators on a cone, existence criteria are developed for triple positive solutions of the boundary value problem. We also include examples to illustrate the usefulness of the results obtained. 相似文献
20.
J.-P. Croisille 《Computing》2006,78(4):329-353
We introduce a new box-scheme, called ``hermitian box-scheme' on the model of the one-dimensional Poisson problem. The scheme
combines features of the box-scheme of Keller, [20], [13], with the hermitian approximation of the gradient on a compact stencil,
which is characteristic of compact schemes, [9], [21]. The resulting scheme is proved to be 4th order accurate for the primitive
unknown u and its gradient p. The proved convergence rate is 1.5 for (u,p) in the discrete L
2 norm. The connection with a non standard mixed finite element method is given. Finally, numerical results are displayed on
pertinent 1-D elliptic problems with high contrasts in the ellipticity, showing in practice convergence rates ranging from
1 to 2.5 in the discrete H
1 norm.
This work has been performed with the support of the GDR MOMAS, (ANDRA, CEA, EDF, BRGM and CNRS): Modélisation pour le stockage des déchets radioactifs. The author thanks especially A. Bourgeat for his encouragements and his interest in this work. 相似文献