Rewriting sin(x^(1/2)) for Laplace Transform

  • Thread starter twotaileddemon
  • Start date
In summary: Note that ea/(4s) is a constant, so you can pull it out of the...I think you should. ea/(4s) is a constant, so you can just use it in the...I think you should. ea/(4s) is a constant, so you can just use it in the...In summary, the student is trying to solve for sin(x^(1/2)) using the binomial theorem but 1/2 is not an even or odd number, so they can't use it. They are also trying to solve for sin(x^(1/2)) using the Laplace transform but they are having trouble with the 1/2 term. The
  • #1
twotaileddemon
260
0

Homework Statement



I would like to know how to rewrite sin[x^(1/2)]? I have to do a laplace transform of a certain equation but I need to get it into a more usable form...


Homework Equations





The Attempt at a Solution



I was thinking of using the binomial theorem but 1/2 is neither a even or odd number so I can't use it?
 
Physics news on Phys.org
  • #2
I don't think there is a way to get in a more usable form. Was your equation to solve directly involving sin(√x)?
 
  • #3
rock.freak667 said:
I don't think there is a way to get in a more usable form. Was your equation to solve directly involving sin(√x)?

essentially, yes, where x = at. and

f(t) = sin(√at)

I didn't know how to solve it because there was a t^(1/2) term X_X. Do you have any ideas?

Can I just use the formula L(f) = F(s) = [tex]\int[/tex]f(t)e-st dt where the integration runs from 0 to infinity, even though I have a t^(1/2) term?
 
Last edited:
  • #4
twotaileddemon said:
Can I just use the formula L(f) = F(s) = [tex]\int[/tex]f(t)e-st dt where the integration runs from 0 to infinity, even though I have a t^(1/2) term?

Of course. That is the definition of the Laplace transform after all (you are missing your integration limits though)
 
  • #5
gabbagabbahey said:
Of course. That is the definition of the Laplace transform after all (you are missing your integration limits though)

can you have a Fourier transform with a complex number in it? The answer I got has a complex number in it...
 
  • #6
If [itex]a[/itex] is real, your answer should be too. When you take the Laplace transform of a real-valued function, you get a real-valued function.

Maybe you should show us your calculation...
 
  • #7
You could try 'u' substitution
 
  • #8
gabbagabbahey said:
If [itex]a[/itex] is real, your answer should be too. When you take the Laplace transform of a real-valued function, you get a real-valued function.

Maybe you should show us your calculation...

well I'm not sure how to do all the symbols and stuff very well, but basically, I did

[tex]\int[/tex]sin[(at)1/2]e-stdt

substituted 1/(2i) * ( exp[i(at)1/2] - exp[-i(at)1/2] ) for the sine term, integrated with limits between 0 and infinity, and eventually got

-i / s

as my answer, where i is an imaginary number.
 
  • #9
twotaileddemon said:
well I'm not sure how to do all the symbols and stuff very well

See my sig :wink:

[tex]\int[/tex]sin[(at)1/2]e-stdt

substituted 1/(2i) * ( exp[i(at)1/2] - exp[-i(at)1/2] ) for the sine term, integrated with limits between 0 and infinity, and eventually got

-i / s

as my answer, where i is an imaginary number.

Your method looks fine, but your result doesn't...show the rest of your calculations
 
  • #10
gabbagabbahey said:
See my sig :wink:



Your method looks fine, but your result doesn't...show the rest of your calculations

Okay!

1/(2i) * [ [tex]\int[/tex] exp[i(at)1/2] - st] dt - [tex]\int[/tex] exp[-i(at)1/2] - st] dt ] =

1/(2i) * [ 1 / [ i(a / 2(at)1/2) - s] - 1 / [ -i(a / 2(at)1/2) - s] ]

*note at this point I plugged in 0 for t for the lower limit and infinity for the upper limit of t from integration*
= 1/(2i) * [1/s + 1/s]
= 1/(2i) * [2/s]
= 1/is = -i/s
 
  • #11
You seem to be claiming that

