Chapter 6
The R1-Fractional Trigonometry

6.1 R1-Trigonometric Functions

The R1-fractional trigonometry [73–75] is defined by taking the parameter a in the argument of the R-function to be imaginary. This is in contrast to Chapter 5 in which both the a and t parameters were taken to be real. Then, the R1-trigonometry is based on c06-math-0001, where R-function is expanded into even and odd series, that is, series with terms that are even and odd powers of a. This yields

We define c06-math-0005 to be the even part of c06-math-0006, that is, the terms in equation (6.1) with even powers of a, and c06-math-0007 to be the odd part of c06-math-0008, that is, the terms in equation (6.1) with odd powers of a. Then, we have c06-math-0009 and c06-math-0010 from equation (6.3). Therefore, c06-math-0011-function is defined as the real part of equation (6.1) and the c06-math-0012-function as the imaginary part. Because t is raised to a fractional power in these series, the series are multivalued or indexed. To elucidate the indexed behavior of these functions, we consider the roots of t. Then, from equation (3.122),

equation

where for c06-math-0013 the exponent of t is

equation

where c06-math-0014 and c06-math-0015 are assumed rational and irreducible, c06-math-0016, and M/D is rational and in minimal form. From Lorenzo and Hartley [74], with permission of Springer:

Then, the multivalued, or indexed, definition for the c06-math-0017 is given by

with c06-math-0019, c06-math-0020, with common factors in M/D removed. The c06-math-0021 and c06-math-0022 are defined in a similar manner

with c06-math-0024, and

