Superellipse (Lamé Curve)
Total Page:16
File Type:pdf, Size:1020Kb
6 Superellipse (Lamé curve) 6.1 Equations of a superellipse A superellipse (horizontally long) is expressed as follows. Implicit Equation x n y n + = 10<b a , n > 0 (1.1) a b Explicit Equation 1 n x n y = b 1 - 0<b a , n > 0 (1.1') a When a =3, b=2 , the superellipses for n = 9.9 , 2, 1 and 2/3 are drawn as follows. Fig1 These are called an ellipse when n= 2, are called a diamond when n = 1, and are called an asteroid when n = 2/3. These are known well. parametric representation of an ellipse In order to ask for the area and the arc length of a super-ellipse, it is necessary to calculus the equations. However, it is difficult for (1.1) and (1.1'). Then, in order to make these easy, parametric representation of (1.1) is often used. As far as the 1st quadrant, (1.1) is as follows . xn yn 0xa + = 1 a n b n 0yb Let us compare this with the following trigonometric equation. sin 2 + cos 2 = 1 Then, we obtain xn yn = sin 2 , = cos 2 a n b n i.e. - 1 - 2 2 x = a sin n ,y = bcos n The variable and the domain are changed as follows. x: 0a : 0/2 By this representation, the area and the arc length of a superellipse can be calculated clockwise from the y-axis. cf. If we represent this as usual as follows, 2 2 x = a cos n ,y = bsin n the area and the arc length of a superellipse can be calculated counterclockwise from the x-axis. However, from this representation, we cannot obtain the Legendre forms of elliptic integrals - 2 - 6.2 Area of a superellipse Formula 6.2.1 When n,a,b ()b a are positive numbers respectively and ()z is the gamma function, the area S of the ellipse of degree n is given by the following expression. 1 2 1+ n S = 4ab 2 1+ n Proof In order to obtain the area of a superellipse, we integrate with the following equation in the 1st quadrant, and should just multiply it by 4. 1 n xn y = b 1- 0<b a , n > 0 a n That is, 1 n a xn S = 4b 1- n dx 0 a However, it is known that this integration cannot be expressed with elementary functions. Then, we use the parametric representation mentioned in the previous section. That is, 2 2 x = a sin n ,y = bcos n The variable and the domain are changed as follows. x: 0a : 0/2 And dx is 2 2a -1 dx = sin n cos d n Then, 2 2 S a 2 2a -1 = ydx = bcos n ( sin n cos)d 4 0 0 n 2 2 2 2 n -1 n +1 = ab sin cos d n 0 According to "岩波 数学公式Ⅰ" p 243 , 2 1 +1 +1 sin cos d = , , > -1 0 2 2 2 Using this, 2 2 2 2 -1+1 +1+1 2 -1 +1 1 n n 1 1 1 sin n cos n d = , = , +1 0 2 2 2 2 n n - 3 - 1 1 1 2 +1 1+ n n n n = = 1 1 2 2 2 + +1 1+ n n n Substituting this for the avove, 1 2 1 2 1+ 1+ S 2 n n n = ab = ab 4 n 2 2 2 1+ 1+ n n Quadrupling both sides, we obtain the desired expression. Note Especially when a =b =1 , 1 2 2 = 4 1+ / 1+ n n n n n n is the area of the unit circle ||x +||y =1 n > 0 , and should be called constant of the circle of degree n . That is, the area of an ellipse of degree n is given by abn , where a is the major axis and b is the minor axis. In fact, when n =2 , 1 2 1 1 2 1+ 2 2 2 1 2 = 4 = 4 = 4 = 2 2 1! 2 1+ 2 2 2 This is the area of the circle x + y = 1 , and the area of an ellipse of degree2 is given by ab . Example: Area of an ellipse of degree 2.3, a =3, b =2 As the result of calculating this area in numerical integral and the formula, both were completely consistent. - 4 - 6.3 A part of area of a superellipse In this section, we calculate the area of the light-blue portion of the following figure. Formula 6.3.1 When n,a,b ()b a are positive numbers respectively, the area s()x from 0 to x of the ellipse of degree n in the 1st quadrant is given by the following expression. r nr+1 1/n ()-1 x (1.1) s()x = bΣ nr r=0 r a nr+1 Proof The area s()x from 0 to x of the superellipse in the 1st quadrant is given by the following integral. 1 n x xn s()x = b 1- n dx 0< x a 0 a This integral is not expressed with elementary functions except for a special case. So, we try the termwise integration. Applying generalized binomial theorem to this integrand, 1 n nr n x r 1/n x 1- n = Σ()-1 a r=0 r a Integrating the right side with respect to x from 0 to x term by term, 1 n nr x n x x r 1/n x 1- n dx = Σ()-1 dx 0 a r=0 r 0 a r nr+1 1/n ()-1 x = Σ nr r=0 r a nr+1 Multiplying this by b , we obtain (1.1) . Example: Area of the ellipse of degree 2.3, a =3, b =2 on x=02 - 5 - We calculate the area of the light-blue portion of the above figure. As the result of calculating by numerical integral and (1.1) , both were completely consistent. By-products Especially integrating this to a , 1 n n r nr+1 r a x 1/n ()-1 a 1/n ()-1 b 1- n dx = bΣ nr = abΣ 0 a r=0 r a nr+1 r=0 r nr+1 Then, the area S of the superellipse is r 1/n ()-1 S = 4abΣ (1.2) r=0 r nr+1 Therefore also, 2 r 1 2 1/n ()-1 1+ / 1+ = Σ (1.3) n n r=0 r nr+1 - 6 - 6.4 Arc length of a superellipse 6.4.1 Arc length of an oblong superellipse In this sub-section, we calculate the length bP of the following oblong superellipse. 1 n xn y = b 1- 0<ba , n > 0 (1.0) a n Formula 6.4.1 n 1 - b n-1 n Let n,a,b ()b a are positive numbers, x be a number s.t. 0< x a 1 + a and ()z s be Pochhammer's symbol. Then, the arc length l()x from 0 to x of (1.0) is given by the followings. 2()n-1 r -1+s r n r ns 1/2 n b 2 x2(-1 ) + +1 (1.1) l()x = ΣΣ nr ns r=0 s=0 r 2()n -1 r a 2 + 2()n -1 r+ns+1 -1 n 2r n r ns 1/2 2()n -1 r 1 b x2(-1 ) + +1 (1.1') = ΣΣ 2nr+ns r=0 s=0 r n s s! a 2()n -1 r+ns+1 2()r-s n r s ns r 1/2 2()n -1 ()r-s 1 b x2()-1 ()- + +1 (1.1") = ΣΣ 2n()r-s +ns r=0 s=0 r-s n s s! a 2()n -1 ()r-s +ns+1 Proof Length l()x of the curve y = f()x x on a plane is given by the following formula. x dy 2 l()x = 1+ dx x dx Differentiating both sides of (1.0) with respect to x, 1 n-1 n -1 - n dy b xn nxn-1 b xn-1 xn = 1- - = - 1- dx n a n a n a a n-1a n n-1 n b xn/a n = - a 1-xn/a n - 7 - Then, 2()n-1 1 n 2 x b 2 xn/a n (1) l()x = 1+ 2 n n dx 0 a 1-x /a This integral is not expressed with elementary functions except for a special case. So, we expand this integrand to series. n 1 2(n-1 ) b - b 2 x n/a n 0< x a 1 + n-1 n 0 < n 1 a a 2 1-x n/a n According to genralized binomial thorem, 2()n-1 1 2()n-1 r 2 2 n n n 2r n n n b x /a 1/2 b x /a 1+ 2 n n = Σ n n a 1-x /a r=0 r a 1-x /a Then, 2()n-1 r 2r n n n 1/2 b x x /a (2) l()x = Σ n n dx r=0 r a 0 1-x /a Since this integral is also a non-elementary function, we try the termwise integration. 1 s+0 = 1 + x1 + x2 + = Σ xs 1-x s=0 0 Differentiating both sides of this with respect to x and dividing it by factorial one by one, 1 1 s+1 0 1 2 3 s 2 = 1x + 2x + 3x + 4x + = Σ x ()1-x 1! s=0 1 1 1 s+2 0 1 2 3 s 3 = 1x + 3x + 6x + 10x + = Σ x ()1-x 2! s=0 2 1 s+m-1 s m = Σ x ()1-x s=0 m-1 m Multiplying both sides by x , m x s+m-1 m+s m = Σ x ()1-x s=0 m-1 n n Replacing x with x /a , n n m n()m+s x /a s+m-1 x n n = Σ 1-x /a s=0 m-1 a Integrating both sides with respect to x from 0 to x, n n m n()m+s x x /a s+m-1 x x n n dx = Σ dx 0 1-x /a s=0 m-1 0 a i.e.