arxiv:math/ v3 [math.pr] 26 Sep 2005

Similar documents
Quadratic harnesses, q-commutations, and orthogonal martingale polynomials. Quadratic harnesses. Previous research. Notation and Motivation

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

X n = c n + c n,k Y k, (4.2)

Stochastic process for macro

Continued fractions for complex numbers and values of binary quadratic forms

A matrix over a field F is a rectangular array of elements from F. The symbol

4 Sums of Independent Random Variables

An almost sure invariance principle for additive functionals of Markov chains

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION

ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS

1 Last time: least-squares problems

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan

Gaussian Processes. 1. Basic Notions

Linearization coefficients for orthogonal polynomials. Michael Anshelevich

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

MATH 205C: STATIONARY PHASE LEMMA

arxiv: v2 [math.pr] 4 Sep 2017

ELEMENTARY LINEAR ALGEBRA

P (A G) dp G P (A G)

Boolean Inner-Product Spaces and Boolean Matrices

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Lecture 12: Detailed balance and Eigenfunction methods

LARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*

A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra

Proofs for Large Sample Properties of Generalized Method of Moments Estimators

Notes on Quadratic Extension Fields

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.

arxiv: v4 [math.pr] 18 May 2017

Eötvös Loránd University, Budapest. 13 May 2005

1 Linear Difference Equations

Lecture 7. Econ August 18

2. Matrix Algebra and Random Vectors

. Find E(V ) and var(v ).

SOMEWHAT STOCHASTIC MATRICES

Supermodular ordering of Poisson arrays

Lecture Notes 1: Vector spaces

I forgot to mention last time: in the Ito formula for two standard processes, putting

Kernel families of probability measures. Saskatoon, October 21, 2011

Generalized Gaussian Bridges of Prediction-Invertible Processes

CAAM 454/554: Stationary Iterative Methods

Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments

Notes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009

Other properties of M M 1

A Concise Course on Stochastic Partial Differential Equations

ELEMENTARY LINEAR ALGEBRA

A CLT FOR MULTI-DIMENSIONAL MARTINGALE DIFFERENCES IN A LEXICOGRAPHIC ORDER GUY COHEN. Dedicated to the memory of Mikhail Gordin

For a stochastic process {Y t : t = 0, ±1, ±2, ±3, }, the mean function is defined by (2.2.1) ± 2..., γ t,

RATE-OPTIMAL GRAPHON ESTIMATION. By Chao Gao, Yu Lu and Harrison H. Zhou Yale University

MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5

Stein s Method and Characteristic Functions

A PRIMER ON SESQUILINEAR FORMS

The properties of L p -GMM estimators

Harmonic Analysis. 1. Hermite Polynomials in Dimension One. Recall that if L 2 ([0 2π] ), then we can write as

Kaczmarz algorithm in Hilbert space

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Stochastic Realization of Binary Exchangeable Processes

On Least Squares Linear Regression Without Second Moment

We simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =

SAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

1. Introduction to commutative rings and fields

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

UC Berkeley Department of Electrical Engineering and Computer Sciences. EECS 126: Probability and Random Processes

ON KRONECKER PRODUCTS OF CHARACTERS OF THE SYMMETRIC GROUPS WITH FEW COMPONENTS

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =

7 Curvature of a connection

IEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion

Large Sample Properties of Estimators in the Classical Linear Regression Model

Notes on Linear Algebra and Matrix Theory

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0}

Chebyshev coordinates and Salem numbers

Foundations of Matrix Analysis

Detailed Proof of The PerronFrobenius Theorem

Independence of some multiple Poisson stochastic integrals with variable-sign kernels

A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University

ELEMENTARY LINEAR ALGEBRA

CHUN-HUA GUO. Key words. matrix equations, minimal nonnegative solution, Markov chains, cyclic reduction, iterative methods, convergence rate

This operation is - associative A + (B + C) = (A + B) + C; - commutative A + B = B + A; - has a neutral element O + A = A, here O is the null matrix

SMSTC (2007/08) Probability.

Elementary linear algebra

Analysis-3 lecture schemes

Mathematical Methods wk 2: Linear Operators

A NOTE ON FREE QUANTUM GROUPS

Characterizations of free Meixner distributions

Some functional (Hölderian) limit theorems and their applications (II)

STAT 7032 Probability Spring Wlodek Bryc

arxiv: v1 [math.pr] 1 Jan 2016

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models

ACO Comprehensive Exam October 14 and 15, 2013

Lecture 13: Simple Linear Regression in Matrix Format. 1 Expectations and Variances with Vectors and Matrices

SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA

ELEMENTARY LINEAR ALGEBRA

Meixner matrix ensembles

Spectral inequalities and equalities involving products of matrices

STAT 100C: Linear models

7 Poisson random measures

Unit Roots in White Noise?!

Transcription:

