Singular Event Detection
|
|
- Dortha Stokes
- 5 years ago
- Views:
Transcription
1 Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate Mentor: Aaron Ames Summer Undergraduate Program in Engineering at Berkeley (SUPERB) 2004 Department of Eletial Engineering & Computer Sienes College of Engineering University of California, Berkeley
2 Singular Event Detetion Rafael S. Garía Abstrat The main objetive in the study of event loation is to find when the solution to an Ordinary Differential Equation (ODE) rosses a swithing surfae. Due to the fat that it is almost always impossible to find the exat solution to an arbitrary ODE, the event must be loated with a numerial approximation to the exat solution. Sine approximate solutions are by nature inexat, the question is whether the real solution displays the same qualitative behavior as its numerial approximation. This raises the need for a number that gauges the validity of this approximation, i.e., a ondition number. The purpose of this researh is to show how numerial approximations of the solutions of simple systems linear systems with simple swithing surfaes hyperplanes an lead to inorret event detetion. Motivated by these examples, a andidate ondition number will be proposed in an attempt to quantify this failure. I. INTRODUCTION Consider an Initial Value Problem (IVP), for a set of ODEs suh that: { f ẋ = + (x) if g(x) 0 f (x) if g(x) < 0 where tɛ[a, b], x(a) = x 0, and x R 2. An event is said to our at time t if g(ˆx + (t )) = 0; the point at whih the vetor field hanges from f + to f. Due to the fat that the majority of the Ordinary Differential Equations (ODEs) do not have an exat solution, it is imperative to utilize numerial methods to approximate the solution of the IVP stated above. As a result, a level of unertainty is introdued in the event loation, whih leads to the need for a guarantee of its reliability. Thus reating a neessity of a ondition number that ould be apable of quantifying the level of orretness of the event loation, i.e, the orretness the approximated solution with respet to the atual solution. In vetor algebra, the ondition number, κ, of a square matrix, A, is defined as: κ(a) = A A 1 where, is a valid vetor norm. Alternatively, if the square matrix A is singular deomposed as A = UΣV T, the ondition number of a matrix an also be stated as: κ(a) = σ max σ min where, σ max and σ min are, respetively, the maximum and minimum salar values of the diagonal matrix Σ. The ondition number of a matrix provides a measure of how sensitive is a matrix to numerial operations and errors in the data. Additionally, it indiates the auray of matrix inversion and the auray of the solution of the lineal equation. If the ondition number of a matrix is approximately 1, the matrix is said to be well-onditioned. In ontrast, if the ondition number of the matrix is muh greater than one, it is said to be ill-onditioned. Ill-onditioned matries do not yield preise numerial solutions [2].This study seeks to propose a ondition number for event loation. This number will serve as a mean to numerially judge its reliability. In fat, the proposed ondition number is: ond#(f) = f + (ˆx(t )) f (ˆx(t )) +g(ˆx(t )) g(ˆx(t )) t p P sign(g(ˆx(t p ))) g(ˆx(t p )) h ˆx(t p ) ˆx(t ) t p P
3 For a given subset, P, suh that: P = {t p R n d : +g(ˆx(t p )) = 0 dt +g(ˆx(t p)) > 0} Further details of this proposed number will be desribed later on. II. NUMERICAL ODE SOLVERS Already knowing that is impossible to find an exat solution to all of the existing ODEs, it is imperative to find numerial ways to approximate the solutions to these equations. In this study, it was desirable to find a numerial solution to an IVP in a given interval a t b. The studied ODEs were of the form: ẋ = ẋ 1 ẋ 2. ẋ n = f(t, x) with given initial values x(a). These equations are of speial interest, for they are utilized to desribe the dynamis and behavior of swithing systems. From [7], it is known that for a given tolerane, τ, the numerial solver an approximate a solution, ˆx, suh that is aurate to the order of r(τ). Furthermore, it an be assumed that the solution is bounded by x(t) ˆx(t) r(h), where h is step size of the solver and r(h) 0 as h 0. Following is a desription of several ODEs that were used to solve the ODEs of the event detetion algorithms performed in setion 3 and 5. All of the methods utilized were either one step methods or multi step methods. The one step methods inluded: Forward Euler s, 2nd order Taylor s, and 4th order Runge-Kotta. Conversely, the multistep methods inluded 3rd order Preditor Correted Adams-Moulton (AM3) and Bakward Differene Formula (BDF). The equation ẋ = f(t, x) was onsidered 2 dimensional, and solved with these methods. The final system of equations is omposed of equations (1) and (2) as follows: A. Forward Euler s Method ẋ 1 = F(t, x 1, x 2 ) (1) ẋ 2 = G(t, x 1, x 2 ) (2) The forward Euler method is attrative beause it is easy to implement and visualize. Although, it has some drawbaks as is the ase of major error results with inreasing step size numbers. This method is based on the one-term Taylor series. From [6] (pg. 142), it an be found that the numerial solution of a system of two variables an be approximated using: x 1 (j + 1) = x 1 (j) + hf(t(j), x 1 (j), x 2 (j)) x 2 (j + 1) = x 2 (j) + hg(t(j), x 1 (j), x 2 (j)) where h is the time interval, h = t j+1 t j, and j=0,1,2,...n. N represents the number of steps in the interval studied. In this method the numerial solution is found by reusing the previously omputed value to approximate the subsequent value of the solution. This proess ontinues iteratively until the desired interval is overed.
4 B. Seond-Order Taylor Method The seond-order Taylor method offers a better loal error than Euler s method; although, it has also showed very unpleasing results with inreasing step size numbers. In [6](pg. 143), it is shown that the numerial solution to an ODE using a seond-order Taylor expansion is given by: x 1 (j + 1) = x 1 (j) + hf (t(j), x 1 (j), x 2 (j)) + h2 2 [ F t (t(j), x 1(j), x 2 (j)) + F t (t(j), x 1(j), x 2 (j))f (t(j), x 1 (j), x 2 (j))] x 2 (j + 1) = x 2 (j) + hg(t(j), x 1 (j), x 2 (j)) + h2 2 [ G t (t(j), x 1(j), x 2 (j)) + G t (t(j), x 1(j), x 2 (j))g(t(j), x 1 (j), x 2 (j))] where h is the time interval, h = t j+1 t j. As in the Forward Euler s method, the numerial solution is found by reusing the previously omputed value to approximate the subsequent value of the solution. C. Runge-Kotta Method The Runge-Kotta method offers another alternative to approximate the solution of an ODE. This method is widely used due to its higher auray in omparison to other one-step methods. Additionally, it offers a better loal error than Euler s and Taylor s, showing good approximations with inreasing step size numbers. From [1], it is found a 4th order Runge-Kotta numerial solution for the system of equations (1) and (2) as follows: where, ˆx 1 (j + 1) = ˆx 1 (j) + h 6 [k 1 + 2k 2 + 2k 3 + k 4 ] ˆx 2 (j + 1) = ˆx 2 (j) + h 6 [m 1 + 2m 2 + 2m 3 + m 4 ] k 1 = F (t(j), ˆx 1 (j), ˆx 2 (j)) m 1 = G(t(j), ˆx 1 (j), ˆx 2 (j)) k 2 = F (t(j) + h 2, ˆx 1(j) + h 2 k 1, ˆx 2 (j) + h 2 m 1) m 2 = G(t(j) + h 2, ˆx 1(j) + h 2 k 1, ˆx 2 (j) + h 2 m 1) k 3 = F (t(j) + h 2, ˆx 1(j) + h 2 k 2, ˆx 2 (j) + h 2 m 2) m 3 = G(t(j) + h 2, ˆx 1(j) + h 2 k 2, ˆx 2 (j) + h 2 m 2) k 4 = F (t(j) + h, ˆx 1 (j) + hk 3, ˆx 2 (j) + hm 3 ) m 4 = G(t(j) + h, ˆx 1 (j) + hk 3, ˆx 2 (j) + hm 3 ) In this method, h is one again the time interval, h = t(j + 1) t(j). The numerial solution is found by reusing the previously omputed value to approximate the subsequent point of the solution. D. Adam-Moulton Method As opposite to the one-step methods developed by Euler, Taylor, and Runge-Kutta; the Adam-Moulton method has the attrative of utilizing previously alulated points to find the subsequent point in the iteration. Methods that follow this desription are often referred as memory methods. Although, it is not self starting and it has to be given previously alulated values to start its omputations depending on the order of its implementation. In order to ahieve a better auray, it is good to utilize a good self-starting one-step method. For this implementation, Runge-Kutta was hosen to perform this task for the three first points of the solution.
5 From [3], it is found that one those 3 initial points are reahed the algorithm used for the 3rd order Preditor Correted Adams-Moulton (AM3) is the following: x 1 (j + 1) = ˆx 1 (j) + h 12 [23F (t(j), ˆx 1(j), ˆx 2 (j)) 16F (t(j 1), ˆx 1 (j 1), ˆx 2 (j 1)) +5F (t(j 2), ˆx 1 (j 2), ˆx 2 (j 2))] x 2 (j + 1) = ˆx 2 (j) + h 12 [23G(t(j), ˆx 1(j), ˆx 2 (j)) 16G(t(j 1), ˆx 1 (j 1), ˆx 2 (j 1)) +5G(t(j 2), ˆx 1 (j 2), ˆx 2 (j 2))] ˆx 1 (j + 1) = h 12 [5F (t(j + 1), x 1(j + 1), x 2 (j + 1)) + 8F (t(j), ˆx 1 (j), ˆx 2 (j)) F (t(j 1), ˆx 1 (j 1), ˆx 2 (j 1))] ˆx 2 (j + 1) = h 12 [5G(t(j + 1), x 1(j + 1), x 2 (j + 1)) + 8G(t(j), ˆx 1 (j), ˆx 2 (j)) F (t(j 1), ˆx 1 (j 1), ˆx 2 (j 1))] As in the previous methods, h represents the time interval, h = t j+1 t j. In this equation the preditor is used to aquire a better approximation. Finally, the solution is found by reusing the previously omputed value to alulate the subsequent point of the solution. E. Bakward Differene Formula Method The BDF method was another memory method implemented beause of its attrativeness of utilizing previously alulated points to find the subsequent point in the iteration. As in the AM3 method, the memory harateristi yields more aurate results than one-step methods suh as similar to the AM3 method, the BDF method is not self starting, it has to be given previously alulated values to start its omputations. For auray issues, a fourth order Runge-Kutta was used to alulate the three first points of the solution. From [6] (pag. 183), it is found that one those 3 initial points are reahed, the algorithm to solve the system is the following: ˆx 1 (j + 1) = 6 11 [6F (t(j + 1), ˆx 1(j + 1), ˆx 2 (j + 1)) + 3ˆx 1 (j) 3 2 ˆx 1(j 1) ˆx 1(j 2)] ˆx 2 (j + 1) = 6 11 [6G(t(j + 1), ˆx 1(j + 1), ˆx 2 (j + 1)) + 3ˆx 2 (j) 3 2 ˆx 2(j 1) ˆx 2(j 2)] As in the previous methods, h also represents the time interval, h = t j+1 t j. In this equation the preditor is used to aquire a better approximation. Finally, the solution is found by reusing the previously omputed value to alulate the subsequent point of the solution. III. EVENT DETECTION TECHNIQUES The next task is to find and apply event detetion tehniques (EDTs), to the already stated IVP given by: { f ẋ = f(x) = + (x) if g(x) 0 f (3) (x) if g(x) < 0 for the given interval tɛ[a, b], and initial ondition x = x 0, in whih x R 2. In order to simplify further explanations, an event detetion as (3) will be referred as f = (f +, f, g) EDTs utilize different algorithms with the objetive of finding the most aurate way to loate events. Two different algorithms were simulated utilizing MATLAB. The first method will be referred as the zero-loation method and is similar to that of [4], while the seond method is referred as the minimum-neighborhood method, whih is similar to the interval bisetion algorithm found in [5], but with some slight hanges.
6 A. Zero-Loation Method The zero-loation method (ZLM) starts by numerially solving the ODE, ẋ = f + (x). The approximated solution is denoted by ˆx + (t) with initial ondition ˆx + (0) = x 0. Subsequently, at eah time step, t n of the interval, it is heked if g(x + (t n )) > 0 and g(x + (t n+1 )) < 0. If this ondition holds, an event has ourred. Consequently, a linear interpolation between the points ˆx + (t n ) and ˆx + (t n+1 ) is performed, and denoted by ˆx + (t). Then, it is numerially solved for t suh that g(ˆx + (t )) = 0 After the ourrene of an event, the same proess is started again, exept that ẋ = f (x) is solved with ˆx (0) = ˆx + (t ). Finally, this algorithm is repeated (with its respetive sign hanges) if multiple events are to be deteted. B. Minimum-Neighborhood Method The minimum-neighborhood method (MNM)follows the first ouple of steps as the zero-loation method. It starts by numerially solving ẋ = f + (x), where the approximated solution is denoted by ˆx + (t) with initial ondition ˆx + (0) = x 0. Subsequently, at eah time step t n of the interval, it is heked if g(x + (t n )) > 0 and g(x + (t n+1 )) < 0. If this ondition holds, an event has ourred and the next step is to subdivide the interval [t n, t n+1 ]. This is performed by seleting times τ 0 = t n < τ 1 < < τ k = t n+1 where k is input by the user. Subsequently, values for i are found suh that g(ˆx(τ i )) > 0 and g(ˆx(τ i+1 )) < 0 holds. One this ondition is true, the interval is deremented again for the points in whih it holds, until values are found suh that g(ˆx(τ i )) and g(ˆx(τ i+1 )) are approximately zero. Consequently, t is set to τ i. After loating the event, the same proess is started again, exept that ẋ = f (x) with initial ondition ˆx (0) = ˆx + (t ). C. Example The following IVP was implemented to illustrate the solutions of both algorithms: ( ) ( ) 1 1 x1 ẋ = 1 1 with an h = 0.01 and x + (a) = ( 2, 2) T. By looking at the graphs, it is notieable that both methods follow x vs. X 1 Euler s INT Runge Kotta INT Euler s MNM Runge Kotta MNM Event at X 1 = X 1 Fig. 1. Event detetion for both methods a very similar behavior. Although, they are slightly different in terms of omputational issues. For instane, the MNM offers greater preision but at the expense of more omputational time. In our ase, the omputational
7 9.148 vs. X 1 Euler s INT Runge Kotta INT Euler s MNM Runge Kotta MNM Event at X 1 = X 1 Fig. 2. Euler s proximity for both methods vs. X 1 Euler s INT Runge Kotta INT Euler s MNM Runge Kotta MNM Event at X 1 = X 1 Fig. 3. Runge-Kotta proximity for both methods time was an issue for the alulation of the ondition number. As a result,the ZLM was used to performed the event detetion tasks. IV. SOLUTION BEHAVIOR Up to this point, the importane of the existene of a ondition number has been pointed out but not properly defined in terms of good or bad behavior of an IVP as given in (3). In order to do this, it is important to define what onstitutes a good senario and what does not. A good senario would be one in whih { r x(t) ˆx(t) + (h) if t [a, t ] r (h) if t [t, b] when x + (t) ˆx + (t) r + (h) and x (t) ˆx (t) r (h) where ˆx ± (t) are the approximate solutions of x ±. This definition of a good senario an be simplified by taking r(h) = max{r + (h), r (h)} yielding the following definition:
8 Definition 4.1: The solution of the IVP f = (f +, f, g) is said to be well-posed if x + (t) ˆx + (t) r(h) and x (t) ˆx (t) r(h) for some r(h) implies that x(t) ˆx(t) r(h) for all t [a, b]. This definition of a well-posed event an be used to give the definition of an ill-posed event Definition 4.2: The solution of the IVP f = (f +, f, g) is said to be ill-posed if x + (t) ˆx + (t) r(h) and x (t) ˆx (t) r(h) for some r(h), but x(t) ˆx(t) > r(h) for some t [a, b]. Alternatively, an f is ill-posed if and r(h) > r(h) for all h. x(t) ˆx(t) r(h) An ill-posed solution reates a very unstable behavior in whih the aurate detetion of an event annot be guaranteed. The above definitions involve knowing the atual solution in order to detet wether an IVP f is ill-posed. Sine it annot be assumed that the atual solution is known, the auses some diffiulty. To irumvent this, we introdution the following lemma: Lemma 4.1: If two different approximate solutions ˆx(t) and, ŷ(t) of an IVP f = (f +, f, g), are ill-posed, i.e., if x ± (t) ˆx ± (t) r(h) and x ± (t) ŷ ± (t) r(h) but, ˆx(t) ŷ(t) r(h) then the atual solution of the IVP is ill-posed: x(t) ˆx(t) r(h) and x(t) ŷ(t) r(h). V. CONDITION NUMBER Following is a detail explanation of the set of speifiations utilized to propose a ondition number for event detetion tehniques. Consequently, the proposed onditioned number will be expressed and thoroughly explained. Finally, two examples are going to be presented with the objetive of providing an idea on the behavior of the ondition number in event detetion appliations.
9 A. Speifiations The following are a set of speifiations reated with the objetive of narrowing the behavior of any prospetive andidate for the ondition number: 1) For the perfet ase, i.e., a ase in whih the solution hits the surfae orthogonally, the ondition number should be one. 2) As the ondition number inreases, the auray of the approximate solution dereases; beoming progressively ill-posed. 3) The ondition number is defined to be infinite for a singular event, i.e., when the vetor field is tangent to the surfae at the event. In order to understand these speifiations, it is important to know their impliation. First of all, these onditions imply that the ondition number will depend on a funtion that is apable of desribing the way in whih the approximate solution intersets a swithing surfae. As a result, the Lie Derivative of g with respet to f, referred as g, and given by g(x) = g(x) x f(x) was targeted to perform this task. This funtion relates the hange of the funtion g(x) in the diretion of f(x). Therefore, in the ase in whih g = 0 at an event, it means that the swithing surfae G 0 = {x g(x) = 0} is tangent to the vetor field at this point, i.e., a singular event loation. Conversely if g = 1, it means that the vetor field is orthogonal to the surfae G 0 and the approximate solution is ompletely reliable, i.e, well-posed. To state ondition (2) of the speifiations in a more preise way, onsider two IVP f 1 = (f1, f 1 +, g 1) and f 2 = (f2, f 2 +, g 2). Beause these are ill-posed, if then x ± 1 (t) ˆx± 1 (t) r 1(h) and x ± 2 (t) ˆx± 2 (t) r 2(h) x 1 (t) ˆx 1 (t) r 1 (h) and x 2 (t) ˆx 2 (t) r 2 (h) where r 1 (h) > r 1 (h) and r 2 (h) > r 2 (h). Now a mesure of how ill-posed these events are an be obtained by onsidering r 1 (h) r 1 (h) and r 2 (h) r 2 (h). Therefore if ond#(f 1 ) denotes the ondition number of f 1 and ond#(f 2 ) denotes the ondition number of f 2, speifiation number (2) an be stated preisly as: if ond#(f 1 ) > ond#(f 2 ) r 1 (h) r 1 (h) > r 2 (h) r 2 (h). (4) This ondition is extremely important in that if a prospetive ondition number an be shown to satisfy this ondition (as well as speifiations number (1) and (3)), then it is a valid ondition number for the event detetion problem. B. Proposed Condition Number Consider the set P given by: P = {t p R n d : +g(ˆx(t p )) = 0 dt +g(ˆx(t p)) > 0} where +g and g are the Lie Derivatives of the event funtion at the right and left side of the event respetively and t p is the point at whih the approximate solution yields an g equal to zero, i.e., the time at whih the approximate solution is tangential to the solution. Subsequently, for a given funtion f = (f +, f, g), in whih g is a hyperplane of the form g(x) = x + α. Having said this, the proposed ondition number will have the form: ond#(f) = γ 1 (ˆx(t )) +g(ˆx(t )) g(ˆx(t )) t p P γ 2 (ˆx(t p )) g(ˆx(t p )) t p P ˆx(t p ) ˆx(t ) γ 3 (ˆx(t p )) (5) where the gamma funtions γ 1, γ 2 and γ 3 will be disussed in the remainder of this setion. This ondition number an be extended to the ase when g(x) is an arbitrary (smooth) funtion beause of a speifi
10 transformation. Namely, for the IVP f = (f, f +, g) setting T (x) = (x, z = g(x)) T we obtain an equivalent IVP f = ( f, f +, g) by setting ( ) ( f x f = ) ( ) ( (x) f z g(x) + x f = + ) (x) g(x, z) = z. z +g(x) Note that in this ase, g has the desired form. Using this transformation, for an IVP f = (f, f +, g) with an arbitrary smooth event funtion g, the ondition number an be defined by ond#(f) = ond#( f). So, from this point on, it an be assumed that g(x) = x + α. In order to disuss the gamma funtions in more detail, some speial ases will be onsidered. In the ase in whih, P =, i.e., there is not a tp that satisfies +g(ˆx(t p )) = 0, the ondition number beomes: ond#(f) = γ 1 (ˆx(t )) +g(ˆx(t )) g(ˆx(t )) Alternatively, when there are no events in the interval [a, b] of the integration, i.e., there does not exists a t suh that g(ˆx(t )) = 0, and P then the ondition number beomes: ond#(f) = γ 2 (ˆx(t p )) g(ˆx(t p )) t p P Finally, if P = and there are no events, the ondition number is not defined. In this ase the IVP is a normal differential equation with no swithing, so a ondition number is not needed. C. The gamma funtions γ 1 (x): The first gamma funtion was taken to be γ 1 (ˆx(t )) = f + (ˆx(t )) f (ˆx(t )) The objetive of this funtion is to sale the ondition number, so that it is insensitive the hanges in the magnitude of the vetor field as it rosses the swithing surfae and sensitive to hanges in the orientation of the vetor field at the swithing surfae. To justify this, onsider the ase when P = but an event ours, i.e., there is not a t p that satisfies +g(ˆx(t p )) = 0, (5) resulting in: ond#(f) = f + (ˆx(t )) f (ˆx(t )) +g(ˆx(t )) g(ˆx(t )). This is the ase for the IVP given by f = (f, f +, g) where ( ) f + (x) = f (x) =, g(x) = x 0 In this ase, without γ 1 (ˆx(t )), for two onstants > > 0, 1 + g(ˆx(t )) g(ˆx(t )) = 1 < 1 1 = + g(ˆx(t ))L f g(ˆx(t )) so the ondition number would sale with magnitude, while it should remain onstant under this saling. By adding the term γ 1 (ˆx(t )), the magnitude of the vetor field is saled down by its magnitude. As a result, the ondition number is no longer sensitive to the magnitude of the vetor field at the swithing surfae, i.e., we have ond#(f ) = f + (ˆx(t )) f (ˆx(t )) + g(ˆx(t )) g(ˆx(t )) = = 1 = = f + (ˆx(t )) f (ˆx(t )) + g(ˆx(t ))L f g(ˆx(t )) = ond#(f )
11 Alternatively, γ 1 (ˆx(t )) also aounts for the orientation of the vetor field at the swithing surfae. For example, onsider the IVP given by f = (f, f +, g) where ( ) f + (x) = f 1 (x) =, g(x) = x without the γ 1 (ˆx(t )) term, the ondition number would not be sensitive to hanges in orientation of vetor field at the swithing surfae beause for two onstants > > 0, 1 + g(ˆx(t )) g(ˆx(t )) = 1 = 1 + g(ˆx(t ))L f g(ˆx(t )) so as the vetor beause more and more tangent to the swithing surfae these values do not hange. Now, with the γ 1 (ˆx(t )) term, we have ond#(f ) = f + (ˆx(t )) f (ˆx(t )) + g(ˆx(t )) g(ˆx(t )) = > 1 + ( ) 2 = f + (ˆx(t )) f (ˆx(t )) + g(ˆx(t ))L f g(ˆx(t )) = ond#(f ) so as the vetor field beomes more tangent to the swithing surfae, the ondition number inreases as desired. γ 2 (x): Another ase for the ondition number is when there are no events our in the interval [a, b], but the set P is not null, i.e., t suh that g(ˆx(t )) = 0, and so P. Assuming for simpliity that P = {p} (5) beomes: ond#(f) = γ 2(ˆx(t p )) g(ˆx(t p )) In this ase we would like ond#(f) to measure how lose the solution gets to the swithing surfae G 0. With this in mind take γ 2 (ˆx(t p )) = sign(g(ˆx(t p ))). Wherein it follows that γ 2 (ˆx(t p )) g(ˆx(t p )) = 1 d(x(t p ), G 0 ) where d(x(t p ), G 0 ) is the distane between x(t p ) and a subset of points G 0 = {x g(x) = 0}. This is due to the fat that the event funtion was hosen to be of the form g(x) = x + α. To see this, note that but beause g(x) = x + α, Now, substituting into (6), as desired. d(x(t p ), G 0 ) = min y G 0 d(x(t p ), y) = min y G 0 x(t p ) y min x(t p ) y = g(x(t p )) = sign(g(x(t p )))g(x(t p )) y G 0 γ 2 (ˆx(t p )) g(ˆx(t p )) = 1 sign(g(x(t p )))g(x(t p )) = 1 d(x(t p ), G 0 ) γ 3 (x): The last ase to analyze is when there is an event in the interval [a, b], but P. In this ase the ondition number will be: ond#(f) = γ 1 (ˆx(t )) +g(ˆx(t )) g(ˆx(t )) t p P γ 2 (ˆx(t p )) g(ˆx(t p )) ˆx(t p ) ˆx(t ) γ 3 (ˆx(t p )) t p P (6)
12 This last equation ontains all of the previously indiated onstants exept from γ 3 (ˆx(t p )). In this ase, γ 3 (ˆx(t p )) has the objetive of aounting for the step size seletion in the ondition number equation. As a result, it was hosen to be equal to h. Additionally, the equation ontains the term ˆx(t p ) ˆx(t ). This term has the objetive of taking into aount how far the point where the event ours, x(t ), and the point tangential to the surfae, x(t p ). This ompensating for big strethes and ompressions of the solution near the surfae. Meaning that the further apart these points are, the more strethed and ill-posed the solution beomes. Conversely, the losest these points are, the more ompressed and well-posed should beome. Therefore, we an summarize by saying that the proposed ondition number is: D. Examples ond#(f) = f + (ˆx(t )) f (ˆx(t )) +g(ˆx(t )) g(ˆx(t )) t p P sign(g(ˆx(t p ))) g(ˆx(t p )) t p P h ˆx(t p ) ˆx(t ) (7) Following are two main examples that were simulated to alulate their respetive ondition numbers. The first example has a linear ODE funtion in f +, while the seond example has a non-linear/polynomial ODE. Reall that the ondition number was simulated with funtions of the form g(x) = bx + α. Also, 4th order Runge-Kotta was utilized to solve the funtion with a time step of ) Linear Funtion: The following IVP was implemented to find its ondition number: ( ) ( ) 1 1 x1 if g(x) x 2 ẋ = ( 0 1 ) if g(x) < 0 with g(x) = x 2 + α were α varies from α = 9 to 10 at 0.01 steps and x + (a) = ( 2, 2) T. Fig. 5, shows the 150 vs. X 1 Euler Runge Kotta x 2 = alpha X 1 Fig. 4. Zoom-in of the ODE solution plot of ondition number as alpha hanges. From this plot, it an be seen that the ondition number inreases and the solution to the IVP beomes progressively ill-posed, i.e, as α Looking at Fig. 4, it an be seen that the two approximate solutions have larger errors than what it should be expeted. This number is lose to where it is supposed to happen exatly when the approximate solution is tangent to the surfae, whih agrees with the large ondition number aquired in Fig. 5.
13 120 Condition Number vs. alpha Condition Number alpha Fig. 5. Condition Number for α = 9 : 0.01 : 10 2) Non-Linear Funtion: The following IVP was also analyzed to find its ondition number: ( ) ẋ = 1 10( 5x 4 1+β 4 ) β 5 ( 0 1 ) if g(x) 0 if g(x) < 0 with g(x) = x 2 + α were α varies from α = 4 to 8 at 0.01 steps and x + (a) = (2.5β, f + 2 (2.5β))T. vs. X Euler Runge Kotta = alpha X 1 Fig. 6. ODE solution In this example, two different things were tested. First of all, the value of alpha was varied to see variations of the ondition number throughout the interval. Seond of all, values of β in f + were varied to see if the ondition number behaved as expeted with hanges in x(t ) and x(t p ). Fig. 6 shows that the event beomes ill-pose at approximately x 2 = Also, Fig. 7, shows the plot for the ondition number as alpha varies and β remains onstant at 1. In this ase, the plot shows that the IVP beomes progressively ill-posed as α Fig.8 shows the effet that β produes on the ondition number. The graph shows that as β inreases, the value of the ondition number inreases as well. This is due to the fat that x(t ) and x(t p ) will be further apart.
14 Condition Number vs. alpha Condition Number alpha Fig. 7. Condition Number for α = 4 : 0.01 : 8 Condition Number vs. alpha Condition Number Fig. 8. Condition number for multiple β This an also be seen in Fig. 9. This figure was saled with the objetive of making all of the magnitudes equal and fousing on the behavior of the graph. As it an be seen, the slope of the bottom part of the surfae inrements as β inrements. As a result, the value of β is taken into aount as expeted. REFERENCES [1] Ordinary Differential equations.. Available Online, June [2] Matrix Condition Number. Available Online, June [3] Ordinary Differential Equations. Computational Siene Eduation Program, [4] W. H. Enrht, K. R. Jakson, S. P. Nørset and P. G. Tompsen. Effetive Soulution of disontinuous IVPs a Rounge-Kotta formula pair with interpolants. Appl. Math. Comp., vol. 27 (1998), pg [5] C. Moler. Are we there yet? Matlab News and Notes. Simulink 2 Speial Edition (1997), pg [6] L. F. Shampine. Numerial Solutions for Ordinary Differential Equations. Chapman and Hall, [7] L. F. Shampine, S. Thompson. Event Loation for Ordinary Differential Equations. Computers with Mathematis and Appliations, vol. 39 (2000), pg
15 Fig. 9. Relation between β,γ, and ondition number
Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationLecture 3 - Lorentz Transformations
Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationTwo Points Hybrid Block Method for Solving First Order Fuzzy Differential Equations
Journal of Soft Computing and Appliations 2016 No.1 (2016) 43-53 Available online at www.ispas.om/jsa Volume 2016, Issue 1, Year 2016 Artile ID jsa-00083, 11 Pages doi:10.5899/2016/jsa-00083 Researh Artile
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationAverage Rate Speed Scaling
Average Rate Speed Saling Nikhil Bansal David P. Bunde Ho-Leung Chan Kirk Pruhs May 2, 2008 Abstrat Speed saling is a power management tehnique that involves dynamially hanging the speed of a proessor.
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationThe gravitational phenomena without the curved spacetime
The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,
More informationSome facts you should know that would be convenient when evaluating a limit:
Some fats you should know that would be onvenient when evaluating a it: When evaluating a it of fration of two funtions, f(x) x a g(x) If f and g are both ontinuous inside an open interval that ontains
More informationRelativistic Dynamics
Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable
More informationA Differential Equation for Specific Catchment Area
Proeedings of Geomorphometry 2009. Zurih, Sitzerland, 3 ugust - 2 September, 2009 Differential Equation for Speifi Cathment rea J. C. Gallant, M. F. Huthinson 2 CSIRO Land and Water, GPO Box 666, Canberra
More information9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry
9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationImprovements in the Modeling of the Self-ignition of Tetrafluoroethylene
Exerpt from the Proeedings of the OMSOL onferene 010 Paris Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene M. Bekmann-Kluge 1 *,. errero 1, V. Shröder 1, A. Aikalin and J. Steinbah
More informationParticle-wave symmetry in Quantum Mechanics And Special Relativity Theory
Partile-wave symmetry in Quantum Mehanis And Speial Relativity Theory Author one: XiaoLin Li,Chongqing,China,hidebrain@hotmail.om Corresponding author: XiaoLin Li, Chongqing,China,hidebrain@hotmail.om
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationThe homopolar generator: an analytical example
The homopolar generator: an analytial example Hendrik van Hees August 7, 214 1 Introdution It is surprising that the homopolar generator, invented in one of Faraday s ingenious experiments in 1831, still
More informationExploring the feasibility of on-site earthquake early warning using close-in records of the 2007 Noto Hanto earthquake
Exploring the feasibility of on-site earthquake early warning using lose-in reords of the 2007 Noto Hanto earthquake Yih-Min Wu 1 and Hiroo Kanamori 2 1. Department of Geosienes, National Taiwan University,
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More informationarxiv: v2 [math.pr] 9 Dec 2016
Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First
More information7 Max-Flow Problems. Business Computing and Operations Research 608
7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates
More informationBreakdown of the Slowly Varying Amplitude Approximation: Generation of Backward Traveling Second Harmonic Light
Claremont Colleges Sholarship @ Claremont All HMC Faulty Publiations and Researh HMC Faulty Sholarship 1-1-003 Breakdown of the Slowly Varying Amplitude Approximation: Generation of Bakward Traveling Seond
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More informationarxiv:gr-qc/ v7 14 Dec 2003
Propagation of light in non-inertial referene frames Vesselin Petkov Siene College, Conordia University 1455 De Maisonneuve Boulevard West Montreal, Quebe, Canada H3G 1M8 vpetkov@alor.onordia.a arxiv:gr-q/9909081v7
More information11.1 Polynomial Least-Squares Curve Fit
11.1 Polynomial Least-Squares Curve Fit A. Purpose This subroutine determines a univariate polynomial that fits a given disrete set of data in the sense of minimizing the weighted sum of squares of residuals.
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More informationAcoustic Waves in a Duct
Aousti Waves in a Dut 1 One-Dimensional Waves The one-dimensional wave approximation is valid when the wavelength λ is muh larger than the diameter of the dut D, λ D. The aousti pressure disturbane p is
More informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More informationEinstein s Three Mistakes in Special Relativity Revealed. Copyright Joseph A. Rybczyk
Einstein s Three Mistakes in Speial Relativity Revealed Copyright Joseph A. Rybzyk Abstrat When the evidene supported priniples of eletromagneti propagation are properly applied, the derived theory is
More informationWord of Mass: The Relationship between Mass Media and Word-of-Mouth
Word of Mass: The Relationship between Mass Media and Word-of-Mouth Roman Chuhay Preliminary version Marh 6, 015 Abstrat This paper studies the optimal priing and advertising strategies of a firm in the
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationAfter the completion of this section the student should recall
Chapter I MTH FUNDMENTLS I. Sets, Numbers, Coordinates, Funtions ugust 30, 08 3 I. SETS, NUMERS, COORDINTES, FUNCTIONS Objetives: fter the ompletion of this setion the student should reall - the definition
More informationVariation Based Online Travel Time Prediction Using Clustered Neural Networks
Variation Based Online Travel Time Predition Using lustered Neural Networks Jie Yu, Gang-Len hang, H.W. Ho and Yue Liu Abstrat-This paper proposes a variation-based online travel time predition approah
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More informationGreen s function for the wave equation
Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More informationVolume 29, Issue 3. On the definition of nonessentiality. Udo Ebert University of Oldenburg
Volume 9, Issue 3 On the definition of nonessentiality Udo Ebert University of Oldenburg Abstrat Nonessentiality of a good is often used in welfare eonomis, ost-benefit analysis and applied work. Various
More informationDO PHYSICS ONLINE. SPECIAL RELATIVITY Frames of Reference
DO PHYSICS ONLINE SPACE SPECIAL RELATIVITY Frames of Referene Spae travel Apollo 11 spaeraft: Earth Moon v ~ 40x10 3 km.h -1 Voyager spaeraft: v ~ 60x10 3 km.h -1 (no sling shot effet) Ulysses spaeraft:
More informationNumerical simulation of a one-dimensional shock tube problem at supercritical fluid conditions
International Journal of Physial Sienes Vol. 3 (1), pp. 314-30, Deember, 008 Available online at http://www.aademijournals.org/ijps ISSN 199-1950 008 Aademi Journals Full ength esearh Paper Numerial simulation
More informationAircraft CAS Design with Input Saturation Using Dynamic Model Inversion
International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This
More informationTowards an Absolute Cosmic Distance Gauge by using Redshift Spectra from Light Fatigue.
Towards an Absolute Cosmi Distane Gauge by using Redshift Spetra from Light Fatigue. Desribed by using the Maxwell Analogy for Gravitation. T. De Mees - thierrydemees @ pandora.be Abstrat Light is an eletromagneti
More informationAssessing the Performance of a BCI: A Task-Oriented Approach
Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,
More informationTheory. Coupled Rooms
Theory of Coupled Rooms For: nternal only Report No.: R/50/TCR Prepared by:. N. taey B.., MO Otober 00 .00 Objet.. The objet of this doument is present the theory alulations to estimate the reverberant
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More informationFinal Review. A Puzzle... Special Relativity. Direction of the Force. Moving at the Speed of Light
Final Review A Puzzle... Diretion of the Fore A point harge q is loated a fixed height h above an infinite horizontal onduting plane. Another point harge q is loated a height z (with z > h) above the plane.
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationOn Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles
Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS
More informationON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS
ON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS MARIO LEFEBVRE and JEAN-LUC GUILBAULT A ontinuous-time and ontinuous-state stohasti proess, denoted by {Xt), t }, is defined from a proess known as
More information22.01 Fall 2015, Problem Set 6 (Normal Version Solutions)
.0 Fall 05, Problem Set 6 (Normal Version Solutions) Due: November, :59PM on Stellar November 4, 05 Complete all the assigned problems, and do make sure to show your intermediate work. Please upload your
More informationMATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II
Journal of Applied Mathematis and Computational Mehanis 2014, 13(2), 141-147 MATHEMATICA AND NUMERICA BAI OF BINARY AOY OIDIFICATION MODE WITH UBTITUTE THERMA CAPACITY. PART II Ewa Węgrzyn-krzypzak 1,
More informationConformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices
AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This
More informationProbabilistic Graphical Models
Probabilisti Graphial Models David Sontag New York University Leture 12, April 19, 2012 Aknowledgement: Partially based on slides by Eri Xing at CMU and Andrew MCallum at UMass Amherst David Sontag (NYU)
More informationarxiv: v1 [physics.gen-ph] 5 Jan 2018
The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the
More informationFrequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach
Frequeny Domain Analysis of Conrete Gravity Dam-Reservoir Systems by Wavenumber Approah V. Lotfi & A. Samii Department of Civil and Environmental Engineering, Amirkabir University of Tehnology, Tehran,
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationEstimating the probability law of the codelength as a function of the approximation error in image compression
Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using
More informationArray Design for Superresolution Direction-Finding Algorithms
Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College
More informationTHEORETICAL PROBLEM No. 3 WHY ARE STARS SO LARGE?
THEORETICAL PROBLEM No. 3 WHY ARE STARS SO LARGE? The stars are spheres of hot gas. Most of them shine beause they are fusing hydrogen into helium in their entral parts. In this problem we use onepts of
More informationCavity flow with surface tension past a flat plate
Proeedings of the 7 th International Symposium on Cavitation CAV9 Paper No. ## August 7-, 9, Ann Arbor, Mihigan, USA Cavity flow with surfae tension past a flat plate Yuriy Savhenko Institute of Hydromehanis
More informationCALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS
International Journal of Modern Physis A Vol. 24, No. 5 (2009) 974 986 World Sientifi Publishing Company CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS PAVEL SNOPOK, MARTIN
More informationLECTURE NOTES FOR , FALL 2004
LECTURE NOTES FOR 18.155, FALL 2004 83 12. Cone support and wavefront set In disussing the singular support of a tempered distibution above, notie that singsupp(u) = only implies that u C (R n ), not as
More informationNUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT
Paper ID ILASS8--7 ILASS 28 Sep. 8-, 28, Como Lake, Italy A44 NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT Anne Bagué, Daniel Fuster, Stéphane Popinet + & Stéphane Zaleski Université
More informationEspen Gaarder Haug Norwegian University of Life Sciences April 4, 2017
The Mass Gap, Kg, the Plank Constant and the Gravity Gap The Plank Constant Is a Composite Constant One kg Is 85465435748 0 36 Collisions per Seond The Mass Gap Is.734 0 5 kg and also m p The Possibility
More informationA multiscale description of failure in granular materials
A multisale desription of failure in granular materials Nejib Hadda, François Niot, Lu Sibille, Farhang Radjai, Antoinette Tordesillas et al. Citation: AIP Conf. Pro. 154, 585 (013); doi: 10.1063/1.4811999
More informationA Numerical Method For Constructing Geo-Location Isograms
A Numerial Method For Construting Geo-Loation Isograms Mike Grabbe The Johns Hopkins University Applied Physis Laboratory Laurel, MD Memo Number GVW--U- June 9, 2 Introdution Geo-loation is often performed
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More information