Abstract Inverse problems for elliptic and parabolic partial dierential equations are considered. It is assumed that a solution of the equation is obs
|
|
- Lorraine Morton
- 6 years ago
- Views:
Transcription
1 Statistical approach to inverse boundary problems for partial dierential equations Golubev G. and Khasminskii R. eierstra{institut fur Angewandte Analysis and Stochastik Mohrenstae 39 D{7, Berlin, Deutchland ayne State University Detroit, MI 48, USA This work was partially supported NSF grant DMS 9645 and ONR grant N Mathematics Subject Classications. 6G5. Key words and phrases. Partial dierential equation, minimax estimator, boundary problem.
2 Abstract Inverse problems for elliptic and parabolic partial dierential equations are considered. It is assumed that a solution of the equation is observed in white Gaussian noise with a small spectral density. The goal is to recover smooth but unknown boundary or initial conditions based on the noisy data. It is shown that the second order minimax estimators are linear as the spectral density of the noise goes to zero. Introduction Ill-posed inverse problems for ordinary and partial dierential equations are very popular in modern mathematics. e refer to the books Hadamard (93), Ivanov & Vasin & Tanana (978), Tihonov & Arsenin (978), avrentiev & Romanov & Shishatskii (98). Usually these problems are associated with recovering of an unknown dierential operator or unknown boundary conditions based on the observations of a solution of the equation in some domain. There are two important questions, which arise naturally in this context. In the rst place one wants to prove the uniqueness of a solution of the inverse boundary problem if it is known that the restored boundary functions belong to a certain functional class. If the solution represents a noisy data, one has to show next that the error in the recovered functions tends to zero when the noise in the data tends to zero too, attes & ions (967). Statistical approach to inverse problems is based on the assumption that solutions are observed in a white noise with a small spectral density. The goal is to nd estimators which are at least consistent. The main advantage of the statistical approach lies in the fact that it helps to compare dierent approaches to inverse problems. Thus we have a mathematical basis for choosing solutions, which are optimal in dierent senses. Sudakov & Khaln (964) and Khaln (978) were among the rst who propose to use the statistical approach to ill-posed problems. In Chow & Khasminskii (997a) and Chow & Ibragimov & Khasminskii (997b) this approach was used for estimation of a source based on solutions of ordinary and partial dierential equations. In the present paper we consider from statistical point ofviewtwo very popular inverse problems i) recovering of initial conditions for parabolic equations based on the observations in a xed time strip, ii) recovering of boundary conditions for elliptic equations based on the observations in an internal domain. Many scientists have paid attention on these problems. For instance, the major part of the monograph attes & ions (967) is devoted to the problem i). From statistical point of view the authors have proposed very interesting estimates and demonstrated that they are consistent. e show in this paper how to construct estimates with optimal rates of convergence when the spectral density of the noise goes to zero ( ). Moreover
3 we indicate these rates up to constants. Here we follow the well-known paper Pinsker (98). Unfortunately these rates are very low and there exist many estimators having the optimal rates. In order to discriminate between the estimators in this situation we nd the second order term in the expansion of the minimax risk and propose second order minimax estimators. In the problem of estimation initial conditions for parabolic equations we propose also exponentially minimax estimators.. Estimation of initial conditions for parabolic equations To simplify technical details we consider the problem in the simplest setting. et u(t x) bea solution of the heat u in the domain D = f x gft>g with the periodic boundary conditions and with the initial condition u(t ) = @u (t ) = (t ) u( x)=(x) (3) Identifying points x = and x =, one could say that the equation () is considered on the surface of the cylinder Cft>g, where C is the circle of the unit length. It is also assumed that (x) x Chas Sobolev smoothness, so where j = ;j and (x) = j=; j=; j exp(ijx) i= p ; (4) j j j jjj (5) The set of all function satisfying to (4) and (5) constitutes the Sobolev class (). Next we consider two models of observations. Model a. Assume that a solution of the equation () is observed at t = T in a white Gaussian noise with the spectral density. Thus we observe the generalized random process Y (x) x [ ], which has the form Y (x) =u(t x)n(x) (6) where n(x) is a white Gaussian noise, t.e. a generalized Gaussian process with the covariance function En(x)n(y) = (x ; y) here () is the Dirac -function. It means that for any ' () (C) we can observe the family of random variables Y = Z C u(t x)' (x)dx (7)
4 where are Gaussian random variables with zero mean and the covariance function E = Thus, our problem is to nd the minimax risk r ( ()) = inf ~ Z C ' (x)' (x)dx sup () E ZC ~(x Y ) ; (x) dx where the inf is taken over all estimators based on observations (6). e want also to nd a minimax estimator i.e. the estimator such that r ( ()) = sup E ( ZC (x Y ) ; (x)) dx () The above problem can be easily reduced to the following one. et u(t x) be a solution of (){(3). Since this function is periodic, it can be expanded into the Fourier series u(t x) = j=; u j (t)exp(ijx) (8) Using the Fourier method (cf. Petrovskii (95), s. 38) one obtains u j (t) = j exp(;4 j t) (9) where j are the Fourier coecients associated with the function (x) (see (4)). Since the Fourier basis is a complete orthonormal system in (C), the problem of estimation (x) is equivalent to estimation of the Fourier coecients j j = based on the observations Y j = j exp(;4 j T ) j () where j are independent N ( ). Or equivalently, we have to estimate j based on the data Z j = j j j () where j =exp(4 j T )and j belong to the following set = 8 < j j=; 9 = j j j jjj () Note the minimax risk admits the following representation in the terms of the Fourier coecients r ( ()) = inf ~ sup E j=; ~j ; j (3) Model b. It is assumed that the solution of () is observed into the strip T t T h inatwo-dimensional white Gaussian noise Y (t x) =u(t x)n(t x) t [T T h] x C 3
5 where n(t x) is a generalized two-dimensional Gaussian eld with the covariance function En(t x )n(t x )=(t ; t x ; x ) By (8) and (9) we get the equivalent model (see also (7)) where ~Z j = j ~ j j (4) ~ j = Z T h T exp(;8 j t)dt ;= = p jjj exp(4 j T ) ; exp(;8 j h) ;= Based on the observation from (4) we have to nd the minimax estimator and to calculate its risk (cf. (3)). Thus we see that the only dierence between models a and b is in the denition of j and ~ j.. Estimation of boundary conditions for elliptic equations Consider the simplest case of the Dirichlet problem for the aplace equation on the circle of radius u = (5) u(cos ' sin ') =f(') It is known (see e.g. Petrovskii (95), s. 9) that the solution of the above problem can be rewritten in the polar coordinates as u(r ')= p p k= where k are the Fourier coecients of f('), so that f(') = p p k= r k ( k cos(k') ;k sin(k')) (6) ( k cos(k') ;k sin(k')) et us assume that a solution of (5) is observed on the circle C of the radius < in a white Gaussian noise. Thus we observe a generalized Gaussian eld Y (r ')=ru(r ') p rn(r ') r ' [ ] (7) where n(r ') is a white Gaussian noise i.e. a generalized Gaussian eld with the covariance function En(r ' )n(r ' )=(r ; r ' ; ' ) 4
6 Substituting (6) in (7), multiplying it by cos(k') sin(k') andintegrating by parts over ', we arrive to the equivalent family of the observations (7) k (r) = k r jkj= n k (r) k = (8) where r and n k (r) are one-dimensional independent (for dierent k) white Gaussian noises. Thus the problem of estimation k based on (8) is equivalent to estimation of k based on the observations Y k = Z r jkj= k (r)dr Z r jkj dr Z ;= = k r dr jkj k where k are independent N ( ). Finally we arrive to the following equivalent problem to estimate k based on the data where k are i.i.d. N ( ) and Y k = k k k k = k = q jkj ;jkj; e see that each problem under consideration is equivalent to the problem of estimation unknown parameters k based on the observations (), while the prior information is provided by (). The type of the considered problem is reected only in the denition of k. For parabolic equations k grow like exp(ck ). In the elliptic case we have a growth of the order exp(ck). This dierence plays an essential role only when we prove that second order minimax estimators are linear for the both models. inear estimation. An upper bound for the minimax risk Consider the following linear estimator ^ j = h j Y j of the parameters j based on the observations (). It is assumed j are subjected to restrictions (). The problem of calculation of the linear minimax risk r ( ())=inf h j sup E has the well-known solution, Pinsker (98). j=; (h j Y j ; j ) 5
7 Proposition r ( ()) = j=; j j=; where is a root of the equation j jjj 4 j ; ; j 3 5 # (9) = () The estimator j = ; j # Y j () is the minimax linear estimator. Proof follows directly from the fact that the functional [h ] = j=; ( ; h j ) j j=; j h j = E j=; (h j Y j ; j ) has a saddle point on l (; ) (see e.g. Pinsker (98)). The components of this saddle point are given by h j = ; j () where [x] = max( x). j = j 4 j # ;. A lower bound for the minimax risk 3 5 (3) To show that the linear estimator dened by (), () is a rst order minimax estimator we will use the following fact. Proposition et ( j ) are dened by (3). Then r ( ()) (j ) ; 4 ( j ) 4 j ; (4) j6= j6= Proof. et j be independent random variables taking values j ; j with the probabilities =. It is clear that r ( ()) j6= inf EEf( ~ j ; j ) j j g (5) ~ j 6
8 It is also evident that the inf in the right-hand side of (5) is attained when ~ j is the Bayesian estimator Y ~ j = j j j th j where Then we have th x = exp(x) ; exp(;x) exp(x) exp(;x) EEf( ~ j ; j ) j j g = Ef( ~ j ; j ) j j = j g (6) 8 < = (j ) E exp Y jj j 8 < ; 9 = ; j = j 9 = ; j = j (j ) E ; exp Y jj j 9 = (j ) 4 ; 3 exp 4( j ) 4 exp ( j ) j j ( j ) 5 ; 3 exp 4( j ) j Since 5 ; 3 exp(x) is a convex function, one easily obtains that for x [ log(5=3)]. Note also that if 5 ; 3exp(x) ; x= log(5=3) (7) 4( j ) j > log(5=3) the right-hand side in (6) is negative. Then (7) and (6) yield EEf( ~ j ; j ) j j g( j ) ; 4( j ) j thus, by virtue of (5), proving the required inequality (4)..3 Asymptotic behaviour of the minimax risk (8) At the rst glance it seems that the lower bound given by Proposition is not very good. This observation is true in the case when j grow not very rapidly. But in our models we have at least an exponential growth of j, and this fact plays an essential role in the proving of the following theorem. 7
9 Theorem et then j j j j j j > j (9) r ( ()) = r ( ()) O (3) where the bandwidth is dened by (). Proof. et us simplify the lower bound for the minimax risk (4). First of all note that j j6=( ) = j ; j # j (3) j6= = j ; j # j ; ; j6= j ; j ; j ; j # j6= j6= In order to estimate from above the last term in the right-hand side of the above equation we use (9) and the inequality ; ( ; x) max( )x, which is valid for x. Then we have j6= j ; j # b c = k= b c;k # 4 ; b c;k ( ) ; b c; k= 3 5 ;(k;) k = C b c; ; (3) On the other hand it is evident that r ( ()) = j ; j j=; b c; ; b c; A O b c; ; (33) Therefore from (3) and (3) we see that j6= ( j ) r ( ()) O (34) Thus to complete the proof it remains to consider the last term in the right-hand side of (4). e have ; (j ) 4 j ; = j ; j # j (35) j6= j6= 8
10 = <jjj= j ; j # <jjj= j j6= j j ; j jjj>= # j ; j # j From (9) we arrive at <jjj= Hence from (3), (33) and (35) it follows that j ( ; ) ; b c; b c ; (j ) 4 j ; O j6= ; r ( ()) This inequality together with (34) completes the proof of the theorem. Our next goal is to specify the asymptotic behaviour of r ( ()) as. Proposition 3 et (9) is fullled, then r ( ()) = where is a root of equation ( ) O (36) ; b c = C (37) here bxc means integer part of x and C is an arbitrary positive constant not depending on. Proof. A simple algebra easily reveals (see (9), ()) that r ( ()) = ; () j j j=; j j=; where is dened by (). At the same time by virtue of (3) and (33) ; j # = O ; r ( ()) # (38) Hence we get by (38) r ( ()) = () O (39) 9
11 Thus it remains to obtain a lower and an upper bound for dened by (). By () = j ; j j jjj j=; # = () j j ; j j b c = () b c;j(b c;j) j= b c;j ; A By Taylor expansion and (9) we get from the above equation that for some constants C C do not depending on () j j ; j j () j j ; j j # # C ; b c C ; b c; Therefore is located between two roots and of the equations ; b c C = ; b c; C = (4) Taking the logarithm in (4) one concludes that j ; j < C. This inequality together with (39) completes the proof. The following theorem easily follows from Proposition 3 and Theorem. Theorem et (9) is fullled and is a root of equation (37). Then r ( ()) = ( ) O et us look how this theorem works in the case of elliptic equation. Remind that in this case k =(jkj )e ;(jkj) log Taking the logarithm of (37) we get where Hence r ( ()) = ; log log(= ) = ~ O() (4) ~ = ; log log ; log log 4 ; log log(= ) log log(= ) log(= ) O log ;;
12 Thus Theorem gives us the expansion of the minimax risk up to the second order term. The linear estimator ~ j = 4 ; j ~ 3 5 Y j (4) with ~ given by (4), is the second order minimax estimator. Remark. The estimator (4) is robust with respect to the constant in the denition of the class () if the true value of is unknown, but it is known that < < then we can use the new bandwidth ~ = ; log log ~ ; log log ~ with an arbitrary ~ [ ], and the estimator (4) is again a second order minimax estimator. In the problem of estimation of the initial condition for the parabolic equation j have the following form (c.f. ()) j = exp(8 j T ) Therefore taking the logarithm of (37), we get = 8 T log ( ; ) ; log log O() (43) Hence, by Theorem r ( ()) = T O log ;;= log(= ) Thus the estimator from (4) with the bandwidth from (43) is a rst order minimax estimator. It is not very dicult to check that the projection estimator ( Zj jjj ^ ^ j = jjj > ^ with ^ = 8 T log ; 8 T log log where >is also the rst order minimax estimator. Remark. It is easy to see that the estimator (4) with the bandwidth dened by (43) and the projection estimator ^ j are also robust w.r.t. in the sense of Remark. Unfortunately Theorem gives us only the rst term in the expansion of the minimax risk. This situation is quite dierent from the elliptic case, where this theorem provides us with the second order term. Therefore, we consider in the next section the parabolic case in more detail.
13 3 Estimation of initial conditions for parabolic equations In this section we assume that j = ;j and log j ; log j = C j o(j) j (44) where C > is some constant. Dene the integer as follows = arg min j j (j) ; ; (45) Denote for brevity by bold letters two dimensional vectors, that is x = (x x ) T, and kxk = x x. The asymptotic behaviour of the minimax risk is given by the following theorem. Proposition 4 et be given by (45) and condition (44) be fullled. Then, as r ( n ()) = inf sup Ek (Y) ; yk (46) ( ) kyk ( ;kyk )) O (exp(;c )) where Y = y s N( E) s = ( ) p here C > is a certain constant and E is the identity matrix. Proof. et us obtain rst a lower bound. Assume that only =( ; ) T and = ( ;;) T are unknown. The others j are assumed to be. Then the observations have the form Z = (47) Z = Based on these observations we have to estimate and provided that k k ( ) k k (( )) (48) Introduce the new variables y = ( ) ;= y = (( )) ;= Then equation (48) can be rewritten as kyk ky k
14 and the observations (47) are equivalent to the following ones Y = y s (49) Y = y s ; ( ) ; where and are two-dimensional independent N ( E). Thus assuming that y = zq ;kyk, where z z are independent random variables taking values ;= p = p with the probabilities =, we have r ( ()) inf sup kyk ky k inf ( ( ) E k (Y Y ) ; yk (5) ) (( )) E k (Y Y ) ; y k sup kyk ( ( ) EE fw ( (Y Y ) ; y) jzg (( )) EE ( (Y Y ) ; z q )) ;kyk z where the new loss function w(x) = minfkxk 4g, and the observations Y Y are given by Y = y s (5) Y = z q ;kyk s ; ( ) ; To simplify the right-hand side of (5) rst of all note that by (8) and (44) ( q EE (Y Y) ) ; z ;kyk z ( ;kyk ) Hence from (5) we have r ( ()) inf sup kyk ; 4( ;kyk ) ;kyk s ( ) O(exp(;C )) ( ( ) EEfw ( (Y Y ) ; y) jzg (5) (( )) ;kyk ) O(exp(;C )) To simplify more the right-hand side of the above equation consider the auxiliary observations Y = s ; ( ) ; (53) It is easy to see that -distance between the corresponding densities of the observations Y Y and Y Y dened by (5) and (53), is suciently small. Indeed 3
15 by Cauchy-Schwartz inequality one obtains Z Z py (x Y x ) ; p (x Y Y x ) dx dx Z Z = py (x )p (x Y ) ; p Y (x )p (x Y ) dx dx Z = py (x ) ; p (x Y ) dx = Z p = Y Z (x ) ; p = p = Y Y (x ) ; p = (x ) Y p = Y (x ) dx = (x )p = Y (x ) dx Hence for any estimator () = p ; exp ; ( ;kyk ) 8s ( ) = exp(;c ) EEfw ( (Y Y ) ; y) j zg = EEfw ( (Y Y ) ; y) j zg O (exp(;c )) Since the inf in (5), which is taken over all mesurable functions (), can be replaced by the inf taken over k k, then using the above equation and (5) we arrive at the following lower bound r ( ()) inf sup ( ) kyk O (exp(;c )) (E k (Y) ; yk ( ) ;kyk ) ( ) In order to prove that this lower bound is attainable consider the following estimator 8 >< Y k jkj < ^ k = k(y k Y ;k ) jkj = > jkj > where k( ) is an arbitrary function R R. Note that by (44) and (45) Therefore ; 5( ) =5( ) ;e (o())c ; j= Denote for brevity = f j r ( (l)) = inf sup = 8 < j ( o()) ; =e ;(o())c P j infsup jjj< jjj= j (j) g. Then we get 8 < jjj= E( j(y j Y ;j ) ; j ) inf sup ( ) kyk j E( j(y j Y ;j ) ; j ) jjj> jjj= jjj= ) (Ek (Y) ; yk ( ;kyk ) ( ) 4 j 9 = O e ;C (54) 9 = O e ;C
16 Together with (54) this inequality completes the proof of the theorem. At the rst glance it seems that Proposition 4, which species the minimax risk expansion up to the exponential term, has very restrictive practical meaning since we have to solve the suciently complicated variational problem to calculate the vector-valued function ( ), at which the inf in (46) is attained. Fortunately, it can be done very easily and second order minimax estimators can be found explicitly. Moreover we will see that there exists a so-called exponentially ecient estimator, whose risk approximates the minimax risk up to terms of the order exp(;c ). Theorem 3 Assume that (44) is fullled and is given by (45), then r ( ()) = ( ) exp(;c ) and the following linear estimator ^ k = is exponentially ecient. 8 >< > Y s ; A (55) jkj < ; Y k jkj = jkj > Proof. According to Proposition 4 it suces to solve the following problem. Assume that we are given the noisy data Y = y s where is N ( E). Based on Y we have to estimate the unknown vector y. More precisely, our goal is to calculate the minimax risk inf sup kyk ( y) up to terms of the order O(exp(;C )), where ( y)=ek (Y) ; yk ( ;kyk ) To getan upper bound consider the linear estimator ^y = hy with h =; Then elementary algebra reveals that inf sup ( y) sup kyk kyk ( ( ; h) kyk h s ( ;kyk )) (56) = s ; 5
17 To get a lower bound assume that y is the Gaussian two-dimensional vector N ( E). Then we have inf sup ( y) inf kyk Pfkyk g E ( y)fkyk g (57) = inf k k Pfkyk g E ( y)fkyk g 5Pfkyk > g inf E ( y) ; Pfkyk g Since y is Gaussian the inf E ( y) is attained when et us chose = s ; (Y) = Noting that s < we have on the one hand Pfkyk > g C exp s y ; C exp ; s On the other hand noting that under such choice of and we obtain that (cf. (56)) inf E ( y) = h = s This inequality together with (56), (57) reveals that ; exp(;c ) s ; inf sup ( y) = s ; exp(;c ) kyk The proof of the theorem follows now from Proposition 4. 6
18 References [] Hadamard, J. (93). e probleme de Cauchy et les equations aux derivees partielles lineares hyperboliques. Paris. Herman et Cie. [] Chow, P.-. Khasminskii, R. (997). Statistical approach to dynamical inverse problems. Comm. Math. Phys. To appear. [3] Chow, P.-. Ibragimov, I. Khasminskii, R. (997). Statistical approach to some inverse problems for linear PDE. Submitted to PTRF [4] Ivanov, B. K. Vasin, B. V. Tanana, V. P. (978). Theory of linear ill-posed problems and its applications. Moscow, Nauka. [5] avrentiev, M. M. Romanov, V. G. and Shishtatskii, S. P. (98). Ill-posed problems of mathematical physics and analysis. Moscow, Nauka. [6] Petrovskii, N. G. (95). ectures on partial dierential equations. Moscow, Gostekhizdat. [7] Khaln,. A. (978). The statistical approach to the solution of ill-posed problems of geophysics. Studies in statistical estimation theory. Zap. Nauchn. Sem. eningrad. Otd. Inst. Steklov. (OMI) 79, 67{8. [8] attes, R. and ions, J.. (967). Methode de quasi-reversibilite et applications. Travoux et Recherche Mathematiques. 5. Dunod, Paris. [9] Pinsker, M. S. (98). Optimal ltering of square integrable signals in Gaussian white noise. Problems of Information Transmission, {33. [] Sullivan, F. O'. (996). A statistical perspective on ill-posed inverse problems. Statist. Science. 5{57. [] Sudakov, V. N. Khaln,. A. (964). Statistical approach to ill-posed problems in mathematical physics. DAN USSR {96. [] Tihonov, A. N. Arsenin, V. Ja. (986). Methods of solutions of ill-posed problems. Moscow, Nauka. [3] alsh, J. B. (984). An introduction to stochastic PDE. ect. Notes in Math. 8, Springer-Varlag, 65{439. 7
Mathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector
On Minimax Filtering over Ellipsoids Eduard N. Belitser and Boris Y. Levit Mathematical Institute, University of Utrecht Budapestlaan 6, 3584 CD Utrecht, The Netherlands The problem of estimating the mean
More informationThe Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov
More informationMarkov Sonin Gaussian rule for singular functions
Journal of Computational and Applied Mathematics 169 (2004) 197 212 www.elsevier.com/locate/cam Markov Sonin Gaussian rule for singular functions G.Mastroianni, D.Occorsio Dipartimento di Matematica, Universita
More information460 HOLGER DETTE AND WILLIAM J STUDDEN order to examine how a given design behaves in the model g` with respect to the D-optimality criterion one uses
Statistica Sinica 5(1995), 459-473 OPTIMAL DESIGNS FOR POLYNOMIAL REGRESSION WHEN THE DEGREE IS NOT KNOWN Holger Dette and William J Studden Technische Universitat Dresden and Purdue University Abstract:
More informationYadernaya Fizika, 56 (1993) Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993)
Start of body part published in Russian in: Yadernaya Fizika, 56 (1993) 45-5 translated in English in: Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993) A METHOD FOR CALCULATING THE HEAT KERNEL
More informationSPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in
SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING A. DELSHAMS, V. GELFREICH, A. JORBA AND T.M. SEARA At the end of the last century, H. Poincare [7] discovered the phenomenon of separatrices splitting,
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationThe Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-19 Wien, Austria The Negative Discrete Spectrum of a Class of Two{Dimentional Schrodinger Operators with Magnetic
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationThe following can also be obtained from this WWW address: the papers [8, 9], more examples, comments on the implementation and a short description of
An algorithm for computing the Weierstrass normal form Mark van Hoeij Department of mathematics University of Nijmegen 6525 ED Nijmegen The Netherlands e-mail: hoeij@sci.kun.nl April 9, 1995 Abstract This
More informationBilinear Stochastic Elliptic Equations
Bilinear Stochastic Elliptic Equations S. V. Lototsky and B. L. Rozovskii Contents 1. Introduction (209. 2. Weighted Chaos Spaces (210. 3. Abstract Elliptic Equations (212. 4. Elliptic SPDEs of the Full
More informationGaussian distributions and processes have long been accepted as useful tools for stochastic
Chapter 3 Alpha-Stable Random Variables and Processes Gaussian distributions and processes have long been accepted as useful tools for stochastic modeling. In this section, we introduce a statistical model
More informationA Stable Finite Dierence Ansatz for Higher Order Dierentiation of Non-Exact. Data. Bob Anderssen and Frank de Hoog,
A Stable Finite Dierence Ansatz for Higher Order Dierentiation of Non-Exact Data Bob Anderssen and Frank de Hoog, CSIRO Division of Mathematics and Statistics, GPO Box 1965, Canberra, ACT 2601, Australia
More informationIntroduction Wavelet shrinage methods have been very successful in nonparametric regression. But so far most of the wavelet regression methods have be
Wavelet Estimation For Samples With Random Uniform Design T. Tony Cai Department of Statistics, Purdue University Lawrence D. Brown Department of Statistics, University of Pennsylvania Abstract We show
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationREPORTRAPPORT. Centrum voor Wiskunde en Informatica. Asymptotics of zeros of incomplete gamma functions. N.M. Temme
Centrum voor Wiskunde en Informatica REPORTRAPPORT Asymptotics of zeros of incomplete gamma functions N.M. Temme Department of Analysis, Algebra and Geometry AM-R9402 994 Asymptotics of Zeros of Incomplete
More informationGarrett: `Bernstein's analytic continuation of complex powers' 2 Let f be a polynomial in x 1 ; : : : ; x n with real coecients. For complex s, let f
1 Bernstein's analytic continuation of complex powers c1995, Paul Garrett, garrettmath.umn.edu version January 27, 1998 Analytic continuation of distributions Statement of the theorems on analytic continuation
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationAN INEQUALITY FOR THE NORM OF A POLYNOMIAL FACTOR IGOR E. PRITSKER. (Communicated by Albert Baernstein II)
PROCDINGS OF TH AMRICAN MATHMATICAL SOCITY Volume 9, Number 8, Pages 83{9 S -9939()588-4 Article electronically published on November 3, AN INQUALITY FOR TH NORM OF A POLYNOMIAL FACTOR IGOR. PRITSKR (Communicated
More informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) YES, Eurandom, 10 October 2011 p. 1/32 Part II 2) Adaptation and oracle inequalities YES, Eurandom, 10 October 2011
More informationComputation Of Asymptotic Distribution. For Semiparametric GMM Estimators. Hidehiko Ichimura. Graduate School of Public Policy
Computation Of Asymptotic Distribution For Semiparametric GMM Estimators Hidehiko Ichimura Graduate School of Public Policy and Graduate School of Economics University of Tokyo A Conference in honor of
More informationThe best expert versus the smartest algorithm
Theoretical Computer Science 34 004 361 380 www.elsevier.com/locate/tcs The best expert versus the smartest algorithm Peter Chen a, Guoli Ding b; a Department of Computer Science, Louisiana State University,
More information106 E.G.Kirjackij i.e. in the circle jzj < 1; and let functions f(z) be normalized by conditions f(0) = 0; f 0 (0) = 1; and f 0 (z) 6= 0 in E: Let als
MMCA{97 Conference, pp. 105{116 R. Ciegis (Ed) c 1997 Vilniaus Gedimino technikos universitetas USING VARIATIONAL FORMULAE FOR SOLVING EXTREMAL PROBLEMS IN LINEARLY-INVARIANT CLASSES E.G.KIRJACKIJ Vilnius
More informationA converse Gaussian Poincare-type inequality for convex functions
Statistics & Proaility Letters 44 999 28 290 www.elsevier.nl/locate/stapro A converse Gaussian Poincare-type inequality for convex functions S.G. Bokov a;, C. Houdre ; ;2 a Department of Mathematics, Syktyvkar
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig Bloch electron in a magnetic eld and the Ising model by Igor V. Krasovsky Preprint no.: 20 2000 Bloch electron in a magnetic eld and
More informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) Yale, May 2 2011 p. 1/35 Introduction There exist many fields where inverse problems appear Astronomy (Hubble satellite).
More informationSources through Sensor. Observations. A. V. Kryazhimskii, V. I. Maksimov, and. Yu. S. Osipov WP August 1996
Working Paper Reconstruction of Boundary Sources through Sensor Observations A. V. Kryazhimskii, V. I. Maksimov, and Yu. S. Osipov WP-96-97 August 1996 IIASA International Institute for Applied Systems
More informationIterative procedure for multidimesional Euler equations Abstracts A numerical iterative scheme is suggested to solve the Euler equations in two and th
Iterative procedure for multidimensional Euler equations W. Dreyer, M. Kunik, K. Sabelfeld, N. Simonov, and K. Wilmanski Weierstra Institute for Applied Analysis and Stochastics Mohrenstra e 39, 07 Berlin,
More informationOn It^o's formula for multidimensional Brownian motion
On It^o's formula for multidimensional Brownian motion by Hans Follmer and Philip Protter Institut fur Mathemati Humboldt-Universitat D-99 Berlin Mathematics and Statistics Department Purdue University
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationCounterexamples to witness conjectures. Notations Let E be the set of admissible constant expressions built up from Z; + ;?;;/; exp and log. Here a co
Counterexamples to witness conjectures Joris van der Hoeven D pt. de Math matiques (b t. 45) Universit Paris-Sud 91405 Orsay CEDEX France August 15, 003 Consider the class of exp-log constants, which is
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationApplications and fundamental results on random Vandermon
Applications and fundamental results on random Vandermonde matrices May 2008 Some important concepts from classical probability Random variables are functions (i.e. they commute w.r.t. multiplication)
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationSupport Vector Regression with Automatic Accuracy Control B. Scholkopf y, P. Bartlett, A. Smola y,r.williamson FEIT/RSISE, Australian National University, Canberra, Australia y GMD FIRST, Rudower Chaussee
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationSupport Detection in Super-resolution
Support Detection in Super-resolution Carlos Fernandez-Granda (Stanford University) 7/2/2013 SampTA 2013 Support Detection in Super-resolution C. Fernandez-Granda 1 / 25 Acknowledgements This work was
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationUpper Bounds for Partitions into k-th Powers Elementary Methods
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 9, 433-438 Upper Bounds for Partitions into -th Powers Elementary Methods Rafael Jaimczu División Matemática, Universidad Nacional de Luján Buenos Aires,
More informationMinimax Risk: Pinsker Bound
Minimax Risk: Pinsker Bound Michael Nussbaum Cornell University From: Encyclopedia of Statistical Sciences, Update Volume (S. Kotz, Ed.), 1999. Wiley, New York. Abstract We give an account of the Pinsker
More informationwith Perfect Discrete Time Observations Campus de Beaulieu for all z 2 R d M z = h?1 (z) = fx 2 R m : h(x) = zg :
Nonlinear Filtering with Perfect Discrete Time Observations Marc Joannides Department of Mathematics Imperial College London SW7 2B, United Kingdom mjo@ma.ic.ac.uk Francois LeGland IRISA / INRIA Campus
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationCompanding of Memoryless Sources. Peter W. Moo and David L. Neuho. Department of Electrical Engineering and Computer Science
Optimal Compressor Functions for Multidimensional Companding of Memoryless Sources Peter W. Moo and David L. Neuho Department of Electrical Engineering and Computer Science University of Michigan, Ann
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More information2 Section 2 However, in order to apply the above idea, we will need to allow non standard intervals ('; ) in the proof. More precisely, ' and may gene
Introduction 1 A dierential intermediate value theorem by Joris van der Hoeven D pt. de Math matiques (B t. 425) Universit Paris-Sud 91405 Orsay Cedex France June 2000 Abstract Let T be the eld of grid-based
More informationOPTIMAL POINTWISE ADAPTIVE METHODS IN NONPARAMETRIC ESTIMATION 1
The Annals of Statistics 1997, Vol. 25, No. 6, 2512 2546 OPTIMAL POINTWISE ADAPTIVE METHODS IN NONPARAMETRIC ESTIMATION 1 By O. V. Lepski and V. G. Spokoiny Humboldt University and Weierstrass Institute
More informationA slow transient diusion in a drifted stable potential
A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationROYAL INSTITUTE OF TECHNOLOGY KUNGL TEKNISKA HÖGSKOLAN. Department of Signals, Sensors & Systems
The Evil of Supereciency P. Stoica B. Ottersten To appear as a Fast Communication in Signal Processing IR-S3-SB-9633 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems Signal Processing
More informationing a minimal set of layers, each describing a coherent motion from the scene, which yield a complete description of motion when superposed. Jepson an
The Structure of Occlusion in Fourier Space S. S. Beauchemin Dept. of Computer Science The University of Western Ontario London, Canada N6A 5B7 e-mail: beau@csd.uwo.ca J. L. Barron Dept. of Computer Science
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationD I S C U S S I O N P A P E R
I N S T I T U T D E S T A T I S T I Q U E B I O S T A T I S T I Q U E E T S C I E N C E S A C T U A R I E L L E S ( I S B A ) UNIVERSITÉ CATHOLIQUE DE LOUVAIN D I S C U S S I O N P A P E R 2014/06 Adaptive
More informationON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay
ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov
More informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More informationMultidimensional Chebyshev Systems (Haar systems) - just a definition
arxiv:0808.2213v2 [math.fa] 17 Apr 2011 Multidimensional Chebyshev Systems (Haar systems) - just a definition Ognyan Kounchev Institute of Mathematics and Informatics Bulgarian Academy of Sciences and
More informationPavel Dvurechensky Alexander Gasnikov Alexander Tiurin. July 26, 2017
Randomized Similar Triangles Method: A Unifying Framework for Accelerated Randomized Optimization Methods Coordinate Descent, Directional Search, Derivative-Free Method) Pavel Dvurechensky Alexander Gasnikov
More informationNumerical Schemes from the Perspective of Consensus
Numerical Schemes from the Perspective of Consensus Exploring Connections between Agreement Problems and PDEs Department of Electrical Engineering and Computer Sciences University of California, Berkeley
More informationQUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS * WEI-CHANG SHANN AND JANN-CHANG YAN Abstract. Scaling equations are used to derive formulae of quadratures involving polynomials and scaling/wavelet
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More information114 EUROPHYSICS LETTERS i) We put the question of the expansion over the set in connection with the Schrodinger operator under consideration (an accur
EUROPHYSICS LETTERS 15 April 1998 Europhys. Lett., 42 (2), pp. 113-117 (1998) Schrodinger operator in an overfull set A. N. Gorban and I. V. Karlin( ) Computing Center RAS - Krasnoyars 660036, Russia (received
More informationSECTION FOR DIGITAL SIGNAL PROCESSING DEPARTMENT OF MATHEMATICAL MODELLING TECHNICAL UNIVERSITY OF DENMARK Course 04362 Digital Signal Processing: Solutions to Problems in Proakis and Manolakis, Digital
More informationConcentration of Measure: Logarithmic Sobolev Inequalities
Concentration o Measure: Logarithmic Sobolev Inequalities Aukosh Jagannath September 15, 13 Abstract In this note we'll describe the basic tools to understand Log Sobolev Inequalities and their relation
More informationSome nonlinear elliptic equations in R N
Nonlinear Analysis 39 000) 837 860 www.elsevier.nl/locate/na Some nonlinear elliptic equations in Monica Musso, Donato Passaseo Dipartimento di Matematica, Universita di Pisa, Via Buonarroti,, 5617 Pisa,
More informationCan we do statistical inference in a non-asymptotic way? 1
Can we do statistical inference in a non-asymptotic way? 1 Guang Cheng 2 Statistics@Purdue www.science.purdue.edu/bigdata/ ONR Review Meeting@Duke Oct 11, 2017 1 Acknowledge NSF, ONR and Simons Foundation.
More informationin the company. Hence, we need to collect a sample which is representative of the entire population. In order for the sample to faithfully represent t
10.001: Data Visualization and Elementary Statistical Analysis R. Sureshkumar January 15, 1997 Statistics deals with the collection and the analysis of data in the presence of variability. Variability
More informationDomain and range of exponential and logarithmic function *
OpenStax-CNX module: m15461 1 Domain and range of exponential and logarithmic function * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License
More informationOn a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI
On a WKB-theoretic approach to Ablowitz-Segur's connection problem for the second Painleve equation Yoshitsugu TAKEI Research Institute for Mathematical Sciences Kyoto University Kyoto, 606-8502, JAPAN
More informationREMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS. Leonid Friedlander
REMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS Leonid Friedlander Abstract. I present a counter-example to the conjecture that the first eigenvalue of the clamped buckling problem
More informationPade approximants and noise: rational functions
Journal of Computational and Applied Mathematics 105 (1999) 285 297 Pade approximants and noise: rational functions Jacek Gilewicz a; a; b;1, Maciej Pindor a Centre de Physique Theorique, Unite Propre
More informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationBayesian Nonparametric Point Estimation Under a Conjugate Prior
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 5-15-2002 Bayesian Nonparametric Point Estimation Under a Conjugate Prior Xuefeng Li University of Pennsylvania Linda
More informationWerner Romisch. Humboldt University Berlin. Abstract. Perturbations of convex chance constrained stochastic programs are considered the underlying
Stability of solutions to chance constrained stochastic programs Rene Henrion Weierstrass Institute for Applied Analysis and Stochastics D-7 Berlin, Germany and Werner Romisch Humboldt University Berlin
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Numerical Algorithms for Inverse and Ill-Posed Problems - A.М. Denisov
NUMERICAL ALGORITHMS FOR INVERSE AND ILL-POSED PROBLEMS Faculty of Computational Mathematics and Cybernetics, Moscow State University, Moscow, Russia Keywords: Direct problem, Inverse problem, Well-posed
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig Uniformly distributed measures in Euclidean spaces by Bernd Kirchheim and David Preiss Preprint-Nr.: 37 1998 Uniformly Distributed
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More informationPart IV Stochastic Image Analysis 1 Contents IV Stochastic Image Analysis 1 7 Introduction to Stochastic Processes 4 7.1 Probability................................... 5 7.1.1 Random events and subjective
More informationIntroduction and mathematical preliminaries
Chapter Introduction and mathematical preliminaries Contents. Motivation..................................2 Finite-digit arithmetic.......................... 2.3 Errors in numerical calculations.....................
More information1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary
Prexes and the Entropy Rate for Long-Range Sources Ioannis Kontoyiannis Information Systems Laboratory, Electrical Engineering, Stanford University. Yurii M. Suhov Statistical Laboratory, Pure Math. &
More informationInverse Gravimetry Problem
Inverse Gravimetry Problem Victor Isakov September 21, 2010 Department of M athematics and Statistics W ichita State U niversity W ichita, KS 67260 0033, U.S.A. e mail : victor.isakov@wichita.edu 1 Formulation.
More informationfunction independent dependent domain range graph of the function The Vertical Line Test
Functions A quantity y is a function of another quantity x if there is some rule (an algebraic equation, a graph, a table, or as an English description) by which a unique value is assigned to y by a corresponding
More informationMinimal Surfaces. December 13, Alex Verzea. MATH 580: Partial Dierential Equations 1. Professor: Gantumur Tsogtgerel
Minimal Surfaces December 13, 2012 Alex Verzea 260324472 MATH 580: Partial Dierential Equations 1 Professor: Gantumur Tsogtgerel 1 Intuitively, a Minimal Surface is a surface that has minimal area, locally.
More informationWavelets and regularization of the Cauchy problem for the Laplace equation
J. Math. Anal. Appl. 338 008440 1447 www.elsevier.com/locate/jmaa Wavelets and regularization of the Cauchy problem for the Laplace equation Chun-Yu Qiu, Chu-Li Fu School of Mathematics and Statistics,
More informationComputation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling
Journal of Computational and Applied Mathematics 144 (2002) 359 373 wwwelseviercom/locate/cam Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling
More informationConcentration inequalities and the entropy method
Concentration inequalities and the entropy method Gábor Lugosi ICREA and Pompeu Fabra University Barcelona what is concentration? We are interested in bounding random fluctuations of functions of many
More informationCreating materials with a desired refraction coefficient: numerical experiments
Creating materials with a desired refraction coefficient: numerical experiments Sapto W. Indratno and Alexander G. Ramm Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More informationFrom Fractional Brownian Motion to Multifractional Brownian Motion
From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December
More informationError Estimates for Trigonometric Interpolation. of Periodic Functions in Lip 1. Knut Petras
Error Estimates for Trigonometric Interpolation of Periodic Functions in Lip Knut Petras Abstract. The Peano kernel method is used in order to calculate the best possible constants c, independent of M,
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More informationNonparametric volatility estimation in scalar diusions: Optimality across observation frequencies.
Nonparametric volatility estimation in scalar diusions: Optimality across observation frequencies. akub Chorowski Institute of Mathematics Humboldt-Universität zu Berlin chorowsj@math.hu-berlin.de Abstract
More information1 Introduction Travelling-wave solutions of parabolic equations on the real line arise in a variety of applications. An important issue is their stabi
Essential instability of pulses, and bifurcations to modulated travelling waves Bjorn Sandstede Department of Mathematics The Ohio State University 231 West 18th Avenue Columbus, OH 4321, USA Arnd Scheel
More informationU U B U P x
Smooth Quantum Hydrodynamic Model Simulation of the Resonant Tunneling Diode Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-184 Abstract Smooth
More informationThe Pacic Institute for the Mathematical Sciences http://www.pims.math.ca pims@pims.math.ca Surprise Maximization D. Borwein Department of Mathematics University of Western Ontario London, Ontario, Canada
More informationConvergence of The Multigrid Method With A Wavelet. Abstract. This new coarse grid operator is constructed using the wavelet
Convergence of The Multigrid Method With A Wavelet Coarse Grid Operator Bjorn Engquist Erding Luo y Abstract The convergence of the two-level multigrid method with a new coarse grid operator is studied.
More informationA theorem on summable families in normed groups. Dedicated to the professors of mathematics. L. Berg, W. Engel, G. Pazderski, and H.- W. Stolle.
Rostock. Math. Kolloq. 49, 51{56 (1995) Subject Classication (AMS) 46B15, 54A20, 54E15 Harry Poppe A theorem on summable families in normed groups Dedicated to the professors of mathematics L. Berg, W.
More information