# Golden Ratio Phenomenon of Random Data Obeying von Karman Spectrum

Authors: Ming Li; Wei Zhao
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

von Karman originally deduced his spectrum of wind speed fluctuation based on the Stokes-Navier equation. Taking into account, the practical issues of measurement and/or computation errors, we suggest that the spectrum can be described from the point of view of the golden ratio. We call it the golden ratio phenomenon of the von Karman spectrum. To depict that phenomenon, we derive the von Karman spectrum based on fractional differential equations, which bridges the golden ratio to the von Karman spectrum and consequently provides a new outlook of random data following the von Karman spectrum in turbulence. In addition, we express the fractal dimension, which is a measure of local self-similarity, using the golden ratio, of random data governed by the von Karman spectrum.


---

## Body

## 1. Instruction

The golden ratio, denoted byφ, is an irrational number given by φ=(1+5)/2 [1]. The paper by Ackermann [2] may likely be the earliest literature on the golden ratio in a mathematics journal in English in 1895, but it attracted and has attracted the interest of scientists and engineers in various fields of sciences and engineering, ranging from chemistry to computer science; see, for example, [1], Benassi [3], Putz [4], Orita et al. [5], Perez [6], Hassaballah et al. [7], Kellerhals [8], Henein et al. [9], Hurtley [10], Coldea et al. [11], Affleck [12], Jones et al. [13], Kaygn et al. [14], Cervantes et al. [15], Chebotarev [16], Benavoli et al. [17], Manikantan et al. [18], Assimakis et al. [19], Good [20], Davis and Jahnke [21], Totland [22], Moufarrège [23], Boeyens [24], Iñiguez et al. [25], Andrews and Zhang [26], Hofri and Rosberg [27], Itai and Rosberg [28], Cassandras and Julka [29], and Tanackov et al. [30], just to mention a few.In the field of random functions, more precisely, turbulence in fluid mechanics, a kind of power spectra density (PSD) function introduced by von Karman [31], known as the von Karman spectra (VKS), has been widely used in the diverse fields, ranging from turbulence to acoustic wave propagation in random media; see, for example, Goedecke et al. [32] and the references therein. Among the von Karman spectra, the spectrum (VKSW for short) expressed in (1) is particularly useful in the field of wind engineering for the modeling of wind speed fluctuation; see, for example, [33–41]. That PSD is in the form
(1)Svon(f)=4uf2bvwf(1+70.8w2)5/6,w=fLuxU,
where f is frequency (Hz), Lux is turbulence integral scale, U is mean speed, uf is friction velocity (ms−1), and bv is friction velocity coefficient such that the variance of wind speed σu2=bvuf2.Note that (1) was conventionally deduced based on the Stokes-Navier equation ([31], Bauer and Zeibig [42], Tropea [43], Monin and Yaglom [44], Xiushu [45]). It does not originally relate to the concept of either the golden ratio or fractal dimension. As a matter of fact, reports regarding turbulence’s fractal dimension derived directly based on the Stokes-Navier equation are rarely seen as Gaoan stated in [46, page 55], letting alone the golden ratio.This paper aims at contributing the following three results. First, we will propose a rigorous but concise derivation of (1). Then, we will generalize (1) such that the generalization may be described from the point of view of the golden ratio. Finally, we will explain the golden ratio phenomenon of the VKSW from the point of view of fractal dimension or local self-similarity. As a result, we achieve the goal of bridging the golden ratio to the VKSW as well as self-similarity of random data, establishing a new outlook of data following the VKSW.The rest of the paper is organized as follows. The preliminaries are briefed in Section2. The results are given in Section 3, which is followed by conclusions.


## 2. Preliminaries

### 2.1. Golden Ratio

One of the conventional ways to deduce the golden ratioφ is to solve the difference equation that produces the Fibonacci sequences. The equation is given by (see, e.g., [1], Jamieson [46], Ranum [47], and Eggar [48])
(2)F(n)=F(n-1)+F(n-2),n∈N(thesetofnaturalnumbers).
Denote by ZF(z) the z-transform of F(n). Then, doing the z-transform on both sides of (2) yields
(3)ZF(z)=ZF(z)(z-1+z-2).
The above expression can be rewritten as
(4)z2-z-1=0.
The solutions to (4) are expressed by
(5)z1,2={1+521-52.
The golden ratio equals z1; that is,
(6)φ=z1=1+52≈1.618.
In addition,
(7)z2=1-52=-1φ≈-0.618.


### 2.2. Fractional Oscillators

There are three types of fractional oscillators. The conventional type, see, for example, Achar et al. [49, 50], is given by
(8)d2-εx(t)dt2-ε+ω02x(t)=e(t),0<ε<1.
The second type was introduced by Lim and Muniandy [51]. It is in the form
(9)(d2dt2+λ)βx(t)=e(t),β>0.
The third type introduced by Lim and Teo [52] is expressed by
(10)(Datα+λ)βxα,β(t)=e(t),0<α<1,β>0.
The symbol Datα is a fractional differential operator; see, for example, Eab and Lim [53, 54], Lim et al. [55], Klafter et al. [56], Machado et al. [57], and Cattani [58]. In what follows, we use (10) in the general sense.We now consider the fractional differential equation with the coefficientA in the form
(11)A(ddt+λ)βyfOU(t)=η(t),β>0,
where η(t) is a white noise. Then, the solution to (11) is the fractional Ornstein-Uhlenbeck (OU) process [59], referring to Coffey et al. [60] for the meaning of OU process.


## 2.1. Golden Ratio

One of the conventional ways to deduce the golden ratioφ is to solve the difference equation that produces the Fibonacci sequences. The equation is given by (see, e.g., [1], Jamieson [46], Ranum [47], and Eggar [48])
(2)F(n)=F(n-1)+F(n-2),n∈N(thesetofnaturalnumbers).
Denote by ZF(z) the z-transform of F(n). Then, doing the z-transform on both sides of (2) yields
(3)ZF(z)=ZF(z)(z-1+z-2).
The above expression can be rewritten as
(4)z2-z-1=0.
The solutions to (4) are expressed by
(5)z1,2={1+521-52.
The golden ratio equals z1; that is,
(6)φ=z1=1+52≈1.618.
In addition,
(7)z2=1-52=-1φ≈-0.618.


## 2.2. Fractional Oscillators

There are three types of fractional oscillators. The conventional type, see, for example, Achar et al. [49, 50], is given by
(8)d2-εx(t)dt2-ε+ω02x(t)=e(t),0<ε<1.
The second type was introduced by Lim and Muniandy [51]. It is in the form
(9)(d2dt2+λ)βx(t)=e(t),β>0.
The third type introduced by Lim and Teo [52] is expressed by
(10)(Datα+λ)βxα,β(t)=e(t),0<α<1,β>0.
The symbol Datα is a fractional differential operator; see, for example, Eab and Lim [53, 54], Lim et al. [55], Klafter et al. [56], Machado et al. [57], and Cattani [58]. In what follows, we use (10) in the general sense.We now consider the fractional differential equation with the coefficientA in the form
(11)A(ddt+λ)βyfOU(t)=η(t),β>0,
where η(t) is a white noise. Then, the solution to (11) is the fractional Ornstein-Uhlenbeck (OU) process [59], referring to Coffey et al. [60] for the meaning of OU process.


## 3. Results

Denote thatgfOU(t) is the impulse response function of (11). Then, it is the solution to the following equation with zero initial conditions
(12)A(ddt+λ)βgfOU(t)=δ(t),
where δ(t) is the Dirac-δ function. Doing the Fourier transforms on both sides of the above equation yields
(13)GfOU(f)=A(λ-j2πf)β,
where GfOU(f) is the Fourier transform of gfOU(t).LetSyfOU(f) be the PSD of yfOU(t). Then, SyfOU(f) is given by
(14)SyfOU(f)=|GfOU(f)|2=A2[λ2+(2πf)2]β.
Thus, we have the theorem below.Theorem 1.
LetXvk(t) be the random function that obeys (1). Then, Xvk(t) is governed by the fractional differential equation that is in the form
(15)Avk(ddt+Bvk)5/6Xvk(t)=η(t).
Its solution in frequency domain is given by
(16)Svk(f)=Avk[(Bvk)2+(2πf)2]5/6.Proof.
ReplacingA, λ, and β in (11) with Avk, Bvk, and 5/6, respectively, yields (15). Substituting A, λ, and β in (14) with Avk, Bvk, and 5/6, respectively, produces (16). Thus, Theorem 1 results.From Theorem1, we obtain the following corollary.Corollary 2 (modified VKSW).
LetXvkφ(t) be the random function that is governed by the fractional differential equation given by
(17)Avk(ddt+Bvk)φ/2Xvkφ(t)=η(t).
Then, its solution in frequency domain is in the form
(18)Svkφ(f)=Avk[(Bvk)2+(2πf)2]φ/2.The proof is straightforward and omitted consequently.From Theorem1 and Corollary 2, we immediately have the remark below.Remark 3.
The VKSW may be approximately expressed by the golden ratio.As a matter of fact,(19)φ2≈56.
Thus, (16) approximately equals (18). This may not be a simple approximation but substantially develops the implication of the VKSW from the point of view of the golden ratio.Note that there are errors in measuring real random data [61–63] and computation errors [64–66]. Thus, from a view of practice, the power of the VKSW may not exactly be the value of 5/6 in most cases in engineering. Rather, it may be in the form
(20)(56)+e,
where e is error. Thus, by using the golden ratio, (18) is quite reasonable to characterize random functions that obey the VKSW.For the purpose of exhibiting the results in time domain, we denote byF and F-1 the operator of the Fourier transform and its inverse, respectively. Then, we get the theorem below.Theorem 4.
The inverse Fourier transform of1/[1+(2πf)2]φ/2 is given by
(21)F-1{1[1+(2πf)2]φ/2}=2πΓ(φ/2)(|τ|2)(φ-1)/2K(φ-1)/2(|τ|),
where Kv(z) is the modified Bessel function of second kind or the MacDonald function and τ is the time lag.Proof.
Because(φ-1)/2>1/2, according to the computation formula in Gelfand and Vilenkin [67, page 188, in Section 2, Chapter 2], (21) holds.Recall that the Fourier transform of|t|α is expressed by [68]
(22)F(|t|α)=-2sin(απ2)Γ(α+1)|ω|-α-1,
where α≠1,3,….Note that(23)Svkφ(f)~1fφ/2forf⟶∞.
Then, denoting rvkφ(τ) is the inverse Fourier transform of Svkφ(f), one has
(24)rvkφ(τ)~|τ|(φ-1)/2forτ⟶0.
The fractal dimension of a process can be determined by its autocorrelation function (ACF) for τ→0 [69]. Thus,
(25)rvk(0)-rvk(τ)~|τ|(φ-1)/2forτ⟶0.
Therefore, with the probability one [69], the fractal dimension of the modified von Karman process based on the golden ratio is given by
(26)Dvkφ=(2-φ-14)=7-φ4.
Approximately, it is expressed by
(27)Dvkφ≈1.346.Note that fractal dimension is a measure of local self-similarity, irregularity, or roughness [70]. High value of fractal dimension of a sample path implies high irregularity of that path. Thus, (26) means that the modified von Karman process with the golden ratio has considerable local irregularity.We would like to call out the work described above as the golden ratio phenomenon of the von Karman process. From the point of view of our work in data science or big data, this research may not be enough. The future work will investigate possible golden ratio phenomena in other topics of data such as those discussed in [71–82], exploring laws associating with the golden ratio in the universe.


## 4. Conclusions

We have given the derivation of the von Karman spectrum based on the fractional differential equation (10). The results suggest that the process obeying VKSW is in the class of fractional OU processes. Moreover, we have explained the reasons why the VKSW may be described from the point of view of the golden ratio. The fractal dimension of random data obeying the VKSW by using the golden ratio has also been discussed.



---

*Source: 130258-2013-10-02.xml*