The principal c06-math-0026-generalized trigonometric functions are real functions of the real variable t. They are obtained by taking c06-math-0027, in equations ((6.4)–(6.6). Thus,

6.7 equation
6.8 equation

and

where the simplified notation introduced earlier is used for the principal, k = 0, functions. These functions, equations (6.4)–(6.9), are generalizations of the circular functions (also known as harmonic functions) of the classical integer-order trigonometry. The indexed behavior seen in equations (6.4)–(6.6) has no parallel in the classical trigonometry. We note, when q = 1 and v = 0, the c06-math-0031-trigonometric functions revert back to the circular functions for t > 0.

The principal c06-math-0032-trigonometric functions are shown in Figures 6.1a, b and 6.2 for various values of q. As the value of q decreases from 1 to 0.1, the oscillatory nature of the c06-math-0033- and c06-math-0034-functions is also seen to decrease. For q < 1, it can also be seen that the curves tend to zero as t increases. The functions grow at an increasing rate, as q increases for values of q > 1 (not shown). For smaller values of q in the c06-math-0035 to c06-math-0036 range, the lack of zero crossings of the c06-math-0037-function strongly affects the behavior of the c06-math-0038-function.

A plot marked (a) with t-Time on the horizontal axis, curves plotted, and q = 0.2,  q = 1.0, and q = 1.2 given in the plotted area.; A plot marked (b) with t-Time on the horizontal axis, curves plotted, and q = 0.2,  q = 1.0, and q = 1.2 given in the plotted area.

Figure 6.1 (a) c06-math-0039 versus t-Time, for q = 0.2–1.2 in steps of 0.1, v = 0, a = 1. (b) c06-math-0040 versus t-Time, for q = 0.2–1.2 in steps of 0.1, v = 0, a = 1. (a, b)

Source: Lorenzo and Hartley 2004b [74]. Adapted with permission of Springer.

A plot with t-Time on the horizontal axis, curves plotted, and q values given in the plotted area.

Figure 6.2 c06-math-0041 versus t-Time, for q = 0.2–1.2 in steps of 0.2, v = 0, a = 1.

Figures 6.36.8 show the effects of the a parameter on the R1Cos-, R1Sin-, and the c06-math-0042-functions, respectively. Increasing the a parameter increases the response rate of the function. In other words, increasing a decreases the apparent period and increases the overshoot. The reader should also consider the effect of the order variable q in these figures. For the c06-math-0043 and the c06-math-0044, it also increases the response rate. For the c06-math-0045, it changes the character of the response.

A plot with t-Time on the horizontal axis, curves plotted, and a = 0.25 and a = 1.0 given in the plotted area.

Figure 6.3 Effect of a on c06-math-0046 for a = 0.25–1.0 in steps of 0.25, with q = 0.25, v = 0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and a = 0.25 and a = 1.0 given in the plotted area.

Figure 6.4 Effect of a on c06-math-0048 for a = 0.25–1.0 in steps of 0.25, with q = 0.75, v = 0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and a = 0.25 and a = 1.0 given in the plotted area.

Figure 6.5 Effect of a on c06-math-0049 for a = 0.25–1.0 in steps of 0.25, with q = 0.25, v = 0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and a = 0.25 and a = 1.0 given in the plotted area.

Figure 6.6 Effect of a on c06-math-0050 for a = 0.25–1.0 in steps of 0.25, with q = 0.75, v = 0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and a values given in the plotted area.

Figure 6.7 Effect of a on c06-math-0051 for a = 0.25–1.0 in steps of 0.25, with q = 0.25, v = 0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and a values given in the plotted area.

Figure 6.8 Effect of a on c06-math-0047 for a = 0.25–1.0 in steps of 0.25, with q = 0.75, v = 0 and k = 0.

Figures 6.96.14 show the effects of secondary order parameter v on the R1Cos-, R1Sin-, and c06-math-0052-functions, respectively. In general, increasing v tends to increase the rate of response of the functions and reduce the period. From the Laplace transforms of the functions (Section 6.5), it is clear that v > 0 fractionally differentiates the v = 0 function, while v < 0 integrates it. The effect of v reduces the period of the c06-math-0053-function; this is shown clearly in Figure 6.14.

A plot with t-Time on the horizontal axis, curves plotted, and v = -0.20 and v = 0.25 values given in the plotted area.

Figure 6.9 Effect of v on c06-math-0054 for v = −0.20–0.20 in steps of 0.10, with q = 0.25, a = 1.0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and v = -0.20 and v = 0.20 values given in the plotted area.

Figure 6.10 Effect of v on c06-math-0056 for v = −0.20–0.20 in steps of 0.10, with q = 0.75, a = 1.0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and v = -0.20 and v = 0.20 values given in the plotted area.

Figure 6.11 Effect of v on c06-math-0057 for v = −0.20–0.20 in steps of 0.10, with q = 0.25, a = 1.0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and v = -0.20 and v = 0.20 values given in the plotted area.

Figure 6.12 Effect of v on c06-math-0058 for v = −0.20–0.20 in steps of 0.10, with q = 0.75, a = 1.0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and v values given in the plotted area.

Figure 6.13 Effect of v on c06-math-0059 for v = −0.20–0.20 in steps of 0.10, with q = 0.25, a = 1.0 and k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and v values given in the plotted area.

Figure 6.14 Effect of v on c06-math-0055 for v = −0.20–0.20 in steps of 0.10, with q = 0.75, a = 1.0 and k = 0.

Figure 6.15 shows the phase plane for variations in the order q, with a = 1.0, v = 0.0, and k = 0. Note that for q < 1 the functions are attracted to the origin from infinity, while for q > 1, the functions start at the origin and spiral out to infinity.

A plot with t-Time on the horizontal axis, curves plotted, and q values given in the plotted area. There are arrows pointing to curves with values.

Figure 6.15 Phase plane c06-math-0060 versus c06-math-0061 for q = 1.1 and q = 0.25–1.75 in steps of 0.25, a = 1.0, v = 0.0.

6.1.1 R1-Trigonometric Properties

From the definitions (6.4) and (6.5), we have that

6.10 equation

and

6.11 equation

These results substituted into equation (6.6) yield

6.12 equation

Using the definitions (6.4) and (6.5) with equation (6.3) gives

which is a fractional Euler equation for the R1-trigonometry. Expanding c06-math-0066 in the same manner as done for equation (6.3) yields

a complementary fractional Euler equation. The addition of equations (6.13) and (6.14) gives

If we subtract equation (6.15) from (6.14), we have

Both equations (6.15) and (6.16) are the fractional generalizations of the classical equations (1.18) and (1.19), respectively. Based on the definitions of the c06-math-0070-trigonometric functions, we have

another R-function expression. Combining this with equation (6.15) gives

Similarly, we have

6.19 equation

taken with equation (6.16) yields

A similar relation for the c06-math-0076 may be obtained from the ratio of equations (6.18) and (6.21), that is,

6.22 equation

The addition of the squares equations (6.15) and (6.16) gives

a fractional-order generalization of the classical identity c06-math-0079. Because the right-hand side of this identity is a real number, this complex R product behaves like the multiplication of complex conjugates.

In a similar manner, the subtraction of the square of equation (6.16) from the square of equation (6.15) gives

Figure 6.16 shows the value of the product for a range of q and t values. This is a generalization of the identity c06-math-0081. Figure 6.17 shows the value of the conjugate sum for a range of q and t values. The substitution of equations (6.18) and (6.21) into (6.13) gives

6.25 equation

the fractional Euler equation in an exponential form. Note that for q = 1/2 and v = 0,

equation

Furthermore, with a = 1, we have

equation
A plot with t-Time on the horizontal axis, curves plotted, and q values given in the plotted area.

Figure 6.16 Conjugate product c06-math-0083, q = 0.2–1.4 in steps of 0.2, k = 0.

A plot with t-Time on the horizontal axis, curves plotted, and q values given in the plotted area.

Figure 6.17 Conjugate sum c06-math-0084, q = 0.2–1.2 in steps of 0.2, k = 0, a = 1.0.

6.2 R1-Trigonometric Function Interrelationship

A relationship between the c06-math-0085-function and the c06-math-0086-function is easily determined. Based on equation (6.4)

Thus,

Conversely,

equation

The substitution c06-math-0089 into the defining series for the fractional trigonometric functions gives

6.28 equation
6.29 equation

and

6.30 equation

demonstrating backward compatibility to common trigonometry.

6.3 Relationships to R1-Hyperbolic Functions

The relationship between the R1-hyperbolic and the R1-trigonometric functions may be determined by appropriate substitutions into the defining series. It is clear from equations (5.21) and (6.17) that

6.31 equation

and conversely

6.32 equation

From equations (5.22) and (6.21), we have that

6.33 equation

and conversely

6.34 equation

These relationships between R1-hyperbolic and R1-trigonometric functions exactly parallel and generalize those between the integer-order trigonometric functions (q = 1, v = 0) and classical hyperbolic functions.

6.4 Fractional Calculus Operations on the R1-Trigonometric Functions

This section determines fractional differintegrals of the principal fractional R-trigonometric functions. This section is from Lorenzo and Hartley [74], with permission of Springer:

The c06-math-0097-order differintegral of c06-math-0098 is determined as follows: Differintegrating term-by-term (see Section 3.16) gives

6.35 equation

Then applying equation (5.37), we have

equation

Thus,

Applying equation (6.27), this may be written as

6.37 equation

For c06-math-0102, this becomes

6.38 equation

The c06-math-0104-order differintegral of c06-math-0105, is determined as follows: We may differintegrate (6.26) term-by-term based on Section 3.16. Then, using equation (5.37)

6.39 equation
equation

Thus,

Applying equation (6.27), this may be written as

6.41 equation

and for c06-math-0109, we have the important case

6.42 equation

6.5 Laplace Transforms of the R1-Trigonometric Functions

Because the series are uniformly convergent (see Chapter 14), the Laplace transforms for the principal c06-math-0111-trigonometric functions are readily determined by transforming the defining infinite series term-by-term.

Laplace Transform of R1Cosq. v(a, k, t)

Then, the transform of equation (6.4) is

6.43 equation

Inverse transforming term-by-term (see Section 3.16),

6.44 equation

where c06-math-0114 and c06-math-0115. Thus,

equation

Replacing cos and sin with exponential forms,

6.45 equation
6.46 equation

Now,

therefore, we have

6.48 equation

where c06-math-0120 and c06-math-0121. For the principal function, c06-math-0122

6.5.2 Laplace Transform of R1Sinq. v(a, k, t)

Determination of this transform proceeds in a similar manner as that for the c06-math-0124; thus,

6.50 equation

From Section 3.16, we may transform term-by-term; thus,

6.51 equation

where c06-math-0127 and c06-math-0128. Thus,

equation

Replacing cos and sin with exponential forms,

6.52 equation
equation

Now applying equation (6.47), we have

6.53 equation

where c06-math-0131 and c06-math-0132. For the principal function, c06-math-0133

The principal c06-math-0135-trigonometric and c06-math-0136-hyperbolic functions, their defining series, Laplace transforms, and R-function relationships are presented in Table 6.1 along with those of the classical circular functions.

Table 6.1 Summary of the principal R1-trigonometric, R1-hyperbolic, and the traditional functions

image

Source: Lorenzo and Hartley 2004b [74]. Adapted with permission of Springer.

6.6 Complexity-Based R1-Trigonometric Functions

In Section 6.1, fractional trigonometric functions that parallel the familiar classical trigonometric functions were derived based on the parity-based series (equation (6.2). In this section, we shall consider complexity-based functions starting from equation (6.1). Equation (6.1) can be written as

6.55 equation
6.56 equation

Furthermore, we have that

6.57 equation

a complex plane version of the R-function. Because the results are R-functions which have been studied earlier, we do not pursue this area further.

6.7 Fractional Differential Equations

Continuing from Ref. [74], with permission of Springer:

The c06-math-0179-trigonometric functions may also be solutions to fractional differential equations. This transform (equation (6.49) indicates that the uninitialized fractional differential equation

has a solution of the form

6.59 equation

when c06-math-0182, that is, an impulse function. As with the c06-math-0183-hyperbolic functions for arbitrary c06-math-0184, the solution is given by the convolution

Assuming composition equation (6.58) may also be written in the form of a fractional integro-differential equation, namely

equation

In a similar manner, equation (6.54) indicates that the uninitialized fractional differential equation

has a solution of the form

6.62 equation

when c06-math-0188, that is, an impulse function. For arbitrary c06-math-0189, the solution is obtained by the convolution similar to equation (6.60). Equation (6.61) may also be written as a fractional integro-differential equation, namely

6.63 equation

Linear combinations of fractional trigonometric functions may also be used to infer solutions for fractional differential forms of noncommensurate order. For example, from the linear combination

or using the relationships (6.17) and (6.20)

equation

we may infer the Laplace transformed equation

6.65 equation

or

equation

From this, it is seen that the uninitialized fractional differential equation

6.66 equation

where

equation

will have a solution of the form given by equation (6.64). Various specializations of the orders and parameters lead to other useful results. Still more varied fractional differential equations may be obtained by considering linear combinations of the both the R1-trigonometric functions and the c06-math-0194-hyperbolic functions, using, for example, forms such as

equation

Chapter 7 explores the effect of imaginary t in the definition of a fractional trigonometry.