ariv:math/0507303v3 [math.pr] 6 Sep 005 q WIENER AND RELATED PROCESSES. A BRYC PROCESSES CONTINUOUS TIME GENERALIZATION PAWEL J. SZABLOWSKI Abstract. We define two Markov processes. The finite dimensional distributions of the first one say = t t 0 depend on one parameter q 1, 1 > and of the second one say Y =Y t t R on two parameters q, α 1, 1 > 0,. The first one resembles Wiener process in the sense that for q = 1 it is Wiener process but also that for q < 1 and n 1 t n/ H n t/ t q, where H n n 0 are so called q Hermite polynomials, are martingales. It does not have however independent increments. The second one resemble Orstein- Ulehnbeck processes. For q = 1 it is a classical OU process. For q < 1 it is stationary with correlation function equal to exp α t s. On the way we deny Wesolowski s martingale characterization of Wiener process 1. Introduction As announced in abstract, we are going to define two continuous time families of Markov processes. One of them will resemble Wiener process and the other Orstein- Uhlenbeck OU process. They will be indexed apart from time parameter by additional parameter q 1, 1 > and in case q = 1 those processes will be classical Wiener and OU process. Of course for OU process there will additional parameter α > 0 responsible for covariance function of the process. For q 1, 1 both processes will assume values in a compact space < / 1 q, / 1 q >. One dimesional probabilities and transitional probabilities of these will be given explicitly. Moreover two families of polynomials orthogonal with respect to these measures will also be presented. Some properties of conditional expectations given the past of these processes including martingale properties will also be exposed. Thus quite detailed knowledge concerning these processes will be presented. The problem is that these processes will be introduced via discrete time one dimensional random fields Bryc processes that have been studied in detail recently see e.g. [1], [], [4], [5], [3]. That is why the first section will be dedicated to definition of Bryc processes and recollection of their basic properties. The interesting thing about Bryc and related processes is that their one dimensional densities belong to so called family of q Gaussian distributions introduced Bożejko et. al. in [10] using purely non-commutative probability arguments including special interpretation of parameter q. Hence these distributions and consequently Bryc processes can be viewed as the mysterious connection between commutative and non-commutative probability theories. Since 1997 there appeared Date: July, 005. 000 Mathematics Subject Classification. Primary 60J5, 60G44; Secondary 05A30, 60K40. Key words and phrases. q Wiener process, q Orstein-Uhlenbeck process, polynomial martingales property, q-gaussian distributions, quadratic harnesses. 1

PAWEL J. SZABLOWSKI couple of papers mostly in journals dedicated to Physics on the properties of q Gaussian distributions. See e.g. [7], [14], [1]. There is also a different path of research followed by W lodek Bryc and Jacek Weso lowski see e.g. [6], [13]. Their starting point is process with continuous time that satisfies several exactly 5 conditions on covariance function and on first and second conditional moments. Under these assumptions they prove that these processes are Markov and several properties of families of polynomials that orthogonalize transitional and one dimensional probabilities. They give several examples illustrating developed theory. One of the processes considered by them is so called q Brownian process. As far as one dimensional probabilities and transitional probabilities are concerned this process is identical with introduced in this paper q Wiener process. The point is that we derive this process under less assumptions on the first and second conditional moments we need only. The construction is also different. Our starting point is discrete time Bryc process. q Wiener process is obtained as continuous time transformation of continuous time generalization of the discrete time process. Besides we list many properties of q Wiener process that justify its name are sort of q analogies of well known martingale properties of Wiener process. Some of these properties can be derived from the definition of q Brownian process presented in [6]. They are not however stated explicitly. Recently W lodek Bryc and Jacek Weso lowski together with Wojtek Matysiak in [8] follow another more general path and considered so called continuous time quadratic harnesses. Vaguely speaking quadratic harnesses Z =Z t t>0 are such L processes with zero mean and covariance function EZ s Z t = min s, t that have respectively linear and quadratic first and second conditional moments of Z t with respect to field Z :, s > < u, for s < t < u. In their paper they show that each quadratic harness is characterized by 5 parameters apart from time parameter. They show also that under mild technical assumptions for each harness one can define a system of orthogonal polynomials {P n x, t} n 1 such that n : P n Z t, t, Z : t t>0 is a martingale. One of the processes that are defined in this paper, namely q Wiener process, although emerging from different context, turns out to be quadratic harnesses with 4 out of 5 parameters assuming value equal to 0, the family of orthogonal polynomials of its one dimensional distribution satisfies 3 term recursive fomulae given in paper of Bryc Matysiak and Weso lowski. They provide also an example of quadratic harness that is not an independent increment process. So far all the concrete examples given by these there authors concerned independent increment processes. Our approach is totally commutative, classical as the one presented in recent papers of Bryc and Weso lowski. q by no means is an operator. It is number parameter. Yet we are touching q series and special functions. The paper is organized as follows. As we mentioned in the first section we recall definition and summarize basic properties of Bryc processes. The next section is still dedicated mainly to certain properties of discrete time Bryc processes. In the third section first introduce continuous version of Bryc processes then we prove their existence, later parts of this section are devoted to presentation of the two processes that were mentioned earlier. We show their connection with emerging quadratic harness theory. We show here that Weso lowski s martingale characterization contained in [9] of Wiener process is not true. The last section contains proofs of the results from previous sections.

GENERALIZATIONS OF BRYC S PROCESSES 3 Although it is still not known if the Bryc process exist with parameter q > 1, we show in this paper that q Wiener and q Ortein-Uhlenbeck processes do not exist for q > 1.. Bryc processes By Bryc processes we mean square integrable random field = { n } n Z indexed by the integers, with non-degenerate covariance matrices and constant first two moments, that satisfy the following two sets of conditions :.1 a, b R; n Z : E n F n = a n 1 + n+1 + b, a.s. and A, B, C R; n Z : E n F n. =.3 A n 1 + n+1 + Bn 1 n+1 + D n 1 + n+1 + C, a.s., where F m := k : k m. Non-singularity of covariance matrices implies that all random variables n are non-degenerate and there is no loss of generality in assuming that E k = 0 and Ek = 1 for all k Z, which implies b = 0. It has been shown in [4] that.1 implies L -stationarity of. Since the case ρ := corr 0, 1 = 0 contains sequences of independent random variables which satisfy.1 and. but can have arbitrary distributions, we shall exclude it from the considerations. Observe that non-singularity of the covariance matrices implies ρ < 1. By Theorem 3.1 from [1] see also Theorems 1 and in [4], corr 0, k = ρ k and one-sided regressions are linear.4 E k..., 1, 0 = ρ k 0. Let us define algebras F m := k : k m, F m := k : k m, F m := k : k m, F m, j := k : k m k j. It turns out that parameters a, ρ, A, B, C are related to one another in the following way:.5 a = ρ 1 + ρ, C = 1 A ρ B, and.6 A In [5] it was shown that ρ + 1ρ + B = 1. Theorem 1. There do not exist standardized random fields = k k Z with nonsingular covariance matrices, satisfying.1 and. and ρ 0, unless D = 0 and either B = 0 or.6 holds with B B 1 B B 3, where { ρ 1 B 1 = ρ 1 + ρ { ρ 1 + ρ B 3 = }, B = 0, ρ 1 + ρ 1 ρ 4 1 + ρ /m ρ /m ρ 4 : m N Furthermore, if D = B = 0 or D = 0 and.6 holds with B B then the relevant fields exist and have the following distributions: ], }.

4 PAWEL J. SZABLOWSKI 1 If B = 0 and A = 1/ then k = RY k, R and Y k are independent, R 0, ER = 1 and all random variables Y k are symmetric with values ±1; If B = 0 and A 1/ then all random variables k are symmetric with values ±1; 3 If.6 holds then for all B 0, ρ / 1 + ρ all one-dimensional distributions of are equal to the uniquely determined symmetric and absolutely continuous distribution depending on B supported on a finite interval; 4 If.6 holds and B = ρ / 1 + ρ then k N 0, 1 for all k Z. Theorem. If.6 holds with D = 0 and B B 1 then there do not exist the relevant random fields satisfying.1-. with ρ 0 and such that sequence k k is uniformly integrable and Pr k 1 = 0 < 1. There also i.e. in [5] the following problem has been put up. Problem 1. Do there exist random fields satisfying.1-. with ρ 0, D = 0, and such that.6 holds with B B 3? Theorem 3 below, strengthens Theorem.3 from [] by stating that is always a Markov chain and relaxing assumptions for the uniquely determined case. Theorem 3. If is a standardized random field with non-singular covariance matrices, satisfying.1 and.1-., ρ 0, and such that D = B = 0 or D = 0 and.6 holds with B B, then is a stationary Markov chain. has uniquely determined distributions unless B = D = 0 and A = 1/. In [1] it has been shown regular B-processes are stationary random chains and that one can redefine parameters by introducing new parameter q < 1, 1 > and define parameters A and B as functions of q and ρ. It turns out that q = 1 is equivalent to the B = 0, case q 1, 1 > is equivalent to B B, the case B B 1 leads to undefined q and finally case B B 3 leads to q { 1/ρ n : n N }. It is known that the case q < 1 implies to nonprobabilistic solutions see e.g.[1]. Following.[1] we can express parameters A, B, C using parameters ρ and q : ρ 1 qρ.7a A = ρ + 11 qρ 4,.7b B = ρ 1 ρ 1 + q ρ + 11 qρ 4. Since we have.5 we can also express C as a function of q and ρ : 1 ρ C = 1 qρ 4. Conversely parameter q can be expressed by parameter B by the relationship: ρ 4 + B ρ + ρ 1 1.8 q =. 1 + ρ B 4 ρ + 1 ρ 1 We can rephrase and summarize the results of [1] and [4] in the following way. Each Bryc process is characterized by two parameters ρ and q. For q outside the

GENERALIZATIONS OF BRYC S PROCESSES 5 set < 1, 1 > { 1/ρ /n : n N } Bryc processes do not exist. For 0 < ρ < 1 and q < 1, 1 > Bryc processes exist and all their finite dimensional distribution are uniquely determined and known. For 0 < ρ < 1 and q { 1/ρ /n : n N } it is not known if Bryc processes exist. Bryc sequences with 0 < ρ < 1 and q 1, 1 > we will denote regular B- processes. To complete recollection of basic properties of regular B-processes let us introduce so called q Hermite polynomials {H n x q} n 1 defined by the following recurrence:.9a.9b H 1 x q = 0, H 0 x q = 1, n 0 : xh n x q = H n+1 x q + [n] q H n 1 x q, where [0] q = 0, [1] q = 1, n [n] q = 1 + q +... + q n 1. W. Bryc in [1] has shown that: n Z, k 1 : EH k n q F n 1 = ρ k H k n 1 q, a.s. and also that {H n x q} n 1 are orthogonal polynomials of distribution of n. The case q = 1 is equivalent to B = 0 and leads to distribution described above. We will concentrate thus on the case q 1, 1 >. It turns out then that this distribution has density for q 1, 1 and x 1 q, 1 q 1 q f H x q = π 4 1 qx 1 + q k 1 qx q k 1 q k+1, k=0 and.10 f H x 1 = 1 exp x, x R. π W. Bryc in [1] also found density of the conditional distribution n n 1 = y and later W. Bryc, W. Matysiak and P. J. Szab lowski in [3] found orthogonal polynomials of this conditional distribution. Namely it turns out that for q 1, 1 and x 1 q, 1 q we have 1 q f x y, q, ρ = π 4 1 qx.11 k=0 k=0 1 ρ q k 1 q k+1 1 + q k 1 qx q k 1 ρ q k 1 qρq k 1 + ρ q k xy + 1 qρ x + y q k, and that polynomials {p n x y, q, ρ} n 1 defined by.1a p 1 x y, q, ρ = 0, p 0 x y, q, ρ = 1,.1b n 0 : p n+1 x y, q = x ρyq n p n x y, q, ρ 1 ρ q n 1 [n] q p n 1 x y q, ρ, are orthogonal with respect to measure defined by the density.11 For q = 1 we have 1.13 f x y, 1, ρ = exp x ρy 1 ρ 1 ρ.

6 PAWEL J. SZABLOWSKI In this paper we are going to consider continuous time generalization of. Namely, more precisely, we are going to study the existence and basic properties of process Y = {Y t } t R, such that for every positive δ, random sequences are regular B-processes. δ n = Y nδ, n Z, 3. Auxiliary properties of Bryc processes The following two properties follow assumptions.1 and. and assumptions concerning parameters ρ and A, B, and C. Lemma 1. If is a Bryc process then: 3.1a 3.1b 3.1c E 4 n = + q, E n m = 1 + ρ n m 1 + q, E n n j n+k = ρ k+j E 4 n = ρ j+k + q. Proposition 1. If is a Bryc process then E n F n k, n+j is a linear function of n k and n+j. More precisely we have 3. E n F n k, n+j = ρk 1 ρ j 1 ρ j+k n k + ρj 1 ρ k 1 ρ n+j. j+k Proposition. If is Bryc process, then for n, j, k N, E n F n j, n+k is a linear function of n j, n+k and n j n+k. In particular 3.3 E n F n j, n+k = A 1 jk n j + A jk n+j + B jk n j n+j + C jk where: 3.4a 3.4b 3.4c 3.4d A 1 jk = ρj 1 ρ k 1 qρ k 1 qρ j+k 1 ρ j+k, A jk = ρk 1 ρ j 1 qρ j 1 qρ j+k 1 ρ j+k B jk = q + 1ρj+k 1 ρ j 1 ρ k 1 qρ j+k 1 ρ j+k, C jk = 1 ρ j 1 ρ k 1 qρ j+k. As a simple corollary of Proposition and its proof we have the following Lemma Lemma. Let { i } i Z Bryc process with parameters ρ and q. Let us define Z k = kj+m for fixed j N, m = 0,...,j 1 and k N. {Z k } k Z is also Bryc process with parameters ρ j and q. 4. Continuous time generalizations of regular B processes 4.1. Existence. From consideration { } of the previous section it follows that for each δ > 0 regular B-process n δ is L stationary with covariance function K n, m = ρ δ n m = E δ 0 δ n m n Z. But we have also for any integer k: ρ δ n = E δ 0 δ n = E δ/k 0 δ/k nk = ρ δ/k nk.

GENERALIZATIONS OF BRYC S PROCESSES 7 From this identity by standard methods we deduce that ρ δ = ρ δ, for certain positive parameter ρ < 1. Another words we deduce that if Y exists, then it is L stationary with covariance function K t, s = K s t = ρ s t. for some ρ 0, 1. Let us introduce new parameter We have then α = log 1 ρ > 0. K s, t = exp α s t. Existence of Y will be shown for two cases separately. Since for q = 1 we have normality of one dimensional and conditional distributions. Thus we deduce that process Y for q = 1 exists. It is in fact well known Orstein-Uhlenbeck process. Now consider fixed q 1, 1. Let S =< / 1 q, / 1 q >. First we will deduce the existence of the process Ỹ = Y t t Q. This follows from Kolmogorov s extension theorem. Since having natural ordering of Q we need only consistency of the family of finite dimensional distributions of Ỹ. This can be however easily shown by the following argument. Let us take finite set of numbers r 1 < r < < r n from Q. Let M denote the smallest common denominator of these numbers. Let us consider regular B process M with q and ˆρ = ρ 1/M or for parameter α : ˆα = α/m. Note that then n M = Y n/m for n N. Then joint distribution of Y r1,...,y rn is in fact a joint distribution of R M 1,...,R M n where numbers R 1,..., R n are defined by the relationships r i = R i /M. Since process M exists we have consistency since if { 1,..., k } {r 1,...,r n } then distribution of Y 1,..., Y k being equal to the distribution of T M 1,...,T M k with Tj defined by j = T j /M for j = 1,...,k is a projection of the distribution of Y r1,..., Y rn. Hence we deduce that process Ỹ with values in compact space S, exist. Now we use separation theorem and view Ỹ as separable modification of process Y itself. Hence Y exists. 4.. q OU processes. Process Y with parameters q, α will be called continuous time q, α OU-process. OU standing for Orstein-Uhlenbeck. Let us summarize what properties of Y can be deduced from the properties of discrete time regular B-processes presented in [1], [], [4], [5] and [3]. Theorem 4. Let Y be a continuous time q, α OU-process. Then: 4.1 4. 1 t R : Y t f H x q, s > t : Y s Y t = y f x y, q, e α s t, Y is a stationary Markov process with f x y, q, e α s t as its transfer function.

8 PAWEL J. SZABLOWSKI 3 Y is time symmetric. Moreover we have for any s R and δ, γ > 0 : EY s F s δ, s+γ = e αδ 1 e αγ Y 1 e αδ+γ s δ + e αγ 1 e αδ Y 1 e αδ+γ s+γ, and for s R and δ, γ > 0 : E Y s F s δ, s+γ = A1 δ, γy s δ + A δ, γy s+γ + B δ, γy s δ Y s+γ + C δ, γ, where: A 1 δ, γ = e αδ 1 e αγ 1 qe αγ 1 qe αδ+γ 1 e αδ+γ, A δ, γ = e αγ 1 e αδ 1 qe αδ 1 qe αδ+γ 1 e αδ+γ B δ, γ = q + 1e αδ+γ 1 e αδ 1 e αγ 1 qe αδ+γ 1 e αδ+γ, C δ, γ = 1 e αδ 1 e αγ 1 qe αδ+γ. 4 the following families of polynomials respectively given by.9 and.1 {H n x q} n 0 and { p x y, q, e α s t }, are orthogonal polynomials of n 0 4.1 and 4.. That is in particular 5 n 1, t R : EH n Y t q = 0 n 1, s > t : E p n Y s Y t, q, e α s t F t = 0 a.s., n 1, s > t : EH n Y s q = e nα s t H n Y t q a.s., 6 Consequently n 1 : H n Y t q t R is a stationary random process with covariance function and K n s, t = e nα s t nα S n ω α = A n ω + n α, as its spectral density, where A n = EH n Y t q. Remark 1. Thus in particular {H n Y t q} t R more and more resembles white noise the greater n we choose. Remark. Note that following Theorem 1 and considerations of the section 4.1, for the existence of q OU process one needed only the following two conditions apart from linear independence for any m Z, n N and d > 0 of random variables from the set { Y md, Y m+1d,...,y m+nd } : 1 4.3 n Z, d > 0 : E Y nd F n 1d, n+1d = e αd 1 + e αd Yn 1d + Y n+1d

GENERALIZATIONS OF BRYC S PROCESSES 9 4.4 n Z, d > 0 : E Y nd F n 1d, n+1d =  Yn 1d + Y n+1d + ˆBY n 1d Y n+1d +Ĉ, where  = A 1 d, d = A d, d, ˆB = B d, d, Ĉ = C d, d. That is we need only symmetric and discrete for all increments d > 0 versions of condition defining E Ys F s δ, s+γ. Remark 3. Note that continuous time q OU process does not exist for q > 1. It is so because for discrete time Bryc process with parameters q, ρ presuming such process existed - compare Theorem 1 and Problem 1 the following relationship between parameters q and ρ must be satisfied for some integer n : ρ q n = 1 or equivalently ratio log q log ρ must be equal to some integer. However if q OU process existed parameter ρ would depend on time parameter t in the following way ρ = exp αt for some fixed positive α and consequently log q αt would have to be integer for all real t which is impossible. 4.3. q Wiener process. Let Y be given q, α OU process. Let us define: 4.5 0 = 0; > 0 : = Y log /α. Process = 0 will be called q Wiener process. Let us also introduce the following filtration: F = t : t = Y t : t log /α = F log /α Following this definition and Theorem 4 of the previous section we have: Theorem 5. Let be a q Wiener process, then it has the following properties: 1, 0 : cov, = min,, > 0 : 1 f H x q 3 = y 1 x + y f y, q, 4 For all n 1 and 0 < we have E n/ H n, q F = n/ 4.6 H n q, a.s E n/ H n q = n/ 4.7 H n q. F Hence n 1 the pair Z n, F 0 is a martingale and the pair 0 U n, where Zn = n/ H n / q,, F, where Un = n/ H n 0 is the reverse martingale. In particular is a martingale and, F is a reversed martingale. >0 q

10 PAWEL J. SZABLOWSKI 4.8a 4.8b 5 We have also: 1 q + γ δ δ E F δ, +γ = γ δ + γ δ + δ δ + γ +γ a.s. E F δγ δ, +γ = δ + γ 1 q + γ + qδ + 1 q + qδ +γ γ + q + 1 δ +γ + δ + γ Corollary 1. Let be a q Wiener process and Y related to it q, α OU process. We have i For q 1, 1 both and Y processes do not allow continuous path modifications. ii 4.9 4.10 F =, a.s. > > 0 : E 4 F = > > 0 : E Hence does not have independent increments. iii δγ var F δ, +γ = 1 q + γ + qδ 1 q + + q + 1 q, a.s 1 1 q +γ δ δ + γ + γ δ δ +γ. δ + γ Thus following [?] we deduce from iii that is a quadratic harness with parameters introduced in[?] θ = η = = sigma = 0 thus q commutation equation is very simple in this case, namely [x.y] = I. Remark 4. From assertion 4. of the Theorem 5 it follows that, F, 0, F are martingales and 0 /, F >0 / 1/, F are >0 reversed martingales. Thus if main result of Weso lowski s paper [9] was true we would deduce that is the Wiener process. We will show that Weso lowski made a mistake in his calculations. Following these martingale properties we have:for all 0 < < : E F = a.s., E F =. In order to simplify calculations we will assume that the following equalities are satisfied with probability 1. Hence E F = + or var F =. Further from reversed martingale property we have : E F =, hence E F = and E 1 F = 1, hence E F = +. Now mimicking way of reasoning of Weso lowski we have: E = E E F = E 4 + E = E E F = E 4 + E. Thus: E 4 + = E4 +..

GENERALIZATIONS OF BRYC S PROCESSES 11 Dividing both sides by we get: E 4 + 1 = E4 + 1. Thus we see that E 4 = c for some constant c. We know that c = + q.see e.g. the proof of Corollary 1 and the fact that H 4 = 4 3 + q + q + 1 + q + q and the fact that EH 4 = 0 and E =. Hence following Weso lowski we get E 3 = E4 = + q, E = + q +, E 3 = E 4 = + q. So far the Weso lowskis and our calculations yield the same partial results. Now: E 4 = E 4 4E3 + 6E 4E 3 + E4 = + q 4 + q + 6 + q + 4 + q + + q = + 3q q = + q + 1 q. In Weso lowski s paper it is wrongly written that E 4 = c + c 3 +, where c = E 4 /. Following Theorem 5 we get: Lemma 3. Let = 0 be a q Wiener process. For s R and 0 < <, we have: E 1 1 q sq k + 1 qs q k 1 F 4.11 = k=0 1 1 q sq k + 1 qs q k 1, a.s. k=0 Remark 5. Bryc and Weso lowski in [6] also introduced q Wiener processes. However the process introduced by them was a particular element of general class of q Meixner processes. To define q Meixner processes they needed 5 assumptions to be satisfied. Hence according to Bryc and Weso lowski q Wiener process is a L process that satisfies 5 the following conditions. Namely: Let t t 0 be L -process with zero expectation and for all 0 s < t < u 1 E t s = min t, s E t F s, u = u t u s s + t s u s t. 3 E t F s, u = A1 s, t, us+a s, t, uu+b s, t, u s u +D s, t, u+ E s, t, u s + F s, t, u s for some nonrandom functions A 1, A, B, C, D, E, F of s, t, u. 4 E t F s = s + t s 5 var t F u = t u u t. The procedure presented in this paper and in particular remark only conditions are needed to define q Wiener process in the class of L processes. Namely, conditions. and 3. for special s = t δ, u = t+γ and assuming A 1 = A and E = F = 0. To see this apply change of time scale defined by 4.5 to conditions 4.3 and 4.4 to obtain conditions 4.8a and 4.8b with γ = expnαdexpαd 1 and δ = exp n 1αdexpαd 1. All other properties mentioned in theorem 5 and assumed in the paper [6] that is conditions 1., 4., 5. follow these two plus

1 PAWEL J. SZABLOWSKI linear independence discussed in remark. In particular properties 4. and 5. follow assertion 4.of the Theorem 5 considered for n = 1 and. Bryc & Weso lowski indicate that in general that is if for these particular δ and γ we would have A 1 A and E + F 0, then conditions 1. 4. 5. could not be deduced from. and 3. One of the examples is discussed by them so called q-poisson process. Above mentioned δ, γ, and A 1 = A when applied to time transform defined by 4.5 imply symmetry with respect to time direction of conditions imposed on first two conditional moments. This symmetry assumption thus turns out extremely important, and crucial. Remark 6. Following Remark 3 we deduce that q Wiener processes do not exist for q > 1. 5. Proofs of the results Proof of Lemma 1. 1. Remembering that H 4 x q = x 4 3+q+q x +1+q+q thus consequently x 4 = H 4 x q + 3 + q + q H x q + + q and we see that En 4 = +q since EH 4 m = EH m = 0. Now using the fact that EH n F m = ρ H m for m < n we get: E n m = E m H n q + 1 = 1 + E meh n q F m = 1 + ρ n m EmH m = 1 + ρ n m E m 4 m To get 3.1c. we have: = 1 + ρ n m + q 1. E n n j n+k = E n n j E n+k F n = ρ k E 3 n n j = ρ k+j E 3 ne n j F n = ρ k+j E 4 n = ρ k+j + q. Proof of Proposition 1. For fixed natural numbers k, j let us denote Z i = E n+i F n k, n+j. Let Z =[Z i ] j 1 i= k+1. Now notice that for each coordinate of vector Z we have: Z i = E i F n k, n+j = EE i F i F n k, n+j = ρ 1 + ρ E i 1 F n k, n+j + E i+1 F n k, n+j ρ 1+ρ Z i 1 + Z i+1 if i k + 1 j 1 ρ = 1+ρ n k + ρ 1+ρ Z k+ if i = k 1. ρ 1+ρ n+j + ρ 1+ρ Z j if i = j 1 Hence we have a vector linear equation Z = AZ+ ρ 1 + ρ n k + n+j, where A = ρ 0 1+ρ... 0 ρ 1+ρ 0... 0......... 0 0... 0 ρ 1+ρ.

