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

  • Thread starter Thread starter twotaileddemon
  • Start date Start date
Click For Summary

Homework Help Overview

The discussion revolves around rewriting the function sin(√x) in preparation for applying the Laplace transform to an equation involving this function. Participants are exploring how to manipulate the expression to make it more usable for integration.

Discussion Character

  • Exploratory, Mathematical reasoning, Problem interpretation

Approaches and Questions Raised

  • Participants discuss the potential use of the binomial theorem and question its applicability due to the nature of the exponent. There are inquiries about the feasibility of using the Laplace transform definition with a square root term. Some participants express uncertainty about handling complex numbers in the context of Fourier transforms and the implications of real versus complex results.

Discussion Status

The conversation is ongoing, with participants providing guidance on the use of substitution methods and integration techniques. There are suggestions to show calculations and clarify steps, indicating a collaborative effort to understand the problem better. Multiple interpretations of the approach are being explored without a clear consensus on the best method.

Contextual Notes

Participants note the presence of complex numbers in their calculations and discuss the implications of integrating functions with square root terms. There is also mention of integration limits and the handling of imaginary components in the context of Laplace transforms.

  • #31
vela said:
Do the signs work out for that? The function ve^{-v^2} is odd, and though the function e^{-v^2} 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.
 
Physics news on Phys.org
  • #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

\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]

where \alpha=\frac{{\rm i}a}{2s}\sqrt{\frac{a}{s}}?

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

\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]

where \alpha=\frac{{\rm i}a}{2s}\sqrt{\frac{a}{s}}?

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:

<br /> \int_{-\infty}^{\infty} e^{-a v^2}dv=\sqrt{\frac{\pi}{a}}<br />
 
  • #36
Cyosis said:
I haven't read through your results but the integral certainly isn't zero. It's the Gaussian integral:

<br /> \int_{-\infty}^{\infty} e^{-a v^2}dv=\sqrt{\frac{\pi}{a}}<br />

isn't this some kind of error function though? from what I understand it's impossible to integrate exp(-x^2)?
 
  • #37
You can evaluate it since the limits are from -infinity to +infinity.
 
  • #38
vela said:
You can evaluate it since the limits are from -infinity to +infinity.

hmm ok. any other limits would fail though, right?
 
  • #39
Because of the symmetry of the Gaussian, you can also figure out the integral from 0 to infinity, but yeah, for arbitrary limits in general, you have to evaluate the integral numerically.
 
  • #40
Any limit would work, but it just so happens that the integral has a neat solution with the limits you have.
 
  • #41
oh ok!
so after combining the two integrals with (ia/2s) in that way, I ended up getting just a pi1/2 for that. the other two with v*(exp(-v2) ended up canceling out, and so I was left with

1/(as)1/2*exp(-a/4s)*(ia/2s)*pi1/2
= (a*pi)1/2 * exp(-a/4s) / s3/2

It looks weird because I guess I'm used to laplace transforms with "neat" answers, but at least there are no imaginary numbers in the answer.
 
  • #42
I think you still have an overall factor of 1/2 missing from when you wrote the sine in terms of exponentials.
 
  • #43
vela said:
I think you still have an overall factor of 1/2 missing from when you wrote the sine in terms of exponentials.

I did bring that out front, but it ended up canceling with the 2 I had after making the u substitution dt = 2u/a du (at the very beginning)
 
  • #44
Your final answer is missing a factor of 1/2, but other than that it looks good.

As an aside, the trick to evaluating \int_{-\infty}^{\infty}e^{-x^2}dx is to realize that it is identical to \int_{-\infty}^{\infty}e^{-y^2}dy and so

\int_{-\infty}^{\infty}e^{-x^2}dx=\left(\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-(x^2+y^2)}dxdy\right)^{1/2}

which you can easily evaluate by swithcing to polar coordinates.
 
  • #45
twotaileddemon said:
I did bring that out front, but it ended up canceling with the 2 I had after making the u substitution dt = 2u/a du (at the very beginning)
Oh, yeah. Well, you still have a two missing somewhere then.
 
  • #46
gabbagabbahey said:
Your final answer is missing a factor of 1/2, but other than that it looks good.

As an aside, the trick to evaluating \int_{-\infty}^{\infty}e^{-x^2}dx is to realize that it is identical to \int_{-\infty}^{\infty}e^{-y^2}dy and so

\int_{-\infty}^{\infty}e^{-x^2}dx=\left(\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-(x^2+y^2)}dxdy\right)^{1/2}

which you can easily evaluate by swithcing to polar coordinates.

that's really quite interesting and good to know - thank you!

I checked my work over, and I don't -think- I'm missing a factor of 1/2. When I converted from sin[(at)^(1/2)] to sin u, I made the substitution u = (at)^(1/2), and so du = (a/2u)dt and dt = (2u/a)du. I brought this 2 out front, and when I rewrote sinu in terms of complex exponential, I brought out the 1/(2i) to the front and the 2's canceled.

EDIT:

vela said:
Oh, yeah. Well, you still have a two missing somewhere then.

actually, yeah, I found it. I had it in my calculations, but when simplifying the terms I forgot to include it in the final answer. Thanks for pointing that out!

Thanks again for your patience and guidance (to everyone involved) =]
 
  • #47
The missing 1/2 came from the factor of ia/2s you got after the v-sub. You had it in your second last line of post#41, but then it inexplicably disappeared in your last line.
 

Similar threads

  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
3
Views
2K
Replies
5
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K