THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
|
|
- Hester Stewart
- 5 years ago
- Views:
Transcription
1 THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT and Ronad J. Stern Department of Mathematics, Concordia University, Montrea, Canada H4X 1J7 and Michae J. Tsatsomeros Department of Mathematics, University of Connecticut, Storrs, CT Research supported in part by US Air Force Research Grant No AFOSR and by NSF Grant No DMS This author woud aso ike to thank NSERC for making it possibe for him to visit Ronad J. Stern in Montrea. Research supported by the Natura Sciences and Engineering Counci of Canada, grant No A4641. Research supported in part by US Air Force Research Grant No AFOSR and by NSF Grant No DMS
2 Abstract Let A be an n n essentiay nonnegative matrix and consider the inear differentia system ẋ(t) = Ax(t), t 0. We show that there exists a constant h(a) > 0 such that the trajectory emanating from x o reaches R+ n at a finite time t o = t(x o ) 0 if and ony if the sequence of points generated by a finite differences approximation from x o, with time step 0 < h < h(a), reaches R+ n at a finite index k o = k(x o ) 0. This generaizes and strengthens earier resuts of two of the authors, where some additiona spectra restrictions were imposed on A. Our proof makes use of the existence of a basis of nonnegative vectors to the Perron eigenspace.
3 1 Introduction Let A be an n n matrix and consider the inear differentia system The soution to (1), given by ẋ(t) = Ax(t), t 0. (1.1) x(t) = e ta x o, t 0, (1.2) wi be referred to as the trajectory of the differentia system emanating from x o. A set Γ R n is caed positivey invariant with respect to A if e ta (Γ) Γ, t 0. This condition has the impication that once a trajectory emanating from x o reaches Γ in some finite time, it remains in Γ thereafter. The set of a initia points of trajectories which reach Γ is caed the reachabiity set for Γ under A and is denoted by X A (Γ). Ceary, X A (Γ) = {x o R n ( t o = t(x o ) 0)( t t o )[e ta x o Γ]} = t 0 e ta Γ. (1.3) In the case when Γ is a proper cone, namey, a cosed, convex, pointed and soid cone, then it was shown in Neumann and Stern [?] that X A (Γ) is a convex cone itsef, which contains Γ, but which is not necessariy cosed or pointed. We sha denote the cosure of X A (Γ) by X A (Γ). In this paper A = (a ij ) wi aways denote an essentiay nonnegative matrix, that is, a matrix whose off diagona entries are nonnegative. This condition is known to be equivaent to the positive invariance of the nonnegative orthant R+ n with respect to A (see Beman [?] and Birkhoff and Varga [?]). In Neumann and Stern [?] and in Berman, Neumann, and Stern [?], expicit formuas were derived for X A (R+) n under additiona spectra assumptions on A. In Neumann and Stern [?], a numerica characterization of X A (R+) n based on a Cauchy Euer approximation to the soution to (1.1) was deveoped. However, there, in addition to A being essentiay nonnegative, it was assumed that A has a rea spectrum or that it possesses a stricty positive generaized eigenvector corresponding to its Perron root. For further description of these reachabiity resuts the reader is referred to 1
4 Chapter 6 of Berman, Neumann, and Stern [?]. In the present work we pose the questions of generaizing and strengthening the aforementioned numerica characterization of X A (R+) n by assuming that A is ony essentiay nonnegative and by numericay characterizing X A (R+) n and not just its cosure. For the sake of competeness we first give one of the main resuts in Neumann and Stern [?]. Define h(a) = sup{h min (1 + ha ii) > 0}. (1.4) 1 i n It is cear that h(a) can be infinite and that for any 0 < h < h(a) the matrix I + ha is nonnegative. For any such h we define the discrete reachabiity set of R+ n to be X A,h (R+) n = {x o R n ( k o = k(x o ) 0)( k k o )[(I + ha) k x o 0]} (1.5) which can be shown to be an invariant cone under I +ha. The sequence of vectors x k = (I + ha) k x o, k = 0, 1,..., wi be referred to as the discrete trajectory emanating from x o. This sequence represents an approximation to the soution of (1.1), in the sense of Cauchy Euer, at a discrete set of times. Next, it is known by the cassica Perron Frobenius theory that the Perron root λ 1 := max{reµ µ σ(a)} σ(a), (1.6) where σ(a) denotes the spectrum of A. Reca now that the Perron (generaized) eigenspace of A is defined to be the A invariant subspace W A = N((λ 1 I A) p ), (1.7) where p is the mutipicity of λ 1 in the minima poynomia of A. The resut of Neumann and Stern which we wish to quote here can now be stated as foows. Theorem 1.1 ([?], Theorem 2.2) Suppose A is an n n essentiay nonnegative matrix. If W A int(r n +), (1.8) then for any 0 < h < h(a) which satisfies we have that 1 + hλ 1 > 1 + hµ, µ σ(a) \ {λ 1 }, (1.9) X A,h (R n +) = X A (R n +). (1.10) 2
5 We shoud emphasize the fact that the time step h in (1.9) depends ony on A and is not necessariy sma. In other words, Theorem 1.1 provides us with a way of testing whether a point ies in X A (R+) n regardess of how much the continuous and discrete trajectories emanating from x o diverge from each other. The main objectives of this paper are to extend and improve the resut of Theorem 1.1. This wi be achieved in two principa stages. In the first stage, Theorem 3.1, we sha show that the concusion of Theorem 1.1 hods for any essentiay nonnegative matrix. Our proof of this resut wi be aided by a bock trianguar form, due to Hartwig, Neumann, and Rose [?], to which, actuay, any nonnegative matrix can be symmetricay permuted. For the sake of competeness, this form is dispayed in Lemma 2.1. In the second stage, Theorem 3.3, we sha show that with a possibe sight additiona restriction on the supremum of h, the equaity in (1.10) hods even when the cosure signs are removed. The proof of Theorem 3.3 requires that we further expose geometrica properties of the discrete reachabiity cone X A,h (R+). n We do this in Lemma Preiminaries Suppose that B is an n n nonnegative matrix and et W B = N((ρ(B)I B) p ) be the generaized Perron eigenspace of B corresponding to its spectra radius ρ(b), where p is the degree of ρ(b) in the minima poynomia of B. The emma to foow assigns to B a certain bock upper trianguar form which wi be essentia to the generaization of Theorem 1.1. It is based on Lemma 3.1 and Coroaries of Hartwig, Neumann, and Rose [?]. We comment that this form coud aso be derived from proofs of the existence of a nonnegative basis for W B, which are based on the Frobenius norma form for B, as given in Rothbum [?] and Richman and Schneider [?]. However, the proof of our main resut in Section 3 does not require that we consider the Frobenius norma form of a matrix. For a nonnegative vector w, we sha et ν(w) denote the number of positive entries in w. Lemma 2.1 Any n n nonnegative matrix B is permutationay simiar 3
6 to a bock upper trianguar matrix B 11 B12... B1p B B2p.... (2.1) 0 Bpp with the foowing properties. If for each j {1, 2,..., p} we set B j = B jj... Bjp Bpp R q,q +, (2.2) then each diagona bock B jj in (2.1) is of size k j k j, where k j = max{ν(w) w W Bj R q +}, q = n k 1... k j 1 and it possesses a stricty positive generaized eigenvector corresponding to ρ( B jj ). Moreover, and u j W Bjj intr k j + (2.3) ρ( B j ) = ρ( B jj ), ρ( B j ) > ρ( B j+1 ), 1 j p 1, (2.4) where m j = index ρ( Bjj ) B jj. W Bj = {u R q u = [ū T, 0] T, ū W Bjj } (2.5) ( B jj ρ( B jj )I) m j 1 u j R k j + \ {0}, (2.6) We concude this section by recaing that R+ n admits an aternative representation as the intersection of n cosed hafspaces (e.g. Rockafear [?]) as foows : Let.,. denote the usua inner product in R n and suppose that {ν i i = 1, 2,..., n} is the set of a outward unit normas to R+ n (that is, ν i = [0,..., 0, 1, 0,..., 0] T, where the nonzero entry occurs in the i-th position). Then, n R+ n = {z R n ν i, z 0}. (2.7) i=1 4
7 3 The Main Resuts Let A be an n n essentiay nonnegative matrix and et h(a) be as defined in (1.4). In what foows we sha assume, without oss of generaity, that the matrix B = B(h) = I + ha is aready in the upper trianguar bock form given in Lemma 2.1. Otherwise, our considerations appy to a permutation simiarity of A. We are now ready to state the first main resut of this paper. Theorem 3.1 Let A be an n n essentiay nonnegative matrix. Then for a 0 < h < h(a), X A,h (R n +) = X A (R n +). (3.1) Proof We wi first show that X A (R+) n X A,h (R+). n Suppose that there exists z X A (R+) n such that Define the vector z X A,h (R n +). (3.2) u = [u T 1, u T 2,..., u T p ] T intr n +, (3.3) where the vectors u j, j = 1, 2,..., p are defined as in Lemma 2.1. Then there exists a stricty decreasing sequence of positive numbers {ɛ m } m=1, with ɛ m 0, such that z + ɛ m u X A,h (R n +), m 1. (3.4) For an arbitrary but fixed m, the excusion in (3.4) means by (2.7) that for each i 1 there exists an outward unit norma ν (m,i) to R n + such that ν (m,i), B i (z + ɛ m u) > 0. (3.5) Since there are ony a finite number of such normas there exists a norma, say ν (m), and a sequence {i (m) k } k=1 such that ν (m), B i(m) k (z + ɛ m u) > 0, k 1. (3.6) Suppose that z and ν (m) are partitioned in conformity with (2.1). Assume further that the nonzero entry of ν (m) occurs in the th bock, namey, z = [z T 1,..., z T,..., z T p ] T and ν (m) = [0,..., 0, (ν (m) ) T, 0,..., 0] T. (3.7) 5
8 Consider now the traiing submatrix (viz. Lemma 2.1) B = I + ha R q,q, q = n k 1... k 1 (3.8) and construct the foowing q vectors: ν (m) := [(ν (m) ) T, 0,..., 0] T, z (m) := [z T,..., z T p ] T, ū (m) := [u T, 0,..., 0] T R q +, û (m) := [0, u T +1,..., u T p ] T R q +. (3.9) Then, k 1, ν (m), B i(m) k ( z (m) + ɛ m ū (m) ) ν (m), B i(m) k ( z (m) + ɛ m ū (m) ) + Thus by (3.5) and (3.10), + ν (m), ɛ m B i(m) k û (m) = ν (m), B i(m) k ( z (m) + ɛ m (ū (m) + û (m) )) = ν (m), B i(m) k (z + ɛ m u). (3.10) ν (m), B i(m) k ( z (m) + ɛ m ū (m) ) > 0, k 1. (3.11) Next, as e ta is upper trianguar t 0, a necessary condition for the assumption z X A (R n +) to hod true is that z X A (R q +). (3.12) Reca now that u intr k +. It was shown in Neumann and Stern [?] that intx A (R k +)=X A (intr k +) and, furthermore, once a trajectory has entered intr k + it cannot in finite time reach the boundary of R k +. Consequenty ν (m), ɛ m e taū (m) < 0, t 0. (3.13) Then, (3.12) and (3.13) have the impication that there exists a sufficienty arge t m 0 such that ν (m), e ta ( z (m) + ɛ m ū (m) ) < 0, t t m. (3.14) We sha next show that (3.11) and (3.14) are incompatibe as we vary m. But first, et λ := max{reµ µ σ(a )} (3.15) 6
9 and write B = (1 + hλ )I + h(a λ I). (3.16) Observe then that, by Lemma 2.1, ū (m) W A R+ q and consider the resoution of z (m) + ɛ m ū (m) into z (m) + ɛ m ū (m) = ( w (m) + ɛ m ū (m) ) + r (m), (3.17) where r (m) is the projection of z (m) (and hence of z (m) + ɛ m ū (m) ) onto U A, the join of a eigenspaces of A corresponding to eigenvaues µ λ aong W A. We caim that, independenty of m, ν (m), w (m) + ɛ m ū (m) = 0. (3.18) Suppose to the contrary. As A λ I is nipotent on W A, et p 1 0 be the argest integer such that and ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ) = 0 (3.19) ν (m), (A λ I) j ( w (m) + ɛ m ū (m) ) = 0, j > p 1. (3.20) If A is a nonnegative matrix, so that h(a) =, then for any h (0, ) the eigenvaues of I + ha satisfy 1 + hλ > 1 + hµ, µ σ(a ) \ {λ } (3.21) or ese (3.15) is vioated. Simiary, if h(a) <, (3.21) hods true for a 0 < h < h(a). Consequenty, the restriction of B /(1 + hλ ) to U A is a convergent matrix, namey, [ im k B 1 + hλ ] i (m) k r (m) = 0. (3.22) Then for i (m) k > p 1 by (3.19) and (3.20) we can write that = + ν (m), B i(m) k ( z (m) + ɛ m ū (m) ) = ( (m) ) i k (1 + hλ ) i(m) p k 1 h p 1 ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ) + p1 p 1 1 j=0 ( i (m) k j ) (1 + hλ ) i(m) k j h j ν (m), (A λ I) j ( w (m) + ɛ m ū (m) ) + + ν (m), B i(m) k r (m). (3.23) 7
10 Note that as k and for a 0 j < p 1, Aso, by (3.22), ( (m) ) / ( (m) ) i k i k j p1 0. (3.24) im k (1 + hλ ) p 1 ( i (m) k p1 ) h p 1 ν (m), [ B 1 + hλ ] i (m) k r (m) = 0. (3.25) Now, since 1 + hλ > 0, upon taking k sufficienty arge, reations (3.11) and (3.23) (3.25) have the impication that ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ) > 0. (3.26) On the other hand, the restriction of A λ I to U A is a stabiity matrix and so im t et (A λ I) r (m) = 0. (3.27) Now, by (3.19) and (3.20) one obtains ν (m), e ta ( z (m) + ɛ m ū (m) = = e tλ ν (m), e t(a λ I) ( w (m) + ɛ m ū (m) ) + e tλ ν (m), e t(a λ I) r (m) = tp 1 e tλ ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ) p 1! + e tλ p 1 1 j=0 t j j! ν(m), (A λ I) ( w (m) + ɛ m ū (m) ) + e tλ ν (m), e t(a λ I) r (m). (3.28) But then, for t sufficienty arge, (3.14), (3.27) and (3.28) give us that ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ) 0, (3.29) a contradiction to (3.26) showing that (3.18) is vaid. On setting w (m) := [0, w (m) ] T R n, by (3.17) we obtain that ν (m), w (m) + ɛ m ν (m), u = 0. (3.30) 8
11 Now vary m. Then, since the number of outward normas is finite, there must exist indices m 2 > m 1 1 such that ν := ν (m1) = ν (m2), and hence w := w (m1) = w (m2), so that which, as ɛ m1 < ɛ m2, is ony possibe if ν, w + ɛ mi ν, u = 0, i = 1, 2 (3.31) ν, u = 0, (3.32) which is not possibe as u intr n +. This shows that X A (R+) n X A,h (R+) n (3.33) and competes the first part of the theorem. We sha next prove the reverse containment, namey, that X A,h (R+) n X A (R+). n (3.34) Let z X A,h (R+) n such that z X A (R+). n (3.35) Suppose that u and {ɛ m } m=1 are chosen as before so that z + ɛ m u X A (R n +), m 1. (3.36) Consider m arbitrary but fixed. The excusion in (3.36) means that for any unbounded sequence of stricty increasing positive times {t (m) j } j=1, there exist ν (m), ν (m), z (m), ū (m), defined in a manner simiar to the first part of the theorem (viz. equations (3.3) (3.11)), so that ν (m), e t(m) j A ( z (m) + ɛ m ū (m) > 0, j 1. (3.37) Next, by equations (2.4) and (2.6) of Lemma 2.1 we can deduce that B j u intr k +, j 0. Observe now that the assumption z X A,h (R+) n has the impication that there exists a sufficienty arge exponent j m so that ν (m), B j ( z(m) + ɛ m ū (m) ) < 0, j j m. (3.38) 9
12 We sha now show that, as we vary m, equations (3.37) and (3.38) are incompatibe. For this, et us resove z (m) + ɛ m ū (m) as in (3.17). Then, an anaysis simiar to the one for equations (2.23) and (2.28) shows that or ese, for some p 1 0, ν (m), w (m) + ɛ m ū (m) = 0 (3.39) sgn ν (m), e t(m) j A ( z (m) +ɛ m ū (m) ) = sgn ν (m), (A λ I) p 1 ( w (m) +ɛ m ū (m) ) (3.40) and sgn ν (m), B j ( z(m) + ɛ m ū (m) ) = sgn ν (m), (A λ I) p 1 ( w (m) + ɛ m ū (m) ), (3.41) which contradict (3.37) and (3.38). Then, in a simiar fashion to (3.30) and (3.31), upon varying m, it foows from (3.39) that for some outward norma ν to R+, n ν, u = 0. This is a contradiction to (3.3) which competes the proof of the theorem As a consequence of the fact that in genera the reachabiity cone is not cosed, in order to achieve our goa of providing a numerica characterization for the eements of X A (R+) n we must improve the resut of Theorem 3.1 by distinguishing between those boundary points which reach R+ n and those which do not. The set of a such points was termed in [?] the effective part of the boundary of X A (R+). n In the remainder of this section we sha utiize some geometrica properties of the continuous and discrete reachabiity cones to show that for a but a finite number of vaues of h in (0, h(a)) we can strengthen the resut of Theorem 3.1, given in equation (3.1), by proving that X A (R+) n = X A,h (R+). n (3.42) For this purpose we require the foowing emma, some of whose causes have aready been estabished in the iterature. For any set Γ R n and A R n,n we define core A (Γ) to be the maxima subset of Γ which is positivey invariant with respect to A and core A,h (Γ) to be the maxima subset of Γ which is invariant with respect to B = B(h) = I + ha. As usua, we denote the boundary of Γ by Γ. Lemma 3.2 Let A be an n n essentiay nonnegative matrix and suppose that 0 < h < h(a) is chosen so that B = B(h) = I + ha is invertibe. Then the foowing hod. 10
13 (i) X A (intr n +) = intx A (R n +). (ii) X A (R n +) is positivey invariant with respect to A. (iii) X A,h (intr n +) = intx A,h (R n +). (iv) X A,h (R n +) is invariant with respect to B. (v) X A (R n +) X A (R n +) = X A (core A ( R n +)). (vi) X A,h (R n +) X A,h (R n +) = X A,h (core A,h ( R n +)). (vii) core A ( R n +) is the union of a the positivey invariant faces of R n +. (viii) core A,h ( R n +) = core A ( R n +). Proof Caims (i), (ii), (v) and (vii) are due to Neumann and Stern [?] (see aso Berman, Neumann, and Stern [?]). We comment that it is evident from the proof of (i) in [?] that once a trajectory emanating from a point in intx A (R+) n has entered intr+, n it remains in intr+ n for a finite time. Now et h be as prescribed. Then the proofs of (iii), (iv), and (vi) are very simiar to their respective counterparts in the continuous case and therefore they are omitted. To show (viii) consider a point z core A ( R+). n By (vii) z beongs to some positivey invariant face F of R+. n Thus, since F is a poyhedra proper cone in the span of F, the restriction of A to F must be essentiay F nonnegative, namey for some α 0, (A + αi)f F (3.43) (see Schneider and Vidyasagar [?] and Stern [?]). Then, for any x F R+ n we have that y:=(a + h 1 I)x R+ n and by (3.43) there exists some β 0 such that y + βx F and so, by the definition of a face, we have the impication 0 y y + βx F y F. (3.44) Thus, hy = Bx F, that is, for any choice of h (0, h(a)), F is B=B(h) invariant. This means that z core A,h ( R n +) and hence core A ( R n +) core A,h ( R n +). (3.45) 11
14 Conversey, suppose that z core A,h ( R+) n and note that, in particuar, z R+. n Thus if z core A ( R+), n then we must have that z X A (intr+). n It foows from (3.1) that intx A (R n +) = intx A,h (R n +) (3.46) which, combined with caims (i) and (iii), has the impication that X A (intr n +) = X A,h (intr n +). (3.47) But then z X A,h (intr+) n contradicting the fact that z core A,h ( R+). n This competes the proof of the emma We can now prove our second main resut. Theorem 3.3 Let A be an n n essentiay nonnegative matrix. Then for a 0 < h < h(a) such that B = B(h) = I + ha is invertibe we have X A (R n +) = X A,h (R n +). (3.48) Proof Let h be as prescribed in the statement of the theorem. First we wi show that X A (R n +) X A,h (R n +). (3.49) By Theorem 3.1 the interior points of the two cones coincide so it suffices to show that X A (R n +) X A (R n +) X A,h (R n +). (3.50) Suppose that z X A (R+) X n A (R+). n Then, by (v) and (vii) of Lemma 3.2, z X A (F 1 ), where F 1 is some positivey invariant and hence, according to the expanation invoving (3.43) and (3.44), a B invariant face of R+. n Without oss of generaity assume that F 1 = {x = (x 1, x 2,..., x n ) T x j = 0, j > k for some k 1} (3.51) in which case, as BF 1 F 1, B has necessariy the upper trianguar form [ ] [ ] B11 B 12 A11 A = I + h 12, (3.52) 0 B 22 0 A 22 where A 11 R n k,n k and A 22 R k,k. Otherwise, our considerations appy to a permutation simiarity of A. Then z has the form z = [z T 1, 0] T, z 1 X A11 (R n k + ). (3.53) 12
15 Thus, since h(a) h(a 11 ), Theorem 3.1 appied to A 11 yieds that and so, by (3.53), z 1 X A11,h(R n k + ) = X A11 (R n k + ) (3.54) z X A,h (F 1 ) = X A (F 1 ). (3.55) If now z 1 intx A11 (R+ n k )=intx A11,h (R+ n k ), then z X A,h (F 1 ) which impies that z X A,h (R+) n (3.56) and we are done. Suppose then that z 1 X A11 (R n k + ) X A11 (R n k + ). (3.57) We may then appy a simiar anaysis to z 1 to show that there exists a B invariant face F 2 of F 1 such that z X A,h (F 2 ) = X A (F 2 ). (3.58) Then, either z X A,h (F 2 ) or we continue the reduction. The process must terminate when the face under consideration is of dimension 0 (the origin) or of dimension 1 (the span of a Perron eigenvector). In either case the probem is trivia and (3.49) is shown. The reverse containment foows simiary by appying (vi) and (viii) of Lemma
Problem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationThe Group Structure on a Smooth Tropical Cubic
The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,
More informationCONGRUENCES. 1. History
CONGRUENCES HAO BILLY LEE Abstract. These are notes I created for a seminar tak, foowing the papers of On the -adic Representations and Congruences for Coefficients of Moduar Forms by Swinnerton-Dyer and
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION
J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,
More informationPricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications
Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationPartial permutation decoding for MacDonald codes
Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationNOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS
NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs
More informationu(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0
Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,
More informationInvestigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l
Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA
More informationUniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete
Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity
More informationFRIEZE GROUPS IN R 2
FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the
More informationMat 1501 lecture notes, penultimate installment
Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationHomogeneity properties of subadditive functions
Annaes Mathematicae et Informaticae 32 2005 pp. 89 20. Homogeneity properties of subadditive functions Pá Burai and Árpád Száz Institute of Mathematics, University of Debrecen e-mai: buraip@math.kte.hu
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationReichenbachian Common Cause Systems
Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,
More informationResearch Article Numerical Range of Two Operators in Semi-Inner Product Spaces
Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda
More informationBASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set
BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a
More informationBALANCING REGULAR MATRIX PENCILS
BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationSTABLE GRAPHS BENJAMIN OYE
STABLE GRAPHS BENJAMIN OYE Abstract. In Reguarity Lemmas for Stabe Graphs [1] Maiaris and Sheah appy toos from mode theory to obtain stronger forms of Ramsey's theorem and Szemeredi's reguarity emma for
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationLECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS
More information#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG
#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG Guixin Deng Schoo of Mathematica Sciences, Guangxi Teachers Education University, Nanning, P.R.China dengguixin@ive.com Pingzhi Yuan
More informationK a,k minors in graphs of bounded tree-width *
K a,k minors in graphs of bounded tree-width * Thomas Böhme Institut für Mathematik Technische Universität Imenau Imenau, Germany E-mai: tboehme@theoinf.tu-imenau.de and John Maharry Department of Mathematics
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationPrimal and dual active-set methods for convex quadratic programming
Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:
More informationTheory of Generalized k-difference Operator and Its Application in Number Theory
Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationThe arc is the only chainable continuum admitting a mean
The arc is the ony chainabe continuum admitting a mean Aejandro Ianes and Hugo Vianueva September 4, 26 Abstract Let X be a metric continuum. A mean on X is a continuous function : X X! X such that for
More informationK p q k(x) K n(x) x X p
oc 5. Lecture 5 5.1. Quien s ocaization theorem and Boch s formua. Our next topic is a sketch of Quien s proof of Boch s formua, which is aso a a brief discussion of aspects of Quien s remarkabe paper
More informationThe Construction of a Pfaff System with Arbitrary Piecewise Continuous Characteristic Power-Law Functions
Differentia Equations, Vo. 41, No. 2, 2005, pp. 184 194. Transated from Differentsia nye Uravneniya, Vo. 41, No. 2, 2005, pp. 177 185. Origina Russian Text Copyright c 2005 by Izobov, Krupchik. ORDINARY
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationThe Partition Function and Ramanujan Congruences
The Partition Function and Ramanujan Congruences Eric Bucher Apri 7, 010 Chapter 1 Introduction The partition function, p(n), for a positive integer n is the number of non-increasing sequences of positive
More informationA NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC
(January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationPowers of Ideals: Primary Decompositions, Artin-Rees Lemma and Regularity
Powers of Ideas: Primary Decompositions, Artin-Rees Lemma and Reguarity Irena Swanson Department of Mathematica Sciences, New Mexico State University, Las Cruces, NM 88003-8001 (e-mai: iswanson@nmsu.edu)
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationPattern Frequency Sequences and Internal Zeros
Advances in Appied Mathematics 28, 395 420 (2002 doi:10.1006/aama.2001.0789, avaiabe onine at http://www.ideaibrary.com on Pattern Frequency Sequences and Interna Zeros Mikós Bóna Department of Mathematics,
More informationMIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI
MIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI KLAUS SCHMIDT AND TOM WARD Abstract. We prove that every mixing Z d -action by automorphisms of a compact, connected, abeian group is
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationJENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF SEVERAL VARIABLES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Voume 128, Number 7, Pages 2075 2084 S 0002-99390005371-5 Artice eectronicay pubished on February 16, 2000 JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More informationLecture 11. Fourier transform
Lecture. Fourier transform Definition and main resuts Let f L 2 (R). The Fourier transform of a function f is a function f(α) = f(x)t iαx dx () The normaized Fourier transform of f is a function R ˆf =
More informationMaejo International Journal of Science and Technology
Fu Paper Maejo Internationa Journa of Science and Technoogy ISSN 1905-7873 Avaiabe onine at www.mijst.mju.ac.th A study on Lucas difference sequence spaces (, ) (, ) and Murat Karakas * and Ayse Metin
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationUNIFORM CONVERGENCE OF MULTIPLIER CONVERGENT SERIES
royecciones Vo. 26, N o 1, pp. 27-35, May 2007. Universidad Catóica de Norte Antofagasta - Chie UNIFORM CONVERGENCE OF MULTILIER CONVERGENT SERIES CHARLES SWARTZ NEW MEXICO STATE UNIVERSITY Received :
More informationYET ANOTHER PROPERTY OF THE SORGENFREY PLANE
Voume 6, 1981 Pages 31 43 http://topoogy.auburn.edu/tp/ YET ANOTHER PROPERTY OF THE SORGENFREY PLANE by Peter de Caux Topoogy Proceedings Web: http://topoogy.auburn.edu/tp/ Mai: Topoogy Proceedings Department
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More informationHomework 5 Solutions
Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationKing Fahd University of Petroleum & Minerals
King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161
More informationOn nil-mccoy rings relative to a monoid
PURE MATHEMATICS RESEARCH ARTICLE On ni-mccoy rings reative to a monoid Vahid Aghapouramin 1 * and Mohammad Javad Nikmehr 2 Received: 24 October 2017 Accepted: 29 December 2017 First Pubished: 25 January
More informationOn Non-Optimally Expanding Sets in Grassmann Graphs
ectronic Cooquium on Computationa Compexity, Report No. 94 (07) On Non-Optimay xpanding Sets in Grassmann Graphs Irit Dinur Subhash Khot Guy Kinder Dor Minzer Mui Safra Abstract The paper investigates
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized
More informationarxiv: v3 [math.ac] 7 Dec 2009
SUPERTROPICAL MATRIX ALGEBRA II: SOLVING TROPICAL EQUATIONS ZUR IZHAKIAN AND LOUIS ROWEN arxiv:0902.2159v3 [math.ac] 7 Dec 2009 Abstract. We continue the study of matrices over a supertropica agebra, proving
More information(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].
PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform
More informationPREPUBLICACIONES DEL DEPARTAMENTO DE ÁLGEBRA DE LA UNIVERSIDAD DE SEVILLA
EUBLICACIONES DEL DEATAMENTO DE ÁLGEBA DE LA UNIVESIDAD DE SEVILLA Impicit ideas of a vauation centered in a oca domain F. J. Herrera Govantes, M. A. Oaa Acosta, M. Spivakovsky, B. Teissier repubicación
More informationGeometric Mapping Properties of Semipositive Matrices
Geometric Mapping Properties of Semipositive Matrices M. J. Tsatsomeros Mathematics Department Washington State University Pullman, WA 99164 (tsat@wsu.edu) July 14, 2015 Abstract Semipositive matrices
More informationThe ordered set of principal congruences of a countable lattice
The ordered set of principa congruences of a countabe attice Gábor Czédi To the memory of András P. Huhn Abstract. For a attice L, et Princ(L) denote the ordered set of principa congruences of L. In a
More informationReflection principles and kernels in R n _+ for the biharmonic and Stokes operators. Solutions in a large class of weighted Sobolev spaces
Refection principes and kernes in R n _+ for the biharmonic and Stokes operators. Soutions in a arge cass of weighted Soboev spaces Chérif Amrouche, Yves Raudin To cite this version: Chérif Amrouche, Yves
More informationCnoidal wave solutions to Boussinesq systems
IOP PUBLISHING Noninearity 0 007) 443 46 NONLINEARITY doi:0.088/095-775/0/6/007 Cnoida wave soutions to Boussinesq systems Hongqiu Chen, Min Chen and Nghiem V Nguyen Department of Mathematica Sciences,
More informationWeak lumpability of finite Markov chains and positive invariance of cones
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Weak umpabiity of finite Markov chains and positive invariance of cones James Ledoux N 2801 Fevrier 1996 PROGRAMME 1 apport de recherche
More informationOn Bounded Nondeterminism and Alternation
On Bounded Nondeterminism and Aternation Mathias Hauptmann May 4, 2016 Abstract We continue our work on the combination of variants of McCreight and Meyer s Union Theorem with separation resuts aong the
More informationON CERTAIN SUMS INVOLVING THE LEGENDRE SYMBOL. Borislav Karaivanov Sigma Space Inc., Lanham, Maryland
#A14 INTEGERS 16 (2016) ON CERTAIN SUMS INVOLVING THE LEGENDRE SYMBOL Borisav Karaivanov Sigma Sace Inc., Lanham, Maryand borisav.karaivanov@sigmasace.com Tzvetain S. Vassiev Deartment of Comuter Science
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationFinite element method for structural dynamic and stability analyses
Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department
More informationFFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary
More informationScott Cohen. November 10, Abstract. The method of Block Cyclic Reduction (BCR) is described in the context of
ycic Reduction Scott ohen November, 99 bstract The method of ock ycic Reduction (R) is described in the context of soving Poisson's equation with Dirichet boundary conditions The numerica instabiityof
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationChapter 5. Wave equation. 5.1 Physical derivation
Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing
More informationAlgorithms to solve massively under-defined systems of multivariate quadratic equations
Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations
More informationSmall generators of function fields
Journa de Théorie des Nombres de Bordeaux 00 (XXXX), 000 000 Sma generators of function fieds par Martin Widmer Résumé. Soit K/k une extension finie d un corps goba, donc K contient un éément primitif
More informationThe Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation
The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric
More informationResearch Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation
Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi
More informationCompletion. is dense in H. If V is complete, then U(V) = H.
Competion Theorem 1 (Competion) If ( V V ) is any inner product space then there exists a Hibert space ( H H ) and a map U : V H such that (i) U is 1 1 (ii) U is inear (iii) UxUy H xy V for a xy V (iv)
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationCourse 2BA1, Section 11: Periodic Functions and Fourier Series
Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........
More informationOn prime divisors of remarkable sequences
Annaes Mathematicae et Informaticae 33 (2006 pp. 45 56 http://www.ektf.hu/tanszek/matematika/ami On prime divisors of remarkabe sequences Ferdinánd Fiip a, Kámán Liptai b1, János T. Tóth c2 a Department
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
SIAM J. NUMER. ANAL. Vo. 0, No. 0, pp. 000 000 c 200X Society for Industria and Appied Mathematics VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW
More informationarxiv: v1 [math.co] 17 Dec 2018
On the Extrema Maximum Agreement Subtree Probem arxiv:1812.06951v1 [math.o] 17 Dec 2018 Aexey Markin Department of omputer Science, Iowa State University, USA amarkin@iastate.edu Abstract Given two phyogenetic
More informationNumerical methods for elliptic partial differential equations Arnold Reusken
Numerica methods for eiptic partia differentia equations Arnod Reusken Preface This is a book on the numerica approximation of partia differentia equations. On the next page we give an overview of the
More informationODE Homework 2. Since M y N x, the equation is not exact. 2. Determine whether the following equation is exact. If it is exact, M y N x 1 x.
ODE Homework.6. Exact Equations and Integrating Factors 1. Determine whether the foowing equation is exact. If it is exact, find the soution pe x sin qdx p3x e x sin qd 0 [.6 #8] So. Let Mpx, q e x sin,
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More informationCo-degree density of hypergraphs
Co-degree density of hypergraphs Dhruv Mubayi Department of Mathematics, Statistics, and Computer Science University of Iinois at Chicago Chicago, IL 60607 Yi Zhao Department of Mathematics and Statistics
More informationb n n=1 a n cos nx (3) n=1
Fourier Anaysis The Fourier series First some terminoogy: a function f(x) is periodic if f(x ) = f(x) for a x for some, if is the smaest such number, it is caed the period of f(x). It is even if f( x)
More informationBourgain s Theorem. Computational and Metric Geometry. Instructor: Yury Makarychev. d(s 1, s 2 ).
Bourgain s Theorem Computationa and Metric Geometry Instructor: Yury Makarychev 1 Notation Given a metric space (X, d) and S X, the distance from x X to S equas d(x, S) = inf d(x, s). s S The distance
More information