
In this paper, the Levenberg-Marquardt backpropagation scheme is used to develop a neural network model for the examination of the fluid flow on a magnetized flat surface with slip boundaries. The tangent hyperbolic fluid is considered along with heat generation, velocity, and thermal slip effects at the surface. The problem is modelled in terms of a non-linear differential system and Lie symmetry is used to get the scaling group of transformation. The order reduction of differential equations is done by using Lie transformation. The reduced system is solved by the shooting method. The surface quantity, namely skin friction, is evaluated at the surface for the absence and presence of an externally applied magnetic field. A total of 88 sample values are estimated for developing an artificial neural network model to predict skin friction coefficient (SFC). Weissenberg number, magnetic field parameter, and power law index are considered three inputs in the first layer, while 10 neurons are taken in the hidden layer. 62 (70%), 13 (15%), and 13 (15%) samples are used for training, validation, and testing, respectively. The Levenberg-Marquardt backpropagation is used to train the network by entertaining the random 62 sample values. Both mean square error and regression analysis are used to check the performance of the developed neural networking model. The SFC is noticed to be high at a magnetized surface for power law index and Weissenberg number.
Citation: Khalil Ur Rehman, Wasfi Shatanawi, Zead Mustafa. Artificial intelligence (AI) based neural networks for a magnetized surface subject to tangent hyperbolic fluid flow with multiple slip boundary conditions[J]. AIMS Mathematics, 2024, 9(2): 4707-4728. doi: 10.3934/math.2024227
[1] | Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Neimark-Sacker bifurcation, chaos, and local stability of a discrete Hepatitis C virus model. AIMS Mathematics, 2024, 9(11): 31985-32013. doi: 10.3934/math.20241537 |
[2] | Yao Shi, Zhenyu Wang . Bifurcation analysis and chaos control of a discrete fractional-order Leslie-Gower model with fear factor. AIMS Mathematics, 2024, 9(11): 30298-30319. doi: 10.3934/math.20241462 |
[3] | Ibraheem M. Alsulami . On the stability, chaos and bifurcation analysis of a discrete-time chemostat model using the piecewise constant argument method. AIMS Mathematics, 2024, 9(12): 33861-33878. doi: 10.3934/math.20241615 |
[4] | Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Discrete Hepatitis C virus model with local dynamics, chaos and bifurcations. AIMS Mathematics, 2024, 9(10): 28643-28670. doi: 10.3934/math.20241390 |
[5] | A. Q. Khan, Ibraheem M. Alsulami, S. K. A. Hamdani . Controlling the chaos and bifurcations of a discrete prey-predator model. AIMS Mathematics, 2024, 9(1): 1783-1818. doi: 10.3934/math.2024087 |
[6] | Kamel Guedri, Rahat Zarin, Ashfaq Khan, Amir Khan, Basim M. Makhdoum, Hatoon A. Niyazi . Modeling hepatitis B transmission dynamics with spatial diffusion and disability potential in the chronic stage. AIMS Mathematics, 2025, 10(1): 1322-1349. doi: 10.3934/math.2025061 |
[7] | Yulong Li, Long Zhou, Fengjie Geng . Dynamics on semi-discrete Mackey-Glass model. AIMS Mathematics, 2025, 10(2): 2771-2807. doi: 10.3934/math.2025130 |
[8] | Xiongxiong Du, Xiaoling Han, Ceyu Lei . Dynamics of a nonlinear discrete predator-prey system with fear effect. AIMS Mathematics, 2023, 8(10): 23953-23973. doi: 10.3934/math.20231221 |
[9] | Muhammad Farman, Ali Akgül, J. Alberto Conejero, Aamir Shehzad, Kottakkaran Sooppy Nisar, Dumitru Baleanu . Analytical study of a Hepatitis B epidemic model using a discrete generalized nonsingular kernel. AIMS Mathematics, 2024, 9(7): 16966-16997. doi: 10.3934/math.2024824 |
[10] | Ali Al Khabyah, Rizwan Ahmed, Muhammad Saeed Akram, Shehraz Akhtar . Stability, bifurcation, and chaos control in a discrete predator-prey model with strong Allee effect. AIMS Mathematics, 2023, 8(4): 8060-8081. doi: 10.3934/math.2023408 |
In this paper, the Levenberg-Marquardt backpropagation scheme is used to develop a neural network model for the examination of the fluid flow on a magnetized flat surface with slip boundaries. The tangent hyperbolic fluid is considered along with heat generation, velocity, and thermal slip effects at the surface. The problem is modelled in terms of a non-linear differential system and Lie symmetry is used to get the scaling group of transformation. The order reduction of differential equations is done by using Lie transformation. The reduced system is solved by the shooting method. The surface quantity, namely skin friction, is evaluated at the surface for the absence and presence of an externally applied magnetic field. A total of 88 sample values are estimated for developing an artificial neural network model to predict skin friction coefficient (SFC). Weissenberg number, magnetic field parameter, and power law index are considered three inputs in the first layer, while 10 neurons are taken in the hidden layer. 62 (70%), 13 (15%), and 13 (15%) samples are used for training, validation, and testing, respectively. The Levenberg-Marquardt backpropagation is used to train the network by entertaining the random 62 sample values. Both mean square error and regression analysis are used to check the performance of the developed neural networking model. The SFC is noticed to be high at a magnetized surface for power law index and Weissenberg number.
Schoenberg introduced the trigonometric spline functions defined by divided differences in [18]. And the trigonometric spline functions have been shown to possess many B-spline-like properties. In view of this, many scholars call the trigonometric splines the trigonometric B-splines, in [4,8]. It is well-known that the trigonometric B-splines are piecewise functions corresponding to the spaces
T2n+1:=span{1,cost,sint,cos2t,sin2t,⋯,cosnt,sinnt}, |
for odd-order, and
H2n:=span{cost2,sint2,cos3t2,sin3t2,⋯,cos(2n−1)t2,sin(2n−1)t2}, |
for even-order. Odd-order trigonometric B-splines in T2n+1 form a partition of unity, a desirable property for curve design. However, even-order B-splines in H2n lack this partition of unity, creating challenges in certain applications where this property is critical. As a result, extensive research has focused on the odd-order basis, particularly on its normalization. The author of [8] established the recurrence relation for the trigonometric B-splines of arbitrary order and derived trigonometric Marsden identity. The author of [20] utilized the trigonometric Marsden identity to derive the normalized odd-order trigonometric B-splines over uniform knots. Building on the trigonometric Marsden identity expansion introduced by [7], [11] explicitly derived the normalization coefficients required to ensure that the trigonometric B-spline basis functions are properly normalized. [15] provided the p-Bézier basis functions in the space T2n+1, which is defined over any interval of length <π. Those basis functions are a subcase of the normalized trigonometric B-splines. [2] presented the normalized Bernstein-like basis functions in the space T2n+1, which is defined over the interval [0,π/2]. [21] established the C-B spline basis. In this work, the C-B spline basis of order 3 is just the normalized trigonometric B-spline basis corresponding to the space T3. In a more general context, the normalized trigonometric B-splines are also considered in [16] as a special case.
Curve and surface design is an important area in Computer Aided Geometric Design (CAGD), where trigonometric B-splines are a foundational tool. Normalized trigonometric B-splines enable enhanced construction and control of curves and surfaces, suggesting their potential applications in aircraft design [9]. Additionally, trigonometric B-splines show promising applications in other fields, such as physical simulations (see [5,12,19]).
The purpose of this paper is to introduce the integral formula for odd-order trigonometric B-splines. Given the significant applications and theoretical importance of integral formulas in Chebyshev systems for fields such as numerical analysis, signal processing, and function approximation, it is notable that the integral properties of sin and cos in trigonometric B-spline bases render them incapable of being directly derived like other bases in Chebyshev systems. The aim of this study is to provide a similar integral formula for trigonometric B-spline bases. To achieve this, this paper first constructs a novel set of even-order trigonometric B-spline curve basis functions and, through integration of these functions, successfully derives the traditional odd-order trigonometric B-spline basis functions, thereby establishing the integral formula for odd-order trigonometric B-spline bases. During this derivation process, a determinant of even order with structural symmetry is obtained. Furthermore, this study refines the conditions for knot sequences to ensure that the corresponding normalized trigonometric B-spline bases possess nonnegativity.
Our main contributions are
● A set of trigonometric spline bases corresponding to the DC component-free space is provided, along with the integral representation of normalized trigonometric B-spline bases corresponding to space T2n+1.
● Adjusting the conditions imposed on the knot sequence to guarantee the nonnegativity of the normalized trigonometric B-spline basis functions.
● A structurally symmetric even-order determinant is presented.
The remainder of this article is organized as follows: Section 2 reviews the related concepts and properties. The improvements to knot sequences are discussed in Section 3. In Section 4, the trigonometric spline basis corresponding to the DC component-free space and the integral formula for the normalized trigonometric B-spline basis are presented. The final conclusions are drawn in the last section.
In this section, we will recall some established concepts and conclusions.
[1] has demonstrated that the knot sequences for B-spline basis functions can be finite, infinite, or bi-infinite. Analogous to the case of traditional B-spline basis functions, this paper focuses on the study of trigonometric B-spline basis functions corresponding to bi-infinite knot sequences.
The normalized trigonometric B-spline basis functions are defined in a manner analogous to the de Boor-Cox formula [4,7,8,11].
Definition 2.1. (Normalized trigonometric B-spline basis functions) Given a knot sequence T={ti}+∞i=−∞, such that
ti≤ti+1,0<ti+2n+1−ti<2π,i∈Z,n∈Z+, | (2.1) |
the normalized trigonometric B-spline basis functions of order 2n+1 are defined as follows:
Ki,2n+1(t)=ζi,2n+1Ni,2n+1(t), | (2.2) |
where
Ni,1(t)={1, if ti≤t<ti+1,0,otherwise. | (2.3) |
Ni,2n(t)=sin(t−ti2)sin(ti+2n−1−ti2)Ni,2n−1(t)+sin(ti+2n−t2)sin(ti+2n−ti+12)Ni+1,2n−1(t), | (2.4) |
Ni,2n+1(t)=sin(t−ti2)sin(ti+2n−ti2)Ni,2n(t)+sin(ti+2n+1−t2)sin(ti+2n+1−ti+12)Ni+1,2n(t), | (2.5) |
ζi,2n+1=1(2n)!∑μn∏j=1costi+μ(2j)−ti+μ(2j−1)2, | (2.6) |
and the sum is taken over all permutations μ:{1,2,⋯,2n}→{1,2,⋯,2n}.
Remark 2.1 In [4,8,11], the knot sequence satisfies the condition 0<ti+2n+1−ti<2π. However, the knot sequence described in [7] satisfies a slightly different condition, 0<ti+2n−ti<2π. In the subsequent section, the conditions of the knot sequence are reiterated.
If ti−1<ti=ti+1=⋯=ti+mi−1<ti+mi, where 1≤mi≤2n, the knots tj, where j=i,i+1,⋯,i+mi−1, are referred to as knots of multiplicity mi. Especially, we set 00=0. The space of trigonometric spline basis is defined by
Γ2n+1[T]:={Ni,2n+1(t)|Ni,2n+1(t)|t∈[ti,ti+1)∈T2n+1, and N(l)i,2n+1(ti−)=N(l)i,2n+1(ti+),0≤l≤2n−mi,i∈Z}. |
The trigonometric B-spline basis functions possess many B-spline-like properties [3,6,8,14].
Property 2.1. (Properties of the trigonometric B-spline basis functions) The trigonometric B-spline basis functions Ni,2n+1(t) defined in Eq (2.5) possess the following properties:
(1) (Local support) For any i∈Z and n∈Z+, there exists
Ni,2n+1(t)={>0,t∈(ti,ti+2n+1),=0,t∉[ti,ti+2n+1]. | (2.7) |
(2) (Continuity) The continuous order of Ni,2n+1(t) at tj (where i≤j≤i+2n+1), denoted as kji,2n+1, can be described as
kji,2n+1={2n−ξ,i≤j≤i+ξ−1,2n−mj, i+ξ≤j≤i+2n+1−η,2n−η, i+2n+2−η≤j≤i+2n+1, | (2.8) |
if ti=ti+1=⋯=ti+ξ−1<ti+ξ≤ti+ξ+1≤⋯≤ti+2n+1−η<ti+2n+2−η=⋯=ti+2n+1.
In this section, we will adjust the conditions of the knot sequence to ensure the nonnegativity of the normalized trigonometric B-spline basis functions Ki,2n+1(t) for i∈Z and n∈Z+.
First, the condition (2.1) can be relaxed.
The nonnegativity of the basis function Ni,2n+1(t) is ensured by the condition 0<ti+2n+1−ti<2π as stated in (2.1). Since the length of the support intervals for Ni,2n(t) and Ni+1,2n(t) in Eq (2.4) is 2n, replacing 0<ti+2n+1−ti<2π with the condition 0<ti+2n−ti<2π still guarantees the nonnegativity of Ni,2n+1(t). Therefore, the condition (2.1) can be relaxed to
ti≤ti+1,0<ti+2n−ti<2π,i∈Z,n∈Z+. | (3.1) |
Second, the condition to ensure the positivity of the normalized coefficients ζi,2n+1 in Eq (2.6) is presented. According to the representation in Eq (2.6), we obtain the condition to ensure the positivity of the normalized coefficient ζi,2n+1, that is
0≤ti+2n−1−ti<π, i∈Z,n∈Z+. | (3.2) |
The Bézier-like basis defined in the space T2n+1 is a subcase of the trigonometric B-spline basis. In [13], it was noted that the space T2n+1 does not have an NTP basis when the domain of T2n+1 is [0,π]. This implies the non-existence of normalized coefficients. Consequently, we conclude that ti+2n−1−ti≠π. In [15], it proved that there exist NTP bases provided that the domain of T2n+1 is any interval of length <π, specifically referred to as p-Bézier bases. This implies that the normalized coefficient satisfies the condition ti+2n−1−ti<π, in the special cases of Bézier.
Figure 1a presents an example of the basis functions Ni,2n+1(t) for i∈Z, whose corresponding knot sequence T satisfies condition (3.1) but fails to satisfy condition (2.1). Here n=2 and the knot sequence T={ti}16i=1={1,1,1,1,1,3,4,5.3,6.5,7.5,8.8,10,10,10,10,10}. Figure 1b shows the trigonometric B-spline basis functions Ni,3(t) and Ki,3(t) corresponding to the knot sequence that satisfies condition (3.1) but does not satisfy condition (3.2), while Figure 1c presents the functions Ni,3(t) and Ki,3(t) corresponding to the knot sequence that satisfies both conditions (3.1) and (3.2).
Although the coefficients of the normalized trigonometric B-spline basis have already been defined in Eq (2.6), for the convenience of deriving and proving the integral formula of the normalized trigonometric B-spline basis, this paper introduces a simplified normalized coefficient expression that is equivalent to Eq (2.6).
Lemma 3.1. Let
Ci,2n+1=1(2n−1)!!∑γ(1,2,…,2n)n∏r=1costi+m2r−ti+m2r−12,i∈Z,n∈Z+, | (3.3) |
where the sum is taken over all permutations γ(1,2,…,2n):{m1,m2,…,m2n}→{1,2,…,2n}, with m1,m2,…,m2n being a permutation of 1,2,…,2n that satisfies m1<m3<…<m2n−1 and m2r−1<m2r for each r=1,…,n. Then there exists
Ci,2n+1=ζi,2n+1, |
where ζi,2n+1 defined in Eq (2.6).
Third, a new definition for the normalized trigonometric B-spline basis that guarantees nonnegativity is presented.
Definition 3.1. (Nonnegative normalized trigonometric B-spline basis functions) Given a knot sequence T={ti}+∞i=−∞, such that
ti≤ti+1, 0<ti+2n−ti<2π and 0≤ti+2n−1−ti<π, i∈Z,n∈Z+, | (3.4) |
the nonnegative normalized trigonometric B-spline basis functions of order 2n+1 are defined as follows:
Ki,1(t)={1, if ti≤t<ti+1,0,otherwise. | (3.5) |
Ki,2n+1(t)=Ci,2n+1Ni,2n+1(t), | (3.6) |
where Ni,2n+1(t), Ci,2n+1 is defined in Eqs (2.5) and (3.3), respectively.
In this section, we generally set the knot sequence to be T={ti}+∞i=−∞ that satisfies condition (3.4), where the multiplicity of ti is mi with 1≤mi≤2n, and kji,2n+1 is defined in Eq (2.8) without further explanation. Let us define the DC component-free space as T2n:=span{cost,sint,⋯,cosnt,sinnt}. Its corresponding piecewise trigonometric polynomial space is
Γ2n[T]:={Fi,2n(t)|Fi,2n(t)|t∈[ti,ti+1)∈T2n, and F(l)i,2n(ti−)=F(l)i,2n(ti+),0≤l≤2n−mi−1, 1≤mi≤2n−1,i∈Z}. |
Clearly, Γ2n[T] is a linear space. We can derive the subspace of Γ2n[T] as follows:
Γ2n[ti,ti+2n]={Fi,2n(t)∈Γ2n[T]|Fi,2n(t)=0,t∉[ti,ti+2n] and Fi,2n(t)≠0,t∈(ti,ti+2n),F(l)i,2n(tj−)=F(l)i,2n(tj+),0≤l≤kji,2n, i≤j≤i+2n}. |
It can be shown that a function in the space \(\Gamma_{2n}[t_i, t_{i + 2n}]\) exhibits local support and a specific order of continuity.
The following determinant and its accompanying proof are presented to facilitate future derivations.
Lemma 4.1. For any n∈Z+, then the following identity holds.
D(t1,t2,⋯,t2n):=|cost1cost2…cost2nsint1sint2…sint2ncos2t1cos2t2…cos2t2nsin2t1sin2t2…sin2t2n⋮⋮⋮⋮cosnt1cosnt2…cosnt2nsinnt1sinnt2…sinnt2n|=22n2−nn!∏1≤l<j≤nsintj−tl2∑γ(1,2,…,2n)n∏r=1costm2r−tm2r−12, | (4.1) |
where γ(1,2,…,2n) defined in Lemma 3.1.
Proof. According to Euler's formula and the identity eitj−eits=2iei2(tj+ts)sintj−ts2, it follows that
D(t1,t2,⋯,t2n)=|eit1+e−it12eit2+e−it22…eit2n+e−it2n2eit1−e−it12ieit2−e−it22i…eit2n−e−it2n2i⋮⋮⋮⋮enit1+e−nit12enit2+e−nit22…enit2n+e−nit2n2enit1−e−nit12ienit2−e−nit22i…enit2n−e−nit2n2i|=(−1)n2+n(2i)ne−nit1…e−nit2n|11…1eit1eit2…eit2n⋮⋮⋮⋮e(n−1)it1e(n−1)it2…e(n−1)it2ne(n+1)it1e(n+1)it2…e(n+1)it2n⋮⋮⋮⋮e2nit1e2nit2…e2nit2n|=22n2−n2ne−nit1…e−nit2n∑1≤m1<m2<…<mn≤2neitm1…eitmn∏1≤s<j≤2nei2(tj+ts)sintj−ts2. |
For the sake of convenience, let
C(t1,t2,⋯,t2n)=∑γ(1,2,…,2n)n∏r=1costm2r−tm2r−12. | (4.2) |
Therefore, it suffices to prove that
C(t1,t2,⋯,t2n)=n!2ne−nit1…e−nit2n∑1≤m1<m2<…<mn≤2neitm1…eitmn∏1≤s<j≤2nei2(tj+ts)=n!2n∑1≤m1<m2<…<mn≤2nexp(i2mn∑q=m1tq−i2(2n∑h=1th−mn∑q=m1tq)). | (4.3) |
We establish the inductive hypothesis for n. For n=1, the result is straightforward. Assume that the conclusion holds for n≤p−1, where \(p\) is any positive integer. Specifically, we have
C(tn1,tn2,⋯,tn2p−2)=(p−1)!2(p−1)∑εexp(i2mp−1∑q=m1tq−i2(n2p−2∑h=n1th−mp−1∑q=m1tq)), |
where the sum is taken over all permutations ε:{m1,m2,⋯,mp−1}→{n1,n2,⋯,n2p−2}. Here, n1,n2,⋯,n2p−2 represent a permutation of 1,2,⋯,2n such that n1<n2<⋯<n2p−2, while m1,m2,⋯,mp−1 is a subset of n1,n2,⋯,n2p−2 satisfying m1<m2<⋯<mp−1.
Next, consider the case where n=p. Based on this assumption, we have
C(t1,t2,⋯,t2p)=cost2p−t2p−12∑γ(1,2,⋯,2p−2)p−1∏r=1costm2r−tm2r−12+cost2p−t2p−22∑γ(1,2,⋯,2p−3,2p−1)p−1∏r=1costm2r−tm2r−12+⋯+cost2p−t22∑γ(1,3,4,⋯,2p−1)p−1∏r=1costm2r−tm2r−12+cost2p−t12∑γ(2,3,⋯,2p−1)p−1∏r=1costm2r−tm2r−12=12(exp(i2(t2p−t2p−1))+exp(−i2(t2p−t2p−1)))C(t1,t2,⋯,t2p−2)+12(exp(i2(t2p−t2p−2))+exp(−i2(t2p−t2p−2)))C(t1,t2,⋯,t2p−3,t2p−1)+⋯+12(exp(i2(t2p−t2))+exp(−i2(t2p−t2)))C(t1,t3,t4,⋯,t2p−1)+12(exp(i2(t2p−t1))+exp(−i2(t2p−t1)))C(t2,t3,⋯,t2p−1)=p!2p∑1≤m1<m2<…<mp≤2pexp(i2mp∑q=m1tq−i2(2p∑h=1th−mp∑q=m1tq)), |
where the sum is taken over all permutations γ(n1,n2,⋯,n2p−2):{m1,m2,⋯,m2p−2}→{n1,n2,⋯,n2p−2}, with m1,m2,⋯,m2p−2 being a permutation of n1,n2,⋯,n2p−2 that satisfies m1<m3<⋯<m2p−3 and m2r−1<m2r for each r=1,⋯,p−1. Here, the sequence {n1,n2,⋯,n2p−2} represents, in order, the sequences {1,2,⋯,2p−2}, {1,2,⋯,2p−3,2p−1}, ⋯, {1,3,4,⋯,2p−1}, and \(\{2, 3, \cdots, 2p-1\}\). Eq (4.3) is valid for any positive integer n. Thus, the lemma is proved.
This subsection defines a set of truncated functions and demonstrates that any function in the space Γ2n[ti,ti+2n] can be expressed as a linear combination of these functions.
Let g2n(t):=sintsin2n−2(t2)(n∈Z+), the truncated functions Gi,2n(t), where i∈Z and n∈Z+, are defined as follows: If ti=ti+1=⋯=ti+ξ−1<ti+ξ, then
Gi,2n(t):={0,t<ti,g(ξ−1)2n(t−ti),t≥ti. | (4.4) |
Figure 2a–2c illustrate examples of the function g4(t) and the truncated functions Gi,4(t), where i∈Z, over single and multiple knots, respectively.
To prove that any function in the space Γ2n[ti,ti+2n] can be represented as a linear combination of {Gi,2n(t)}i∈Z, the following lemmas will be utilized.
Lemma 4.2. There exist n+1 real numbers d0,d1,⋯,dn such that
sin2n(t2)=d0+d1cost+d2cos2t+⋯+dncosnt, |
where d0=(2n−1)!!(2n)!!.
Proof. It is well known that
sin2n(t2)∈span{1,cost,cos2t,⋯,cosnt}, |
since
sin2n(t2)=(sin2(t2))n=(1−cost2)n. |
Thus, there are n+1 real numbers d0,d1,⋯,dn such that
sin2n(t2)=d0+d1cost+d2cos2t+⋯+dncosnt. |
We deduce that, based on Euler's formula,
(eit2−e−it22i)2n=d0+d1eit+e−it2+d2e2it+e−2it2+⋯+dnenit+e−nit2. | (4.5) |
Expanding the left side of Eq (4.5) using the binomial theorem results in d0=(2n−1)!!(2n)!!.
Lemma 4.3. The function Gi,2n(t) defined in Eq (4.4) lies in Γ2n[T] for any i∈Z,n∈Z+.
Proof. First, according to Lemma 4.2, we obtain
sin2n−2(t2)∈span{1,cost,cos2t,⋯,cos(n−1)t}, |
Thus, we have
sintsin2n−2(t2)∈span{sint,sin2t,⋯,sinnt}, |
and
g2n(t−ti)∈T2n for i∈Z. |
Second, the order of continuity of g2n(t−ti) at ti is 2n−2, because that
g(l)2n(t−ti)|t=ti=0,l=0,1,⋯,2n−2, and g(2n−1)2n(t−ti)|t=ti≠0. |
Third, it is easy to see that the order of continuity of Gi,2n(t) is 2n−ξ−1 at ti and is ∞ at the other knots.
In conclusion, Gi,2n(t)∈Γ2n[T], for i∈Z.
Lemma 4.4. The functions Gi,2n(t) for i∈Z in Eq (4.4) are linearly independent.
Proof. If ti−1<ti=ti+1=⋯=ti+mi−1<ti+mi, the mi functions Gj,2n(t) for i≤j≤i+mi−1 are linearly independent since they have different continuous orders at ti.
For the sake of simplicity, let r=i+mi and tr−1<tr=tr+1=⋯=tr+mr−1<tr+mr. We can then similarly conclude that the functions Gj,2n(t) for r≤j≤r+mr−1 are linearly independent. In addition, it is straightforward to derive that
span{Gi,2n(t),Gi+1,2n(t),⋯,Gi+mi−1,2n(t)}∩span{Gr,2n(t),Gr+1,2n(t),⋯,Gr+mr−1,2n(t)}={0}. |
Thus, Gi,2n(t),Gi+1,2n(t),⋯,Gi+mi−1,2n(t),Gr,2n(t),Gr+1,2n(t),⋯,Gr+mr−1,2n(t) are linearly independent. Consequently, the functions Gi,2n(t) for i∈Z are linearly independent.
From the above lemmas, we conclude that any function in the space Γ2n[ti,ti+2n] can be expressed as a linear combination of Gi,2n(t) for i∈Z.
Theorem 4.1. For any function Fi,2n(t)∈Γ2n[ti,ti+2n], there exist 2n−η+1 real numbers νi,νi+1,⋯,νi+2n−η such that
Fi,2n(t)=i+2n−η∑j=iνjGj,2n(t), t∈[ti,ti+2n), |
where η is the multiplicity of ti+2n in the interval [ti,ti+2n] and the functions Gj,2n(t), where i≤j≤i+2n−η, are defined in Eq (4.4).
Proof. Since Fi,2n(t) in space Γ2n[ti,ti+2n] is a piecewise function, it can first be linearly represented by the functions Gj,2n(t) for i≤j≤i+2n−η in Eq (4.4) over a non-zero interval within its support interval.
Suppose ti=ti+1=⋯=ti+ξ−1<ti+ξ≤⋯≤ti+2n−η<ti+2n−η+1=⋯=ti+2n. We consider that the function Fi,2n(t) can be linearly represented in the interval [ti,ti+ξ). Thus, we prove that there exist ξ real numbers κi,κi+1,⋯,κi+ξ−1 such that
Fi,2n(t)|[ti,ti+ξ)=i+ξ−1∑j=iκjGj,2n(t). |
For simplicity, assume that
f1(t)=Fi,2n(t)|[ti,ti+ξ), t∈[ti,ti+ξ). |
Based on the definition of the space Γ2n[ti,ti+2n], we know that f1(t)∈T2n, which means that there exist 2n real numbers x1,x2,⋯,x2n such that
f1(t)=x1cost+x2sint+⋯+x2n−1cosnt+x2nsinnt, |
and
f(s)1(ti−)=f(s)1(ti)=f(s)1(ti+), 0≤s≤2n−ξ−1, |
which implies that
f(s)1(ti)=0, 0≤s≤2n−ξ−1. |
Let ρ=2n−ξ−1, Ωρ:={f1(t)|f1(t)∈T2n,f(s)1(ti)=0,0≤s≤ρ}, and Ψρ:=span{g2n(t−ti),g′2n(t−ti+1),⋯,g(2n−ρ−2)2n(t−ti+ξ−1)}=span{g2n(t−ti),g′2n(t−ti),⋯,g(2n−ρ−2)2n(t−ti)}.
Since g2n(t−ti),g′2n(t−ti),⋯,g(2n−ρ−2)2n(t−ti) are linearly independent, we conclude that the dimension of Ψρ is 2n−ρ−1. According to the definition of the function g2n(t), we obtain that g(l)2n(t−ti)∈Ωρ, l=0,1,⋯,2n−ρ−2. Hence, it follows naturally that Ψρ is the subspace of Ωρ.
The dimension of the space Ωρ is equal to the dimension of the solution space corresponding to the following linear equations.
{f1(ti)=0,f′1(ti)=0,⋮f(ρ)1(ti)=0. | (4.6) |
Thus, the following linear equations holds:
{x1costi+x2sinti+x3cos2ti+x4sin2ti+⋯+x2n−1cosnti+x2nsinnti=0,−x1sinti+x2costi−x32sin2ti+x42cos2ti+⋯−x2n−1nsinnti+x2nncosnti=0,−x1costi−x2sinti−x322cos2ti−x422sin2ti+⋯−x2n−1n2cosnti−x2nn2sinnti=0,x1sinti−x2costi+x323sin2ti−x423cos2ti+⋯+x2n−1n3sinnti−x2nn3cosnti=0,⋮x1cos(π2ρ+ti)+x2sin(π2ρ+ti)+⋯+x2n−1nρcos(π2ρ+nti)+x2nnρsin(π2ρ+nti)=0, | (4.7) |
where the corresponding coefficient matrix is given as
(costisinticos2tisin2ti…cosntisinnti−sinticosti−2sin2ti2cos2ti…−nsinntincosnti−costi−sinti−22cos2ti−22sin2ti…−n2cosnti−n2sinntisinti−costi23sin2ti−23cos2ti…n3sinnti−n3cosnti⋮⋮⋮⋮⋮⋮cos(π2ρ+ti)sin(π2ρ+ti)2ρcos(π2ρ+2ti)2ρsin(π2ρ+2ti)…nρcos(π2ρ+nti)nρsin(π2ρ+nti)). | (4.8) |
By performing elementary row operations, it is demonstrated that the matrix (4.8) maintains full row rank, independent of the parity of ρ. This result establishes that the dimension of the solution space for the linear system (4.7) is
2n−(ρ+1)=2n−ρ−1. |
The dimension of Ωρ is determined to be 2n−ρ−1. Consequently, it follows that Ωρ=Ψρ. In other words, we know that there exist 2n−ρ−1=2n−(2n−ξ−1)−1=ξ real numbers κi,κi+1,⋯,κi+ξ−1 such that
f1(t)=i+ξ−1∑j=iκjg(j−i)2n(t−tj). |
According to the definition of Gi,2n(t), the function \(f_1(t)\) can be expressed as
f1(t)=i+ξ−1∑j=iκjGj,2n(t). |
Consider the non-zero interval [ti+ξ,ti+ξ+mi+ξ). We define
f2(t)=Fi,2n(t)|[ti+ξ,ti+ξ+mi+ξ),t∈[ti+ξ,ti+ξ+mi+ξ). |
According to the continuous order of Fi,2n(t) at ti+ξ, we deduce that
f(s)1(ti+ξ)=f(s)1(ti+ξ−)=f(s)2(ti+ξ+)=f(s)2(ti+ξ),0≤s≤2n−mi+ξ−1. |
This implies that
f(s)2(ti+ξ)−f(s)1(ti+ξ)=0,0≤s≤2n−mi+ξ−1. |
Thus, from the above analysis, it follows that there are mi+ξ real numbers κi+ξ,κi+ξ+1,⋯,κi+ξ+mi+ξ−1 such that
f2(t)−f1(t)=i+ξ+mi+ξ−1∑j=i+ξκjGj,2n(t). |
Consequently, we have
f2(t)=i+ξ+mi+ξ−1∑j=iκjGj,2n(t). |
For every non-zero subinterval of the support interval of Fi,2n(t), we consider it this way. It can be concluded that there exist 2n−η+1 real numbers νi,νi+1,⋯,νi+2n−η such that
Fi,2n(t)=i+2n−η∑j=iνjGj,2n(t), t∈[ti,ti+2n−η+1)=[ti,ti+2n), |
where νj for i≤j≤i+2n−η is expressed as a linear combination of κi,κi+1,⋯,κi+2n−η.
The subsection demonstrates that the dimension of Γ2n[ti,ti+2n] is 1. To support this, we first require the following lemma.
Lemma 4.5. Given an integer i and positive integers n,η such that η≤2n−1, then the determinant
|Gi,2n(ti+2n)Gi+1,2n(ti+2n)…Gi+2n−η−1,2n(ti+2n)G′i,2n(ti+2n)G′i+1,2n(ti+2n)…G′i+2n−η−1,2n(ti+2n)⋮⋮⋮⋮G(2n−η−1)i,2n(ti+2n)G(2n−η−1)i+1,2n(ti+2n)…G(2n−η−1)i+2n−η−1,2n(ti+2n)|≠0, |
where the functions Gj,2n(t) for i≤j≤i+2n−η−1 are defined in Eq (4.4).
Proof. Use reduction to absurdity. Assume that the determinant is equal to zero. Then, there are 2n−η real numbers zi,zi+1,⋯,zi+2n−η−1, which are not all equal to zero, satisfying
i+2n−η−1∑l=izl(Gl,2n(ti+2n)G′l,2n(ti+2n)⋮G(2n−η−1)l,2n(ti+2n))=0. |
Let Y(t)=i+2n−η−1∑l=izlGl,2n(t). Then, we find that the first 2n−η−1 derivatives of Y(t) at ti+2n are all zero. According to Theorem 4.1, there exist mi+2n real numbers, such that
Y(t)=i+2n+mi+2n−η−1∑u=i+2n−ηzuGu,2n(t). |
Therefore,
Y(t)=i+2n−η−1∑l=izlGl,2n(t)=i+2n+mi+2n−η−1∑u=i+2n−ηzuGu,2n(t), |
and
i+2n−η−1∑l=izlGl,2n(t)−i+2n+mi+2n−η−1∑u=i+2n−ηzuGu,2n(t)=0. |
According to Lemma 4.4, we obtain that zi=⋯=zi+2n−η−1=0. This conflicts with the assumption. So, the lemma is proved.
From Theorem 4.1 and Lemma 4.5, we obtain the dimension of Γ2n[ti,ti+2n].
Theorem 4.2. The dimension of the linear space Γ2n[ti,ti+2n] is 1.
Proof. Suppose that u(t) is an arbitrary function in Γ2n[ti,ti+2n]. Thus, according to Theorem 4.1, there are 2n−η+1 real numbers νi,νi+1,⋯,νi+2n−η such that
u(t)=i+2n−η∑j=iνjGj,2n(t),t∈[ti,ti+2n)=[ti,ti+2n−η+1), |
where η denotes the multiplicity of ti+2n in the interval [ti,ti+2n]. Consider the continuous order of function u(t) at ti+2n. We have
u(l)(ti+2n−)=u(l)(ti+2n)=u(l)(ti+2n+)=0,l=0,1,2,⋯,2n−η−1, |
which implies that
u(l)(ti+2n)=0,l=0,1,2,⋯,2n−η−1. |
Thus, the equation representing the continuous order of u(t)s at ti+2n is given by
{u(ti+2n)=0,u′(ti+2n)=0,⋮u(2n−η−1)(ti+2n)=0. | (4.9) |
which can be expressed in matrix form as follows:
(Gi,2n(ti+2n)Gi+1,2n(ti+2n)…Gi+2n−η,2n(ti+2n)G′i,2n(ti+2n)G′i+1,2n(ti+2n)…G′i+2n−η,2n(ti+2n)⋮⋮⋮⋮G(2n−η−1)i,2n(ti+2n)G(2n−η−1)i+1,2n(ti+2n)…G(2n−η−1)i+2n−η,2n(ti+2n))(νiνi+1⋮νi+2n−η)=0. | (4.10) |
This system consists of linear equations with νi,νi+1,⋯,νi+2n−η as variables. Based on Lemma 4.5, the coefficient matrix of these linear equations is full row rank, with a rank of 2n−η+1−(2n−η)=1. This indicates that the dimension of Γ2n[ti,ti+2n] is 1.
In this subsection, we consider single knots and assume that ti<ti+1 for any i∈Z. Then the even-order trigonometric spline basis functions corresponding to the DC component-free space T2n, and the integral expression for the trigonometric B-spline basis Ki,2n+1(t) are presented.
According to Theorem 4.2 and the definition of Γ2n[ti,ti+2n], it is established that Fi,2n(t)∈Γ2n[ti,ti+2n] and the dimension of Γ2n[ti,ti+2n] is 1. If we find a function H(t)∈Γ2n[ti,ti+2n], then it follows that Fi,2n(t)=αH(t),t∈[ti,ti+2n) for some real number α. Thus, the following theorem is provided.
Theorem 4.3. (The function expression in the space Γ2n[ti,ti+2n] over single knots) For any function Fi,2n(t)∈Γ2n[ti,ti+2n], there exists a real number α such that
Fi,2n(t)=αH(t),t∈[ti,ti+2n), |
where
H(t)=|Gi,2n(t)Gi+1,2n(t)…Gi+2n−1,2n(t)Gi+2n,2n(t)costicosti+1…costi+2n−1costi+2nsintisinti+1…sinti+2n−1sinti+2ncos2ticos2ti+1…cos2ti+2n−1cos2ti+2nsin2tisin2ti+1…sin2ti+2n−1sin2ti+2n⋮⋮⋮⋮⋮cosnticosnti+1…cosnti+2n−1cosnti+2nsinntisinnti+1…sinnti+2n−1sinnti+2n|, |
and α=−1D(ti,ti+1,⋯,ti+2n−1). Here D(ti,ti+1,⋯,ti+2n−1) and the functions Gj,2n(t) for i≤j≤i+2n are defined in Eqs (4.1) and (4.4), respectively.
Proof. According to Theorem 4.1, there exist 2n real numbers αi,αi+1,⋯,αi+2n−1 such that
Fi,2n(t)=i+2n−1∑j=iαjGj,2n(t),t∈[ti,ti+2n), | (4.11) |
By the continuous order of Fi,2n(t) at ti+2n, it follows that
F(d)i,2n(ti+2n)=i+2n−1∑j=iαjG(d)j,2n(ti+2n)=0,t∈[ti,ti+2n),d=0,1,2,⋯,2n−2. | (4.12) |
This can be expressed equivalently as a system of equations
{i+2n−1∑j=iαjGj,2n(ti+2n)=0,i+2n−1∑j=iαjG′j,2n(ti+2n)=0,⋮i+2n−1∑j=iαjG(2n−2)j,2n(ti+2n)=0. | (4.13) |
In addition, we have
sin(t−ti+2n)sin2n−2(t−ti+2n2)=i+2n−1∑j=iβjsin(t−tj)sin2n−2(t−tj2), | (4.14) |
where βi,βi+1,⋯,βi+2n−1 are real numbers. Since Eq (4.14) when t=ti+2n is equivalent to Eq (4.13). We obtain that αj=βj,j=i,⋯,i+2n−1. Thus, we only focus on βj,i≤j≤i+2n−1. By proving Lemma 4.3, it can be concluded that there exist n numbers η1,η2,⋯,ηn such that
sin(t−tj)sin2n−2(t−tj2)=n∑l=1ηlsinl(t−tj),j=i,i+1,⋯,i+2n. | (4.15) |
So we can rewrite (4.14) as follows:
(costsint…cosntsinnt)(−η1sinti+2nη1costi+2n−η2sin2ti+2nη2cos2ti+2n⋮−ηnsinnti+2nηncosnti+2n)=(costsint…cosntsinnt)(−η1sinti…−η1sinti+2n−1η1costi…η1costi+2n−1−η2sin2ti…−η2sin2ti+2n−1η2cos2ti…η2cos2ti+2n−1⋮⋮⋮−ηnsinnti…−ηnsinnti+2n−1ηncosnti…ηncosnti+2n−1)(βiβi+1⋮βi+2n−2βi+2n−1). | (4.16) |
Based on the properties of matrix operations, we deduce that
(−η1sinti…−η1sinti+2n−1η1costi…η1costi+2n−1−η2sin2ti…−η2sin2ti+2n−1η2cos2ti…η2cos2ti+2n−1⋮⋮⋮−ηnsinnti…−ηnsinnti+2n−1ηncosnti…ηncosnti+2n−1)(βiβi+1⋮βi+2n−2βi+2n−1)=(−η1sinti+2nη1costi+2n−η2sin2ti+2nη2cos2ti+2n⋮−ηnsinnti+2nηncosnti+2n). | (4.17) |
The coefficient matrix of (4.17) is non-zero due to the linear independence of the functions cost,sint,⋯, cosnt,sinnt. According to Cramer's Rule, we have
βi+j−1=D(ti,⋯,ti+j−2,ti+2n,ti+j,⋯,ti+2n−1)D(ti,ti+1,⋯,ti+2n−1), for j=1,2,⋯,2n. |
We apply this result to Eq (4.11) to yield
Fi,2n(t)=2n∑j=1αi+j−1Gi+j−1,2n(t)=2n∑j=1D(ti,⋯,ti+j−2,ti+2n,ti+j,⋯,ti+2n−1)D(ti,ti+1,⋯,ti+2n−1)Gi+j−1,2n(t)=D(ti+2n,ti+1,⋯,ti+2n−1)D(ti,ti+1,⋯,ti+2n−1)Gi,2n(t)+D(ti,ti+2n,ti+2,⋯,ti+2n−1)D(ti,ti+1,⋯,ti+2n−1)Gi+1,2n(t)+⋯+D(ti,⋯,ti+2n−2,ti+2n)D(ti,ti+1,⋯,ti+2n−1)Gi+2n−1,2n(t)=(−1)2n−1D(ti+1,⋯,ti+2n)D(ti,ti+1,⋯,ti+2n−1)Gi,2n(t)+(−1)2n−2D(ti,ti+2,⋯,ti+2n)D(ti,ti+1,⋯,ti+2n−1)Gi+1,2n(t)+⋯+D(ti,⋯,ti+2n−2,ti+2n)D(ti,ti+1,⋯,ti+2n−1)Gi+2n−1,2n(t)+D(ti,⋯,ti+2n−2,ti+2n−1)D(ti,ti+1,⋯,ti+2n−1)Gi+2n,2n(t)=−1D(ti,ti+1,⋯,ti+2n−1)H(t), t∈[ti,ti+2n), |
where Gi+2n,2n(t)=0,t∈[ti,ti+2n).
The linear space Γ2n[ti,ti+2n] corresponds to space T2n. Similarly, the space Γ2n+1[ti,ti+2n+1], which corresponds to space T2n+1, can be defined as follows:
Γ2n+1[ti,ti+2n+1]:={Mi,2n+1(t)∈Γ2n+1[T]|Mi,2n+1(t)=0,t∉[ti,ti+2n+1] and Mi,2n+1(t)≠0,t∈(ti,ti+2n+1), M(l)i,2n+1(tj−)=M(l)i,2n+1(tj+),l=0,1,⋯,kji,2n+1, i≤j≤i+2n+1.} |
We can derive the following theorem analogously.
Theorem 4.4. The dimension of the linear space Γ2n+1[ti,ti+2n+1] is 1.
The proof is similar to the proof of Theorem 4.2.
Inspired by the method for constructing normalized B-basis in extended Chebyshev space presented in [10], we use the functions Fi,2n(t), where i∈Z, to construct the normalized function Mi,2n+1(t). Therefore, the following theorem holds.
Theorem 4.5. (The normalized function) Suppose that
Mi,2n+1(t)=∫t−∞Fi,2n(s)ds∫+∞−∞Fi,2n(t)dt−∫t−∞Fi+1,2n(s)ds∫+∞−∞Fi+1,2n(t)dt, | (4.18) |
where Fi,2n(t)∈Γ2n[ti,ti+2n] and Fi+1,2n(t)∈Γ2n[ti+1,ti+2n+1] are defined in Theorem 4.3, with i∈Z and n∈Z+. Then we have
+∞∑i=−∞Mi,2n+1(t)≡1, t∈[tj,tj+1). |
Proof. For t∈[tj,tj+1), there exists
+∞∑i=−∞Mi,2n+1(t)=+∞∑i=−∞(∫t−∞Fi,2n(s)ds∫+∞−∞Fi,2n(t)dt−∫t−∞Fi+1,2n(s)ds∫+∞−∞Fi+1,2n(t)dt)=j∑i=j−2n(∫t−∞Fi,2n(s)ds∫+∞−∞Fi,2n(t)dt−∫t−∞Fi+1,2n(s)ds∫+∞−∞Fi+1,2n(t)dt)=j∑i=j−2n(∫ttiFi,2n(s)ds∫ti+2ntiFi,2n(t)dt−∫tti+1Fi+1,2n(s)ds∫ti+2n+1ti+1Fi+1,2n(t)dt)=∫ttj−2nFj−2n,2n(s)ds∫tjtj−2nFj−2n,2n(t)dt−∫ttj+1Fj+1,2n(s)ds∫tj+2n+1tj+1Fj+1,2n(t)dt=∫tjtj−2nFj−2n,2n(s)ds+∫ttjFj−2n,2n(s)ds∫tjtj−2nFj−2n,2n(t)dt−0=1+∫ttjFj−2n⋅2n(s)ds∫tjtj−2nFj−2n⋅2n(t)dt−0=1. |
Lemma 4.6. The function Mi,2n+1(t) in Eq (4.18) lies in the space Γ2n+1[ti,ti+2n+1].
Proof. First, Eq (4.18) indicates that the support interval of the function Mi,2n+1(t) is the union of the interval of Fi,2n(t) and Fi+1,2n(t), denoted as [ti,ti+2n+1).
Second, since the integral operator increases the continuous order by 1, the continuous order of Mi,2n+1(t) at tj is greater than or equal to
kji,2n+1=kji,2n+1,i≤j≤i+2n+1. |
So, it is natural that
Mi,2n+1(t)∈Γ2n+1[ti,ti+2n+1]. |
According to Theorem 4.4 and Lemma 4.6, the normalized function Mi,2n+1(t) in Eq (4.18) must be equal to the trigonometric B-spline function Ni,2n+1(t) in Eq (2.5) multiplied by a constant. Therefore, the following theorem is established.
Theorem 4.6. (Integral representation of the normalized trigonometric B-spline basis) Given a knot sequence T={ti}+∞i=−∞ satisfying
ti<ti+1, 0<ti+2n−ti<2π and 0≤ti+2n−1−ti<π, i∈Z,n∈Z+, |
then there holds
Ki,2n+1(t)=Mi,2n+1(t), t∈[ti,ti+2n+1), |
where Ki,2n+1(t) and Mi,2n+1(t) are separately defined in Definition 3.1 and Eq (4.18).
Proof. We will demonstrate that the expression of Mi,2n+1(t), as defined in Eq (4.18), is identical to that of Ki,2n+1(t) defined in Definition 3.1. The notations C(ti+1,ti+2,⋯,ti+2n) in Lemma 4.1 and U(ti,ti+1,⋯,ti+2n) (see [17]) are used in the following proof, where
U(ti,ti+1,⋯,ti+2n)=|11…11costicosti+1…costi+2n−1costi+2nsintisinti+1…sinti+2n−1sinti+2ncos2ticos2ti+1…cos2ti+2n−1cos2ti+2nsin2tisin2ti+1…sin2ti+2n−1sin2ti+2n⋮⋮⋮⋮⋮cosnticosnti+1…cosnti+2n−1cosnti+2nsinntisinnti+1…sinnti+2n−1sinnti+2n|=22n2∏i≤l<j≤i+2nsintj−tl2. | (4.19) |
Based on Definition 2.1, it suffices to prove the explicit expression of the function Mi,2n+1(t) over a non-zero subinterval within its support interval. Therefore, according to Lemmas 3.1, 4.1, and 4.2, and Theorems 4.3 and 4.5, we deduce that
Mi,2n+1(t)=∫ttiFi,2n(s)ds∫ti+2ntiFi,2n(t)dt=D(ti+1,ti+2,⋯,ti+2n)(2n−1)!!(2n)!!U(ti,ti+1,⋯,ti+2n)sin2n(t−ti2)=22n2−nn!∏i+1≤l<j≤i+2nsintj−tl2C(ti+1,ti+2,⋯,ti+2n)(2n−1)!!(2n)!!22n2∏i≤l<j≤i+2nsintj−tl2sin2n(t−ti2)=C(ti+1,ti+2,⋯,ti+2n)(2n−1)!!∏j=1,2,⋯,2nsinti+j−ti2sin2n(t−ti2)=Ci,2n+1∏j=1,2,⋯,2nsinti+j−ti2sin2n(t−ti2), t∈[ti,ti+1). |
Since
Ni,2n+1(t)=sin2n(t−ti2)∏j=1,2,⋯,2nsinti+j−ti2, t∈[ti,ti+1), |
we have
Mi,2n+1(t)=Ci,2n+1Ni,2n+1(t), t∈[ti,ti+1), |
consequently
Mi,2n+1(t)=Ci,2n+1Ni,2n+1(t)=Ki,2n+1(t), t∈[ti,ti+2n+1). |
In this subsection, we consider multiple knots and assume that the multiplicity of the knot ti in the interval [ti,ti+2n) is ξ, while the multiplicity of the knot ti+2n in the same interval is η. Similarly to the single knot case, there exist the following theorems.
Theorem 4.7. (The function expression in the space Γ2n[ti,ti+2n] over multiple knots) Let
Au,v=((cos(t+π2u))|t=tv, (sin(t+π2u))|t=tv, ⋯,nu(cos(nt+π2u))|t=tv, nu(sin(nt+π2u))|t=tv )T, |
Bu,v=((−1)u−1Au−1,v, (−1)u−2Au−2,v, ⋯,(−1)1A1,v, (−1)0A0,v ), |
Eu,v=(G(u−1)v,2n(t), G(u−2)v,2n(t), ⋯,G′v,2n(t), Gv,2n(t) ). |
For any function Fi,2n(t)∈Γ2n[ti,ti+2n], there exists a real number α such that
Fi,2n(t)=αH(t),t∈[ti,ti+2n), |
where
α=(−1)η|Bξ,iBmi+ξ,i+ξBmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Bmi+2n−η,i+2n−ηBη−1,i+2n|, |
and
H(t)=|Eξ,iEmi+ξ,i+ξEmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Emi+2n−η,i+2n−ηEη,i+2nBξ,iBmi+ξ,i+ξBmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Bmi+2n−η,i+2n−ηBη,i+2n|. |
Here the functions Gj,2n(t) for i≤j≤i+2n are defined in Eq (4.4).
Proof. To simplify the notation, we define
φu,v=(sin(t−tv)sin2n−2(t−tv2))(u). |
Thus, we conclude that
φη−1,i+2n=βiφξ−1,i+βi+1φξ−2,i+⋯+βi+ξ−1φ0,i+βi+ξφmi+ξ−1,i+ξ+⋯+βi+ξ+mi+ξ−1φ0,i+ξ+⋯+βi+2n+1−η−mi+2n−ηφmi+2n−η−1,i+2n−η+⋯+βi+2n−ηφ0,i+2n−η+βi+2n−η+1φ0,i+2n+⋯+βi+2n−1φη−2,i+2n. | (4.20) |
Similar to Theorem 4.3, there exist
Fi,2n(t)=i+2n−η∑j=iαjGj,2n(t),t∈[ti,ti+2n)=[ti,ti+2n+1−η), | (4.21) |
F(d)i,2n(ti+2n)=i+2n−η∑j=iαjG(d)j,2n(ti+2n)=0,t∈[ti,ti+2n+1−η),d=0,1,2,⋯,2n−1−η. | (4.22) |
Additionally, it follows that Eq (4.20) is equivalent to Eq (4.22) when t=ti+2n. Thus, by applying the proof strategy from Theorem 4.3, we derive that
Fi,2n(t)=(−1)η|Eξ,iEmi+ξ,i+ξEmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Emi+2n−η,i+2n−ηEη,i+2nBξ,iBmi+ξ,i+ξBmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Bmi+2n−η,i+2n−ηBη,i+2n||Bξ,iBmi+ξ,i+ξBmi+ξ+mi+ξ,i+ξ+mi+ξ⋯Bmi+2n−η,i+2n−ηBη−1,i+2n|. |
Theorem 4.8. (The normalized function over generalized knots) Suppose that Eq (4.18) still holds, where Fi,2n(t)∈Γ2n[ti,ti+2n] and Fi+1,2n(t)∈Γ2n[ti+1,ti+2n+1] are defined in Theorem 4.7, with i∈Z and n∈Z+, then
+∞∑i=−∞Mi,2n+1(t)≡1, t∈[tj,tj+1). |
In addition, when Fi,2n(t)=0, we set
∫t−∞(∫+∞−∞Fi,2n(s)ds)−1Fi,2n(s)ds={0t<ti,1t≥ti. |
Theorem 4.9. (Integral representation of the normalized trigonometric B-spline basis over generalized knots) Given a knot sequence T={ti}+∞i=−∞ satisfying condition (3.4), then there holds
Ki,2n+1(t)=Mi,2n+1(t), t∈[ti,ti+2n+1), |
where Ki,2n+1(t) and Mi,2n+1(t) are separately defined in Definition 3.1 and Theorem 4.8.
This subsection presents examples of curve modeling to demonstrate that a curve possesses the convex hull property when the knot sequence satisfies condition (3.4).
Figures 3 and 4 illustrate examples of open and closed curves over different knot sequences, respectively. The red curves correspond to the knot sequences that satisfy condition (3.4) as defined in Definition 3.1, while the orange curves correspond to the knot sequences that satisfy condition (2.1) as defined in Definition 2.1. Additionally, the knot sequences in Figure 3 are T={0,2,3,3,3.5,5.3,6.1,6.6,8.4,9.1,9.5,11.2,12.2,12.6,14.2,15}, and T={0,0.8,3,3,3.5,5.3,6.6,8.6,9.1,9.1,10,12.2,12.2,13,14.2,15}, while those in Figure 4 are T={−8,−8,−8,−8,−8,−7,−6.3,−5.6,−4.9,−4.5,−1.2,−1.2,−1.2,0,1.2,1.2, 1.2,4.5,4.9,5.6,6.3,7,8,8,8,8,8}, and T={−8,−8,−8,−8,−8,−7,−6.3,−5.3,−4,−3.3,−1.2,−1.2, −1.2,0,1.2,1.2,1.2,3.3,4,5.3,6.3,7,8,8,8,8,8}. Clearly, if the knot sequence only satisfies condition (2.1), the convex hull property of the curve cannot be guaranteed.
In the Chebyshev system, due to the integral properties of sine and cosine in the trigonometric B-spline basis, the trigonometric B-spline basis cannot be directly derived from lower-order bases through integration, which leads to an unnormal Chebyshev system. This paper successfully derives the integral formula for the normalized odd-order trigonometric B-spline basis by constructing a new set of even-order trigonometric B-spline bases. This integral formula allows for the transition from even-order to odd-order bases but cannot be obtained through stepwise integration, indicating that it is only similar to a segment of the integral formula in the Chebyshev system. Although we are currently unable to provide a direct recursive formula for integrating from lower-order trigonometric spline bases to higher-order ones, we hope to use this as a foundation for further exploration of this issue in future work.
Mei Li: Investigation, methodology, software, validation, writing—original draft preparation, writing—review and editing; Wanqiang Shen: Conceptualization, methodology, software, writing—original draft preparation, writing—review and editing. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research is supported in part by the National Natural Science Foundation of China (Grant No. 61772013).
The authors declare that there are no conflicts of interest regarding the publication of this paper.
[1] |
B. Manvi, J. Tawade, M. Biradar, S. Noeiaghdam, U. Fernandez-Gamiz, V. Govindan, The effects of MHD radiating and non-uniform heat source/sink with heating on the momentum and heat transfer of Eyring-Powell fluid over a stretching, Results Eng., 14 (2022), 100435. https://doi.org/10.1016/j.rineng.2022.100435 doi: 10.1016/j.rineng.2022.100435
![]() |
[2] |
H. Shahzad, Q. U. Ain, A. A. Pasha, K. Irshad, I. A. Shah, A. Ghaffari, et al., Double-diffusive natural convection energy transfer in magnetically influenced Casson fluid flow in trapezoidal enclosure with fillets, Int. Commun. Heat Mass, 137 (2022), 106236. https://doi.org/10.1016/j.icheatmasstransfer.2022.106236 doi: 10.1016/j.icheatmasstransfer.2022.106236
![]() |
[3] |
A. A. Pasha, M. M. Alam, T. Tayebi, S. Kasim, A. S. Dogonchi, K. Irshad, et al., Heat transfer and irreversibility evaluation of non-Newtonian nanofluid density-driven convection within a hexagonal-shaped domain influenced by an inclined magnetic field, Case Stud. Therm. Eng., 41 (2023), 102588. https://doi.org/10.1016/j.csite.2022.102588 doi: 10.1016/j.csite.2022.102588
![]() |
[4] |
C. N. Guled, J. V. Tawade, P. Kumam, S. Noeiaghdam, I. Maharudrappa, S. M. Chithra, et al., The heat transfer effects of MHD slip flow with suction and injection and radiation over a shrinking sheet by optimal homotopy analysis method, Results Eng., 18 (2023), 101173. https://doi.org/10.1016/j.rineng.2023.101173 doi: 10.1016/j.rineng.2023.101173
![]() |
[5] |
G. Dharmaiah, J. L. R. Prasad, K. S. Balamurugan, I. Nurhidayat, U. Fernandez-Gamiz, S. Noeiaghdam, Performance of magnetic dipole contribution on ferromagnetic non-Newtonian radiative MHD blood flow: An application of biotechnology and medical sciences, Heliyon, 2 (2023). https://doi.org/10.1016/j.heliyon.2023.e13369 doi: 10.1016/j.heliyon.2023.e13369
![]() |
[6] |
T. Sajid, W. Jamshed, M. R. Eid, S. Algarni, T. Alqahtani, R. W. Ibrahim, et al., Thermal case examination of inconstant heat source (sink) on viscous radiative sutterby nanofluid flowing via a penetrable rotative cone, Case Stud. Therm. Eng., 48 (2023), 103102. https://doi.org/10.1016/j.csite.2023.103102 doi: 10.1016/j.csite.2023.103102
![]() |
[7] |
M. Waqas, Y. J. Xu, M. Nasir, M. M. Alam, A. A. Pasha, K. Irshad, et al., Darcy-Forchheimer mangetized flow based on differential type nanoliquid capturing Ohmic dissipation effects, Propuls. Power Res., 12 (2023), 443−455. https://doi.org/10.1016/j.jppr.2023.08.003 doi: 10.1016/j.jppr.2023.08.003
![]() |
[8] |
N. S. Akbar, S. Nadeem, T. Hayat, A. A. Hendi, Effects of heat and mass transfer on the peristaltic flow of hyperbolic tangent fluid in an annulus, Int. J. Heat Mass., 54 (2011), 4360−4369. https://doi.org/10.1016/j.ijheatmasstransfer.2011.03.064 doi: 10.1016/j.ijheatmasstransfer.2011.03.064
![]() |
[9] |
S. Akram, S. Nadeem, Consequence of nanofluid on peristaltic transport of a hyperbolic tangent fluid model in the occurrence of apt (tending) magnetic field, J. Magn. Magn. Mater., 358 (2014), 183−191. https://doi.org/10.1016/j.jmmm.2014.01.052 doi: 10.1016/j.jmmm.2014.01.052
![]() |
[10] |
M. Naseer, M. Y. Malik, S. Nadeem, A. Rehman, The boundary layer flow of hyperbolic tangent fluid over a vertical exponentially stretching cylinder, Alex. Eng. J., 53 (2014), 747−750. https://doi.org/10.1016/j.aej.2014.05.001 doi: 10.1016/j.aej.2014.05.001
![]() |
[11] |
S. A. Gaffar, V. R. Prasad, O. A. Bég, Numerical study of flow and heat transfer of non-Newtonian tangent hyperbolic fluid from a sphere with Biot number effects, Alex. Eng. J., 54 (2015), 829−841. https://doi.org/10.1016/j.aej.2015.07.001 doi: 10.1016/j.aej.2015.07.001
![]() |
[12] |
T. Hayat, M. Shafique, A. Tanveer, A. Alsaedi, Magnetohydrodynamic effects on peristaltic flow of hyperbolic tangent nanofluid with slip conditions and Joule heating in an inclined channel, Int. J. Heat Mass, 102 (2016), 54−63. https://doi.org/10.1016/j.ijheatmasstransfer.2016.05.105 doi: 10.1016/j.ijheatmasstransfer.2016.05.105
![]() |
[13] |
T. Hayat, S. Qayyum, A. Alsaedi, S. A. Shehzad, Nonlinear thermal radiation aspects in stagnation point flow of tangent hyperbolic nanofluid with double diffusive convection, J. Mol. Liq., 223 (2016), 969−978. https://doi.org/10.1016/j.molliq.2016.08.102 doi: 10.1016/j.molliq.2016.08.102
![]() |
[14] |
K. U. Rehman, A. A. Malik, M. Y. Malik, N. U. Saba, Mutual effects of thermal radiations and thermal stratification on tangent hyperbolic fluid flow yields by both cylindrical and flat surfaces, Case Stud. Therm. Eng., 10 (2017), 244−254. https://doi.org/10.1016/j.csite.2017.07.003 doi: 10.1016/j.csite.2017.07.003
![]() |
[15] |
K. G. Kumar, B. J. Gireesha, M. R. Krishanamurthy, N. G. Rudraswamy, An unsteady squeezed flow of a tangent hyperbolic fluid over a sensor surface in the presence of variable thermal conductivity, Results Phys., 7 (2017), 3031−3036. https://doi.org/10.1016/j.rinp.2017.08.021 doi: 10.1016/j.rinp.2017.08.021
![]() |
[16] |
V. Nagendramma, A. Leelarathnam, C. S. K. Raju, S. A. Shehzad, T. Hussain, Doubly stratified MHD tangent hyperbolic nanofluid flow due to permeable stretched cylinder, Results Phys., 9 (2018), 23−32. https://doi.org/10.1016/j.rinp.2018.02.019 doi: 10.1016/j.rinp.2018.02.019
![]() |
[17] |
S. M. Atif, S. Hussain, M. Sagheer, Heat and mass transfer analysis of time-dependent tangent hyperbolic nanofluid flow past a wedge, Phys. Lett. A, 383 (2019), 1187−1198. https://doi.org/10.1016/j.physleta.2019.01.003 doi: 10.1016/j.physleta.2019.01.003
![]() |
[18] |
Z. Ullah, G. Zaman, A. Ishak, Magnetohydrodynamic tangent hyperbolic fluid flow past a stretching sheet, Chinese J. Phys., 66 (2020), 258−268. https://doi.org/10.1016/j.cjph.2020.04.011 doi: 10.1016/j.cjph.2020.04.011
![]() |
[19] |
W. Khan, I. A. Badruddin, A. Ghaffari, H. M. Ali, Heat transfer in steady slip flow of tangent hyperbolic fluid over the lubricated surface of a stretchable rotatory disk, Case Studi. Therm. Eng., 24 (2021), 100825. https://doi.org/10.1016/j.csite.2020.100825 doi: 10.1016/j.csite.2020.100825
![]() |
[20] |
S. Sindhu, B. J. Gireesha, Scrutinization of unsteady non-Newtonian fluid flow considering buoyancy effect and thermal radiation: Tangent hyperbolic model, Int. Commun. Heat Mass, 135 (2022), 106062. https://doi.org/10.1016/j.icheatmasstransfer.2022.106062 doi: 10.1016/j.icheatmasstransfer.2022.106062
![]() |
[21] |
A. Hussain, N. Farooq, A. Ahmad, L. Sarwar, Impact of double diffusivity on the hyperbolic tangent model conveying nano fluid flow over the wedge, Int. Commun. Heat Mass, 145 (2023), 106849. https://doi.org/10.1016/j.icheatmasstransfer.2023.106849 doi: 10.1016/j.icheatmasstransfer.2023.106849
![]() |
[22] |
M. A. Elogail, Peristaltic flow of a hyperbolic tangent fluid with variable parameters, Results Eng., 17 (2023), 100955. https://doi.org/10.1016/j.rineng.2023.100955 doi: 10.1016/j.rineng.2023.100955
![]() |
[23] |
N. S. Akbar, S. Nadeem, R. U. Haq, Z. H. Khan, Numerical solutions of magnetohydrodynamic boundary layer flow of tangent hyperbolic fluid towards a stretching sheet, Indian J. Phys., 87 (2013), 1121−1124. https://doi.org/10.1007/s12648-013-0339-8 doi: 10.1007/s12648-013-0339-8
![]() |
[24] |
K. U. Rehman, N. U. Saba, M. Y. Malik, A. A. Malik, Encountering heat and mass transfer mechanisms simultaneously in Powell-Erying fluid through Lie symmetry approach, Case Stud. Therm. Eng., 10 (2017), 541−549. https://doi.org/10.1016/j.csite.2017.10.011 doi: 10.1016/j.csite.2017.10.011
![]() |
[25] |
K. U. Rehman, M. Y. Malik, I. Zehra, M. S. Alqarni, Group theoretical analysis for MHD flow fields: A numerical result, J. Braz. Soc. Mech. Sci., 41 (2019), 1−9. https://doi.org/10.1007/s40430-019-1662-6 doi: 10.1007/s40430-019-1662-6
![]() |
[26] |
Z. Ullah, G. Zaman, A. Ishak, Magnetohydrodynamic tangent hyperbolic fluid flow past a stretching sheet, Chinese J. Phys., 66 (2020), 258−268. https://doi.org/10.1016/j.cjph.2020.04.011 doi: 10.1016/j.cjph.2020.04.011
![]() |
[27] |
H. Sadaf, Z. Asghar, N. Iftikhar, Cilia-driven flow analysis of cross fluid model in a horizontal channel, Comput. Part. Mech., 9 (2022), 1−8. https://doi.org/10.1007/s40571-022-00539-w doi: 10.1007/s40571-022-00539-w
![]() |
[28] |
S. Arulmozhi, K. Sukkiramathi, S. S. Santra, R. Edwan, U. Fernandez-Gamiz, S. Noeiaghdam, Heat and mass transfer analysis of radiative and chemical reactive effects on MHD nanofluid over an infinite moving vertical plate, Results Eng., 14 (2022), 100394. https://doi.org/10.1016/j.rineng.2022.100394 doi: 10.1016/j.rineng.2022.100394
![]() |
[29] |
Y. Nawaz, M. S. Arif, K. Abodayeh, M. Mansoor, Finite difference schemes for MHD mixed convective Darcy-Forchheimer flow of Non-Newtonian fluid over oscillatory sheet: A computational study, Front. Phys., 11 (2023), 16. https://doi.org/10.3389/fphy.2023.1072296 doi: 10.3389/fphy.2023.1072296
![]() |
[30] |
Y. Nawaz, M. S. Arif, K. Aboda, Predictor-corrector scheme for electrical magnetohydrodynamic (MHD) Casson nanofluid flow: A computational study, Appl. Sci., 13 (2023), 1209. https://doi.org/10.3390/app13021209 doi: 10.3390/app13021209
![]() |
[31] |
C. S. Liu, The Lie-group shooting method for boundary-layer problems with suction/injection/reverse flow conditions for power-law fluids, Int. J. Nonlin. Mech., 46 (2011), 1001−1008. https://doi.org/10.1016/j.ijnonlinmec.2011.04.016 doi: 10.1016/j.ijnonlinmec.2011.04.016
![]() |
[32] | W. M. K. A. D. Zaimi, B. Bidin, N. A. A. Bakar, R. A. Hamid, Applications of Runge-Kutta-Fehlberg method and shooting technique for solving classical Blasius equation, World Appl. Sci. J., 17 (2012), 10−15. |
[33] |
K. U. Rehman, W. Shatanawi, U. Firdous, A comparative thermal case study on thermophysical aspects in thermally magnetized flow regime with variable thermal conductivity, Case Stud. Therm. Eng., 44 (2023), 102839. https://doi.org/10.1016/j.csite.2023.102839 doi: 10.1016/j.csite.2023.102839
![]() |
[34] | A. Usman, M. Rafiq, M. Saeed, A. Nauman, A. Almqvist, M. Liwicki, Machine learning computational fluid dynamics, In: 2021 Swedish Artificial Intelligence Society Workshop (SAIS), IEEE, Sweden, 2021. https://doi.org/10.1109/SAIS53221.2021.9483997 |
[35] |
D. Drikakis, F. Sofos, Can artificial intelligence accelerate fluid mechanics research?, Fluids, 8 (2023), 212. https://doi.org/10.3390/fluids8070212 doi: 10.3390/fluids8070212
![]() |
[36] |
Z. Said, P. Sharma, R. M. Elavarasan, A. K. Tiwari, M. K. Rathod, Exploring the specific heat capacity of water-based hybrid nanofluids for solar energy applications: A comparative evaluation of modern ensemble machine learning techniques, J. Energy Storage, 54 (2022), 105230. https://doi.org/10.1016/j.est.2022.105230 doi: 10.1016/j.est.2022.105230
![]() |
[37] |
A. Shafiq, A. B. Çolak, T. N. Sindhu, T. Muhammad, Optimization of Darcy-Forchheimer squeezing flow in nonlinear stratified fluid under convective conditions with artificial neural network, Heat Transf. Res., 53 (2022). https://doi.org/10.1615/HeatTransRes.2021041018 doi: 10.1615/HeatTransRes.2021041018
![]() |
[38] |
K. U. Rehman, W. Shatanawi, Non-Newtonian mixed convection magnetized flow with heat generation and viscous dissipation effects: A prediction application of artificial intelligence, Processes, 11 (2023), 986. https://doi.org/10.3390/pr11040986 doi: 10.3390/pr11040986
![]() |
[39] |
M. Adamu, A. B. Çolak, Y. E. Ibrahim, S. I. Haruna, M. F. Hamza, Prediction of mechanical properties of rubberized concrete incorporating fly ash and nano silica by artificial neural network technique, Axioms, 12 (2023), 81. https://doi.org/10.3390/axioms12010081 doi: 10.3390/axioms12010081
![]() |
[40] |
M. Fathizadeh, M. Madani, Y. Khan, N. Faraz, A. Yıldırım, S. Tutkun, An effective modification of the homotopy perturbation method for MHD viscous flow over a stretching sheet, J. King Saud Univ.-Sci., 25 (2013), 107−113. https://doi.org/10.1016/j.jksus.2011.08.003 doi: 10.1016/j.jksus.2011.08.003
![]() |