Research article Special Issues

Fractal approximation of chaos game representations using recurrent iterated function systems

  • Received: 01 October 2018 Accepted: 19 March 2019 Published: 20 January 2019
  • MSC : 37H99, 92D20

  • We demonstrate that chaos game representations of Cannabis sativa may be approximated by the chaos game approximation of a recurrent iterated function system attractor. Via numerical experiments, we then study the fractal scaling properties of both objects and apply a wavelet decomposition in order to investigate scale-invariant patterns. We show that the attractor of a recurrent iterated function system scales similarly to the chaos game representation and has a similar wavelet multiresolution analysis profile.

    Citation: Martin Do Pham. Fractal approximation of chaos game representations using recurrent iterated function systems[J]. AIMS Mathematics, 2019, 5(6): 1824-1840. doi: 10.3934/math.2019.6.1824

    Related Papers:

    [1] Mehmet Kunt, Artion Kashuri, Tingsong Du, Abdul Wakil Baidar . Quantum Montgomery identity and quantum estimates of Ostrowski type inequalities. AIMS Mathematics, 2020, 5(6): 5439-5457. doi: 10.3934/math.2020349
    [2] Marwa M. Tharwat, Marwa M. Ahmed, Ammara Nosheen, Khuram Ali Khan, Iram Shahzadi, Dumitru Baleanu, Ahmed A. El-Deeb . Dynamic inequalities of Grüss, Ostrowski and Trapezoid type via diamond-α integrals and Montgomery identity. AIMS Mathematics, 2024, 9(5): 12778-12799. doi: 10.3934/math.2024624
    [3] Atiqe Ur Rahman, Khuram Ali Khan, Ammara Nosheen, Muhammad Saeed, Thongchai Botmart, Nehad Ali Shah . Weighted Ostrowski type inequalities via Montgomery identity involving double integrals on time scales. AIMS Mathematics, 2022, 7(9): 16657-16672. doi: 10.3934/math.2022913
    [4] Bandar Bin-Mohsin, Muhammad Uzair Awan, Muhammad Zakria Javed, Sadia Talib, Hüseyin Budak, Muhammad Aslam Noor, Khalida Inayat Noor . On some classical integral inequalities in the setting of new post quantum integrals. AIMS Mathematics, 2023, 8(1): 1995-2017. doi: 10.3934/math.2023103
    [5] Humaira Kalsoom, Muhammad Amer Latif, Muhammad Idrees, Muhammad Arif, Zabidin Salleh . Quantum Hermite-Hadamard type inequalities for generalized strongly preinvex functions. AIMS Mathematics, 2021, 6(12): 13291-13310. doi: 10.3934/math.2021769
    [6] Da Shi, Ghulam Farid, Abd Elmotaleb A. M. A. Elamin, Wajida Akram, Abdullah A. Alahmari, B. A. Younis . Generalizations of some q-integral inequalities of Hölder, Ostrowski and Grüss type. AIMS Mathematics, 2023, 8(10): 23459-23471. doi: 10.3934/math.20231192
    [7] Humaira Kalsoom, Muhammad Idrees, Artion Kashuri, Muhammad Uzair Awan, Yu-Ming Chu . Some New (p1p2,q1q2)-Estimates of Ostrowski-type integral inequalities via n-polynomials s-type convexity. AIMS Mathematics, 2020, 5(6): 7122-7144. doi: 10.3934/math.2020456
    [8] Muhammad Amer Latif, Mehmet Kunt, Sever Silvestru Dragomir, İmdat İşcan . Post-quantum trapezoid type inequalities. AIMS Mathematics, 2020, 5(4): 4011-4026. doi: 10.3934/math.2020258
    [9] Mustafa Gürbüz, Yakup Taşdan, Erhan Set . Ostrowski type inequalities via the Katugampola fractional integrals. AIMS Mathematics, 2020, 5(1): 42-53. doi: 10.3934/math.2020004
    [10] Miguel Vivas-Cortez, Muhammad Aamir Ali, Hüseyin Budak, Ifra Bashir Sial . Post-quantum Ostrowski type integral inequalities for functions of two variables. AIMS Mathematics, 2022, 7(5): 8035-8063. doi: 10.3934/math.2022448
  • We demonstrate that chaos game representations of Cannabis sativa may be approximated by the chaos game approximation of a recurrent iterated function system attractor. Via numerical experiments, we then study the fractal scaling properties of both objects and apply a wavelet decomposition in order to investigate scale-invariant patterns. We show that the attractor of a recurrent iterated function system scales similarly to the chaos game representation and has a similar wavelet multiresolution analysis profile.


    In [4] the authors obtained the following generalization of Montgomery identity for quantum calculus.

    Lemma 1. [4] (Quantum Montgomery identity) Let f:[a,b]R, be an arbitrary function with daqf quantum integrable on [a,b], then the following quantum identity holds:

    f(x)1babaf(t)daqt=(ba)10Kq,x(t)Daqf(tb+(1t)a)d0qt (1.1)

    where Kq,x(t) is defined by

    Kq,x(t)={qt,0txaba,qt1,xaba<t1. (1.2)

    Using this identity, the authors have obtained two Ostrowski type inequalities for quantum integrals and applied it in several special cases.

    Unfortunately, in the proof of this lemma an error is made when calculating the integrals involving the kernel Kq,x(t) on the interval [xaba,1]. Also, in the proofs of Theorem 3 and Theorem 4 a small mistake related to the convexity of |Daqf|r is made.

    In the present paper we prove that the identity (1.1) and, thus, all of the consequent results are incorrect and provide corrections for these results.

    The q-derivative of a function f:[a,b]R for q0,1 (see [5] or [2] for a=0) is given by

    Daqf(x)=f(x)f(a+q(xa))(1q)(xa),forxa,b]Daqf(a)=limxaDaqf(x)

    We say that f:[a,b]R is q-differentiable if limxaDaqf(x) exists. The q-derivative is a discretization of the ordinary derivative and if f is a differentiable function then ([1,3])

    limq1 Daqf(x)=f(x).

    Further, the q-integral of f is defined by

    xaf(t)daqt=(1q)(xa)k=0qkf(a+qk(xa)), x[a,b].

    If the series on the right hand-side is convergent, then the q-integral xaf(t)daqt exists and f:[a,b]R is said to be q-integrable on [a,x]. If f is continuous on [a,b] the series (1q)(xa)k=0qkf(a+qk(xa)) tends to the Riemann integral of f as q1 ([1], [3])

    limq1xaf(t)daqt=xaf(t)dt.

    If ca,x the q-integral is defined by

    xcf(t)daqt=xaf(t)daqtcaf(t)daqt.

    Obviously, the q-integral depends on the values of f at the points outside the interval of integration and an important difference between the definite q-integral and Riemann integral is that even if we are integrating a function over the interval [c,x], a<c<x<b, for q-integral we have to take into account its behavior at t=a as well as its values on [a,x]. This is the main reason for mistakes made in [4] since in the proof of Lemma 1 the following error was made:

    1xabaKq,x(t)Daqf(tb+(1t)a)d0qt=10(qt1)Daqf(tb+(1t)a)d0qtxaba0(qt1)Daqf(tb+(1t)a)d0qt.

    But Kq,x(t)(qt1) for t[0,1] or for t[0,xaba], so the equality does not hold.

    Now, we give a proof that the quantum Montgomery identity (1.1) is not correct, since it does not hold for all x[a,b]. As we shall see, the identity (1.1) is valid only if x=a+qm+1(ba) for some mN{0}. We have

    (ba)10Kq,x(t)Daqf(tb+(1t)a)d0qt=(ba)(1q)k=0qkKq,x(qk)Daqf(a+qk(ba)).

    For q0,1 let mN{0} be such that

    qm+1xaba<qm,

    in other words

    m=logqxaba1.

    Then

    Kq,x(qk)={qk+11,km,qk+1,km+1,

    and

    (ba)(1q)k=0qkKq,x(qk)Daqf(a+qk(ba))=(ba)(1q)(mk=0qk(qk+11)f(a+qk(ba))f(a+qk+1(ba))(1q)qk(ba)+k=m+1qk(qk+1)f(a+qk(ba))f(a+qk+1(ba))(1q)qk(ba))=mk=0(f(a+qk(ba))f(a+qk+1(ba)))+k=0qk+1(f(a+qk(ba))f(a+qk+1(ba)))=f(a+qm+1(ba))f(b)+k=0qk+1(f(a+qk(ba))f(a+qk+1(ba))).

    If we put S=k=0qkf(a+qk(ba))=1(1q)(ba)baf(t)daqt, we have

    k=0(qk+1)(f(a+qk(ba))f(a+qk+1(ba)))=qS(Sf(b))

    and

    1babaf(t)daqt+(ba)10Kq,x(t)Daqf(tb+(1t)a)d0qt=1babaf(t)daqt+(f(a+qm+1(ba))f(b))+qS(Sf(b))=(1q)S+f(a+qm+1(ba))f(b)+qSS+f(b)=f(a+qm+1(ba))

    which is obviously not equal to f(x), unless x=a+qm+1(ba).

    This is no surprise since Jackson integral takes into account only f(a+qk(xa)) for kN{0}. Thus, we have proved the next lemma which is a corrected version of Lemma 1 from [4].

    Lemma 2. (Quantum Montgomery identity) Let f:[a,b]R, be an arbitrary function with Daqf quantum integrable on [a,b], then for all xa,b the following quantum identity holds:

    f(a+qlogqxaba(ba))1babaf(t)daqt=(ba)10Kq,x(t)Daqf(tb+(1t)a)d0qt

    where Kq,x(t) is defined by

    Kq,x(t)={qt,0txaba,qt1,xaba<t1.

    In Theorem 3 and Theorem 4 from [4] the authors have used the identity (1.1) to derive Ostrowski type inequalities for functions f for which Daqf is quantum integrable on [a,b] and |Daqf|r, r1 is a convex function. Since these inequalities depends on the validity of Lemma 1, our discussion invalidates all the results from [4].

    More precisely, in all the inequalities an additional assumption x=a+qm(ba) for some mN{0} should be added. In Theorems 3 and 4 |Daqf(a)|r and |Daqf(b)|r should be swapped, since in the proofs of Theorem 3 and Theorem 4, when applying the convexity of |Daqf|r the following mistake was made

    |Daqf(tb+(1t)a)|rt|Daqf(a)|r+(1t)|Daqf(b)|r.

    Lastly, the integral K4(a,b,x,q) is incorrectly computed and should read:

    K4(a,b,x,q)=1q1+q(bxba)+q1+q(bxba)2.

    The main goal of this paper was to point out that some results in [4] are not correct. We have concentrated on Lemma 3 (Quantum Montgomery identity). The statement of that Lemma is not correct as we have shown. We also found and analyzed the mistake in the proof of Lemma 3.

    However, we went one step further and stated and proved the correct version of Lemma 3 (it is Lemma 2 in our paper). We have also explained how can all inequalities derived from Quantum Montgomery identity be corrected.

    Domagoj Kovačević was supported by the QuantiXLie Centre of Excellence, a project co financed by the Croatian Government and European Union through the European Regional Development Fund-the Competitiveness and Cohesion Operational Programme (Grant KK.01.1.1.01.0004).

    The authors declare that they have no competing interests.

    After our Correction was accepted we were contacted by the first author of [4], Professor Kunt, who suggested an alternate way to correct the results of [4].

    The incorrect version of Montgomery identity from [4]

    f(x)1babaf(t)daqt=(ba)10Kq,x(t)Daqf(tb+(1t)a)d0qt

    can be fixed in two ways: either by changing the left hand side or by changing the right hand side of this equation. In Lemma 2 we showed how to fix the identity by correcting the left hand side. This makes it easier to salvage the rest of results in [4], as all the results remain valid with the added assumption that x=a+qm(ba) for some mN{0}.

    Professor Kunt suggested correcting the right hand side of this equation to obtain the identity:

    f(x)1babaf(t)daqt=(ba)[xaba0qtDaqf(tb+(1t)a)d0qt+1xaba(qt1)Daqf(tb+(1t)a)d0qt].  ()

    By doing so, the proofs of all the remaining results have to be corrected as the bound used

    |1xaba(qt1)Daqf(tb+(1t)a)d0qt|1xaba|(qt1)Daqf(tb+(1t)a)|d0qt

    does not hold for q-integrals in general. This is discussed, for example, on page 12 in [1,Section 1.3.1,Remark (ii)].

    When xaba=qm or equivalently x=a+qm(ba) for some mN{0} the bound above does hold, which is why there is no need to change the rest of the results in [4] if one takes our approach. Nevertheless, we list below the results that can be obtained using identity (3.1). The results below are due to Professor Kunt.

    Theorem 3 in [4] should be as follows:

    Theorem 3. Let f:[a,b]R be an arbitrary function with Daqf is quantum integrable on [a,b]. If |Daqf|r, r1 is a convex function, then the following quantum integral inequality holds:

    |f(x)1babaf(t)daqt|(ba)[(11+q)11r[|Daqf(b)|r1(1+q)(1+q+q2)+|Daqf(a)|rq1+q+q2]1r+(xaba)[|Daqf(b)|r(xaba)11+q+|Daqf(a)|r(1(xaba)11+q)]1r] (3.1)

    for all x[a,b].

    Proof. Using convexity of |Daqf|r, we have that

    |Daqf(tb+(1t)a)|rt|Daqf(b)|r+(1t)|Daqf(a)|r. (3.2)

    By using (), quantum power mean inequality and (3.2), we have that

    |f(x)1babaf(t)daqt| (3.3)
    =(ba)|xaba0qtDaqf(tb+(1t)a)d0qt+1xaba(qt1)Daqf(tb+(1t)a)d0qt|=(ba)|10(qt1)Daqf(tb+(1t)a)d0qt+xaba0Daqf(tb+(1t)a)d0qt|(ba)|10(qt1)Daqf(tb+(1t)a)d0qt|+|xaba0Daqf(tb+(1t)a)d0qt|
    (ba)[10(1qt)|Daqf(tb+(1t)a)|d0qt+xaba0|Daqf(tb+(1t)a)|d0qt]
    (ba)[(101qtd0qt)11r(10(1qt)|Daqf(tb+(1t)a)|rd0qt)1r+(xaba0d0qt)11r(xaba0|Daqf(tb+(1t)a)|rd0qt)1r](ba)[(10(1qt)d0qt)11r×(|Daqf(b)|r10(1qt)td0q+|Daqf(a)|r10(1qt)(1t)d0qt)1r+(xaba0d0qt)11r×(|Daqf(b)|rxaba0td0qt+|Daqf(a)|rxaba0(1t)d0qt)1r]

    On the other hand, calculating the following quantum integrals we have

    10(1qt)d0qt=(1q)n=0qn(1qn+1)=(1q)[11qq1q2]=11+q, (3.4)
    10(1qt)td0qt=(1q)n=0qn[(1qn+1)qn]=(1q)[11q2q1q3]=11+qq1+q+q2=1(1+q)(1+q+q2), (3.5)
    10(1qt)(1t)d0qt=101qtd0qt10(1qt)td0qt=11+q1(1+q)(1+q+q2)=q1+q+q2, (3.6)
    xaba0d0qt=(1q)(xaba)n=0qn=xaba, (3.7)
    xaba0td0qt=(1q)(xaba)n=0qn(qn(xaba))=(xaba)211+q, (3.8)
    xaba0(1t)d0qt=xaba0d0qtxaba0td0qt=xaba(xaba)211+q=(xaba)[1(xaba)11+q]. (3.9)

    Using (3.4)–(3.9) in (3.3), we have (3.1).

    Theorem 4 in [4] should be as follows:

    Theorem 4. Let f:[a,b]R be an arbitrary function with Daqf is quantum integrable on [a,b]. If |Daqf|r, r>1 and 1r+1p=1 is convex function, then the following quantum integral inequality holds:

    |f(x)1babaf(t)daqt| (3.10)
    (ba)[(10(1qt)pd0qt)1p(|Daqf(b)|r11+q+|Daqf(a)|rq1+q)1r+(xaba)[|Daqf(b)|r(xaba)11+q+|Daqf(a)|r(1(xaba)11+q)]1r]

    for all x[a,b].

    Proof. By using () and quantum Hölder inequality, we have

    |f(x)1babaf(t)daqt|(ba)|10(qt1)Daqf(tb+(1t)a)d0qt|+|xaba0Daqf(tb+(1t)a)d0qt|(ba)[10(1qt)|Daqf(tb+(1t)a)|d0qt+xaba0|Daqf(tb+(1t)a)|d0qt](ba)[(10(1qt)pd0qt)1p(10|Daqf(tb+(1t)a)|rd0qt)1r+(xaba0d0qt)1p(xaba0|Daqf(tb+(1t)a)|rd0qt)1r](ba)[(10(1qt)pd0qt)1p(10[t|Daqf(b)|r+(1t)|Daqf(a)|r]d0qt)1r+(xaba0d0qt)1p(xaba0[t|Daqf(b)|r+(1t)|Daqf(a)|r]d0qt)1r](ba)[(10(1qt)pd0qt)1p(|Daqf(b)|r10td0qt+|Daqf(a)|r10(1t)d0qt)1r+(xaba0d0qt)1p(|Daqf(b)|rxaba0td0qt+|Daqf(a)|rxaba0(1t)d0qt)1r]=(ba)[(10(1qt)pd0qt)1p(|Daqf(b)|r11+q+|Daqf(a)|rq1+q)1r+(xaba)[|Daqf(b)|r(xaba)11+q+|Daqf(a)|r(1(xaba)11+q)]1r].

    We conclude this section by noting that the bounds obtained in the original paper [4] which, as we have previously shown, do hold with the added assumption x=a+qm(ba) for some mN{0}, are tighter than the bounds obtained above by Professor Kunt. Professor Kunt's bounds, however, hold for all x[a,b].



    [1] J. S. Almeida, J. A. Carrico, A. Maretzek, et al. Analysis of genomic sequences by Chaos Game Representation, Bioinformatics, 17 (2001), 429-437. doi: 10.1093/bioinformatics/17.5.429
    [2] M. F. Barnsley, Superfractals, 1st edition, Cambridge University Press, Cambridge, 2006.
    [3] M. F. Barnsley, J. H. Elton and D. P. Hardin, Recurrent iterated function systems, Constr. Approx., 5 (1989), 3-31. doi: 10.1007/BF01889596
    [4] M. F. Barnsley and S. Demko, Iterated function systems and the global construction of fractals, Proceedings of the Royal Society of London A, 399 (1985), 243-275. doi: 10.1098/rspa.1985.0057
    [5] P. J. Deschavanne, A. Giron, J. Vilain, et al. Genomic signature: characterization and classification of species assessed by chaos game representation of sequences, Mol. Biol. Evol., 16 (1999), 1391-1399. doi: 10.1093/oxfordjournals.molbev.a026048
    [6] A. Fiser, G. E. Tusnady, I. Simon, Chaos game representation of protein structures, J. Mol. Graph. Model., 12 (1994), 302-304. doi: 10.1016/0263-7855(94)80109-6
    [7] J. M. Gutierrez, M. A. Rodriguez, G. Abramson, Multifractal analysis of DNA sequences using a novel chaos-game representation, Physica A: Statistical Mechanics and its Applications, 300 (2001), 271-284. doi: 10.1016/S0378-4371(01)00333-8
    [8] J. C. Hart, Fractal Image Compression and Recurrent Iterated Function Systems, IEEE Comput. Graph., 16 (1996), 25-33. doi: 10.1109/38.511849
    [9] H. J. Jeffrey, Chaos game representation of gene structure, Nucleic Acids Research, 18 (1990), 2163-2170. doi: 10.1093/nar/18.8.2163
    [10] H. Jia-Jing and F. Wei-Juan, Wavelet-based multifractal analysis of DNA sequences by using chaosgame representation, Chinese Phys. B, 19 (2010), 10205.
    [11] L. Kari, K. A. Hill, A. S. Sayem, et al. Mapping the space of genomic signatures, PLOS ONE, 10 (2015), 119815.
    [12] S. G. Mallat, A theory for multiresolution signal decomposition: the wavelet representation, IEEE transactions on pattern analysis and machine intelligence, 11 (1989), 674-693. doi: 10.1109/34.192463
    [13] P. Mayukha, B. Satish, K. Srinivas, et al. Multifractal detrended cross-correlation analysis of coding and non-coding DNA sequences through chaos-game representation, Physica A: Statistical Mechanics and its Applications, 436 (2015), 596-603. doi: 10.1016/j.physa.2015.05.018
    [14] H. Oh, B. Seo, S. Lee, et al. Two complete chloroplast genome sequences of Cannabis sativa varieties, Mitochondrial DNA Part A: DNA mapping, sequencing, and analysis, 27 (2016), 2835-2837.
    [15] D. Vergara, K. H. White, K. G. Keepers, et al. The complete chloroplast genomes of Cannabis sativa and Humulus lupulus, Mitochondrial DNA Part A: DNA mapping, sequencing, and analysis, 27 (2016), 3793-3794.
    [16] Y. Wang, K. Hill, S. Singh, et al. The spectrum of genomic signatures: from dinucleotides to chaos game representation, Gene, 346 (2005), 173-185. doi: 10.1016/j.gene.2004.10.021
    [17] J-Y. Yang, Z-L. Peng, Y. Zu-Guo, et al. Prediction of protein structural classes by recurrence quantification analysis based on chaos game representation, J. Theor. Biol., 257 (2009), 618-626. doi: 10.1016/j.jtbi.2008.12.027
    [18] Y. Zu-Guo, V. Anh, K-S. Lau, Chaos game representation of protein sequences based on the detailed HP model and their multifractal and correlation analyses, J. Theor. Biol., 226 (2004), 341-348. doi: 10.1016/j.jtbi.2003.09.009
    [19] Y. Zu-Guo, X. Qian-Jun, S. Long, et al. Chaos game representation of functional protein sequences, and simulation and multifractal analysis of induced measures, Chinese Phys. B, 19 (2010), 68701.
  • This article has been cited by:

    1. Andrea Aglić Aljinović, Lana Horvat Dmitrović , Ana Žgaljić Keko , Cauchy-Schwarz inequality for shifted quantum integral operator, 2024, 51, 12236934, 106, 10.52846/ami.v51i1.1749
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5183) PDF downloads(332) Cited by(1)

Figures and Tables

Figures(12)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog