Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

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

Search the Archive

matrix switching Markov sequence

  • To: mathgroup at smc.vnet.net
  • Subject: [mg51609] matrix switching Markov sequence
  • From: Roger Bagula <tftn at earthlink.net>
  • Date: Wed, 27 Oct 2004 01:54:20 -0400 (EDT)
  • Reply-to: tftn at earthlink.net
  • Sender: owner-wri-mathgroup at wolfram.com

These four matrices give minimal Pisot type eigenvalues.
I call them pseudopermutations because they act
very much like permutations and are Determinant one.
They don't really form a matrix multiplication group, but are an
SL(3,R) set of matrices.
The second version below is the vector array version set up to give a 
single
output integer.
I had considered them as a marriage group type that would take a
{1/4,1/2,1/4} type gene array of four types of marriage types and 
permutate it
so that lethal crosses wouldn't occur. Most such marriage groups are 
based on 4by4
marriage group permutations not 3by3. I've exchanged the matrices for 
the types
and made the operations on the genetic types.
I have no way of proving
they will work better than the ones the primitive societies used!

(* all real 3by3 Markov sequence of  a four element pseudopermutation *)
digits=25
(* each of these matrices gives a Minimal Pisot type of eigenvalue*)
s1={{0,1,0},{1,0,1},{1,0,0}}
s2={{0,0,1},{1,0,1},{0,1,0}}
s3={{0,1,0},{0,0,1},{1,1,0}}
s4={{0,1,1},{1,0,0},{0,1,0}}
s={s1,s2,s3,s4}
A[n_]:=s[[1+Mod[n,4]]].A[n-1];
A[0]:=s[[1]];
(* sequence of 3by3 matrices made with  a four element pseudopermutation *)
b=Flatten[Table[A[n],{n,0,digits}]]

{0,1,0,1,0,1,1,0,0,1,0,0,1,1,0,1,0,1,1,1,0,1,0,1,2,1,0,3,1,1,1,1,0,1,0,1,1,1,
  
0,4,1,2,3,1,1,3,1,1,4,2,1,4,1,2,4,2,1,4,1,2,7,3,2,11,4,4,4,2,1,4,1,2,4,2,1,
  
15,5,6,11,4,4,11,4,4,15,6,5,15,5,6,15,6,5,15,5,6,26,10,9,41,15,15,15,6,5,15,
  5,6,15,6,5,56,20,21,41,15,15,41,15,15,56,21,20,56,20,21,56,21,20,56,20,21,
  
97,36,35,153,56,56,56,21,20,56,20,21,56,21,20,209,76,77,153,56,56,153,56,56,
  209,77,76,209,76,77,209,77,76,209,76,77,362,133,132,571,209,209,209,77,76,
  
209,76,77,209,77,76,780,285,286,571,209,209,571,209,209,780,286,285,780,285,
  286,780,286,285,780,285,286,1351,495,494,2131,780,780,780,286,285,780,285,
  286,780,286,285,2911,1065,1066,2131,780,780,2131,780,780,2911,1066,1065,
  2911,1065,1066}

(* all real 3by3 Markov sequence of  a four element pseudopermutation *)
digits=50
(* each of these matrices gives a Minimal Pisot type of eigenvalue*)
s1={{0,1,0},{1,0,1},{1,0,0}}
s2={{0,0,1},{1,0,1},{0,1,0}}
s3={{0,1,0},{0,0,1},{1,1,0}}
s4={{0,1,1},{1,0,0},{0,1,0}}
s={s1,s2,s3,s4}
A[n_]:=s[[1+Mod[n,4]]].A[n-1];
A[0]:={0,1,1};
(* sequence of 3 arrays made with  a four element pseudopermutation *
b=Table[A[n][[1]],{n,0,digits}]

{{0,1,1,3,1,3,4,11,4,11,15,41,15,41,56,153,56,153,209,571,209,571,780,2131,780,
  
2131,2911,7953,2911,7953,10864,29681,10864,29681,40545,110771,40545,110771,
  151316,413403,151316,413403,564719,1542841,564719,1542841,2107560,5757961,
  2107560,5757961,7865521}
Respectfully, Roger L. Bagula

tftn at earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
alternative email: rlbtftn at netscape.net
URL :  http://home.earthlink.net/~tftn


  • Prev by Date: Re: Integration of UnitStep has bugs!? help!
  • Next by Date: Re: Fitting question
  • Previous by thread: Elliptic curve over reals
  • Next by thread: MovieParametericPlot-3D