The Computable Differential Equation Lecture ... - Bruce E. Shapiro
The Computable Differential Equation Lecture ... - Bruce E. Shapiro The Computable Differential Equation Lecture ... - Bruce E. Shapiro
170 CHAPTER 9. DIFFERENTIAL ALGEBRAIC EQUATIONS Definition 9.2. The function φ(t) is a solution of 9.1 on an interval I if F(t, φ(t), φ ′ (t)) = 0 (9.43) for all t ∈ I. We call the DAE solvable on I if a solution (or family of solutions that depend on a parameter) exists on I. While this construction described earlier suggests that one can solve a DAE by converting it to an ODE, this is not always the best method, since differentiation can be expensive and knowing how to differentiate to get the desired result is a difficult problem. However, the concepts of solvability and the definition of the Index of the DAE are closely related. 9.2 Linear Differential Algebraic Equations with Constant Coeffcients Our focus will be primarily on linear equations of the form A(t)y ′ + B(t)y = f(t) (9.44) Definition 9.3. A linear system is said to be in standard canonical form if it can be partitioned into the form ( ) ( ) ( ( ) I 0 x ′ C 0 x 0 N y ′ + = 0 I) y ( ) f(t) g(t) (9.45) where N is a strictly 3 triangular matrix. For linear time-varying systems, both N and C depend on t. Expanding, x ′ + Cx = f(t) (9.46) Ny ′ + y = g(t) (9.47) We begin with the linear constant coefficient DAE (9.37). Definition 9.4 (Matrix Pencil). Let A, B be square matrices and λ ∈ C. Then the matrix pencil of A and B is given by P is said to be a regular pencil if det λA + B ≠ 0. P = λA + B (9.48) Definition 9.5 (Nilpotent Matrix). A matrix N is said have a Nilpotency k if k is the smallest integer such that N k = 0 (9.49) 3 A strictly triangular matrix is a triangular matrix with all zeros on the diagonal. Math 582B, Spring 2007 California State University Northridge c○2007, B.E.Shapiro Last revised: May 23, 2007
CHAPTER 9. DIFFERENTIAL ALGEBRAIC EQUATIONS 171 Theorem 9.1. Let Ay ′ + By = f be a linear constant coefficient DAE with regular pencil λA + B. Then there exists nonsingular matrices P, Q, such that where u = Q −1 y, P AQ = P AQu ′ + P BQu = P f (9.50) ( ) I 0 ; P BQ = 0 N ( ) C 0 0 I (9.51) and N is a matrix of nilpotency k, for some integer k ≥ 0. If N = 0 define k = 1; if A is nonsignular, P AQ = I and k = 0. If det(λA + B) is constant then P AQ = n, P BQ = I. Definition 9.6. The index of the pencil λA + B is the degree of nilpotency in equation 9.51. Theorem 9.2. Let Ay ′ + By = f be a linear constant coefficient DAE with regular pencil λA + B. Then the index of the pencil equals the index of the DAE. This theorem says that there exist nonsingular matrices P , Q, such that we can write ( ) ( ) ( ( ) ( I 0 u ′ C 0 u g 0 N v ′ + = (9.52) 0 I) v h) where Qy = ( ( u g , P f = v) h) (9.53) Then u ′ + Cu = g (9.54) Nv ′ + v = h (9.55) and N is k-nilpotent, where k is the index of the DAE. Equation 9.54 has a unique solution that is fully determined by its initial conditions (this is the existence theorem for ordinary differential equations). Equation 9.55 also has a unique solution, but does not require initial conditions, because the matrix N is nilpotent (N k =0). To see this let D be the derivative operator; then 9.55 becomes But since NDv + v = h (9.56) (ND + I)v = h (9.57) v = (ND + I) −1 h (9.58) (I + ND) −1 = I − ND + N 2 D 2 − N 3 D 3 + N 4 D 4 + · · · + N k−1 D k−1 (9.59) c○2007, B.E.Shapiro Last revised: May 23, 2007 Math 582B, Spring 2007 California State University Northridge
- Page 125 and 126: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 127 and 128: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 129 and 130: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 131 and 132: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 133 and 134: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 135 and 136: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 137 and 138: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 139 and 140: CHAPTER 6. LINEAR MULTISTEP METHODS
- Page 141 and 142: Chapter 7 Delay Differential Equati
- Page 143 and 144: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 145 and 146: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 147 and 148: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 149 and 150: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 151 and 152: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 153 and 154: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 155 and 156: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 157 and 158: CHAPTER 7. DELAY DIFFERENTIAL EQUAT
- Page 159 and 160: Chapter 8 Boundary Value Problems 8
- Page 161 and 162: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 163 and 164: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 165 and 166: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 167 and 168: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 169 and 170: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 171 and 172: CHAPTER 8. BOUNDARY VALUE PROBLEMS
- Page 173 and 174: Chapter 9 Differential Algebraic Eq
- Page 175: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 179 and 180: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 181 and 182: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 183 and 184: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 185 and 186: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 187 and 188: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 189 and 190: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 191 and 192: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 193 and 194: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 195 and 196: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 197 and 198: CHAPTER 9. DIFFERENTIAL ALGEBRAIC E
- Page 199 and 200: Chapter 10 Appendix on Analytic Met
- Page 201 and 202: CHAPTER 10. APPENDIX ON ANALYTIC ME
- Page 203 and 204: CHAPTER 10. APPENDIX ON ANALYTIC ME
- Page 205 and 206: Bibliography [1] Ascher, Uri M., Ma
170 CHAPTER 9. DIFFERENTIAL ALGEBRAIC EQUATIONS<br />
Definition 9.2. <strong>The</strong> function φ(t) is a solution of 9.1 on an interval I if<br />
F(t, φ(t), φ ′ (t)) = 0 (9.43)<br />
for all t ∈ I. We call the DAE solvable on I if a solution (or family of solutions<br />
that depend on a parameter) exists on I.<br />
While this construction described earlier suggests that one can solve a DAE by<br />
converting it to an ODE, this is not always the best method, since differentiation can<br />
be expensive and knowing how to differentiate to get the desired result is a difficult<br />
problem. However, the concepts of solvability and the definition of the Index of the<br />
DAE are closely related.<br />
9.2 Linear <strong>Differential</strong> Algebraic <strong>Equation</strong>s with Constant<br />
Coeffcients<br />
Our focus will be primarily on linear equations of the form<br />
A(t)y ′ + B(t)y = f(t) (9.44)<br />
Definition 9.3. A linear system is said to be in standard canonical form if it<br />
can be partitioned into the form<br />
( ) ( ) ( ( )<br />
I 0 x<br />
′ C 0 x<br />
0 N y ′ +<br />
=<br />
0 I)<br />
y<br />
( ) f(t)<br />
g(t)<br />
(9.45)<br />
where N is a strictly 3 triangular matrix. For linear time-varying systems, both N<br />
and C depend on t. Expanding,<br />
x ′ + Cx = f(t) (9.46)<br />
Ny ′ + y = g(t) (9.47)<br />
We begin with the linear constant coefficient DAE (9.37).<br />
Definition 9.4 (Matrix Pencil). Let A, B be square matrices and λ ∈ C. <strong>The</strong>n<br />
the matrix pencil of A and B is given by<br />
P is said to be a regular pencil if det λA + B ≠ 0.<br />
P = λA + B (9.48)<br />
Definition 9.5 (Nilpotent Matrix). A matrix N is said have a Nilpotency k<br />
if k is the smallest integer such that<br />
N k = 0 (9.49)<br />
3 A strictly triangular matrix is a triangular matrix with all zeros on the diagonal.<br />
Math 582B, Spring 2007<br />
California State University Northridge<br />
c○2007, B.E.<strong>Shapiro</strong><br />
Last revised: May 23, 2007