Loading [MathJax]/jax/element/mml/optable/Latin1Supplement.js
 

Modelling chemistry and biology after implantation of a drug-eluting stent. Part Ⅱ: Cell proliferation

  • The aim of a drug eluting stent is to prevent restenosis of arteries following percutaneous balloon angioplasty. A long term goal of research in this area is to use modelling to optimise the design of these stents to maximise their efficiency. A key obstacle to implementing this is the lack of a mathematical model of the biology of restenosis. Here we investigate whether mathematical models of cancer biology can be adapted to model the biology of restenosis and the effect of drug elution. We show that relatively simple, rate kinetic models give a good description of available data of restenosis in animal experiments, and its modification by drug elution.

    Citation: Adam Peddle, William Lee, Tuoi Vo. Modelling chemistry and biology after implantation of a drug-eluting stent. Part Ⅱ: Cell proliferation[J]. Mathematical Biosciences and Engineering, 2018, 15(5): 1117-1135. doi: 10.3934/mbe.2018050

    Related Papers:

    [1] Yiqin Bao, Qiang Zhao, Jie Sun, Wenbin Xu, Hongbing Lu . An edge cloud and Fibonacci-Diffie-Hellman encryption scheme for secure printer data transmission. Mathematical Biosciences and Engineering, 2024, 21(1): 96-115. doi: 10.3934/mbe.2024005
    [2] Takao Komatsu, Haotian Ying . The p-Frobenius and p-Sylvester numbers for Fibonacci and Lucas triplets. Mathematical Biosciences and Engineering, 2023, 20(2): 3455-3481. doi: 10.3934/mbe.2023162
    [3] Zhongxue Yang, Yiqin Bao, Yuan Liu, Qiang Zhao, Hao Zheng, Wenbin Xu . Lightweight blockchain fuzzy decision scheme through MQTT and Fibonacci for sustainable transport. Mathematical Biosciences and Engineering, 2022, 19(12): 11935-11956. doi: 10.3934/mbe.2022556
    [4] Yan Wang, Guichen Lu, Jiang Du . Calibration and prediction for the inexact SIR model. Mathematical Biosciences and Engineering, 2022, 19(3): 2800-2818. doi: 10.3934/mbe.2022128
    [5] Anh Quang Tran, Tien-Anh Nguyen, Van Tu Duong, Quang-Huy Tran, Duc Nghia Tran, Duc-Tan Tran . MRI Simulation-based evaluation of an efficient under-sampling approach. Mathematical Biosciences and Engineering, 2020, 17(4): 4048-4063. doi: 10.3934/mbe.2020224
    [6] Yong Tao, Fan Ren, Youdong Chen, Tianmiao Wang, Yu Zou, Chaoyong Chen, Shan Jiang . A method for robotic grasping based on improved Gaussian mixture model. Mathematical Biosciences and Engineering, 2020, 17(2): 1495-1510. doi: 10.3934/mbe.2020077
    [7] Zhipeng Zhao, Zhenyu Hu, Zhiyu Zhao, Xiaoyu Du, Tianfei Chen, Lijun Sun . Fault-tolerant Hamiltonian cycle strategy for fast node fault diagnosis based on PMC in data center networks. Mathematical Biosciences and Engineering, 2024, 21(2): 2121-2136. doi: 10.3934/mbe.2024093
    [8] Auwalu Saleh Mubarak, Zubaida Said Ameen, Fadi Al-Turjman . Effect of Gaussian filtered images on Mask RCNN in detection and segmentation of potholes in smart cities. Mathematical Biosciences and Engineering, 2023, 20(1): 283-295. doi: 10.3934/mbe.2023013
    [9] Sergio Bermudo, Robinson A. Higuita, Juan Rada . k-domination and total k-domination numbers in catacondensed hexagonal systems. Mathematical Biosciences and Engineering, 2022, 19(7): 7138-7155. doi: 10.3934/mbe.2022337
    [10] Lal Hussain, Wajid Aziz, Ishtiaq Rasool Khan, Monagi H. Alkinani, Jalal S. Alowibdi . Machine learning based congestive heart failure detection using feature importance ranking of multimodal features. Mathematical Biosciences and Engineering, 2021, 18(1): 69-91. doi: 10.3934/mbe.2021004
  • The aim of a drug eluting stent is to prevent restenosis of arteries following percutaneous balloon angioplasty. A long term goal of research in this area is to use modelling to optimise the design of these stents to maximise their efficiency. A key obstacle to implementing this is the lack of a mathematical model of the biology of restenosis. Here we investigate whether mathematical models of cancer biology can be adapted to model the biology of restenosis and the effect of drug elution. We show that relatively simple, rate kinetic models give a good description of available data of restenosis in animal experiments, and its modification by drug elution.


    The Fibonacci sequence is defined Fn=Fn1+Fn2 for n2 with the initial conditions F0=0,F1=1. The Fibonacci sequence and the golden ratio are used in many fields such as cryptology, coding theory and quantum physics in [1,2,3,4,5,6,7,8,9,10,11]. Horadam defined the Gaussian Fibonacci numbers in [12,13] and gave some general identities about the Gaussian Fibonacci numbers. Jordan generalized the Gaussian Fibonacci numbers with a similar definition in [14]. The Gaussian Fibonacci sequence is defined as GFn=GFn1+GFn2 with n>1 where GF0=i,GF1=1 in [14]. It can also be easily seen Fn that GFn=Fn+iFn1 where Fn is the nth Fibonacci number.

    Asci and Gurel generalized these studies and defined the k-order Gaussian Fibonacci numbers in [15] by the following recurrence relation

    GF(k)n=kj=1GF(k)nj,forn>0andk2

    with boundary conditions for 1kn0,

    GF(k)n={1i,k=1ni,k=2n0,otherwise.

    It can also be seen that GF(k)n=F(k)n+iF(k)n1 where F(k)n is the nth k-order Fibonacci number.

    In [16], Asci and Aydinyuz defined the k-order Gaussian Fibonacci polynomials and gave some important results. The k-order Gaussian Fibonacci polynomials {GFn(k)}n=0 are

    GF(k)n(x)=kj=1xkjGF(k)nj(x)

    for n>0 and k2 with the initial conditions for 1kn0,

    GF(k)n(x)={1ix,k=1ni,k=2n0,otherwise.

    It can be seen that

    GF(k)n(x)=F(k)n(x)+iF(k)n1(x)

    where F(k)n(x) is the nth k-order Fibonacci polynomial.

    In order to ensure information security in terms of data transfer over communication channels, a lot of work has been done on this subject and continues to be done. Therefore, coding/decoding algorithms play an important role to ensure information security. Especially, Fibonacci coding theory is one of the most preferred in this field. We can see examples of these in many studies. For example, in [17], Stakhov gave a new approach to a coding theory using the generalization of the Cassini formula for Fibonacci p-numbers and Qp matrices in 2006. In 2009, Basu and Prasad in [18] presented the generalized relations among the code elements for the Fibonacci coding theory. Also, Basu and Das introduced a new coding theory for Tribonacci matrices in [19] and in 2014 they defined the coding theory for Fibonacci n-step numbers by generalizing the coding theory for Tribonacci numbers in [20]. Esmaeili in [21] also described a Fibonacci polynomial-based coding method with error detection and correction.

    In this paper, we remember the coding theory for k-order Gaussian Fibonacci polynomials given in [16] by taking x=1 for k-order Gaussian Fibonacci numbers. We describe the coding theory for the k-order Gaussian Fibonacci numbers. We give illustrative examples. The interesting relation between the elements of the code-message matrix is analyzed and derived. The most important of this paper is that we deal with error detection and error correction methods and examples of error probability is given.

    When we take x=1 in the coding theory defined in [16], we obtain the k-order Gaussian Fibonacci coding method. First, we define the Qk, Rk and En(k) matrices, which plays an important role in this coding theory. Qk, Rk and En(k) are defined the k×k matrices as the following:

    Qk=[111111000001000000000010]k×k,Rk=[GF(k)k1GF(k)k2GF(k)k3GF(k)2GF(k)10GF(k)k2GF(k)k3GF(k)k4GF(k)100GF(k)k3GF(k)k4GF(k)k5000GF(k)2GF(k)1000GF(k)10000i0000i1i]k×k

    and

    E(k)n=[GF(k)n+k1GF(k)n+k2GF(k)n+1GF(k)nGF(k)n+k2GF(k)n+k3GF(k)nGF(k)n1GF(k)n+1GF(k)nGF(k)nk+3GF(k)nk+2GF(k)nGF(k)n1GF(k)nk+2GF(k)nk+1]k×k

    where GF(k)n is the nth k-order Gaussian Fibonacci number.

    Theorem 2.1. For n1, we get in [15] as follow:

    QnkRk=E(k)n. (2.1)

    Corollary 2.1. For k=2, we get

    QnR=[1110]n[1ii1i]
    =[GFn+1GFnGFnGFn1]

    where GFn is the nth usual Gaussian Fibonacci number in [15].

    Corollary 2.2. For k=3, we get

    Qn3R3=[111100010]n[1+i1010i0i1i]
    =[GTn+2GTn+1GTnGTn+1GTnGTn1GTnGTn1GTn2]

    where GTn is the nth Gaussian Tribonacci number in [23].

    In this section, we redefine k-order Gaussian Fibonacci coding theory using k-order Gaussian Fibonacci numbers and Qk, Rk and En(k) play very important role in the construction of k-order Gaussian Fibonacci coding theory.

    We now obtain the matrix En(k) using Qk, Rk matrices for the k=2,k=3 values and examine the inverse.

    • For k=2, introducing the square matrix Q2 of order 2 as:

    Q2=[1110]

    and the square matrix R2 of order 2 as:

    R2=[1ii1i]

    for n=1, we can use (2.1)

    E(2)1=Q2.R2
    =[1+i11i]
    =[GF(2)2GF(2)1GF(2)1GF(2)0]

    such that

    detE(2)1=det(Q2.R2)
    =(1).(2i)
    =2+i.

    The inverse of E(2)1 is as:

    (E(2)1)1=[1525i25+15i25+15i1535i]
    =1detE(2)1[GF0GF1GF1GF2]
    =12+i[GF0GF1GF1GF2]

    such that

    det(E(2)1)1=12+i
    =2515i.

    Also, by (2.1) for n=2, we can get E(2)2 as follows:

    E(2)2=Q22.R2
    =[2+i1+i1+i1]
    =[GF(2)3GF(2)2GF(2)2GF(2)1]

    such that

    detE(2)2=det(Q22.R2)
    =(1)2.(2i)
    =2i.

    The inverse of E(2)2 is as:

    (E(2)2)1=[25+15i1535i1535i35+45i]
    =1detE(2)2[GF(2)1GF(2)2GF(2)2GF(2)3]
    =12i[GF(2)1GF(2)2GF(2)2GF(2)3]

    such that

    det(E(2)2)1=12i
    =25+15i.

    Theorem 2.1.1. E(2)n=[GF(2)n+1GF(2)nGF(2)nGF(2)n1] where E(2)1=[1+i11i].

    Theorem 2.1.2. (E(2)n)1=1detE(2)n[GF(2)n1GF(2)nGF(2)nGF(2)n+1]

    =1(detQ2)n.detR2[GF(2)n1GF(2)nGF(2)nGF(2)n+1].

    • For k=3, introducing the square matrix Q3 of order 3 as:

    Q3=[111100010]

    and the square matrix R3 of order 3 as:

    R3=[1+i1010i0i1i]

    for n=1, we can use (2.1)

    E(3)1=Q3.R3
    =[2+i1+i11+i1010i]
    =[GF(3)3GF(3)2GF(3)1GF(3)2GF(3)1GF(3)0GF(3)1GF(3)0GF(3)1]

    such that

    detE(3)1=det(Q3.R3)
    =1.2i
    =2i.

    The inverse of E(3)1 is as:

    (E(3)1)1=[121212i12i1212i1+i1212i12i1212i12i]
    =1detE(3)1[GF(3)1GF1(3)(GF(3)0)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)0GF2(3)(GF(3)1)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)3GF1(3)(GF(3)1)2GF(3)1GF2(3)GF(3)0GF3(3)(GF(3)1)2GF(3)0GF2(3)GF(3)1GF2(3)GF(3)0GF3(3)GF(3)1GF3(3)(GF(3)2)2]
    =1detR3[GF(3)1GF1(3)(GF(3)0)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)0GF2(3)(GF(3)1)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)3GF1(3)(GF(3)1)2GF(3)1GF2(3)GF(3)0GF3(3)(GF(3)1)2GF(3)0GF2(3)GF(3)1GF2(3)GF(3)0GF3(3)GF(3)1GF3(3)(GF(3)2)2]

    such that

    det(E(3)1)1=12i
    =i2.

    Also, by (2.1) for n=2, we can get E(3)2 as follows:

    E(3)2=Q23.R3
    =[4+2i2+i1+i2+i1+i11+i10]
    =[GF(3)4GF(3)3GF(3)2GF(3)3GF(3)2GF(3)1GF(3)2GF(3)1GF(3)0]

    such that

    detE(3)2=det(Q23.R3)
    =12.2i
    =2i.

    The inverse of E(3)2 is as:

    (E(3)2)1=[12i1212i12i1212i112+32i12i12+32i1+12i]
    =1detE(3)2[GF(3)0GF(3)2(GF(3)1)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)1GF(3)3(GF(3)2)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)0GF(3)4(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)1GF(3)3(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)2GF(3)4(GF(3)3)2]
    =1detR3[GF(3)0GF(3)2(GF(3)1)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)1GF(3)3(GF(3)2)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)0GF(3)4(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)1GF(3)3(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)2GF(3)4(GF(3)3)2].

    Theorem 2.1.3. E(3)n=[GF(3)n+2GF(3)n+1GF(3)nGF(3)n+1GF(3)nGF(3)n1GF(3)nGF(3)n1GF(3)n2] where E(3)1=[2+i1+i11+i1010i].

    Theorem 2.1.4.

    (E(3)n)1=1detR3[GF(3)n2GF(3)n(GF(3)n1)2GF(3)n1GF(3)nGF(3)n2GF(3)n+1GF(3)n1GF(3)n+1(GF(3)n)2GF(3)n1GF(3)nGF(3)n2GF(3)n+1GF(3)n2GF(3)n+2(GF(3)n)2GF(3)nGF(3)n+1GF(3)n1GF(3)n+2GF(3)n1GF(3)n+1(GF(3)n)2GF(3)nGF(3)n+1GF(3)n1GF(3)n+2GF(3)nGF(3)n+2(GF(3)n+1)2].

    For arbitrary k-positive integers, the square matrix E(k)n of order k and inverses can be found similarly.

    In this section, we describe a new k-order Gaussian Fibonacci coding theory. We put our message in a matrix of M and let us represent the initial message in the form of the square matrix M of order k. We take the k-order Gaussian Fibonacci matrix E(k)n as a coding matrix and its inverse matrix (E(k)n)1 as a decoding matrix for an arbitrary positive integer k. The transformation M×E(k)n=C is called k-order Gaussian Fibonacci coding and we name the transformation C×(E(k)n)1=M as k-order Gaussian Fibonacci decoding. We define C as a code matrix.

    The given example is solved using the alphabet table below.

    Using the arbitrary value of s0, we write the following alphabet according to mod30. We can extend the characters in the table according to our wishes. We begin the “s” for the first character in Table 1.

    Table 1.  Alphabet table.
    A B C D E F G H I J
    s s+1 s+2 s+3 s+4 s+5 s+6 s+7 s+8 s+9
    K L M N O P Q R S T
    s+10 s+11 s+12 s+13 s+14 s+15 s+16 s+17 s+18 s+19
    U V W X Y Z 0 ! ? .
    s+20 s+21 s+22 s+23 s+24 s+25 s+26 s+27 s+28 s+29

     | Show Table
    DownLoad: CSV

    Example 2.3.1. Let us consider the message matrix for the following message text:

    {\rm{“CODE”}} (2.2)

    Step 1: Let’s create the message matrix using the message text:

    M = {\left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right]_{2 \times 2}}

    Step 2: Let’s write the message matrix M according to the alphabet table for the arbitrary value “s” we choose. For s = 2 ;

    M = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]

    Step 3: For k = 2 , n = 2 , we use (2.1);

    E_2^{(2)} = Q_2^2.{R_2}
    = \left[ {\begin{array}{*{20}{c}} {2 + i}&{1 + i} \\ {1 + i}&1 \end{array}} \right]

    Step 4: The code message is:

    C = M \times E_2^{(2)}
    = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {2 + i}&{1 + i} \\ {1 + i}&1 \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right]

    Step 5: The decode message is:

    M = C \times {\left( {E_2^{(2)}} \right)^{ - 1}}
    = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\frac{2}{5} + \frac{1}{5}i}&{ - \frac{1}{5} - \frac{3}{5}i} \\ { - \frac{1}{5} - \frac{3}{5}i}&{\frac{3}{5} + \frac{4}{5}i} \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right] .

    Example 2.3.2. Let us consider the message matrix for the following message text:

    {\rm{"PUBLIC \;KEY"}}

    Step 1: Let’s create the message matrix using the message text:

    M = {\left[ {\begin{array}{*{20}{c}} P&U&B \\ L&I&C \\ K&E&Y \end{array}} \right]_{3 \times 3}}

    Step 2: Let’s write the message matrix M according to the alphabet table for the arbitrary value “s” we choose. For s = 5 :

    M = \left[ {\begin{array}{*{20}{c}} {20}&{25}&6 \\ {16}&{13}&7 \\ {15}&9&{29} \end{array}} \right]

    Step 3: For k = 3 , n = 6 , we use (2.1):

    E_6^{(3)} = Q_3^6.{R_3}
    = \left[ {\begin{array}{*{20}{c}} {44 + 24i}&{24 + 13i}&{13 + 7i} \\ {24 + 13i}&{13 + 7i}&{7 + 4i} \\ {13 + 7i}&{7 + 4i}&{4 + 2i} \end{array}} \right]

    Step 4: The code message is:

    C = M \times E_6^{(3)}
    = \left[ {\begin{array}{*{20}{c}} {1558 + 847i}&{847 + 459i}&{459 + 252i} \\ {1107 + 602i}&{602 + 327i}&{327 + 178i} \\ {1253 + 680i}&{680 + 374i}&{374 + 199i} \end{array}} \right]

    Step 5: The decode message is:

    M = C \times {\left( {E_3^{(6)}} \right)^{ - 1}}
    = \left[ {\begin{array}{*{20}{c}} P&U&B \\ L&I&C \\ K&E&Y \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} {20}&{25}&6 \\ {16}&{13}&7 \\ {15}&9&{29} \end{array}} \right] .

    In this section, we consider the k-order Gaussian Fibonacci coding/decoding method for k = 2 . We have an interesting relation among the elements of a code matrix C that has a crucial role in the error-correction process, outlined as follows:

    C = M \times E_n^{(2)} = \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right] \times \left[ {\begin{array}{*{20}{c}} {GF_{n + 1}^{(2)}}&{GF_n^{(2)}} \\ {GF_n^{(2)}}&{GF_{n - 1}^{(2)}} \end{array}} \right] (3.1)
    = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] (3.2)

    and

    M = C \times {\left( {E_n^{(2)}} \right)^{ - 1}} (3.3)
    = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] \times \left( {\frac{1}{{{{\left( {\det \left( {{Q_2}} \right)} \right)}^n}\det {R_2}}}\left[ {\begin{array}{*{20}{c}} {GF_{n - 1}^{(2)}}&{ - GF_n^{(2)}} \\ { - GF_n^{(2)}}&{GF_{n + 1}^{(2)}} \end{array}} \right]} \right) (3.4)
    = \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right]

    For the case for an even integer n = 2m , we obtain the following equation

    \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] \times \left( {\frac{1}{{2 - i}}\left[ {\begin{array}{*{20}{c}} {GF_{n - 1}^{(2)}}&{ - GF_n^{(2)}} \\ { - GF_n^{(2)}}&{GF_{n + 1}^{(2)}} \end{array}} \right]} \right) (3.5)

    since the \det \left({{Q_2}} \right) = 1 and \det \left({{R_2}} \right) = 2 - i .

    It follows from (3.5) that the elements of the matrix M can be obtained according to the following formulas:

    {m_1} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_1} - GF_n^{(2)}{c_2}} \right) (3.6)
    {m_2} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_1} + GF_{n + 1}^{(2)}{c_2}} \right) (3.7)
    {m_3} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_3} - GF_n^{(2)}{c_4}} \right) (3.8)
    {m_4} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_3} + GF_{n + 1}^{(2)}{c_4}} \right) (3.9)

    Since the s \geqslant 0 , the elements of the matrix M are

    {m_1} \geqslant 0 , {m_2} \geqslant 0 , {m_3} \geqslant 0 \;{\rm{ and}}\; {m_4} \geqslant 0 (3.10)

    Because of the condition (3.10) we can write the equalities (3.6)–(3.9) as:

    {m_1} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_1} - GF_n^{(2)}{c_2}} \right) \geqslant 0 (3.11)
    {m_2} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_1} + GF_{n + 1}^{(2)}{c_2}} \right) \geqslant 0 (3.12)
    {m_3} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_3} - GF_n^{(2)}{c_4}} \right) \geqslant 0 (3.13)
    {m_4} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_3} + GF_{n + 1}^{(2)}{c_4}} \right) \geqslant 0 (3.14)

    From the Eqs (3.11) and (3.12), we can get

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}}{c_2} \leqslant {c_1} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}}{c_2} (3.15)

    or

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}} \leqslant \frac{{{c_1}}}{{{c_2}}} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}} (3.16)

    Similarly, from (3.13) and (3.14), we can obtain:

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}}{c_4} \leqslant {c_3} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}}{c_4} (3.17)

    or

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}} \leqslant \frac{{{c_3}}}{{{c_4}}} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}} (3.18)

    Since two consecutive Gaussian Fibonacci numbers approach the golden ratio, we also obtain the following equations in (3.16) and (3.18) that connect the elements of the code matrix in (3.1):

    {c_1} \approx \tau {c_2} (3.19)
    {c_3} \approx \tau {c_4} (3.20)

    where \tau = \frac{{1 + \sqrt 5 }}{2} is the “golden ratio”.

    Similarly, for an odd integer n = 2m + 1 , we can write the same approximate equalities (3.19) and (3.20) that connect the elements {c_1}, {c_2}, {c_3} and {c_4} of the code matrix in (3.1).

    Thus, we have generated some important identities that connect the elements of the code matrix in (3.1) for the case k = 2 .

    If we consider the coding/decoding method for k = 3 , we can obtain interesting identities as in k = 2 . Since two consecutive Gaussian Tribonacci numbers approach \Phi where \Phi = 1.8393. However, for general case of k we can find mathematical identities that connect the code matrix elements similar to the (3.19) and (3.20).

    Example 3.1. Suppose k = 2 . Thus, we have \tau = \frac{{1 + \sqrt 5 }}{2} \cong 1.618 . Assume that the following message-matrix is to be transmitted in (2.2):

    For s = 2 , we can get

    M = \left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]

    • If n = 2 then

    C = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right]

    In this case, we have the following numbers rounded off to their first five digits:

    \frac{{{c_1}}}{{{c_2}}} = 1.3462 + 0.7307i , \frac{{{c_3}}}{{{c_4}}} = 1.5822 + 0.2808i

    • For n = 5 , we have

    C = \left[ {\begin{array}{*{20}{c}} {112 + 68i}&{68 + 44i} \\ {70 + 43i}&{43 + 27i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.6171 - 0.0463i , \frac{{{c_3}}}{{{c_4}}} = 1.6179 - 0.0159i

    • For n = 10 , we have

    C = \left[ {\begin{array}{*{20}{c}} {1236 + 764i}&{764 + 472i} \\ {775 + 479i}&{479 + 296i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.618 + 0.0003i , \frac{{{c_3}}}{{{c_4}}} = 1.618 + 0.0001i

    • For n = 15 , we have

    C = \left[ {\begin{array}{*{20}{c}} {13708 + 8472i}&{8472 + 5236i} \\ {8595 + 5312i}&{5312 + 3283i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.618 , \frac{{{c_3}}}{{{c_4}}} = 1.618

    These show that for n = 15 the relation (3.19) and (3.20) holds very well.

    Therefore, we select the value of n large enough, as seen in Examples 2.3.1 we get a good approximation of \tau . If we apply similar operation in Example 2.3.2, when n is chosen large enough we get a good approximation of \Phi where \Phi = 1.8393 and \Phi is the ratio between two Gaussian Tribonacci numbers.

    The k-order Gaussian Fibonacci coding/decoding method considered above provides an interesting possibility to detect and correct “errors” in code message C . Error detection and correction for k = 2 depends on the states of the determinant of the matrix C and the connections with the code matrix elements at (3.19) and (3.20). These mathematical relations given by (3.19), (3.20) and the determinant of the code matrix C play a role in checking relations of the k-order Gaussian Fibonacci coding/decoding method. The error-correction algorithm for Fibonacci coding has been explained in [10] and [13]. This defined algorithm is valid and applicable to the k-order Gaussian Fibonacci numbers we have given in this study.

    First, let’s calculate the determinant of our message matrix M and send it to the communication channel right after the code matrix elements. Here, the determinant of M will be the control element for the C code matrix received from the communication channel. After taking the determinant of the code matrix C and its control element M , we can calculate the determinant of the code matrix C and compare the determinant of M with the determinant of C with respect to the control relation. If the operations we have done after this check are exactly appropriate, we can conclude that the elements of the code matrix C are transmitted on the communication channel without errors. If the opposite happens, we can say that there are errors in the elements of the code matrix C or the determinant of the control element M .

    Suppose the elements of the C code matrix errors. This matrix can have a one-fold error, double-fold error or {k^2} - fold error. To explain how to correct these errors, we consider the 2 \times 2 type matrix obtained for k = 2 . These cases are considered.

    Our first hypothesis is that we have the case of a “single” error in the code matrix C received from the communication channel. It is clearly seen that there are four types of single-fold errors in the C code matrix:

    (a) \left[ {\begin{array}{*{20}{c}} x&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] (b) \left[ {\begin{array}{*{20}{c}} {{c_1}}&y \\ {{c_3}}&{{c_4}} \end{array}} \right] (c) \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ z&{{c_4}} \end{array}} \right] (d) \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&t \end{array}} \right] (4.1)

    where x, y, z, t are possible “destroyed” elements.

    Since it is C = M \times E_n^{(2)} , if we take the determinant of both sides, we get the following equations:

    x{c_4} - {c_2}{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a\; possible \;“single \;error” \;is \;in \;the \;element}}\; {c_1} ) (4.2)
    {c_1}{c_4} - y{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a\; possible \;“single \;error” \;is \;in \;the\; element}}\; {c_2} ) (4.3)
    {c_1}{c_4} - {c_2}z = \det (M){( - 1)^n}(2 - i) ({\rm{a \;possible \;“single \;error”\; is \;in \;the \;element}}\; {c_3} ) (4.4)
    {c_1}t - {c_2}{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a \;possible \;“single \;error” \;is \;in \;the \;element}}\; {c_4} ) (4.5)

    or equivalently

    x = \frac{{\det (M){{( - 1)}^n}(2 - i) + {c_2}{c_3}}}{{{c_4}}} (a) (4.6)
    y = \frac{{ - \det (M){{( - 1)}^n}(2 - i) + {c_1}{c_4}}}{{{c_3}}} (b) (4.7)
    z = \frac{{ - \det (M){{( - 1)}^n}(2 - i) + {c_1}{c_4}}}{{{c_2}}} (c) (4.8)
    t = \frac{{\det (M){{( - 1)}^n}(2 - i) + {c_2}{c_3}}}{{{c_1}}} (d) (4.9)

    The formulas (4.6)–(4.9) give four possible states of a “single error”. However, we need to choose the correct case from among x, y, z, t integer solutions. We have to choice such solutions, which satisfies to the additional “checking relations” (3.19) and (3.20). If calculations by formulas (4.6)–(4.9) do not give an integer result, we can conclude that our “single error” hypothesis is incorrect or that the determinant of control element M is “error”. In the second case, we can use approximate Eqs (3.19) and (3.20) to check the correctness of the C code matrix.

    Similarly, we can easily check “double errors” of the C code matrix. Suppose that the code matrix C has elements x and y error as shown below:

    \left[ {\begin{array}{*{20}{c}} x&y \\ {{c_3}}&{{c_4}} \end{array}} \right] (4.10)

    Then,

    x{c_4} - y{c_3} = \det (M){( - 1)^n}(2 - i) (4.11)

    However, according to the (3.19) there is the following relation between x and y

    x \approx \tau y (4.12)

    Again, only integer solutions are acceptable. If no integer solution is obtained, then two-fold errors have not occurred.

    Hence, we can show that with this type of approach, there is the possibility to correct all possible triple-fold errors in the C code matrix.

    \left[ {\begin{array}{*{20}{c}} x&y \\ z&{{c_4}} \end{array}} \right] .

    Therefore, our method of error correction depends on the verification of different hypotheses regarding errors in the code matrix using the determinant of the control relationship C matrix, (3.19) and (3.20), and the elements of the code matrix as integers. If all our solutions do not reach integer solutions, it means that the determinant of the control item M is wrong or there is a four-fold error in the C code matrix, and we should reject it. The code matrix C is defective and not correctable.

    As a result, there are 15 error conditions in the C code matrix. According to the method given in [17], it means that the correctable probability of the method is equal since 14 cases between them can be corrected and

    {S_{cor}} = \frac{{14}}{{15}} = 0.9333 = 93.33\%

    If we generalize this equation as in [20], since only {k^2} - fold errors cannot be corrected with this method, the error correction capacity of the method is

    \frac{{{2^{{k^2}}} - 1}}{{{2^{{k^2}}}}}

    where k is the order of the message-matrix. So far sufficiently large values of k , the probability of decoding error is almost zero.

    Consequently, for a sufficiently large value of k , in the case of n = k , the correct ability of this method is

    \frac{{{2^{{k^2}}} - 2}}{{{2^{{k^2}}} - 1}} .

    Therefore, for large value of k , the correct possibility of the method is

    \frac{{{2^{{k^2}}} - 2}}{{{2^{{k^2}}} - 1}} \approx 1 = 100\% .

    In this paper, we obtained the coding theory for k-order Gaussian Fibonacci polynomials given in [16] by taking x = 1 for k-order Gaussian Fibonacci numbers and gave illustrative examples. This coding method differs from classical algebraic coding methods in this respect. We can show this difference with the following features:

    1)      Since the n, k and s \geqslant 0 values are arbitrarily chosen in the k-order Gaussian Fibonacci coding theory, it is difficult to estimate the information transmission between the two channels by a third channel, which increases the reliability of the information transmission. In addition, this method depends on matrix multiplication and can be performed quickly and easily by today’s computers.

    2)      With k-order Gaussian Fibonacci coding theory, we can encrypt and send messages of the desired length by enlarging the k value sufficiently.

    3)      The main practical feature of this method is that matrix elements are error detection and correction objects.

    4)      The Gaussian Fibonacci coding method, which is the simplest coding method obtained for the k = 2 value, has been redefined by handling the 2 \times 2 type matrices with “single-fold, double-fold and triple-fold errors”.

    5)      In the simplest case, for k = 2 , the correct capability of the method is essentially equal to 93.33%, exceeding all well-known correction codes.

    6)      This error correction method given for k = 2 is generalized and the correct ability of the errors of this method increases as k increases, and the correct ability approaches 100% for a sufficiently large value of k .

    7)      This article is just a brief outline of a new coding theory based on the k-order Gaussian Fibonacci matrices of the articles in [17,18,20,21] and [22].

    The authors declare there is no conflict of interest.

    [1] [ World Health Organisation, Global Status Report on Noncommunicable Diseases 2010, Report of the World Health Organisation, 2011. Available from: http://www.who.int/nmh/publications/ncd_report_full_en.pdf.
    [2] [ A. Bertuzzi, A. Fasano, A. Gandolfi and C. Sinisgalli, Tumour cords and their response to anticancer agents, in Modelling and Simulation in Science, Engineering, and Technology (ed. N. Bellomo), Birkhauser, (2008), 183–206.
    [3] [ F. Bozsak,J. M. Chomaz,A. I. Barakat, Modeling the transport of drugs eluted from stents: Physical phenomena driving drug distribution in the arterial wall, Biomechanics and Modeling in Mechanobiology, 13 (2014): 327-347.
    [4] [ F. Bozsak,D. Gonzalez-Rodriguez,Z. Sternberger,P. Belitz,T. Bewley,J. M. Chomaz,A. I. Barakat, Optimization of drug delivery by drug-eluting stents, PLoS One, 10 (2015): e0130182.
    [5] [ A. Caiazzo,D. Evans,J.-L. Falcone,J. Hegewald,E. Lorenz,B. Stahl,D. Wang,J. Bernsdorf,B. Chopard,J. Gunn,R. Hose,M. Krafczyk,P. Lawford,R. Smallwood,D. Walker,A. Hoekstra, A Complex Automata approach for in-stent restenosis: Two-dimensional multiscale modelling and simulations, Journal of Computational Science, 2 (2011): 9-17.
    [6] [ Y. S. Chatzizisis,A. U. Coskun,M. Jonas,E. R. Edelman,C. L. Feldman,P. H. Stone, Role of endothelial shear stress in the natural history of coronary atherosclerosis and vascular remodeling: molecular, cellular, and vascular behavior, Journal of the American College of Cardiology, 49 (2007): 2379-2393.
    [7] [ D. E. Drachman,E. R. Edelman,P. Seifert,A. R. Groothuis,D. A. Bornstein,K. R. Kamath,M. Palasis,D. Yang,S. H. Nott,C. Rogers, Neointimal thickening after stent delivery of paclitaxel: Change in composition and arrest of growth over six months, Journal of the American College of Cardiology, 36 (2000): 2325-2332.
    [8] [ J. A. Ferreira,J. Naghipoor,P. de Oliveira, Analytical and numerical study of a coupled cardiovascular drug delivery model, Journal of Computational and Applied Mathematics, 275 (2015): 433-446.
    [9] [ F. J. H. Gijsen,F. Migliavacca,S. Schievano,L. Socci,L. Petrini,A. Thury,J. J. Wentzel,A. F. W. van der Steen,P. W. S. Serruys,G. Dubini, Simulation of stent deployment in a realistic human coronary artery, BioMedical Engineering OnLine, 7 (2008): 1-23.
    [10] [ W. Khan,S. Farah,A. Nyska,A. J. Domb, Carrier free rapamycin loaded drug eluting stent: In vitro and in vivo evaluation, Journal of Controlled Release, 168 (2013): 70-76.
    [11] [ W. Khan,S. Farah,A. J. Domb, Drug eluting stents: developments and current status, Journal of Controlled Release, 161 (2012): 703-712.
    [12] [ B. M. Mazzag,J. S. Tamaresis,A. I. Barakat, A model for shear stress sensing and transduction in vascular endothelial cells, Biophysical Journal, 84 (2003): 4087-4101.
    [13] [ S. McGinty,S. McKee,R. Wadsworth,C. McCormick, Modeling arterial wall drug concentrations following the insertion of a drug-eluting stent, SIAM Journal on Applied Mathematics, 73 (2013): 2004-2028.
    [14] [ S. McGinty,T. T. N. Vo,M. Meere,S. McKee,C. McCormick, Some design considerations for polymer-free drug-eluting stents: A mathematical approach, Acta Biomaterialia, 18 (2015): 213-225.
    [15] [ S. McGinty,M. Wheel,S. McKee,C. McCormick, Does anisotropy promote spatial uniformity of stent-delivered drug distribution in arterial tissue?, International Journal of Heat and Mass Transfer, 90 (2015): 266-279.
    [16] [ J. C. Panetta, A mathematical model of breast and ovarian cancer treated with paclitaxel, Mathematical Biosciences, 146 (1997): 89-113.
    [17] [ S. Pant,N. W. Bressloff,A. I. Forrester,N. Curzen, The influence of strut-connectors in stented vessels: A comparison of pulsatile flow through five coronary stents, Annals of Biomedical Engineering, 38 (2010): 1893-1907.
    [18] [ G. Pontrelli,F. de Monte, Modelling of mass convection-diffusion in stent-based drug delivery, XXV Congresso Nazionale UIT sulla Trasmissione del Calore, null (2007).
    [19] [ H. Tahir,I. Niculescu,C. Bona-Casas,R. M. H. Merks,A. G. Hoekstra, An in silico study on the role of smooth muscle cell migration in neointimal formation after coronary stenting, Journal of the Royal Society Interface, 12 (2015): 20150358.
    [20] [ A. R. Tzafriri,A. Groothuis,G. S. Price,E. R. Edelman, Stent elution rate determines drug deposition and receptor-mediated effects, Journal of Controlled Release, 161 (2012): 918-926.
    [21] [ T. T. N. Vo,A. Peddle,W. Lee, Modelling chemistry and biology after implantation of a drug-eluting stent. Part Ⅰ: Drug transport, Mathematical Biosciences and Engineering, 14 (2017): 491-509.
    [22] [ T. T. N. Vo,R. Yang,Y. Rochev,M. Meere, A mathematical model for drug delivery, Progress in Industrial Mathematics at ECMI 2010, Mathematics in Industry, 17 (2012): 521-528.
    [23] [ J. Watt,S. Kennedy,C. McCormick,E. O. Agbani,A. McPhaden,A. Mullen,P. Czudaj,B. Behnisch,R. M. Wadsworth,K. G. Oldroyd, Soccinobucal-eluting stents increase neointimal thickening and peri-strut inflammation in a porcine coronary model, Catheterization and Cardiovascular Interventions, 81 (2013): 698-708.
    [24] [ T. K. Yeung,C. Germond,X. Chen,Z. Wang, The mode of action of taxol: Apoptosis at low concentration and necrosis at high concentration, Biochemical and Biophysical Research Communications, 263 (1999): 398-404.
    [25] [ P. Zunino,C. D'Angelo,L. Petrini,C. Vergara,C. Capelli,F. Migliavacca, Numerical simulation of drug eluting coronary stents: Mechanics, fluid dynamics, and drug release, Computer Methods in Applied Mechanics and Engineering, 198 (2009): 3633-3644.
  • This article has been cited by:

    1. Öznur ÖZTUNÇ KAYMAK, Coding theory for h(x)-Fibonacci polynomials, 2024, 26, 1301-7985, 226, 10.25092/baunfbed.1347379
    2. Bahar Demirtürk, Combinatorial Analysis of k-Oresme and k-Oresme–Lucas Sequences, 2025, 17, 2073-8994, 697, 10.3390/sym17050697
  • Reader Comments
  • © 2018 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(4181) PDF downloads(682) Cited by(2)

Figures and Tables

Figures(8)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog