Differential Equation with Summation

jrcdude

1. The problem statement, all variables and given/known data

[tex]y''+0.1y'+y=1+2\sum_{k=1}^{n}(-1)^{k}u_{k\pi}(t)[/tex]
and quiescent initial conditions.

2. Relevant equations

None.

3. The attempt at a solution
[tex](s^{2}+0.1s+1)Y(s)=\mathcal{L}\{1\}+2\sum_{k=1}^{n}(-1)^{k}\mathcal{L}\big\{ u_{k\pi}(t)\big\}[/tex]

I'm not sure if this step was correct, and how to proceed since the result of that is quite nasty.

Any help would be appreciated, thanks.

LCKurtz

We aren't mind readers here. What is ## u_{k\pi}(t)##? And what do you get for the transform of the right side?

jrcdude

I'm sorry, I thought it was common notation that [itex]u_{k\pi}(t)[/itex] is the indicator function denoting the unit step function:
[tex]u_{k\pi}(t)=u_{k}(t)-u_{\pi}(t)=\begin{cases}
0, & t<k\quad\text{or}\quad t\ge\pi\\
1, & k\le t<\pi
\end{cases}
[/tex]
As for what I got for the Laplace transform of the right, it is:
[tex]
\frac{1}{s}+2\sum_{k=1}^{n}(-1)^{k}\cdot\frac{e^{-ks}-e^{-\pi s}}{s}
[/tex]

LCKurtz

For example, if ##k=5## your definition says ##u_{5\pi} = 0\hbox{ if } t<5\hbox{ or }t>\pi##, which is identically zero. Same for any ##k\ge 4##. Something looks fishy to me.

jrcdude

Courtesy of Brannan,


The only other possibility is that [itex]c=k\pi[/itex].


Edit: Yes I think thats actually the problem faced here.. I believe it works if I let
[tex]u_{k\pi}(t)=u(t-k\pi)[/tex]

LCKurtz

It shows what a horrible notation that subscripted u is. If you just write u(t-c) when you want it, there is no ambiguity. So apparently you are saying your equation is$$
y''+0.1y'+y=1+2\sum_{k=1}^{n}(-1)^{k}u(t-k\pi)$$Is that correct? And when you say "it works", do you mean you see how to finish the question?

jrcdude

Yes, that subscript notation is definitely horrible as it caused me a lot of grief haha.

And yes, by "it works" I meant I solved the DE resulting in the solution of,

[tex]
y(t)=h(t)+2\sum_{k=1}^{n}(-1)^{k}h(t-k\pi)
[/tex]

where
[tex]
h(t)=1-e^{-0.05t}\cos(\sqrt{0.9975}t)-\frac{0.05e^{-0.05t}}{\sqrt{0.9975}}\sin(\sqrt{0.9975}t)
[/tex]

Thanks for the 'indirect' help!

LCKurtz

Is that your solution for a particular value of ##n##? Otherwise wouldn't your solution depend on ##n##? Out of curiosity have you graphed your forcing function and solution? What you have is a system being driven for a few cycles by a square wave whose frequency is very near the natural frequency of the system.

[Edit]Never mind about the n. I see where it is now.

jrcdude

The next part of the question was actually to graph the forcing function and solution on the same set of coordinates. So you are interpreting it correctly.

LCKurtz

OK. Have fun with it. I let Maple do the grunt work to look at the graphs. It's a nice example.

jrcdude

Yeah, we too let Maple do the grunt work. Funny thing is Maple was developed at my university so they push us to use it.