[tex]\int e^{h(t)}dt=\frac{e^{h(t)}}{h'(t)}[/tex]

Is that really true? (Use the product rule to differentiate [tex]\frac{e^{h(t)}}{h'(t)}[/tex] and see if you get [itex]e^{h(t)}[/tex])

Instead, start with the substitution [itex]u=\sqrt{at}[/itex] and then complete the square on the exponent...
 
  • #12
gabbagabbahey said:
You seem to be claiming that

[tex]\int e^{h(t)}dt=\frac{e^{h(t)}}{h'(t)}[/tex]

Is that really true? (Use the product rule to differentiate [tex]\frac{e^{h(t)}}{h'(t)}[/tex] and see if you get [itex]e^{h(t)}[/tex])

Instead, start with the substitution [itex]u=\sqrt{at}[/itex] and then complete the square on the exponent...

I did the substitution and now I have an exponent equal to

exp(iu - su2/a)

and another equal to

exp(-iu - su2/a)

can I complete the square with this and solve it?
 
  • #13
Yes, start by completing the square...
 
  • #14
gabbagabbahey said:
Yes, start by completing the square...

ok well..

exp(iu - su2/a)

= exp[-(su2/a - iu)]
= exp[-(s/a)(u2 - i(a/s)u)]
= exp[-(s/a)(u-(ia/2s)2 - (i2/(4(s/a))]
= exp[-(s/a)(u-(ia/2s)2 + a /(4s)]
...?
 
  • #15
twotaileddemon said:
ok well..

exp(iu - su2/a)

= exp[-(su2/a - iu)]
= exp[-(s/a)(u2 - i(a/s)u)]
Up to here it looks okay
= exp[-(s/a)(u-(ia/2s)2 - (i2/(4(s/a))]
= exp[-(s/a)(u-(ia/2s)2 + a /(4s)]
...?
Your parentheses don't match up here, but what you actually have is probably right.

Note that ea/(4s) is a constant, so you can pull it out of the integral.
 
  • #16
Right, now [tex]e^{-\frac{s}{a}\left(u-\frac{ia}{2s}\right)^2+\frac{a}{4s}}=e^{\frac{a}{4s}}e^{-\frac{s}{a}\left(u-\frac{ia}{2s}\right)^2}[/tex] and [tex]e^{\frac{a}{4s}}[/tex] is a constant and can be pulled out of the integral.

Now, use the substitution [itex]v=\sqrt{\frac{s}{a}}\left(u-\frac{ia}{2s}\right)[/tex]...
 
  • #17
gabbagabbahey said:
Right, now [tex]e^{-\frac{s}{a}\left(u-\frac{ia}{2s}\right)^2+\frac{a}{4s}}=e^{\frac{a}{4s}}e^{-\frac{s}{a}\left(u-\frac{ia}{2s}\right)^2}[/tex] and [tex]e^{\frac{a}{4s}}[/tex] is a constant and can be pulled out of the integral.

Now, use the substitution [itex]v=\sqrt{\frac{s}{a}}\left(u-\frac{ia}{2s}\right)[/tex]...

also.. I have a "u term" from what I did

u = (at)1/2
t = u2/a
du = (a/2u)*dt
dt = (2u/a)*du

so I actually have a 1/ai in front of the integral and another u in it.
That means I have to do integration by parts in addition to what you said, right?
 
  • #18
You won't need to use IBP...what does your integral become after the v-sub?
 
  • #19
gabbagabbahey said:
You won't need to use IBP...what does your integral become after the v-sub?

a big mess!
the limits of integration are from -(i/2)(a/s)1/2 to infinity now, and the equation under the integral is
(a/s)1/2 * [(a/s)1/2v + (ia/2s)]exp(-v2) dv

note: there was a u term before the whole exponential term, so I accounted for that. and then made a substitution for du as well
 
Last edited:
  • #20
Well, I'm sure you know how to integrate [itex]\int ve^{-v^2}dv[/itex] right? And if you use the substitution [itex]v=\sqrt{\frac{s}{a}}\left(u+\frac{ia}{2s}\right)[/itex] on the [tex]e^{-\frac{s}{a}\left(u+\frac{ia}{2s}\right)^2[/tex]...I think you will get some fortuitous cancellation
 
  • #21
gabbagabbahey said:
Well, I'm sure you know how to integrate [itex]\int ve^{-v^2}dv[/itex] right? And if you use the substitution [itex]v=\sqrt{\frac{s}{a}}\left(u+\frac{ia}{2s}\right)[/itex] on the [tex]e^{-\frac{s}{a}\left(u+\frac{ia}{2s}\right)^2[/tex]...I think you will get some fortuitous cancellation

I do know how to integrate that, but how can you make a substitution twice using the same variable? you set v equal to two different things in two different equations?
 
  • #22
Separate into two integrals first...your substitution variable is then a dummy variable in each integral
 
  • #23
gabbagabbahey said:
Separate into two integrals first...your substitution variable is then a dummy variable in each integral

I did separate into two integrals, using different substitutions for v, namely the one with the +(ia/2s) in the equation and the other with -(ia/2s) in the equation.

but then the limits of integration are different for each integral and you can't combine and cancel.

hence.. I was thinking of splitting up each integral further, and having to end up using integration by parts at least once.
 
  • #24
The integrands don't have any poles, so you can shift the contours without changing the values of the integral. In particular, you can just drop the imaginary part in the limits, so the integrals both go from 0 to ∞.
 
  • #25
You'll want to use the fact that if [itex]g(v)[/itex] is even,

[tex]\begin{aligned}\int_{-\alpha}^{\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv &= -\int_{\alpha}^{-\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv \\ &= \int^{\alpha}_{-\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv \\ &= \int_{-\infty}^{\infty} g(v)dv\end{aligned}[/tex]
 
  • #26
vela said:
The integrands don't have any poles, so you can shift the contours without changing the values of the integral. In particular, you can just drop the imaginary part in the limits, so the integrals both go from 0 to ∞.

so if the limits of the first integral were -i(1/2)(a/s)1/2 to infinity and the second were +i(1/2)(a/s)1/2 to infinity, I can just make it 0 to infinity in both cases by ignoring the imaginary parts?

I do see how there are no poles, as we covered this a while ago. I never knew this could affect limits of integration though.. thanks for letting me know!

gabbagabbahey said:
You'll want to use the fact that if [itex]g(v)[/itex] is even,

[tex]\begin{aligned}\int_{-\alpha}^{\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv &= -\int_{\alpha}^{-\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv \\ &= \int^{\alpha}_{-\infty} g(v)dv+\int_{\alpha}^{\infty} g(v)dv \\ &= \int_{-\infty}^{\infty} g(v)dv\end{aligned}[/tex]

we don't have the same function under both integrals though,

under one integral I have (a/s)1/2 * [(a/s)1/2v + ia/(2s)]exp(-v2) and the other (a/s)1/2 * [(a/s)1/2v - ia/(2s)]exp(-v2)

where the first has limits from -i(1/2)(a/s)1/2 to infinity and the other has +i(1/2)(a/s)1/2 to infinity.
 
Last edited:
  • #27
Break up the integrals first...[tex]\int\left[\sqrt{\frac{a}{s}}v + \frac{ia}{2s}\right]e^{-v^2}dv=\sqrt{\frac{a}{s}}\int ve^{-v^2}dv +\frac{ia}{2s}\int e^{-v^2}dv[/tex]
 
  • #28
gabbagabbahey said:
Break up the integrals first...[tex]\int\left[\sqrt{\frac{a}{s}}v + \frac{ia}{2s}\right]e^{-v^2}dv=\sqrt{\frac{a}{s}}\int ve^{-v^2}dv +\frac{ia}{2s}\int e^{-v^2}dv[/tex]
[tex]\int\left[\sqrt{\frac{a}{s}}v + \frac{ia}{2s}\right]e^{-v^2}dv=\sqrt{\frac{a}{s}}\int ve^{-v^2}dv +\frac{ia}{2s}\int e^{-v^2}dv[/tex]

---

[tex]\int\left[\sqrt{\frac{a}{s}}v - \frac{ia}{2s}\right]e^{-v^2}dv=\sqrt{\frac{a}{s}}\int ve^{-v^2}dv -\frac{ia}{2s}\int e^{-v^2}dv[/tex]

then I can combine the two integrals with exp(-v2) from this using the equation you gave me? okay, let me try that...

BTW - thank you both for all of your help so far :) it has been invaluable

so , I ended up getting:

[tex]\begin{aligned}\int_{-\alpha}^{\infty} g(v)dv-\int_{\alpha}^{\infty} g(v)dv &= -\int_{\alpha}^{-\infty} g(v)dv-\int_{\alpha}^{\infty} g(v)dv \\ &= \int^{\alpha}_{-\infty} g(v)dv-\int_{\alpha}^{\infty} g(v)dv \\ &= 0\end{aligned}[/tex]

because it's even... right?
 
Last edited:
  • #29
Do the signs work out for that? The function [itex]ve^{-v^2}[/itex] is odd, and though the function [itex]e^{-v^2}[/itex] is even, the integrals with it as the integrand have the wrong overall signs in the beginning, right?

(I'll just mention I just scribbled some stuff on paper, not being particularly careful, so I could very well have messed up the signs myself.)
 
  • #30
vela said:
Do the signs work out for that? The function [itex]ve^{-v^2}[/itex] is odd, and though the function [itex]e^{-v^2}[/itex] is even, the integrals with it as the integrand have the wrong overall signs in the beginning, right?

(I'll just mention I just scribbled some stuff on paper, not being particularly careful, so I could very well have messed up the signs myself.)

I don't even know, to be honest.

I've been working on this problem for a long, long time today and it's really confusing me =|. I got to a point where I had to integrate exp(-x2) which is apparently not integrable, and wolfram suggested to use the error function and gave me some weird calculation. My answer now not only involves complex numbers but the error function as well.

I'm sure I did the substitution, completing the square, and distributing the integrals while keeping the limits of integration correct, but I'm truly at a loss at this point! I think I may just go for extra help on tuesday unless I can figure out something tomorrow.

I'll consider everything that has been said though and think about it some more. Thanks for your help!
 
  • #31
vela said:
Do the signs work out for that? The function [itex]ve^{-v^2}[/itex] is odd, and though the function [itex]e^{-v^2}[/itex] is even, the integrals with it as the integrand have the wrong overall signs in the beginning, right?

(I'll just mention I just scribbled some stuff on paper, not being particularly careful, so I could very well have messed up the signs myself.)

Works out fine for me (I even verified with Mathematica, and I get the correct results)...you end up subtracting the two odd integrals (which leaves you an integral from -something to +something of the odd function, yielding zero) and adding the two even integrals.
 
  • #32
twotaileddemon said:
I don't even know, to be honest.

I've been working on this problem for a long, long time today and it's really confusing me =|. I got to a point where I had to integrate exp(-x2) which is apparently not integrable, and wolfram suggested to use the error function and gave me some weird calculation. My answer now not only involves complex numbers but the error function as well.

I'm sure I did the substitution, completing the square, and distributing the integrals while keeping the limits of integration correct, but I'm truly at a loss at this point! I think I may just go for extra help on tuesday unless I can figure out something tomorrow.

I'll consider everything that has been said though and think about it some more. Thanks for your help!

Are you alright getting to the point

[tex]\mathcal{L}\left[\sin(\sqrt{at}\right]=\frac{{\rm e}^{-\frac{a}{4s}}}{{\rm i}\sqrt{as}}\left[\int_{-\alpha}^{\infty}\left(\sqrt{\frac{a}{s}}v+\frac{{\rm i}a}{2s}\right){\rm e}^{-v^2}dv-\int_{\alpha}^{\infty}\left(\sqrt{\frac{a}{s}}v-\frac{{\rm i}a}{2s}\right){\rm e}^{-v^2}dv\right][/tex]

where [itex]\alpha=\frac{{\rm i}a}{2s}\sqrt{\frac{a}{s}}[/itex]?

If so, why not post your work from then on...
 
  • #33
gabbagabbahey said:
Works out fine for me (I even verified with Mathematica, and I get the correct results)...you end up subtracting the two odd integrals (which leaves you an integral from -something to +something of the odd function, yielding zero) and adding the two even integrals.
Oh, OK, I was looking at just one integral, assuming the problem would work out differently, and wasn't getting what I thought I would need.
 
  • #34
gabbagabbahey said:
Are you alright getting to the point

[tex]\mathcal{L}\left[\sin(\sqrt{at}\right]=\frac{{\rm e}^{-\frac{a}{4s}}}{{\rm i}\sqrt{as}}\left[\int_{-\alpha}^{\infty}\left(\sqrt{\frac{a}{s}}v+\frac{{\rm i}a}{2s}\right){\rm e}^{-v^2}dv-\int_{\alpha}^{\infty}\left(\sqrt{\frac{a}{s}}v-\frac{{\rm i}a}{2s}\right){\rm e}^{-v^2}dv\right][/tex]

where [itex]\alpha=\frac{{\rm i}a}{2s}\sqrt{\frac{a}{s}}[/itex]?

If so, why not post your work from then on...

I have some work, but I cannot understand how to do the latex for it (i did read your signature, but I easily get lost in all the latex text) But basically, you distribute the integrals, then for the one that has -(ia/2s)*exp(-v^2), change the limits of the integral from -alpha to -infinity because it is even and add a plus sign in front. then switch the limits from (-alpha to -infinity) to (-infinity to -alpha) by adding a minus sign in front again. so then you have

(integral of -alpha to infinity) of (ia/2s)*exp(-v^2) + (integral of -infinity to -alpha) of (-ia/2s)*exp(-v^2),

which is equal to integral from - infinity to infinity of (ia/2s)*exp(-v^2) ... which is zero.
is this right?
 
  • #35
I haven't read through your results but the integral certainly isn't zero. It's the Gaussian integral:

[tex]
\int_{-\infty}^{\infty} e^{-a v^2}dv=\sqrt{\frac{\pi}{a}}
[/tex]
 

Similar threads

Back
Top