Stationary solution to differential equation

In summary, you need to solve the differential equation to find the amplitude to the sine wave. You can do this by finding the particular solution or using the method of undetermined coefficients.
  • #1
Twinflower
109
0

Homework Statement



What is the stationary amplitude to the sine wave?

[tex]4 \frac{dy}{dt} + 3y = sin(t) , y(0) = 0[/tex]


2. The attempt at a solution

I figure that the stationary solution to this eq. is found when the derivate is 0.

[tex]0+3y = sin(x)[/tex]
[tex]y = \frac{sin(x)}{3}[/tex]

The amplitude to this function is [itex]\frac{1}{3}[/itex]. The correct answer is supposed to be 0.2

So where did I take the wrong turn?
 
Physics news on Phys.org
  • #2
Hi Twinflower,

I'm afraid that a stationary amplitude to the sine wave does not imply that you have a stationary point in the solution for your DE.

Basically you need to solve the DE.
Do you know how to find the particular solution?
 
  • #3
To be honest, I thougt that the stationary amplitude was given by the particular solution.

But I'll get back to this after dinner :) Hopefully with some more calculations
 
  • #4
Twinflower said:
To be honest, I thougt that the stationary amplitude was given by the particular solution.

But I'll get back to this after dinner :) Hopefully with some more calculations

It is! :)
That is, you can deduce the stationary amplitude from the particular solution.
It's just that sin(t)/3 is not a particular solution.
 
  • #5
Ok, I am a little confused now.

In my earlier exercises, (with DC voltage), the stationary result was always the particular one.
And the particular one was always the "original" equation where the derivative was 0.

Is something different here because of the sine function, or have I misunderstood something?
 
  • #6
Yes, I guess something is different because of the sine function.

In other cases you probably had a nice constant result for y which would indeed be an amplitude.
In this case you find that the amplitude of y equals sin(t)/3, but this is not a constant expression, and you can't tell yet for which value of t you get the amplitude.
The only thing you can tell, is that the amplitude is at most 1/3.
 
  • #7
Yes, I indeed have had nice constants :)
Made my life easy and uptil now I felt that I had somewhat in control of differential equations.

Suppose this sine function messed it up a bit for me.

What exactly is a constant expression? Do I have to derivate the particular one?
 
  • #8
"1" is a constant expression.
Or "C", assuming that C is a constant (meaning it does not depend on t or y).

What you need to find is a solution to your DE.
Since the RHS is a sine, it seems likely that a combination of a cosine and a sine with unknown amplitudes would do the trick...
 
  • #9
Ah, Ok. So I have to solve the homogenous (or something) part and use the "expected solution format", and then find the unknown new constant using the given start condition? (y(0) = 0).

So, the expected solution for the homogeous part is A sin(t) + B cos (t) or something?
(We never learn to find those expected solutions, unfortunately)
 
  • #10
Yep!

(It's called an inhomogeneous DE, and a single solution for it is called a particular solution.
This method to solve it is called the method of undetermined coefficients.)
 
  • #11
The solution of the homogeneous part of this ode tends to zero with time, leaving the partial solution only. The partial solution is a general sine function with the same angular frequency as the sine function on the right-hand side, but with unknown amplitude and phase Y=Asin(ωt+θ) (ω=1 in the example). This oscillatory solution is stationary in the sense that the amplitude and phase constant does not change with time. The problem asks the amplitude A.

ehild
 
  • #12
OK, I feel I'm getting a few steps in the right direction, but I'm not completely sure where I'm headed.

This is what I think so far:

Particular solution:
[tex]y_p = \frac{sin(t)}{3}[/tex]Homogeneous solution
[tex]y_h = 4 \frac{dy}{dx} + 3y = 0[/tex]

Converting the derivate io [itex]\lambda[/itex] and the variable (y) into 1:
[tex]y_h = 4 \lambda +3 = 0[/tex]
[tex]\lambda = -\frac{3}{4}[/tex]

