MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Simplifying polynomial and rounding problems

  • To: mathgroup at smc.vnet.net
  • Subject: [mg82913] Re: [mg82801] Simplifying polynomial and rounding problems
  • From: DrMajorBob <drmajorbob at bigfoot.com>
  • Date: Sat, 3 Nov 2007 03:25:00 -0500 (EST)
  • References: <13454456.1194022919313.JavaMail.root@m35>
  • Reply-to: drmajorbob at bigfoot.com

You can approximate the polynomial with, for instance:

a1 = 1.0000000000000002 - 0.6721560269293043 s^2 -
    0.32784397307069685 s^4;
a2 = Rationalize[a1, 10^-14]
a2 // Factor

1 - (8292319 s^2)/12336896 - (4044577 s^4)/12336896

-(((-1 + s) (1 + s) (12336896 + 4044577 s^2))/12336896)

1 is an exact root of a2, and the approximation error is

a1 - a2 // InputForm

2.220446049250313*^-16 - 2.220446049250313*^-15*
   s^2 + 9.992007221626409*^-16*s^4

Depending on your application and the value of s, that may -- or may not 
-- be small enough (in theory).

In practice, I doubt that you know the coefficients to more than six  
digits in the first place.

Bobby

On Wed, 31 Oct 2007 06:14:15 -0500, Isaac Martinez G.  
<isaac.martinez at sbcglobal.net> wrote:

> I am having problems with Mathematica precision/accuracy (or whatever 
> you call it)
> I have the following expression
> A1=1.0000000000000002-0.6721560269293043s^2 - 0.32784397307069685s^4
> Now when I factorize:
> There are a couple of (-1+s)(1+s) factors or something very close to 1.  
> Like 0.999999999999999997+s
> I want to elliminate those factors but due to the fact that it is not 
> really 1, but 0.999999999999999997 Mathematica does not always cancel 
> those factors neither with A1/((-1.+s)(1.+s)) nor A1/((-1+s)(1+s)).
> How can I avoid the problem of 1 versus 0.999999999999999997
> I am doing directional coupler design and I do need to elliminate some 
> factors from an ABCD matrix where 5 to 6 decimal places is just fine in  
> the last step.
> I have tried
> Factor[]
> Cancel[]
> FactorList[]
> but I always have this problem and it screws up my program.
> Sorry, I am  new to Mathematica
>
> I am new to Mathematica.....
>
>



-- 
DrMajorBob at bigfoot.com


  • Prev by Date: Re: Problem with supplying custom tick marks
  • Next by Date: Re: Dynamic palette of open notebooks
  • Previous by thread: Re: Simplifying polynomial and rounding problems
  • Next by thread: Re: Which symbolic method for this Integral ?