- #1

asdf1

- 734

- 0

I've tried shifting the coordinate axis, but it doesn't work!

It's also not exact...

You are using an out of date browser. It may not display this or other websites correctly.

You should upgrade or use an alternative browser.

You should upgrade or use an alternative browser.

- Thread starter asdf1
- Start date

- #1

asdf1

- 734

- 0

I've tried shifting the coordinate axis, but it doesn't work!

It's also not exact...

- #2

TD

Homework Helper

- 1,022

- 0

You can write this as

[tex]y' = \frac{{ - x + 2y - 3}}{{2x - 4y + 5}}[/tex]

Now you can see nominator and denominator as two lines. Normally, you'd shift to the intersection but these lines are parallel (since [itex]\frac{{ - 1}}{2} = \frac{2}{{ - 4}}[/itex])

Now, write it as

[tex]y' = \frac{{ - x + 2y - 3}}{{ - 2\left( { - x + 2y} \right) + 5}}[/tex]

And do the substitution [itex]u = - x + 2y[/itex]

Your DE should become separable now.

[tex]y' = \frac{{ - x + 2y - 3}}{{2x - 4y + 5}}[/tex]

Now you can see nominator and denominator as two lines. Normally, you'd shift to the intersection but these lines are parallel (since [itex]\frac{{ - 1}}{2} = \frac{2}{{ - 4}}[/itex])

Now, write it as

[tex]y' = \frac{{ - x + 2y - 3}}{{ - 2\left( { - x + 2y} \right) + 5}}[/tex]

And do the substitution [itex]u = - x + 2y[/itex]

Your DE should become separable now.

Last edited:

- #3

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

[tex]y=\frac{3}{4}ln\left[11-4(2y-x)\right]+\frac{1}{4}(2y-x)+\frac{11}{16}ln\left[11-4(2y-x)\right]+c[/tex]

See what I mean. Same boat as the other one. So I'd say treat it parametrically, that is:

[tex]\frac{dy}{dt}=-x+2y-3[/tex]

[tex]\frac{dx}{dt}=2x-4y+5[/tex]

This time, let's solve it with differential operators:

Rearranging:

[tex]

\begin{align*}

y^{'}-2y+x&=-3 \\

\\

4y+x^{'}-2x&=5

\end{align}

[/tex]

Expressing these in terms of differential operators:

[tex]

\begin{align*}

(D-2)y&+ x &= -3 \\

\\

4y &+ (D-2)x &= 5

\end{align}

[/tex]

Now just "operate" on them to eliminate y and x and end up with two second-order (non-homogeneous) ODEs for x(t) and y(t).

You familiar with this method?

- #4

GCT

Science Advisor

Homework Helper

- 1,742

- 0

are you sure its not exact?

- #5

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

GCT said:are you sure its not exact?

For:

[tex](x-2y+3)dx+(2x-4y+5)dy=0[/tex]

[tex]\frac{\partial}{\partial y}(x-2y+3)=-2[/tex]

[tex]\frac{\partial}{\partial x}(2x-4y+5)=2[/tex]

- #6

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

Educated guessing yields a simpler solution than last time: if you try y = mx + b, you get [itex]y = x/2 + 11/8[/itex]. (The other linear solution is singular! How amusing to have a singular singular solution. At least I feel that these linear solutions are singular in some sense)

I notice that this one is halfway between the two lines...

I notice that this one is halfway between the two lines...

Last edited:

- #7

asdf1

- 734

- 0

i did think trying to parameterize it, except that i had a hard time on it~

- #8

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Hurkyl said:Educated guessing yields a simpler solution than last time: if you try y = mx + b, you get [itex]y = x/2 + 11/8[/itex]. (The other linear solution is singular! How amusing to have a singular singular solution. At least I feel that these linear solutions are singular in some sense)

I notice that this one is halfway between the two lines...

Edit: The linear solutions are NOT singular solutions. See discussion below

Thanks Hurkyl. I had no idea. Your solution y=x/2+11/8 is indeed a singular solution by virtue that it is not a particulr case of the general solution and is an envelop of the general solution as per the attached plot which shows 7 solutions of the 1-parameter general solution all of which are tangent to the singular solution. That it, the singular solution is an envelop of the general solution. Can you answer this for me:

1. How does one know if an ODE will have a singular solution and how does one calculate it analytically (that is, without guessing)?

2. What do you mean by the "other linear solution"?

2. What do you mean by "singular singular solution"?

Added plots:

Plot 2: typical solution below singular solution

Plot 3: Solutions showing vertical tangent when denominator of ODE is 0 that is, whenever y=x/2+5/4

Last edited:

- #9

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

When I plug in y = mx + b and solve, I get two equations. The first is quadratic in *m*, and the second is linear in *b* (with coefficients depending on *m*). This system ought to have two solutions, but when I use the second solution for *m*, *m = -1/2*, it kills the *b* in the other equation, rendering it inconsistent.

So the second linear "solution" is y = -x/2 + ∞.

Note that all of the nonsingular solutions are tangent to this one too!

I don't know, in general, how to spot the singular solutions.

However, in this case, we can find some information by looking asymptotically. When*x* and *y* are very large, we have:

(2x - 4y) y' + (x - 2y) = 0

(x - 2y) (2y' + 1) = 0

So there are two options for asymptotic behavior:

(1) y' ~ -1/2

(2) x - 2y ~ 0

There ought to be singular solutions that follow only one of the asymptotes. In this case, we're lucky, because those solutions are the asymptotes themselves:

(1) y = (-1/2)x + b

(2) y = (1/2)x + b

Ooh, cool idea. I decided this suggests using:

y(x) = (1/2) x + e(x)

Solving this for the general solution gives:

y = C - (1/2) x - (1/8) ln |y - (1/2) x - (11/8)|

If I assume that e(x) is not the constant (11/8).

Which suggests a nice closed form of the type y = f(x) is impossible. However, we can read off the behavior from this!

The function generally either looks like:

y = C - (1/2) x

or

y = (1/2) x + (11/8)

which is exactly what the asymptotic analysis told us. We could pin down further details too... for instance, the y = (1/2) x + (11/8) behavior asymptotically occurs in the direction of positive x.

So the second linear "solution" is y = -x/2 + ∞.

Note that all of the nonsingular solutions are tangent to this one too!

I don't know, in general, how to spot the singular solutions.

However, in this case, we can find some information by looking asymptotically. When

(2x - 4y) y' + (x - 2y) = 0

(x - 2y) (2y' + 1) = 0

So there are two options for asymptotic behavior:

(1) y' ~ -1/2

(2) x - 2y ~ 0

There ought to be singular solutions that follow only one of the asymptotes. In this case, we're lucky, because those solutions are the asymptotes themselves:

(1) y = (-1/2)x + b

(2) y = (1/2)x + b

Ooh, cool idea. I decided this suggests using:

y(x) = (1/2) x + e(x)

Solving this for the general solution gives:

y = C - (1/2) x - (1/8) ln |y - (1/2) x - (11/8)|

If I assume that e(x) is not the constant (11/8).

Which suggests a nice closed form of the type y = f(x) is impossible. However, we can read off the behavior from this!

The function generally either looks like:

y = C - (1/2) x

or

y = (1/2) x + (11/8)

which is exactly what the asymptotic analysis told us. We could pin down further details too... for instance, the y = (1/2) x + (11/8) behavior asymptotically occurs in the direction of positive x.

Last edited:

- #10

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

I find this so confusing:

Solving for the parametric solutions (which I'll post the details later), I obtain:

[tex]x(t)=c_1-1/2t+c_2e^{4t}[/tex]

[tex]y(t)=d_1-1/4t+d_2e^{4t}[/tex]

Where the constants are all functions of a single initial condition: [itex]y(x_0)=y_0[/itex]

Forming the quotient of the respective derivatives I obtain:

[tex]\frac{dy}{dx}=\frac{-1/4+4d_2e^{4t}}{-1/2+4c_2e^{4t}}[/tex]

Now this is the definition of a singular solution as per Rainville and Bedient:

(a) is not a special case of the general solution

(b) is, at each of its points, tangent to some element of the one-parameter family that is the general solution.

But the derivative above is NEVER equal to the derivative of the singular solutions y=x/2+11/8 and y=-x/2+[itex]\infty[/itex] except at t=[itex]\pm \infty[/itex]. That is:

[tex]\mathop\lim\limit_{t\to\pm \infty}\frac{-1/4+4d_2e^{4t}}{-1/2+4c_2e^{4t}}=\mp 1/2[/tex]

The limit as [itex]t\to+\infty[/itex] is due to the algebraic relation between the constants. For example, the case y(0)=a, the constants are:

[tex]c_2=11/8-a[/tex]

[tex]d_2=\frac{a}{2}-\frac{11}{16}[/tex]

I don't understand how then could the linear solutions be truely singular solutions when the only point of tangency is at infinity.

Solving for the parametric solutions (which I'll post the details later), I obtain:

[tex]x(t)=c_1-1/2t+c_2e^{4t}[/tex]

[tex]y(t)=d_1-1/4t+d_2e^{4t}[/tex]

Where the constants are all functions of a single initial condition: [itex]y(x_0)=y_0[/itex]

Forming the quotient of the respective derivatives I obtain:

[tex]\frac{dy}{dx}=\frac{-1/4+4d_2e^{4t}}{-1/2+4c_2e^{4t}}[/tex]

Now this is the definition of a singular solution as per Rainville and Bedient:

(a) is not a special case of the general solution

(b) is, at each of its points, tangent to some element of the one-parameter family that is the general solution.

But the derivative above is NEVER equal to the derivative of the singular solutions y=x/2+11/8 and y=-x/2+[itex]\infty[/itex] except at t=[itex]\pm \infty[/itex]. That is:

[tex]\mathop\lim\limit_{t\to\pm \infty}\frac{-1/4+4d_2e^{4t}}{-1/2+4c_2e^{4t}}=\mp 1/2[/tex]

The limit as [itex]t\to+\infty[/itex] is due to the algebraic relation between the constants. For example, the case y(0)=a, the constants are:

[tex]c_2=11/8-a[/tex]

[tex]d_2=\frac{a}{2}-\frac{11}{16}[/tex]

I don't understand how then could the linear solutions be truely singular solutions when the only point of tangency is at infinity.

Last edited:

- #11

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

Then I suppose it's not singular, then.

Doesn't my linear solution arise when c2 = d2 = 0?

Doesn't my linear solution arise when c2 = d2 = 0?

- #12

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Hurkyl said:Then I suppose it's not singular, then.

Doesn't my linear solution arise when c2 = d2 = 0?

Sure looks like it Hurkyl although c2 and d2 being zero affects what c1 and d1 can be which in turn limits what the initial condition can be. I'll work it through. Thanks!

- #13

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

saltydog said:Sure looks like it Hurkyl although c2 and d2 being zero affects what c1 and d1 can be which in turn limits what the initial condition can be. I'll work it through. Thanks!

Alright it works out. Through the parametric analysis, I got for the constants in the case y(0)=a:

[tex]c_1=a-11/8[/tex]

[tex]c_2=-c_1[/tex]

[tex]d_2=a/2-11/16[/tex]

[tex]d_1=a-d_2[/tex]

Which if [itex]d_2=c_2=0[/itex], that leaves a=11/8 which when solving for t in the expression for x(t) and substituting in for y(t) yield the solution:

[tex]y(x)=11/8+\frac{1}{2} x[/tex]

Nice!

Last edited:

- #14

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

That's the word I was looking for: my linear solution is a

Now, I just need to figure out, generally, what being critical means with respect to the algebra.

- #15

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

I'd like to better understand the relation of this system:

[tex]\frac{dx}{dt}=Ax+By+E[/tex]

[tex]\frac{dy}{dt}=Cx+Dy+F[/tex]

to this system:

[tex]\frac{dx}{dt}=ax+by[/tex]

[tex]\frac{dy}{dt}=cx+dy[/tex]

I understand when AD-BC[itex]\ne[/itex]0, the first system can be converted to the second with the eigenvalues of the second system governing its dynamics. But how is this dynamics related to the original system? Is it the same?

I think all of this is related somehow. Suppose in summary I should ask: How can one determine the global behavior of the first system and can the dynamics be catagorized like that for the second system?

Also in the case of the second system, the straight-line solutions are the eigenvectors of the system. How is that related to the straight-line solution:

[tex]y(x)=x/2+11/8[/tex]

for the system studied in this thread?

[tex]\frac{dx}{dt}=Ax+By+E[/tex]

[tex]\frac{dy}{dt}=Cx+Dy+F[/tex]

to this system:

[tex]\frac{dx}{dt}=ax+by[/tex]

[tex]\frac{dy}{dt}=cx+dy[/tex]

I understand when AD-BC[itex]\ne[/itex]0, the first system can be converted to the second with the eigenvalues of the second system governing its dynamics. But how is this dynamics related to the original system? Is it the same?

I think all of this is related somehow. Suppose in summary I should ask: How can one determine the global behavior of the first system and can the dynamics be catagorized like that for the second system?

Also in the case of the second system, the straight-line solutions are the eigenvectors of the system. How is that related to the straight-line solution:

[tex]y(x)=x/2+11/8[/tex]

for the system studied in this thread?

Last edited:

- #16

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Consider the system:

[tex]\frac{dx}{dt}=2x-4y+a[/tex]

[tex]\frac{dy}{dt}=-x+2y+b[/tex]

of which the problem above is a special case of a=5 and b=-3.

If a and b are zero, then this reduces to the second kind of linear system with eigenvalues 0 and 4. Thus the general solution for this case is:

[tex]

\left(\begin{array}{c} x(t) \\

y(t)

\end{array}\right)=k_1

\left(\begin{array}{c} x_1 \\

y_1

\end{array}\right)+k_2e^{4t}

\left(\begin{array}{c} x_2 \\

y_2

\end{array}\right)

[/tex]

Note this is very similar to the solution with the case a=5 and b=-3. In fact, I could write that case in matrix form as:

[tex]

\left(\begin{array}{c} x(t) \\

y(t)

\end{array}\right)=k_1

\left(\begin{array}{c} x_1 \\

y_1

\end{array}\right)+k_2e^{4t}

\left(\begin{array}{c} x_2 \\

y_2

\end{array}\right)+k_3t

\left(\begin{array}{c} x_3 \\

y_3

\end{array}\right)

[/tex]

Thus I suggest investigating (me anyway) if the general solution of:

[tex]\frac{dx}{dt}=Ax+By+E[/tex]

[tex]\frac{dy}{dt}=Cx+Dy+F[/tex]

consist of a "homogeneous" part obtained from the solution of:

[tex]\frac{dx}{dt}=Ax+By[/tex]

[tex]\frac{dy}{dt}=Cx+Dy[/tex]

and a non-homogeneous part determined by the values of a and b and having the form:

[tex]

k_3t

\left(\begin{array}{c} x_3 \\

y_3

\end{array}\right)

[/tex]

- #17

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

Also, when AD-BC=0 as the case above, then we really have no eigenvalues to determine its behavior.

Yes you do.

Let's look at this one again:

x' = -2x + 4y

y' = x - 2y

a.k.a.

[tex]A := \left(

\begin{array}{cc}

-2 & 4 \\

1 & -2

\end{array}

\right)

[/tex]

We compute the characteristic equation of A:

[tex]f(\lambda) = \left|

\begin{array}{cc}

\lambda + 2 & -4 \\

-1 & \lambda + 2

\end{array}

\right|

= (\lambda + 2)^2 - 4 = \lambda^2 + 4\lambda = \lambda(\lambda + 4)

[/tex]

so its two eigenvalues are 0 and -4.

The eigenvector associated to -4 is [2, -1]^T

The eigenvector associated to 0 is [2, 1]^T

So the general solution should be, I suppose,

This looks right, I think. We recover the translational symmetry observed in the solutions.

Of course, that first base solution is unsatisfying. But it comes from an assumption I made in setting up the problem... I assumed that in order for the ODE to be zero, that I had to pick x' and y' cleverly so that the two terms cancelled out. But, I have another option: I can pick x and y cleverly so that both terms are zero!

Because I tried to pick x' and y' cleverly they both zero out at this particular solution, which means that when I try to take "infinitessimal steps" of size [x', y']^T, I'm not going anywhere. But, now that I've spotted the problem, I can just take the solution for x and y that work, and I've already figured that out from the eigenvector:

[x, y] = t [2, 1]

But I can easily read it off of the equation too, by factoring:

(2x - 4y) y' + (x - 2y) x' = (x - 2y) (2y' + x')

So I just need x - 2y = 0, or y = (1/2) x.

Note that the original equation is not linear, so we cannot just add this to the previously determined solution!

Working it out, we have two families:

[x, y] = (At + B) [2, 1]

and

[x, y] = B [2, 1] + C exp(-4 t) [2, -1]

(Because if you plug in, you discover that AC must be zero)

Stupid

(You notice that the equation is of the form x^T M dx = 0 for some matrix M?)

Last edited:

- #18

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Hurkyl said:Working it out, we have two families:

[x, y] = (At + B) [2, 1]

and

[x, y] = B [2, 1] + C exp(-4 t) [2, -1]

(You notice that the equation is of the form x^T M dx for some matrix M?)

Is the negative exponent a typo? Should it not be positive? Also, I'll spend time with your analysis. Thanks!

- #19

Hurkyl

Staff Emeritus

Science Advisor

Gold Member

- 14,967

- 19

(Notice our original ODE's are off by a factor of -1)

- #20

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Hurkyl said:

(Notice our original ODE's are off by a factor of -1)

Thanks for the clairification Hurkyl. I'll spend time with it. Also, I tell you what, this:

[tex]\frac{dy}{dt}=Ax+By+E[/tex]

[tex]\frac{dx}{dt}=Cx+Dy+F[/tex]

should be comprehensively catagorized and I suspect it is somewhere already. That is, what does it do as a function of the parameters? Surely someone has already worked this out. I just don't have the reference.

- #21

lurflurf

Homework Helper

- 2,452

- 148

matrix form is easiersaltydog said:Thanks for the clairification Hurkyl. I'll spend time with it. Also, I tell you what, this:

[tex]\frac{dy}{dt}=Ax+By+E[/tex]

[tex]\frac{dx}{dt}=Cx+Dy+F[/tex]

should be comprehensively catagorized and I suspect it is somewhere already. That is, what does it do as a function of the parameters? Surely someone has already worked this out. I just don't have the reference.

x'=Ax+b

related homogeneous problem

x'=Ax

solution

x=exp(At)x0

variation of parameters

x=exp(At)u

u'=exp(-At)b

[tex]x=e^{At}x_0+e^{At}\int e^{-At}b dt[/tex]

also undetermined coeficents works well

for b constant

b'=0

Last edited:

- #22

asdf1

- 734

- 0

wow!

@@a

this is the first time i've seen matrix form used in solving differential equations...

@@a

this is the first time i've seen matrix form used in solving differential equations...

- #23

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

lurflurf said:matrix form is easier

x'=Ax+b

related homogeneous problem

x'=Ax

solution

x=exp(At)x0

variation of parameters

x=exp(At)u

u'=exp(-At)b

[tex]x=e^{At}x_0+e^{At}\int e^{-At}b dt[/tex]

also undetermined coeficents works well

for b constant

b'=0

Jesus Lurflurf. Why am I not surprised . . . Can you kindly help me with a matter of notation:

If:

[tex]A=\left[\begin{array}{c} ax+by \\

cx+dy

\end{array}\right]

[/tex]

Then what is:

[tex]e^{At}[/tex]

Or if I'm mis-interpreting your notation would you explain a bit further?

- #24

qbert

- 185

- 5

[tex]

\left( \begin{array}{c} x'(t) \\ y'(t) \end{array} \right) =

\left[ \begin{array}{cc} a & b \\ c & d \end{array} \right]

\left( \begin{array}{c} x(t) \\ y(t) \end{array} \right) +

\left( \begin{array}{c} b_1 \\ b_2 \end{array} \right)

[/tex]

gives

[tex] \mathbf{x'} = \mathbf{Ax + b} [/tex]

given some matrix A then exp(A) is

[tex] I + A + \frac{1}{2!}A^2 + \frac{1}{3!}A^3 + \ldots [/tex]

which I always thought was a singularly unhelpful way to define anything.

It's a little easier if A is diagonal. Then the elements of A

are just a

exponentiate along the main diagonal.

example:

[tex]

\exp{\left[ \begin{array}{cc} a & 0 \\ 0 & d \end{array} \right] } =

\left[ \begin{array}{cc} e^a & 0 \\ 0 & e^d \end{array} \right]

[/tex]

If you can diagonalize A. (That is if there is some matrix such that

TDT

A

= (TD

[tex] \exp{A} = I + A + \frac{1}{2!}A^2 + \frac{1}{3!}A^3 + \ldots [/tex]

[tex] = TT^{-1} + TDT^{-1} + \frac{1}{2!}(TDT^{-1})^2 + \frac{1}{3!}(TDT^{-1})^3 + \ldots [/tex]

[tex] = TT^{-1} + TDT^{-1} + (T \frac{1}{2!}D^2T^{-1}) + (T\frac{1}{3!}D^3T^{-1}) + \ldots [/tex]

[tex] = T(I + D + \frac{1}{2!}D^2 + \frac{1}{3!}D^3 + \ldots)

T^{-1} [/tex]

[tex] \exp{A} = Te^DT^{-1} [/tex]

Of course, it just so happens if you make a matrix whose columns

are the eigenvectors of A you get T. (I'll leave the proof to whomever

is interested.)

If A is invertible and b is constant you can solve the

x' = Ax + b system super easy.

note x' = A(x + A

make a linear change of variables

y = x + A

then y' = x' and the system

is y' = Ay

whose solution is y = exp(At) c

where c is a column vector of arbitrary constants

and

x = y - A

Now the questions comes up, under what

conditions can you diagonalize A? what should

you do if you can't? what if A is not invertible?

etc.

- #25

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

qbert said:

[tex]

\left( \begin{array}{c} x'(t) \\ y'(t) \end{array} \right) =

\left[ \begin{array}{cc} a & b \\ c & d \end{array} \right]

\left( \begin{array}{c} x(t) \\ y(t) \end{array} \right) +

\left( \begin{array}{c} b_1 \\ b_2 \end{array} \right)

[/tex]

gives

[tex] \mathbf{x'} = \mathbf{Ax + b} [/tex]

given some matrix A then exp(A) is

[tex] I + A + \frac{1}{2!}A^2 + \frac{1}{3!}A^3 + \ldots [/tex]

which I always thought was a singularly unhelpful way to define anything.

It's a little easier if A is diagonal. Then the elements of A^{n}

are just a_{ii}^{n}. That is, if A is diagonal you just

exponentiate along the main diagonal.

example:

[tex]

\exp{\left[ \begin{array}{cc} a & 0 \\ 0 & d \end{array} \right] } =

\left[ \begin{array}{cc} e^a & 0 \\ 0 & e^d \end{array} \right]

[/tex]

If you can diagonalize A. (That is if there is some matrix such that

TDT^{-1}= A). Then you can get exp(A) fairly easy.

A^{n}= (TDT^{-1}) (TDT^{-1})...(TDT^{-1})

= (TD^{n}T^{-1})

[tex] \exp{A} = I + A + \frac{1}{2!}A^2 + \frac{1}{3!}A^3 + \ldots [/tex]

[tex] = TT^{-1} + TDT^{-1} + \frac{1}{2!}(TDT^{-1})^2 + \frac{1}{3!}(TDT^{-1})^3 + \ldots [/tex]

[tex] = TT^{-1} + TDT^{-1} + (T \frac{1}{2!}D^2T^{-1}) + (T\frac{1}{3!}D^3T^{-1}) + \ldots [/tex]

[tex] = T(I + D + \frac{1}{2!}D^2 + \frac{1}{3!}D^3 + \ldots)

T^{-1} [/tex]

[tex] \exp{A} = Te^DT^{-1} [/tex]

Of course, it just so happens if you make a matrix whose columns

are the eigenvectors of A you get T. (I'll leave the proof to whomever

is interested.)

If A is invertible and b is constant you can solve the

x' = Ax + b system super easy.

note x' = A(x + A^{-1}b)

make a linear change of variables

y = x + A^{-1}b,

then y' = x' and the system

is y' = Ay

whose solution is y = exp(At) c

where c is a column vector of arbitrary constants

and

x = y - A^{-1}b = exp(At) c - A^{-1}b

Now the questions comes up, under what

conditions can you diagonalize A? what should

you do if you can't? what if A is not invertible?

etc.

Thanks a lot Qbert. I'm not supprised neither but I digress. It's all very interesting to me, and I wish to better obtain a handle on the global dynamics of the system as a whole.

- #26

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Salty

- #27

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

[tex]\frac{dy}{dx}=\frac{-x+2y-3}{2x-4y+5},\quad y(x_0)=y_0[/tex]

This is (parametrically) equivalent to:

[tex]\left(\begin{array}{c} x^{'} \\

y^{'}

\end{array}\right)=

\left(\begin{array}{cc} 2 & -4 \\

-1 & 2

\end{array}\right)

\left(\begin{array}{c} x \\

y

\end{array}\right)+

\left(\begin{array}{c} 5 \\

-3

\end{array}\right);

\quad

\left(\begin{array}{c} x_0 \\

y_0

\end{array}\right)

[/tex]

Now from what I've reviewed thus far, if the determinant of the eigenvectors is non-zero, this is easily solved. Wonderful . . . that's just wonderful . . .

Seem to me that this is by far the best approach to solving this ODE. I'm sure that's what you guys meant up there. No doubt I'd be interested in learning what to do if the eigenvalue determinant is zero but first I'll have to work this one through to completion just to get caught up.

- #28

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Asdf, I'm doing this because I want to, it helps me (relearn) it, and I enjoy it and maybe others interested in the details can learn from it. You, well just do whatever you want ok.

[tex]\frac{dy}{dx}=\frac{-x+2y-3}{2x-4y+5};\quad y(0)=1[/tex]

Converting to parametric form:

[tex]\frac{dx}{dt}=2x-4y+5;\quad x(0)=0[/tex]

[tex]\frac{dy}{dt}=-x+2y-3;\quad y(0)=1[/tex]

Expressing this as a non-homogeneous system of equations:

[tex]\left(\begin{array}{c}x \\ y \end{array}\right)^{'}=

\left(\begin{array}{cc} 2 & -4 \\

-1 & 2 \end{array}\right)

\left(\begin{array}{c}x \\ y \end{array}\right)+

\left(\begin{array}{c} 5 \\ -3 \end{array}\right);\quad

\left(\begin{array}{c}0 \\ 1 \end{array}\right)

[/tex]

Solving for the eigenvalues and eigenvectors for the homogeneous system we obtain:

[tex]\left(\begin{array}{c} x \\ y \end{array}\right)_c=

a_1\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

a_2\left(\begin{array}{c} -2 \\ 1 \end{array}\right)e^{4t}

[/tex]

Employing variation of parameter to find a particular solution we let:

[tex]\left(\begin{array}{c} x \\ y \end{array}\right)_p=

a_1(t)\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

a_2(t)\left(\begin{array}{c} -2 \\ 1 \end{array}\right)e^{4t

[/tex]

Substituting this expression into the non-homogeneous system above and noting that the complimentary solution satisfies the system, results in:

[tex]a_1^{'}(t)\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

a_2^{'}(t)\left(\begin{array}{c} -2 \\ 1 \end{array}\right)e^{4t}=

\left(\begin{array}{c} 5 \\ -3 \end{array}\right)

[/tex]

or:

[tex]\left(\begin{array}{cc} 2 & -2 \\

1 & 1 \end{array}\right)

\left(\begin{array}{c} a_1^{'}(t) \\

a_2^{'}(t)e^{4t} \end{array}\right)=

\left(\begin{array}{c} 5 \\ -3 \end{array}\right)

[/tex]

Now, if the determinant is not zero, this is easily solved for [itex]a_1^{'}[/tex] and [itex]a_2^{'}[/tex] via Cramer's rule.

We find:

[tex]a_1^{'}(t)=-1/4[/tex]

[tex]a_2^{'}(t)=-11/4e^{-4t}[/tex]

Integrating:

[tex]a_1(t)=-1/4t[/tex]

[tex]a_2(t)=11/16e^{-4t}[/tex]

Thus the particular solution is:

[tex]\left(\begin{array}{c} x \\ y \end{array}\right)_p=

-\frac{1}{4}t\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

11/16\left(\begin{array}{c} -2 \\ 1 \end{array}\right)

[/tex]

and therefore the general solution is:

[tex]\left(\begin{array}{c} x \\ y \end{array}\right)_g=

-\frac{1}{4}t\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

11/16\left(\begin{array}{c} -2 \\ 1 \end{array}\right)+

a_1\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

a_2\left(\begin{array}{c} -2 \\ 1 \end{array}\right)e^{4t

[/tex]

where [itex]a_1[/itex] and[itex]a_2[/itex] are now arbitrary constants.

Finally, solving for the initial conditions:

[tex]\left(\begin{array}{c}0 \\ 1 \end{array}\right)=

\left(\begin{array}{c} -22/16 \\

11/16 \end{array}\right)+

a_1\left(\begin{array}{c} 2 \\ 1 \end{array}\right)+

a_2\left(\begin{array}{c} -2 \\ 1 \end{array}\right)

[/tex]

Solving for [itex]a_1[/itex] and [itex]a_2[/itex] and substituting into the general solution, we obtain the same answer as above:

[tex]x(t)=-3/8-\frac{1}{2}t+3/8e^{4t}[/tex]

[tex]y(t)=19/16-\frac{1}{4}t-3/16e^{4t}[/tex]

Oh yea, thanks Hurkyl, lurflurf, and qbert .

Last edited:

- #29

asdf1

- 734

- 0

@@a

jeepers!

i've gotta get out my mathematical physics textbook out and review it too~

thanks! :)

jeepers!

i've gotta get out my mathematical physics textbook out and review it too~

thanks! :)

- #30

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

For the system:

[tex]\frac{dx}{dt}=2x-4y+5;\quad x(0)=0[/tex]

[tex]\frac{dy}{dt}=-x+2y-3;\quad y(0)=1[/tex]

Letting:

[tex]\mathcal{L}\{x(t)\}=u(s)[/tex]

[tex]\mathcal{L}\{y(t)\}=v(s)[/tex]

and taking the Laplace Transform of both sides of both equations we obtain the algebraic system of equations:

[tex]su=2u-4v+5/s[/tex]

[tex]sv-1=-u+2v-3/s[/tex]

Solving for u and v:

[tex]u=-\frac{s+2}{(s-4)s^2}[/tex]

[tex]v=\frac{1-5s+s^s}{(s-4)s^2}[/tex]

Upon taking the Inverse Laplace Transforms of these expressions, we once again obtain the same formulas for x(t) and y(t) as above:

[tex]x(t)=-3/8-\frac{1}{2}t+3/8 e^{4t}[/tex]

[tex]y(t)=19/16-\frac{1}{4}t-3/16e^{4t}[/tex]

- #31

saltydog

Science Advisor

Homework Helper

- 1,591

- 3

Code:

```
alist = {u, v} /.
Solve[{s u == 2 u - 4 v + 5/s, s v - 1 == - u + 2 v - 3/s}, {u, v}]
x = InverseLaplaceTransform[alist[[1, 1]], s, t]
y = InverseLaplaceTransform[alist[[1, 2]], s, t]
```

I find that amazing! Granted, in general, I'd have to include two extra lines to first calculate the transform and this doesn't help one learn the math; I would not recommend this to anyone just learning the technique, but once learned, this provides an effective, concise means of approching the global behavior of these systems.

Share:

- Last Post

- Replies
- 2

- Views
- 2K

- Last Post

- Replies
- 2

- Views
- 2K

- Last Post

- Replies
- 1

- Views
- 3K

- Last Post

- Replies
- 2

- Views
- 1K

- Last Post

- Replies
- 2

- Views
- 3K

- Last Post

- Replies
- 9

- Views
- 3K

- Last Post

- Replies
- 3

- Views
- 2K

- Last Post

- Replies
- 2

- Views
- 879

- Replies
- 11

- Views
- 3K

- Last Post

- Replies
- 2

- Views
- 1K