Finding extremums of function with complex derivate



#2

Given is the function

I want to calculate the maximums auf this function using the FCN->EXTR function of the 50g. This fails with the message 'Bad Guess(es)', because the derivate of f is complex. Why is f' complex, whereas f only returns real values?
Adding a ->NUM to the end of the Y-Var program slows down the plot and sends the calculator into an infinite loop when trying to find the maximum of f.
The TI-83 fails at this, too. It finds the maximum at the guessed one.
Has anyone an idea how to find the maximum of f (using the 50g)?


#3

For most values of x the function would approach 4. You can see in this function , can you use the estimation?


Edited: 15 Oct 2011, 8:56 a.m.


#4

I think the URL you pasted didn't work right,how about this?

Yeah, the exponential will be of the form (cos(t)+isin(t)).

So that plus 1 has a maximum absolute value when it equals 2. (when t=0, 2pi, 4pi, etc)

So you need to solve for (x+3.1e-6)/x = an integer.

Or x = 3.1e-6/(n-1), with n an integer.

I'm not sure why it seems like the derivative is complex; it shouldn't be. I'd have to check what the HP50g gives for that later.

Edited: 15 Oct 2011, 10:58 a.m. after one or more responses were posted


#5

My URL didn't work either; the parentheses get lost. Oh, well.

#6

Well, plotting and finding it's maximums did the trick. Why didn't I think of this earlier? :(

f' seems to be complex, because the derivate of |a+bi| is signum(a+bi), which returns a complex value :(


#7

I tried playing around with the HP50g with that expression, and it can give wildly varying forms depending on what you do with it. Some of which do not contain complex numbers at all.

At least, make sure rigorous is on, so it doesn't simply return abs(x) as x.

This is probably a situation for which the TI89 gives a more straight forward result.


Edited: 15 Oct 2011, 7:25 p.m.

#8

Quote:
Well, plotting and finding it's maximums did the trick. Why didn't I think of this earlier? :(

f' seems to be complex, because the derivate of |a+bi| is signum(a+bi), which returns a complex value :(


Are you sure that Re(f) = |f| ? If this is true, then you can graph Re(f) and get your answer that way.

How about using the fact that |z|^2 = z * conj(z) and e^(i*t) = cos(t)+i*sin(t)?


#9

It seems that: Re(z)!=|z|^2
That said, max of f is 4 for x=3.1e-6/k, k being 0,1,2....

My two cents...

Regards


Possibly Related Threads…
Thread Author Replies Views Last Post
  HP50g: Writing a function that returns a function Chris de Castro 2 2,143 12-10-2013, 06:49 PM
Last Post: Han
  HP Prime: complex numbers in CAS. Alberto Candel 1 1,926 12-06-2013, 02:36 PM
Last Post: parisse
  [HP Prime] Plots containing complex numbers bug? Chris Pem10 7 3,675 12-05-2013, 07:40 AM
Last Post: cyrille de Brébisson
  Complex Number Entry on Prime Jeff O. 19 5,254 11-16-2013, 12:34 PM
Last Post: Jeff O.
  Extremums on Prime Richard Berler 2 1,247 10-13-2013, 01:21 PM
Last Post: Richard Berler
  HP Prime complex results Javier Goizueta 0 1,004 10-06-2013, 12:59 PM
Last Post: Javier Goizueta
  HP Prime Solving Nonlinear System of Equations for Complex Results Helge Gabert 11 4,310 09-30-2013, 03:44 AM
Last Post: From Hong Kong
  HP 35s derivatives/extremums Richard Berler 16 4,645 09-14-2013, 04:55 PM
Last Post: Matt Agajanian
  [HP-Prime xcas] operations with complex numbers + BUGs + Request CompSystems 9 3,538 09-08-2013, 10:40 PM
Last Post: CompSystems
  Elliptic integrals of 1st and 2nd kind calculated by complex agm Gjermund Skailand 3 1,492 06-29-2013, 03:39 PM
Last Post: Gjermund Skailand

Forum Jump: