home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #27 / NN_1992_27.iso / spool / sci / math / stat / 2414 < prev    next >
Encoding:
Internet Message Format  |  1992-11-24  |  4.2 KB
  1. Path: sparky!uunet!cis.ohio-state.edu!pacific.mps.ohio-state.edu!linac!uwm.edu!ux1.cso.uiuc.edu!roundup.crhc.uiuc.edu!focus!hougen
  2. From: hougen@focus.csl.uiuc.edu (Darrell Roy Hougen)
  3. Newsgroups: sci.math.stat
  4. Subject: Re: Approximation and Estimation (HELP!!)
  5. Date: 24 Nov 1992 04:41:57 GMT
  6. Organization: Center for Reliable and High-Performance Computing, University of Illinois at Urbana-Champaign
  7. Lines: 76
  8. Message-ID: <1esbqlINN18d@roundup.crhc.uiuc.edu>
  9. References: <24069@galaxy.ucr.edu> <1eribpINNne1@roundup.crhc.uiuc.edu>
  10. NNTP-Posting-Host: focus.csl.uiuc.edu
  11.  
  12. hougen@focus.csl.uiuc.edu (Darrell Roy Hougen) writes:
  13.  
  14. baez@guitar.ucr.edu (john baez) writes:
  15.  
  16. % In article <1ejbebINNber@roundup.crhc.uiuc.edu> hougen@focus.csl.uiuc.edu (Darrell Roy Hougen) writes:
  17. %% I need to estimate a 'hill' function between -pi/2 and pi/2.  It must
  18. %% be 0 at -pi/2 and pi/2 and it must be 1 at 0.  It must increase
  19. %% monotonically from 0 to 1 on [-pi/2,0] and decrease monotonically from
  20. %% 1 to 0 on [0,pi/2].  It must be smooth, ie., differentiable everywhere
  21. %% in (-pi/2,pi/2), and have a derivative of 0 at 0.
  22. %% 
  23. %% My first question is, does there exist a basis for functions of the
  24. %% above form?
  25.  
  26. % Sure, infinitely many.  It's best to say what *kind* of functions  In
  27. % other words, smooth, continuous, or L^2 functions.  But in any of the case
  28. % I just mentioned, there is a countable basis.
  29.     
  30. First, thanks for the information.  Perhaps I wasn't clear, but I
  31. would prefer to have smooth, continuous functions, ie., the first
  32. derivative must exist and be continuous on (-pi/2,pi/2).  Right now I
  33. don't care about the norm too much, as long as it is reasonable.
  34.  
  35. %% One possible basis that I've been thinking about is (cos(x))^k with k
  36. %% ranging from -inf to +inf.  All the members of this set satisfy the
  37. %% conditions given above and it appears that any weighted sum of such
  38. %% functions will also although it would be nice to prove that.
  39.  
  40. % Nope, it's not true, because all the functions you listed are "even"
  41. % (satisfy f(x) = f(-x)), so you can't approximate anything but even 
  42. % functions.    Also for k negative the functions you list are unbounded,
  43. % which is technically a pain in the butt.  
  44.  
  45. Let me simplify things a little by asking for even functions, ie.,
  46. assume that the function to be estimated is even.  Also, I made an
  47. error when I allowed k to be negative.  Instead, I should have defined
  48. my potential basis to be (cos(x))^k for k = 1,2,3,... and
  49. (cos(x))(1/j) for j = 1,2,3,...
  50.  
  51. % The best basis to use is  sin(2kx) and cos(2kx), taken together.  
  52. % This basis is used in the Fourier transform and is well-known to
  53. % be a complete basis.  I recommend this basis because a huge amount
  54. % is known about it and there is a huge amount of software written to
  55. % deal with it.
  56.  
  57. The problem with this basis, a basis of which I am well aware and that
  58. several people have suggested in email, is that the estimate of the
  59. function might be negative in places or have the the wrong slope in
  60. places.  Let me stress that any estimate, g, of the function, f,
  61. **MUST** satisfy the following conditions:
  62.     g(-pi/2) = g(pi/2) = 0
  63.     g(x) > 0  for x in (-pi/2,pi/2)
  64.     g'(x) > 0 for x in (-pi/2,0)
  65.     g'(x) = 0 for x = 0
  66.     g'(x) < 0 for x in (0,pi/2)
  67.  Corresponding to what I said above, it is also exceptable if g(x) =
  68. g(-x).  This reduces the problem to finding a function on [0,pi/2]
  69. that satisfies the above conditions.  
  70.  
  71. It is the above list of conditions that make the problem hard.  In
  72. particular, it appears that every basis function in the set of basis
  73. functions must satisfy the above conditions and that all of the
  74. coefficients must be nonnegative.  If there is a different set of
  75. conditions that yield the same result, it would be nice to hear about
  76. it.  For example, if there is some way to formulate the above problem
  77. as a constrained optimization, it would be nice to know.  My current
  78. proposal for solving the problem is to use my amended set of basis
  79. functions listed above with nonnegative coefficients, ie.,
  80.     g(x) = sum_i a_i * g_i(x) 
  81. where all of the coefficients are forced to satisfy a_i >= 0.
  82.  
  83. Can anyone tell whether this set of basis functions is sufficient for
  84. the problem as I have now stated it, ie., is it a complete basis?  If
  85. not, is there another basis which is?
  86.  
  87. Darrell R. Hougen
  88.