Group Abstract Group Abstract

Message Boards Message Boards

0
|
5.8K Views
|
9 Replies
|
5 Total Likes
View groups...
Share
Share this post:

Problems using InverseFunction

Kaushik Kandakatla
Posted 11 years ago
POSTED BY: Kaushik Kandakatla
9 Replies
Sort By:

Clarke,

Thank you for your timely reply :)

Finally I could get correct answers provided I give the correct value for x to search (i.e., 0.1 in the above code you mentioned).

Can you please give me further tips on how to select the value for x to search for?
POSTED BY: Kaushik Kandakatla

Having many orders of magnitude between the two results is fairly normal for a floating point error. Very often you have to come up with a technique specific to the problem. In this case rephrasing as a root finding problem seems sufficient to get a reasonable answer:

FindRoot[200/3 (-(109/500) Log[3/200 - (203 x)/1000] + (203 Log[x])/1000) - 201, {x, 0.1}]
POSTED BY: Sean Clarke

I think I was not clear in my previous post.

I tried to find the inverse function using 2 methods for the same expression. Directly solving using InverseFunction gave me a result which is 0.0101108 ; while solving it first by using Rationalize and then solving it using N gave me a result which is 4.26674*10^9.

Please note that there are many orders of magnitude difference between the both but not a floating point error.

POSTED BY: Kaushik Kandakatla

The difference in value is probably floating point error. If you're comfortable knowing the result is close to 0, then the result might be fine. How much precision are you look for in your result? You may be able to get it using arbitrary precision arithmetic.

http://reference.wolfram.com/language/tutorial/ArbitraryPrecisionNumbers.html

If InverseFunction isn't working out at all, then some serious numerical work will have to be done to get a good answer. Typically, computing the inverse of a function is viewed as a root finding problem. You may want to experiment with different Methods for FindRoot. The three statements below are essentially equivalent. 

NSolve[f[x] == b,x]

x = InverseFunction[f][b]

FindRoot[f[x]-b, {x,0}]
POSTED BY: Sean Clarke

I would really appreciate if someone could help me deal with this issue.
POSTED BY: Kaushik Kandakatla

Hi Clarke,

I tried what you said.

But doing so I got a completely different answer when I used the method you suggested me. Explanation as follows:

When I tried to compute the following expression directly in float, I got the following:

 InverseFunction[ (0.203*Log[#] - 0.218*Log[0.015 - 0.203*#])/0.015 &][1]
 0.0101108

But when I did the same using the method you suggested me (first using rationalize and then solving it using N ), provided me altogether different solution and the solution doesn't change even when I change the input after the input is 24.

Rationalize[(0.203*Log[#] - 0.218*Log[0.015 - 0.203*#])/0.015]
200/3 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)
N[InverseFunction[ 200/3 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000) &][1]]
4.26674*10^9

Please tell me how to deal with this problem.
POSTED BY: Kaushik Kandakatla

I may have miscopied the code. The suggestion however is to try computing the values symbolically in some way. Or you may want to investigate why this isn't numerically stable. 

Rationalize[(0.203*Log[#] - 0.218*Log[0.015 - 0.203*#])/0.015]
200/3 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)

InverseFunction[200/3 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000) &][201]

This gives back a Root expression, which is just a complicated way of describing a number. You can run N on it to try returning a decimal approximation. You may need to use arbitrary precision computing however.

POSTED BY: Sean Clarke

Hi Clarke,

Thanks for your reply. But I would like to point to you that there is a small mistake in what you have mentioned.

After using rationalize for the whole function in first step of the following,

Rationalize@(0.203*Log[#] - 0.218*Log[0.015 - 0.203*#])/0.015
66.6667 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)

why did you again divide the entire output using 0.015 in the inverse function?

InverseFunction[66.6` (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)/0.015 &][201]

Also, if I put it correctly to solve again, Mathematica couldn't solve and give a proper output
POSTED BY: Kaushik Kandakatla

You are using floating point numbers. In most cases InverseFunction is used with symbolic numbers instead because it is a symbolic function. There are major differences between floating point and symbolic numbers.

Use Rationalize to convert most of your formula to symbolic numbers. 

Rationalize@(0.203*Log[#] - 0.218*Log[0.015 - 0.203*#])/0.015
66.6667 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)

This is enough to allow InverseFunction to work. 

InverseFunction[66.6` (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000)/0.015 &][201]

Remember though that you have to be careful with the precision you have. If we use Rationalize[#,0]& as work with a completely symbolic formula, it's possible to get a situation where there is no precision left in your result. 

InverseFunction[200/3 (-(109/500) Log[3/200 - (203 #1)/1000] + (203 Log[#1])/1000) &][201.]

It is easy to make your calculations meaningless by losing all precision. I'm not sure how many digits of accuracy your coefficients have to begin with, but it's possible that InverseFunction might not be a very stable operation in this case.

POSTED BY: Sean Clarke
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard