Finding basis functions for approximating transcendental function

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
5 replies · 3K views
sauravrt
Messages
15
Reaction score
0
I am working on a problem where I want to approximate a transcendental function of the form
[itex]f(x) = x^Ne^{x}[/itex] for [itex]x \geq 0[/itex] as a linear combination of functions of the form [itex]x^v \text{where} -1 < v < 0[/itex].
How can I find the basis functions of the desired form to represent my transcendental function as a finite linear combination?

If not, what would be approach to obtain finte approximate a transcendental function of the form above ?
 
Mathematics news on Phys.org
Hi mathman, I realize that error. Instead of that I would like my basis to be of the form [itex]\frac{x^v}{1 + x^2} ~\text{where} ~ 0< v < 1[/itex].

What is the general procedure to come up with a basis set?
 
I can't give you an off-hand answer. However your proposal has at least one problem. The functions go to zero for large x, while the functions you are trying to represent become infinite very rapidly.
 
Hi mathman:

I am sorry about the error in the original post. The function I am trying to approximate is [itex]f(x) = x^Ne^{-x}[/itex]. This function goes to zero for large x.
 
sauravrt said:
I am working on a problem where I want to approximate a transcendental function of the form
[itex]f(x) = x^Ne^{x}[/itex] for [itex]x \geq 0[/itex] as a linear combination of functions of the form [itex]x^v \text{where} -1 < v < 0[/itex].
How can I find the basis functions of the desired form to represent my transcendental function as a finite linear combination?

If not, what would be approach to obtain finte approximate a transcendental function of the form above ?

You could use orthogonal polynomials and project your function onto these basis polynomials through an integral transform.

The polynomials themselves will be determined on the basis as well as the interval for projection. If you want to project onto polynomial basis, then there are some texts out there that cover this. If however you want to project onto non-polynomial basis (like discontinuous functions or transendental functions of some sort), then you will need to accommodate for that.

What you can do to derive such polynomials is to use the Gram-Schmidt process to generate the basis functions from first principles. To do this you need to define an inner product, and this is dependent on the interval that you are dealing with.

You should be aware that you will have to use a finite interval: it won't make sense to use an infinite interval as your function is not complete in the L^2 space.

So to sum up:

1) Choose a power for your highest degree for the polynomial
2) Choose the interval that you are approximating for your function
3) Using gram-schmidt and the inner product definition for an L^2 space, generate the basis polynomials
4) Project your function to the basis polynomials to get your coeffecients
5) Write out your approximate function using a linear combination of your basis coeffecients with your basis polynomials.