GENERALIZATIONS OF BRYC S PROCESSES 13 Now notice that matrix I A where I denotes unity matrix is nonsingular. This is so because sum of the absolute values of elements in each row of the matrix is less than 1 which means that eigenvalues of the matrix I A are inside circle at center in 1 and radius less than 1, thus nonzero. Consequently each component of the vector Z is a linear function of n k and n+j. Having linearity of E i F n k, n+j with respect to n k and n+j we get 3.. Proof of Proposition. Let us denote m i,j = E n+i n+j F n l, n+k, for i, j = l, l + 1,...,k 1, k. Notice that using 3. we get m i,j = m j,i and that m l,j = n l E n+j F n l, n+k = ρl+j 1 ρ k j 1 ρ k+l n l + ρk j 1 ρ l+j 1 ρ k+l n l n+k and m k,j = n+k E n+j F n l, n+k = ρl+j 1 ρ k j 1 ρ k+l n l n+k + ρk j 1 ρ l+j 1 ρ k+l n+k m l, l = n 1, m k,k = n+k. Besides we have for i, j = j + 1,...,k 1 and i j m i,j = E n+i n+j F n l, n+k = E n+i E n+j F n+j F n l, n+k = ρ 1 + ρ m i,j 1 + m i,j+1 and if i = j m i,i = E n+i F n l, n+k = E E n+i F n+i Fn l, n+k = Am i 1,i 1 + m i+1,i+1 + Bm i 1,i+1 + C Notice also that we have in fact l + k unknowns and l + k linear equations. Moreover random variables m ij are well defined since conditional expectation is uniquely defined up to set of probability 1. Thus we get the main assertion of the proposition. Now we know that for some A j, B j, C j we have 5.1 E n F n j, n+k = A 1 jk n j + A jk n+k + B jk n j n+k + C jk. First thing to notice is that 1 = A 1 jk + A jk + ρj+k B jk + C jk. Secondly let us multiply 5.1 by n j, n+k and n j n+k and let us take expectation of both sides of obtained in that way equalities. In doing so we apply assertions of Lemma 1. In this way we will get three equations: 1 + ρ j 1 + q = A 1 jk 1 + ρ k 1 + q = A 1 jk + q + A jk + B jk ρ j+k + q + 1 A 1 1 + ρ j+k 1 + q 1 + ρ j+k 1 + q jk A jk ρj+k B jk + A jk + q + B jk ρ j+k + q + 1 A 1 jk A jk ρj+k B jk ρ j+k + q = A 1 jk ρj+k + q + A jk ρj+k + q + B jk 1 + ρ j+k 1 + q + ρ j+k 1 A 1 jk A jk ρj+k B jk. Solution of this system of equations is, as it can be easily, checked 3.4. Secondly consider discrete time random field Z = {Z k } k Z such that Z k = kj. Obviously

14 PAWEL J. SZABLOWSKI we have EZ k F k = E kj F k 1j k+1j = ρj 1 + ρ j k 1j + k+1j = ρ j 1 + ρ j Z k 1 + Z k+1, and E Z k F k = E kj F k 1j k+1j = A j k 1j + k+1j + B j k 1j k+1j + C j = A j Z k 1 + Z k+1 + Bj Z k 1 Z k+1 + C j. where A j = A 1 jj = A jj, B j = B jj. Thus Z is a Bryc random field with different parameters. Notice that one dimensional distributions of processes and Z are the same. Hence parameters q for both processes Z and are the same. On the other hand parameter ρ Z of the process Z is related to parameter ρ of the process by the following relationship ρ Z = EZ 0 Z 1 = E 0 j = ρ j. Thus applying formulae.7 we get 3.4. Proof of Theorem 4. The proof of this theorem is in fact contained in above mentioned papers. For the proof of 5. notice that from 4. it follows that n, m 1, s, t R : EH m Y t q H n Y s q = δ n m e nα s t EHn Y t q. Further we use spectral decomposition theorem. Proof of Theorem 5. 1. Suppose <, cov, = EY log /α Y log /α = e α log /α log /α = e log// = =.. & 3. are just assertion 1 of theorem 4. 4. Notice that / = Y log /α.we have using assertion 5 of Theorem 4: E n/ H n / F = n/ E H n Ylog /α Flog /α = n/ e αnlog /α log /α H n Ylog /α = n/ H n /. a.s. Similarly E n/ H n q = n/ E H n Ylog /α q F log /α = F n/ e αlog /α log /α H n Ylog /α q = n/ n Hn Ylog /α q = n/ H n q. 5. to get first part we proceed as follows: E Ylog /α F log δ/α, log+γ/α = δ/ 1 / + γ / + γ1 δ/ Y log δ/α + Y log+γ/α 1 δ/ + γ 1 δ/ + γ = γ δ + γ δ + δ δ + γ +γ

GENERALIZATIONS OF BRYC S PROCESSES 15. To get second part we have: E F δ,+γ = E Y δ 1 +γ 1 q +γ = Y 1 δ +γ 1 q δ +γ +γ 1 δ 1 q δ + Y 1 δ 1 q δ q + 1 + δ +γ 1 δ +γ log /α F log δ/α, log+γ/α +γ 1 +γ +γ 1 δ 1 q δ +γ Keeping in mind that = Y log /α we get: log δ/α log+γ/α Y log δ/α Y log+γ/α + E F δ,+γ γ 1 q + γ = δ + γ 1 q + γ + qδ δ 1 δ 1 +γ. 1 q δ +γ δ 1 q + qδ + δ + γ 1 q + γ + qδ +γ q + 1δγ + δ + γ 1 q + γ + qδ δ +γ δγ + 1 q + γ + qδ. Proof of Corollary 1. Recall that from well known Levy s characterization of Wiener process it follows that a process = t t 0 with continuous a.s. trajectories is a Wiener process iff condition t, F t, t t, F t are martingales. Recall that H x q = x 1. From assertion 3. of the above Proposition it follows that for any q Wiener process with q < 1, both, F and, F are martingales. Hence iff they had continuos trajectories they would be ordinary Wiener processes i.e. q = 1. A contradiction. Since we obtain q, α OU process from q Wiener process by continuous argument transformation, we deduce that any q, α OU process with q < 1 cannot have continuous a.s. trajectories. To prove ii we calculate: E F = E F + = + =. E 4 F = E 4 F 4 E 3 F + 6 E F 4 3 E F + 4 = E 4 F 4 E 3 F + 6 4 + 6 44 + 4 = E 4 F 4 E 3 F + 3 4 + 6.

16 PAWEL J. SZABLOWSKI Now recalling that H 3 x q = x 3 +qx and H 4 x q = x 4 3+q+q x + 1+q+q we see that 3/ H 3 q = E 3/ H 3 q = E 3 F + qe F and H 4 q = E H t 4 q F F = E 4 F 3 + q + q E F + 1 + q + q. Thus E 3 F = + q + 3 + q = + q + 3 and E 4 F = 3 + q + q + 1+q+q + 4 3+q+q + 1 + q + q = 4 + 3 + q + q + 1 + q + q. So E 4 F = 3 4 + 6 + 4 + 3 + q + q + 3 + q + q 1 + q + q 4 4 4 + q = 1 q + + q 1 + q + q = 1 q + + q + 1 q. If had independent increments then E 4 F would not depend on. Note that for q = 1 E 4 F = 3 = 3E F. iii We have to subtract 4.8b square of 4.8b. var F δ,+γ δγ = 1 q + γ + qδ + + δ γ 1 q + γ δ + γ 1 q + γ + qδ γ δ + γ + +γ δ 1 q + qδ δ δ + γ 1 q + γ + qδ + δ + γ q + 1δγ δ +γ δ + γ 1 q + γ + qδ γδ δ + γ. After some simplifications we get = var F δ,+γ δγ = 1 q + γ + qδ 1 qδγ + γ 1 qδγ δ + γ + δ 1 q + γ + qδ δ + γ + δ 1 q + γ + qδ +γ 1 qγδ δ + γ γ + δ 1 q + γ + qδ δ +γ δγ = 1 q + γ + qδγ + δ γ + δ + 1 q + γ δ + δ +γ δ + γ δ +γ δγ 1 q + γ + qδγ + δ γ + δ + 1 q +γ δ δ +γ + γ δ

GENERALIZATIONS OF BRYC S PROCESSES 17 Proof of Lemma 3. We have P t x = P Y log /α x = P Y log /α x/. Knowing that f H is a density of Y t we get immediately assertion 1. To get assertion we have P x = y = P x + y = y = P Y log /α x + y Y log /α = y = P Y log /α x+y Y log /α = y. Now we recall that f x y, q, e α s t is the density of Y t given Y s = s. For q < 1 to get assertion 3 first we multiply both sides of 4.6 by s n / [n] q and sum over n from 0 to. Changing order of summation and conditional expectation and applying formula 1.4 of [3] giving generating function of q Hermite polynomials, leads to 4.11. For q = 1 we multiply both sides of 4.6 by s n /n! and after summing up with respect to n we apply well known formula for generating functions of Hermite polynomials. References [1] Bryc, W. 001 Stationary random fields with linear regressions. Annals of Probability 9, No. 1, 504-519. [] Bryc, W. 001 Stationary Markov chains with linear regressions. Stochastic Processes and Applications 93, 339-348. [3] Bryc, W., Matysiak, W., Szab lowski, P.J. 003 Probabilistic aspects of Al-Salam Chihara polynomials. Proceedings of the American Mathematical Society 133, 117-1134. [4] Matysiak, W., Szab lowski, P.J. 00 A few remarks on Bryc s paper on random fields with linear regressions. Annals of Probability 30, No. 3, 1486-1491. [5] Matysiak, W., Szab lowski, P. J. 005, Bryc s Random Fields: The Existence and Distributions Analysis, Ariv:math.PR/math/050796 [6] Bryc, W., Weso lowski, J. 005, Conditional Moments of $q$-meixner Processes, Probab. Theory Rel. Fields 131, 415-441 [7] Bryc, W., Bożejko, M. 005, On a Class of Free Levy Laws Related to a Regression Problem, ariv:math.oa/0410601 [8] Bryc, W., Matysiak, W., Weso lowski, J. 005, Quadratic Harnesses, $q$-commutations, and Orthogonal Martingale Polynomials, ariv:math.pr/0504194 [9] Weso lowski, Jacek 1990, A Martingale Characterization of the Wiener Process, Statistics & Probability Letters 10, 13-15 [10] Bożejko, M. Kümmerer, B., Speicher, R. 1997, $q$-gaussian Processes: No-Commutative and Classical Aspects, Comm. Math. Phys. 1851, 19-154 [11] Manuy, R., Yor, M. 005, Levy Bridges and Monsieur Jourdain, Stochastic Processes and Their Applications 115, 39-338 [1] Anshelevich, M. 004, $q$-levy Processes, ariv:math:oa/0309147 [13] Bryc, W., Weso lowski, Jacek, Bi-Poisson Process, arxiv.org:math/pr/040441 [14] Bryc, W. 001, Classical versions of Gaussian processes. Conditional moments and Bell s inequality, Comm. Math. Phys. 19, 59-70 Acknowledgement 1. References Acknowledgement. The Author would like to thank W lodek Bryc and Jacek Weso lowski for their critical remarks concerning first version of the paper. Department of Mathematics and Information Sciences, Warsaw University of Technology, pl. Politechniki 1, 00-661 Warszawa, Poland E-mail address: pszablowski@elka.pw.edu.pl