Cím: Az Arnold-féle diszkrét macska-leképezés
Szerző(k):  Mincsovicsné Sélley Fanni 
Füzet: 2017/december, 514 - 520. oldal  PDF  |  MathML 

A szöveg csak Firefox böngészőben jelenik meg helyesen. Használja a fenti PDF file-ra mutató link-et a letöltésre.

 
1. Bevezetés

Legyenek x0 és y0 a 0 és 1 közötti számok, valamint
x1=x0+y0(mod1),y1=x0+2y0(mod1),(1)
ahol a mod 1 kifejezés törtrész vételét jelenti. Az 1. ábra szemlélteti, hogy mi történik az egységnégyzettel, ha minden pontjára végrehajtjuk ezt a leképezést: egy irányba megnyúlik, egy másik irányba összeszűkül, majd az így kapott parallelogrammát a mod 1 művelet visszadarabolja az egységnégyzetbe.


 

1. ábra. Az egységnégyzet képe az Arnold-féle macska-leképezés hatása alatt
 

Vlagyimir Arnold orosz matematikus után Arnold-féle macska-leképezésnek szokás ezt nevezni, mivel Arnold egy macska képével szemléltette a leképezés hatását. A leképezés érdekessége abból ered, hogy egyszerűsége ellenére erősen kaotikus. Ezalatt azt értjük, hogy ha ismételten végrehajtjuk a leképezést, az egységnégyzet pontjai gyorsan és alaposan megkeverednek. Ennek az az oka, hogy bármely pont egy irányban távolodik az eredeti ,,szomszédaitól'' (amerre nyúlik a négyzet), egy másik irányban pedig közeledik hozzájuk (amerre szűkül a négyzet).
De hogyan is szimulálhatta Arnold ezt a leképezést? A következő eljárás a kézenfekvő: tekintsünk egy N×N pixel méretű képet, és alkalmazzuk az (1) leképezést minden pixelre. A pixelek koordinátáit legkényelmesebb egész számokban megadni, azaz a következő leképezést alkalmazzuk:
xk+1=xk+yk(modN),yk+1=xk+2yk(modN),k=0,1,2,...,(2)
ahol xk és yk0 és N-1 közötti számok. Mivel a mod N kifejezés az N-nel való osztás maradékát veszi, xk+1 és yk+1 szintén 0 és N-1 közé esik.


 

2. ábra. Macskából káoszba (k= iterációk száma, N=400)
 

Kezdetben úgy tűnik, hogy a (2) leképezés teljesen összekeveri a képünk pontjait, ahogyan a folytonos változata is. Viszont némi számítógépes kísérletezés után meglepő módon azt tapasztaljuk, hogy kezdeti képünk előbb-utóbb újra megjelenik. De ha egy kicsit jobban belegondolunk, akkor láthatjuk, hogy maga a visszatérés ténye még nem igazán meglepő.
 
1. feladat. Tegyük fel, hogy a kép minden pixele csak fekete vagy fehér lehet. Mutassuk meg, hogy 2N2 iteráció alatt legalább egy visszatérést tapasztalunk.
 

Ami viszont meglepő, hogy közel sincs szükség ilyen sok iterációra. A 3. ábrán egy olyan esetet láthatunk, ahol kevesebb, mint N/2 iterációra van szükség a visszatéréshez. A következőkben áttekintünk pár egyszerűbb állítást a visszatérési időről.


 

3. ábra. Visszatér a macska (k= iterációk száma, N=241)
 

 
2. Visszatérési idő

Visszatérési időnek fogjuk nevezni azt az iterációszámot, amelyre először visszatér az eredeti képünk. Pontosabban, a visszatérési idő az a legkisebb mN szám, amelyre
xmN=x0,ymN=y0
teljesül az egységnégyzet valamennyi (x0,y0) pontjára.
Fogalmazzuk át ezt a definíciót.
 