OK now. The first thing I am uncertain of is where to use [itex]\lambda[/itex].
So let's just try something:

Since the inhomogeneous part is [itex]k sin (\alpha x)[/itex] the guessed form of the homogeneous part would be [itex]K cos (\alpha x) + M sin (\alpha x)[/itex]
(referring to http://en.wikipedia.org/wiki/Method...ents#Typical_forms_of_the_particular_integral )

Now, since the [itex]\alpha[/itex] is used the same way I am used to use [itex]\lambda[/itex], I am going to try to put lambda where alpha is.[tex]K cos (-\frac{3}{4} t) + M sin (-\frac{3}{4} t)[/tex]And I need both parts here, that is [itex]y(t) = y_p + y_h[/itex]

Thus:
[tex]y(0) = 0 = \frac{sin(t)}{3} + K cos (-\frac{3}{4} t) + M sin (-\frac{3}{4} t)[/tex]
Because t = 0, most parts say poof and disappear.
That is, leaving me with 0 = K because cos(0) = 1

Now it looks like I have to dig another step down to find the M coefficient.
That means that I have to find y'(0), and that's where I get really unsure of what to do..
 
  • #13
The value of lambda=(-3/4) tells you the homogeneous solution is e^(-3t/4). And you aren't following your rule to guess a particular form. Since the inhomogeneous part is sin(t), I would guess the form y=K*cos(t)+M*sin(t). Now find K and M. Put your guessed form into 4*dy/dt+3y=sin(t).
 
  • #14
Twinflower said:
OK, I feel I'm getting a few steps in the right direction, but I'm not completely sure where I'm headed.

I'm impressed that you're diving in! :smile:
Twinflower said:
This is what I think so far:

Particular solution:
[tex]y_p = \frac{sin(t)}{3}[/tex]

A particular solution should be a solution to the (inhomogeneous) DE.
Can you substitute this solution in the DE? See if it fits?
Twinflower said:
Homogeneous solution
[tex]y_h: 4 \frac{dy}{dx} + 3y = 0[/tex]

Converting the derivate io [itex]\lambda[/itex] and the variable (y) into 1:
[tex]y_h: 4 \lambda +3 = 0[/tex]
[tex]\lambda = -\frac{3}{4}[/tex]

OK now. The first thing I am uncertain of is where to use [itex]\lambda[/itex].

This is the homogeneous part of the DE, with the RHS (that is independent of y' and y) replaced by zero.
You have the right method, solving it with λ, but you should know that the resulting homogeneous solution is:
[tex]y_h = C_1 e^{\lambda t}[/tex]
Twinflower said:
So let's just try something:

Since the inhomogeneous part is [itex]k sin (\alpha x)[/itex] the guessed form of the homogeneous part would be [itex]K cos (\alpha x) + M sin (\alpha x)[/itex]
(referring to http://en.wikipedia.org/wiki/Method...ents#Typical_forms_of_the_particular_integral )

In your case the inhomogeneous part is [itex]\sin t[/itex].
You should match this with [itex]k sin (\alpha x)[/itex], deducing k and alpha.
The particular solution would then be of the form:
[tex]y_p = K \cos (\alpha t) + M \sin (\alpha t)[/tex]
 
Last edited:
  • #15
Aha, so I need to "transform" the particular part as well?
I've thought that the method of undetermined coefficients was only used in the homogeneous part.

But am I right in asuming that I can disregard the homogeneous part as the problem ask about the stationary condition and the homogeneous part "disappears" after about 5 time constants?
 
  • #16
Meaning that
[tex]y(t) = y_p + y_h[/tex]
[tex]y(t) = K cos (\alpha t) + M sin (\alpha t) + D e^{\lambda t}[/tex]

But first I have to find the undetermined coefficients K and M.
Det undet. constant D in the homogeneous part is later found when K and M is found, and when inserting the starting condition y(0) = 0.

So basicly I am stuck at the particular part.
Following is a few suggestions/questions:[tex]y_p = K cos (\alpha t) + M sin (\alpha t) = ?[/tex]
What should I put on the RHS to deduce K and M?

Am I right in asuming that I cannot put the derivate to zero because the resulting funtion is a sine where the derivative is constantly changing?
 
  • #17
Twinflower said:
[tex]y_p = K cos (\alpha t) + M sin (\alpha t) = ?[/tex]
What should I put on the RHS to deduce K and M?

You should put the function with the actual alpha value - what is it?

Then substitute yp into the original de and find the constants K and M.

Twinflower said:
Am I right in asuming that I cannot put the derivate to zero because the resulting funtion is a sine where the derivative is constantly changing?

Yes, that is right.

ehild
 
  • #18
Is the function with the actual alpha value just sin(t) ?
And in that case, is alpha 1 ?
 
  • #19
Yes, alpha=1.

ehild
 
  • #20
Ok, thanks ehild.

Still leaves me with 2 unknowns in one equation.
So, if t = 0, I get that K = 0.

Should I derive the Kcos(αt)+Msin(αt) and insert K = 0 to find M ?

I am completely out in the bushes on this one :/
 
  • #21
No, Yp is not equal to sin(t) and K is not zero.

Yp=Kcos(t)+Msin(t) is a solution of the original ode:

[tex]4 \frac{dy}{dt} + 3y = sin(t)[/tex]

Find the derivative of Yp and substitute dYp/dt for dy/dt,
and Yp for y.

ehild
 
  • #22
Start with:
[tex]4 \frac{dy}{dt} + 3y = sin(t)[/tex]
Substitute:
[tex]y=K\cos t + M\sin t[/tex]
What is the derivative of this expression? For you need to substitute that as well.
 
  • #23
I like Serena said:
Start with:
[tex]4 \frac{dy}{dt} + 3y = sin(t)[/tex]
Substitute:
[tex]y=K\cos t + M\sin t[/tex]
What is the derivative of this expression? For you need to substitute that as well.
Ok..

the derivate of [itex]y = K cos(t) + M sin(t)[/itex] equals [itex]M cos(t) - K sin(t)[/itex]
This is now dYp/dt which I am now supposed to insert as dy/dt in the original equation:
(also y is swapped with yp)

[tex]4\times M cos(t) - K sin(t) + 3( K cos(t) + M sin(t)) = sin(t)[/tex]

(still in the bushes picking cherries)
 
  • #24
Twinflower said:
Ok..

the derivate of [itex]y = K cos(t) + M sin(t)[/itex] equals [itex]M cos(t) - K sin(t)[/itex]


[tex]4\times M cos(t) - K sin(t) + 3( K cos(t) + M sin(t)) = sin(t)[/tex]

(still in the bushes picking cherries)

You miss some parentheses.

ehild
 
  • #25
After you have put in the missing parentheses, collect the cos(t) terms and also the sin(t) terms. The result is something like a cos(t)+b sin(t)=0.
What did yo get?
(If you find any cherries send over some...:smile:)

ehild
 
  • #26
ehild said:
After you have put in the missing parentheses, collect the cos(t) terms and also the sin(t) terms. The result is something like a cos(t)+b sin(t)=0.
What did yo get?
(If you find any cherries send over some...:smile:)

ehild


Ok, new parantheses is ordered and installed as requested:


[tex]4\times (M cos(t) - K sin(t)) + 3( K cos(t) + M sin(t)) = sin(t)[/tex]

Moving the RHS to the coolhandside:

[tex]4\times (M cos(t) - K sin(t)) + 3\times (K cos(t) + M sin(t)) - sin(t) = 0[/tex]



(And by the way:
48px-Fruit-cherries.svg.png
 
  • #27
Twinflower said:
Ok, new parantheses is ordered and installed as requested:[tex]4\times (M cos(t) - K sin(t)) + 3( K cos(t) + M sin(t)) = sin(t)[/tex]

Moving the RHS to the coolhandside:

[tex]4\times (M cos(t) - K sin(t)) + 3\times (K cos(t) + M sin(t)) - sin(t) = 0[/tex](And by the way:
48px-Fruit-cherries.svg.png

Thank you for the cherries! Sweet and juicy and red and big...
But what about collecting the terms with sin(t) and cos(t)?

I meant [tex]cos(t)(4M+3K) +sin(t)(3M-4K-1)= 0[/tex]

This equation must be valid for every t. What do you get for the unknow K and M in case t=?
And what comes out if t=pi/2, that is cos(t)=0?ehild
 
Last edited:
  • #28
Ok, so that's going to be something like this:

[tex]7\times M cos(t) -(K+1) sin (t) = 0[/tex]...right?
 
  • #29
Twinflower said:
Ok, so that's going to be something like this:

[tex]7\times M cos(t) -(K+1) sin (t) = 0[/tex]...right?

No... I did it in the previous post already. Read it. I have to leave for a while.

[tex]cos(t)(4M+3K) +sin(t)(3M-4K-1)= 0[/tex]Solve for K and M.
 
  • #30
Cedar+Waxwing+eating+berries.jpg
 
  • #31
I am more on this level of cherry picking:

http://tothevillagesquare.org/images/cherry-picking-numbers.jpg

I also regret to inform you that I have officially given up. This kind of problem is unlikely to show up on my exam tomorrow, and *if* one pops up, I can live with not completing it.
It won't ruin my A (or B or C or D for that matter).

I'm sorry to dissappoint you guys, after all your effort. Hopefully my teacher will cover this topic more thouroughly during the next semester. (which I've heard is almost a dedicated differential equation semester)
 
Last edited by a moderator:
  • #32
Aaaah, you like 'em big! :smile:

Will you let us know how your exam went?
(And if something like this popped up?)
 
  • #33
Of course :)

(and btw, sorry about the huge dimmensions )
 
  • #34
It went very well for my exam, and luckily no differential equations containing those spooky sine functions :)

After studying the solution which was published earlier this evening, I feel like I'm in the right spot for an A :)
 
  • #35
Good! :smile:

And I guess we'll get back to DE's later. ;)
 

1. What is a stationary solution to a differential equation?

A stationary solution to a differential equation is a solution that does not change over time. In other words, the derivative of the solution with respect to time is equal to zero. This means that the solution remains constant and does not depend on any external factors.

2. How is a stationary solution different from a general solution?

A general solution to a differential equation includes all possible solutions, while a stationary solution is a specific type of solution that remains constant over time. A general solution may include both stationary and non-stationary solutions.

3. What are some real-life examples of stationary solutions?

Some real-life examples of stationary solutions include the temperature of a room that is kept at a constant level by a thermostat, the population of a species that is at equilibrium, and the voltage across a capacitor in a steady state circuit.

4. How do you find a stationary solution to a differential equation?

To find a stationary solution to a differential equation, you can set the derivative of the solution with respect to time equal to zero and solve for the constant value. This will give you the specific solution that remains constant over time.

5. What is the significance of a stationary solution in scientific research?

Stationary solutions are important in scientific research as they represent stable and predictable states in systems. They can help scientists understand the behavior of a system over time and make predictions about its future behavior.

Similar threads

  • Calculus and Beyond Homework Help
Replies
18
Views
1K
  • Calculus and Beyond Homework Help
Replies
2
Views
719
  • Calculus and Beyond Homework Help
Replies
25
Views
312
  • Calculus and Beyond Homework Help
Replies
7
Views
254
  • Calculus and Beyond Homework Help
Replies
10
Views
305
  • Calculus and Beyond Homework Help
Replies
7
Views
674
  • Calculus and Beyond Homework Help
Replies
21
Views
826
  • Calculus and Beyond Homework Help
Replies
6
Views
840
  • Calculus and Beyond Homework Help
Replies
1
Views
819
  • Calculus and Beyond Homework Help
Replies
2
Views
252
Back
Top