🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: sci.math.num-analysis
Article #1205

Re: Solutions to Power Series

#1205
From: Brian Smith
Date: Mon, 26 Jul 2004 12:51
123 lines
5035 bytes
Jon <jon8338@peoplepc.com> wrote in message news:<MSHMc.12907$f4.5254@newsread3.news.atl.earthlink.net>...
> a correction is applied below:
> 
> Jon wrote:
> > If T*N+a[0] is the nth degree polynomial where T=(t,t^2,t^3,...,t^n), 
> > N=(a[1],a[2],a[3],...,a[n]), if T is described from two orthogonal 
> > systems, its linear requirement precipitates.  T remains invariant
> > despite the coordinate system used to describe it. Consequently, my 
> > revised web site showing the details,
> > 
> > APPLICATION OF
> > http://mypeoplepc.com/members/jon8338/polynomial/id7.html
> > 
> > EXAMPLE
> > 

Alright Jon, its time for some strong criticism.  Take a look at this:
t^3 - 2t^2 - 2t + 4 = 0
t = 2, sqrt(2), -sqrt(2)

> > (t-2)(t-3)(t-4)=0
> > (t^2-5t+6)(t-4)=0
> >       t-4
> > -------------
> > t^3-5t^2+ 6t
> >    -4t^2+20t-24
> > ----------------
> > t^3-9t^2+26t-24
> > 

N = (1, -2, -2)

> > N=+/-(26,-9,1)
> > N*U[2]=0 at U[2]=(4,13,13)
> >               i    j  k
> > U[3]=NxU[2]= 26   -9  1 =(-130,-334,374)
> >               4   13 13
> > 
> > N/|N|=
> > u[1]=( .944,-.327,.036)=(u[11],u[12],u[13])
> > u[2]=( .212, .691,.691)=(u[21],u[22],u[23])
> > u[3]=(-.249,-.640,.717)=(u[31],u[32],u[33])
> > 

I have decided to try this polynomial with two sets of vectors one is
marked with primes and the other is not.

u[1] = (1/3, -2/3, -2/3)
u[2] = (2/3, -1/3, 2/3)
u[3] = (2/3, 2/3, -1/3)

u'[1] = (1/3, -2/3, -2/3)
u'[2] = (2/sqrt(5), 1/sqrt(5), 0)
u'[3] = (2/sqrt(5), 0, 1/sqrt(5))

> > (-a[0]/|N|)u[1]=(.8717)u[1]
> > Q = +/-(-.823,.285,-.0314)
> > 

Q = -(4)/3 * (1/3, -2/3, -2/3) = (-4/9, 8/9, 8/9)
Q' = -(4)/3 * (1/3, -2/3, -2/3) = (-4/9, 8/9, 8/9)

> >     .212  -.249       1 0 0
> > A = .691  -.640    I= 0 1 0
> >     .691   .717       0 0 1
> > 
    [ 2/3  2/3]
A = [-1/3  2/3]
    [ 2/3 -1/3]

     [2/sqrt(5) 2/sqrt(5)]
A' = [1/sqrt(5) 0        ]
     [0         1/sqrt(5)]

> > T = (AA^T -I)^-1 Q   =(1.085,6.552,6.607)=(t,t^2,t^3)

For the unprimed set of vectors, the expression T = (AA^T -I)^-1 Q
does not exist since (AA^T -I) is not invertible in this case.  This
is an EXTREME problem.  You can't even guarantee the equations
generated to be solvable!

For the primed set of vectors, T' = (4/9, -8/9, -8/9).  Now Jon, since
you did some rounding in your example, how would you round these to
the correct values of 2, 1.41421356237309504880168872421, and
-1.41421356237309504880168872421?

> > 
> > t = 1.085  throw out
> > t = 2.560~3
> > t = 1.876~2
> > 
> > 
> > EXAMPLE
> > 
> > n=8
> > 
> > T=[((AA^T-I)^-1)^T((AA^T-I)^-1)]^-1((AA^T-I)^-1)^T Q
> correct above to,
> T=[(AA^T-I)^T(AA^T-I)]^-1(AA^T-I)^T Q
> 
> > 
> > (Least Squares)
> > 
> >           LEAST SQUARES SOLUTION TO MATRIX               
> > k   a[k]                   SUMMARY           
> > 0  -30     t^k         |t^k|         t        ERROR   
> > 1    9   2.384E-01   2.384E-01   2.384E-01  -2.731E+01   
> > 2    9   2.031E+01   2.031E+01   4.507E+00   1.590E+06   
> > 3    1  -1.902E+01   1.902E+01   2.669E+00   2.557E+04   
> > 4    3  -4.038E+00   4.038E+00   1.418E+00   2.165E+02   
> > 5    6   4.742E+00   4.742E+00   9.485E-01   1.830E+00   
> > 6    1  -3.753E+01   3.753E+01   1.830E+00   1.416E+03   
> > 7    1   3.320E+01   3.320E+01   1.649E+00   6.553E+02   
> > 8    9   3.010E+01   3.010E+01   1.530E+00   3.787E+02   
> > 
This has to be the worst, YOUR ERROR VALUES ARE UP TO HALF A MILLION
TIMES LARGER THAN THE T VALUES!!!!!!!!!!  When an engineer would see
those kinds of error values, he would scream!

> > 
> > Jon Giffen
> > jon8338@peoplepc.com
> >
Technical Details
Message-ID: <12f59340.0407260451.3bbd755@posting.google.com>
Path: arc1-ca.pugleaf.net!g2news1.google.com!postnews2.google.com!not-for-mail
References: <8aGMc.12317$iK.7994@newsread2.news.atl.earthlink.net> <MSHMc.12907$f4.5254@newsread3.news.atl.earthlink.net>

Show Full Headers
Complete Article Headers
X-Google-Thread: 107079,5c2c35d267fdd702
X-Google-Thread: 109d8a,5c2c35d267fdd702
X-Google-Thread: f72df,5c2c35d267fdd702
X-Google-Thread: 103471,5c2c35d267fdd702
X-Google-Thread: 1139f8,5c2c35d267fdd702
X-Google-Thread: 10d67e,5c2c35d267fdd702
X-Google-Thread: 1130d6,5c2c35d267fdd702
X-Google-Thread: f6df2,5c2c35d267fdd702
X-Google-Attributes: gid107079,gid109d8a,gidf72df,gid103471,gid1139f8,gid10d67e,gid1130d6,gidf6df2,public
X-Google-Language: ENGLISH,ASCII-7-bit
Path: arc1-ca.pugleaf.net!g2news1.google.com!postnews2.google.com!not-for-mail
From: brianscsmith@yahoo.com (Brian Smith)
Newsgroups: sci.math.num-analysis,sci.math,japan.sci.math,han.sci.math,fr.sci.maths,aus.mathematics,alt.math.recreational,bln.sci.mathematik
Subject: Re: Solutions to Power Series
Date: 26 Jul 2004 05:51:24 -0700
Organization: http://groups.google.com
Lines: 123
Message-ID: <12f59340.0407260451.3bbd755@posting.google.com>
References: <8aGMc.12317$iK.7994@newsread2.news.atl.earthlink.net> <MSHMc.12907$f4.5254@newsread3.news.atl.earthlink.net>
NNTP-Posting-Host: 207.69.137.10
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: 8bit
X-Trace: posting.google.com 1090846284 12186 127.0.0.1 (26 Jul 2004 12:51:24 GMT)
X-Complaints-To: groups-abuse@google.com
NNTP-Posting-Date: Mon, 26 Jul 2004 12:51:24 +0000 (UTC)