2. feladat. Lássuk be, hogy
xk=u2k-1x0+u2ky0(modN),yk=u2kx0+u2k+1y0(modN),(3)
ahol ui az i-edik Fibonacci szám (azaz u0=0, u1=1, továbbá ui+1=ui+ui-1, i=0,1,...).
 

Tehát a legkisebb olyan k számot keressük, amelyre
u2k-11(modN),u2k0(modN)
teljesül (azaz u2k-1 N-nel való osztás után 1 maradékot ad, u2k pedig nullát). Ez a két feltétel elég, hiszen ezekből u2k-11(modN) következik. Úgy is mondhatjuk, hogy a visszatérési idő egyenlő a Fibonacci számok modulo N periódusának felével. Ezt a periódust szokás Pisano-periódusnak is nevezni. Sajnos zárt formula nem ismert rá, de az idevágó ismereteink jó összefoglalását adja D. D. Wall cikke [9].
A pontos képlet hiányának ellenére léteznek eredmények, amelyek felső korlátot adnak a visszatérési időre. Ezek az állítások egyszerű számelméleti eszközökkel, ám helyenként rendkívül aprólékos munkával bizonyíthatók. Az alábbi állítást Freeman J. Dyson és Harold Falk bizonyította:
 
1. tétel (Dyson‐Falk [5]). Legyen N>2. Ekkor a visszatérési időre fennáll az
mNN22
felső korlát.

 

Mielőtt rátérnénk a bizonyításra, megjegyezzük, hogy Freeman Dyson neves elméleti fizikus és matematikus, aki leginkább a kvantumelektrodinamika elméletének kidolgozásában vállalt szerepe miatt híres. Jelentőségét a matematikában a róla elnevezett Dyson-transzformált bizonyítja, amely az additív számelmélet egyik alapvető eszköze.
Lássuk most Dyson és Falk bizonyítását, amely a N. Vorobiev könyvében [8] leírt módszert követi.
 
Bizonyítás. Legyen Φk az uk Fibonacci-szám N-nel való osztási maradéka. Tekintsük a
Φ0,Φ1,Φ1,Φ2,...,Φk,Φk+1,...(4)
rendezett párokat. Vegyük észre, hogy Φ0,Φ1=0,1. Ebben a sorozatban legfeljebb N2 különböző elem lehet, tehát bármely N2+1 elem között szükségszerűen lesz legalább két megegyező. Most bebizonyítunk egy lemmát, hogy aztán tovább haladhassunk a tétel bizonyításával.
 
1. lemma. Az első ismétlődő páros a 0,1.
 
Bizonyítás. Indirekten fogunk érvelni. Tegyük fel, hogy az első ismétlődő pár Φk,Φk+1 valamely k>0 számra. Tekintsünk ekkor egy olyan Φr,Φr+1, r>k párost, amely megegyezik vele, azaz Φk=Φr és Φk+1=Φr+1. Ekkor a Fibonacci-számok definíciója alapján
Φr-1Φr+1-Φr(modN),Φk-1Φk+1-Φk(modN),
azaz Φr-1=Φk-1. Tehát
Φk-1,Φk=Φr-1,Φr,
azaz Φk-1,Φk is ismétlődik, és korábban van a (4) sorozatban, mint Φk,Φk+1 ‐ ellentmondásra jutottunk. Tehát k=0, és így Φk,Φk+1=Φ0,Φ1=0,1.  
 
3. feladat. Lássuk be az előző lemma segítségével a következőt: tetszőleges N számra igaz, hogy az (u0 utáni) első N2 Fibonacci-szám között lesz legalább egy N-nel osztható.
 
2. lemma. Legyen N>2. Ha uk0(modN) és uk+11(modN), akkor k páros.
 
Bizonyítás. Vezessük be a 
Dk=ukuk+2-uk+12
mennyiséget. D0=-1, valamint Dk+1=-Dk, hiszen
uk+1uk+3-uk+22=(uk+2-uk)(uk+1+uk+2)-uk+22==uk+2uk+1-ukuk+1-ukuk+2==-ukuk+2+(uk+2-uk)uk+1==-ukuk+2+uk+12.
Tehát Dk=-1, ha k páros, és Dk=1, ha páratlan. Amennyiben uk0(modN) és uk+11(modN), akkor Dk-1(modN) ‐ így k páros.  
 

Innen már könnyű a tétel bizonyítása: tekintsük a Φ0,Φ1,... sorozatot. A Φ0,Φ1,...,ΦN2+1 kezdőszeletben ismétlődni fog 0,1 ‐ legyen az első ismétlődés Φt,Φt+1. A 2. lemma alapján t páros, a visszatérési idő definíciója szerint pedig 2mN=t. Azaz mNN2/2.  
 

Egy általánosabb, de kevésbé explicit képletet ad a visszatérési időre Gregory Gaspari:
 
2. tétel (Gaspari [6]). Legyen N prímtényezős felbontása N=p1α1...pkαk, ahol pj prím, és αjN minden j=1,...,k esetében. Ekkor
mN=LKKT{mp1α1,...,mpkαk},
ahol LKKT a legkisebb közös többszöröst jelöli.

 

Megemlítjük, hogy a tétel Fibonacci-számok periódusára vonatkozó megfogalmazása már szerepelt D. D. Wall jóval korábbi cikkében [9]. A Gaspari által adott bizonyítás a következő.
 
Bizonyítás. Kezdjük a bizonyítást egy észrevétellel: azt állítjuk, hogy ha K osztja N-et, akkor mK osztja mN-et. Mivel
u2mN-11(modN),u2mN0(modN),
és K osztja N-et, azért u2mN-1-1K-val osztva is 1 maradékot ad, valamint u2mNK-val is osztható. Azaz
u2mN-11(modK),u2mN0(modK).(5)

Mivel mK a legkisebb szám, amire az (5) kongruencia-rendszer teljesül, mK kisebb (vagy egyenlő) mint mN. Tegyük fel, hogy mN az mK-val osztva nemnulla maradékot ad, azaz mN=qmK+r, ahol 0<r<mK, és q egész. Ekkor qmK iteráció alatt visszatér a képünk q-szor, és mivel mN iteráció alatt visszatér, r iteráció alatt is vissza kell térnie. De r<mK, és mK volt a legkisebb idő, ami alatt visszatér a kép, tehát ellentmondásra jutottunk. Azaz r=0, és ezzel beláttuk az észrevételt.
Térjünk rá a tétel bizonyítására. Tetszőleges j{1,...,k} esetében pjαj osztja N-et, tehát az előző észrevételünk miatt mpjαj osztja mN-et. Ebből következik, hogy mN közös többszöröse az mpjαj, j{1,...,k} számoknak. Legyen M egy közös többszöröse a mpjαj, j{1,...,k} számoknak. Ekkor
u2M-11(modp)jαj,u2M0(modp)jαj.
Mivel a pjαj számok különböző prímek hatványaiként páronként relatív prímek,
u2M-11(modp)1α1...pkαk=N,u2M0(modp)1α1...pkαk=N.
Viszont ez azt jelenti, hogy mN osztja M-et. Mivel mN a legkisebb egész szám, amire a fenti kongruencia teljesül, mNlegkisebb közös többszöröse az mpjαj számoknak.  
 

Tehát elég prímhatványokra tudni a visszatérési időt, ebből már tetszőleges számra kiszámítható. De ez még így sem egyszerű feladat. Gaspari cikkének a függelékében m=1,...,195 periódusokra kigyűjtötte az összes p prímet, amelyre mp=m ‐ ez nyújthat némi segítséget.
 
4. feladat. Lássuk be, hogy
a) m241=120 (3. ábra),
b) m400=300 (4. ábra).


 

4. ábra. N=400, mN=300
 

 
4*. feladat. Tegyük fel, hogy mN páros. Milyen N-ek esetében fordul elő hogy mN/2 iteráció után fejtetőn jelenik meg a macska, ahogyan a 3. ábrán is látható? A 4. ábra mutatja, hogy nem mindig ez a helyzet. Itt úgynevezett szellemképek jelennek meg, a magyarázatért lásd a [3] hivatkozást.
 
