Question: How can I get Maple to find the Fourier transform of a logistic pdf? Also, a question about subs() vs. algsubs().

I'm a Maple novice.  I have two questions.  (1) I'm trying to use Maple to confirm some Fourier transforms of selected probability density functions.  In general I've succeeded, but Maple fails to find the Fourier transform of the pdf of a logistic random variable with mean 0.  Please explain how I can get Maple to carry out this request.  I've attached a Maple file with my work to this question.  (2) At several points in my computations, I wish to substitute 2*Pi*xi for omega.  I have an expression containing two omegas.  If I use algsubs(), only one of the two omegas is replaced.  I have to use subs() to replace both omegas.  Why is this?
 

Use Maple to confirm selected Fourier transform of logistic random variable

with(inttrans)

[addtable, fourier, fouriercos, fouriersin, hankel, hilbert, invfourier, invhilbert, invlaplace, invmellin, laplace, mellin, savetable]

(1)

assume(a > 0)

NULL

"f(x):=((e)^(x/(a)))/(a*(1+(e)^(x/(a)))^(2))"

proc (x) options operator, arrow, function_assign; exp(x/a)/(a*(1+exp(x/a))^2) end proc

(2)

int(f(x), x = -infinity .. infinity)

1

(3)

fourier(f(x), x, omega)

fourier(exp(x/a)/(1+exp(x/a))^2, x, omega)/a

(4)

"f(x) := 1/(4*a*(cosh(x/(2*a)))^(2))"

proc (x) options operator, arrow, function_assign; (1/4)/(a*cosh((1/2)*x/a)^2) end proc

(5)

int(f(x), x = -infinity .. infinity)

1

(6)

fourier(f(x), x, omega)

fourier(1/(exp((1/2)*x/a)+exp(-(1/2)*x/a))^2, x, omega)/a

(7)

Wikipedia's article on "Logistic distribution" gives a characteristic function that implies that the Fourier transform of this pdf should equal Pi*a*omega/sinh(Pi*a*omega).

Unit rectangular function

rect := proc (x) options operator, arrow; Heaviside(x+1/2)-Heaviside(x-1/2) end proc

proc (x) options operator, arrow; Heaviside(x+1/2)-Heaviside(x-1/2) end proc

(8)

fourier(rect(a*x), x, omega)

2*sin((1/2)*omega/a)/omega

(9)

algsubs(omega = 2*Pi*xi, 2*sin((1/2)*omega/a)/omega)

2*sin(Pi*xi/a)/omega

(10)

subs(omega = 2*Pi*xi, 2*sin(Pi*xi/a)/omega)

sin(Pi*xi/a)/(Pi*xi)

(11)

``

NULL


 

Download Logistic_pdf.mw

 

Please Wait...