3. Alkalmazások

Bár első ránézésre Arnold macska-leképezése csak egy matematikai játéknak tűnik, a kaotikusságát kihasználva praktikus alkalmazásai is lehetnek. A következő gyűjtés a [7] hivatkozásra támaszkodik.
A legnyilvánvalóbb egy kép vagy szöveg titkosítása: a kép pixeleire, vagy a szöveg N×N-es blokkba rendezett karaktereire alkalmazzuk a macska-leképezés egy megfelelő hatványát. Így alaposan megkeverednek a képpontok (vagy a betűk), avatatlan szemlélő nem képes az eredeti üzenetet visszafejteni. Tovább bonyolítható a helyzet, ha a titkosító egy általánosabb macska-leképezést használ, például az
xk+1=xk+ayk(modN),yk+1=axk+(a2+1)yk(modN),k=0,1,2,...(6)
vagy az
xk+1=xk+ayk(modN),yk+1=bxk+(ab+1)yk(modN),k=0,1,2,...(7)
iterációt. Így az eredeti üzenet csakis az a,b egész számok ismeretében fejthető vissza. Ezeknek a macska-leképezéseknek a viselkedése teljesen hasonló Arnold macskájához. A visszatérési időkről J. Bao és Q. Yang ír a [2] cikkben.
Egy kicsit izgalmasabb alkalmazás a szteganográfiához köthető. A szteganográfia olyan titkos üzenetek létrehozásának tudománya, amelyek létezéséről a feladón kívül csak a címzett tud ‐ szemben a kriptográfiával, ahol az üzenet léte nem rejtély, csak a tartalma. A két módszer együttes alkalmazása természetesen a leghatékonyabb. Tegyük fel, hogy van egy képünk, amiről később majd meg akarjuk állapítani, hogy valaki manipulálta-e. A következő a módszerünk: alkalmazzuk a képünkre a macska-leképezés k darab iterációját, majd helyezzünk rá egy kis vízjelet. Ezután mN-k iterációval állítsuk vissza az eredeti képet, amelyen a vízjel egy hétköznapi szemlélőnek láthatatlan, hiszen a pixelei alaposan szétszóródtak. Ha később meg akarjuk állapítani, hogy valaki módosította-e a képet, elég a macska-leképezés k iterációját alkalmazni rá: ha a kép nem lett manipulálva, akkor bár a kép káoszba fullad, a vízjel eredeti állapotában megjelenik a sarokban. A részletek a [4] cikkben találhatók.
 
Hivatkozások


[1]V.I. Arnold and A. Avez, Ergodic problems of classical mechanics, W.A. Benjamin, New York (1968).
[2]Jianghong Bao and Qigui Yang, Period of the discrete arnold cat map and general cat map, Nonlinear Dynamics, 70(2) (2012), 1365‐1375.
[3]Ehrhard Behrends, The ghosts of the cat, Ergodic Theory and Dynamical Systems, 18(2) (1998), 321‐330.
[4]Young-Long Chen, Her-Terng Yau, and Guo-Jheng Yang, A maximum entropy-based chaotic time-variant fragile watermarking scheme for image tampering detection, Entropy, 15(8) (2013), 3170‐3185.
[5]Freeman J. Dyson and Harold Falk, Period of a discrete cat mapping, The American Mathematical Monthly, 99(7) (1992), 603‐614.
[6]Gregory Gaspari, The Arnold cat map on prime lattices, Physica D: Nonlinear Phenomena, 73(4) (1994), 352‐372.
[7]Fredrik Svanström, Properties of a generalized Arnold's discrete cat map (2014).
[8]Nicolai N. Vorobiev, Fibonacci numbers, Birkhäuser (2012).
[9]D.D. Wall, Fibonacci series modulo m, The American Mathematical Monthly, 67(6) (1960), 525‐532.