دانشنامه ریاضی و کامپیوتر

سایت جامع در باب کتب و جزوات رشته های ریاضی و کامپیوتر با دانلود مستقیم.

ورود
عضویت




    • مطلبی یافت نشد.
    • مطلبی یافت نشد.
    • مطلبی یافت نشد.



فروشگاه سی شارپ
فروشگاه کدهای php
فروشگاه asp.net

Ellipse
EllipseBipolar Ellipse construction

An ellipse is a curve that is the locus of all points in the plane the sum of whose distances r_1 and r_2 from two fixed points F_1 and F_2 (the foci) separated by a distance of 2c is a given positive constant 2a (Hilbert and Cohn-Vossen 1999, p. 2). This results in the two-center bipolar coordinate equation

 r_1+r_2=2a,
(1)

where a is the semimajor axis and the origin of the coordinate system is at one of the foci. The corresponding parameter b is known as the semiminor axis.

The ellipse is a conic section and a Lissajous curve.

An ellipse can be specified in Mathematica using Circle[{x, y}, {a, b}].

If the endpoints of a segment are moved along two intersecting lines, a fixed point on the segment (or on the line that prolongs it) describes an arc of an ellipse. This is known as the trammel construction of an ellipse (Eves 1965, p. 177).


Elliptic gears

It is possible to construct elliptical gears that rotate smoothly against one another (Brown 1871, pp. 14-15; Reuleaux and Kennedy 1876, p. 70; Clark and Downward 1930; KMODDL).

 

The ellipse was first studied by Menaechmus, investigated by Euclid, and named by Apollonius. The focus and conic section directrix of an ellipse were considered by Pappus. In 1602, Kepler believed that the orbit of Mars was oval; he later discovered that it was an ellipse with the Sun at one focus. In fact, Kepler introduced the word "focus" and published his discovery in 1609. In 1705 Halley showed that the comet now named after him moved in an elliptical orbit around the Sun (MacTutor Archive). An ellipse rotated about its minor axis gives an oblate spheroid, while an ellipse rotated about its major axis gives a prolate spheroid.

 

A ray of light passing through a focus will pass through the other focus after a single bounce (Hilbert and Cohn-Vossen 1999, p. 3). Reflections not passing through a focus will be tangent to a confocal hyperbola or ellipse, depending on whether the ray passes between the foci or not.

 

Let an ellipse lie along the x-axis and find the equation of the figure (1) where F_1 and F_2 are at (-c,0) and (c,0). In Cartesian coordinates,

 

 sqrt((x+c)^2+y^2)+sqrt((x-c)^2+y^2)=2a.
(2)

 

Bring the second term to the right side and square both sides,

 

 (x+c)^2+y^2=4a^2-4asqrt((x-c)^2+y^2)+(x-c)^2+y^2.
(3)

 

Now solve for the square root term and simplify

 

sqrt((x-c)^2+y^2) = -1/(4a)(x^2+2xc+c^2+y^2-4a^2-x^2+2xc-c^2-y^2)
(4)
= -1/(4a)(4xc-4a^2)
(5)
= a-c/ax.
(6)

 

Square one final time to clear the remaining square root,

 

 x^2-2xc+c^2+y^2=a^2-2cx+(c^2)/(a^2)x^2.
(7)

 

Grouping the x terms then gives

 

 x^2(a^2-c^2)/(a^2)+y^2=a^2-c^2,
(8)

 

which can be written in the simple form

 

 (x^2)/(a^2)+(y^2)/(a^2-c^2)=1.
(9)

 

Defining a new constant

 

 b^2=a^2-c^2
(10)

 

puts the equation in the particularly simple form

 

 (x^2)/(a^2)+(y^2)/(b^2)=1.
(11)

 

The parameter b is called the semiminor axis by analogy with the parameter a, which is called the semimajor axis (assuming b<a). The fact that b as defined above is actually the semiminor axis is easily shown by letting r_1 and r_2 be equal. Then two right triangles are produced, each with hypotenuse a, base c, and height b=sqrt(a^2-c^2). Since the largest distance along the minor axis will be achieved at this point, b is indeed the semiminor axis.

 

If, instead of being centered at (0, 0), the center of the ellipse is at (x_0, y_0), equation (◇) becomes

 

 ((x-x_0)^2)/(a^2)+((y-y_0)^2)/(b^2)=1.
(12)

 

As can be seen from the Cartesian equation for the ellipse, the curve can also be given by a simple parametric form analogous to that of a circle, but with the x and y coordinates having different scalings,

 

x = acost
(13)
y = bsint.
(14)

 

The general quadratic curve

 

 ax^2+2bxy+cy^2+2dx+2fy+g=0
(15)

 

is an ellipse when, after defining

 

Delta = |a b d; b c f; d f g|
(16)
J = |a b; b c|
(17)
I = a+c,
(18)

 

Delta!=0, J>0, and Delta/I<0. Also assume the ellipse is nondegenerate (i.e., it is not a circle, so a!=c, and we have already established is not a point, since J=ac-b^2!=0). In that case, the center of the ellipse (x_0,y_0) is given by

 

x_0 = (cd-bf)/(b^2-ac)
(19)
y_0 = (af-bd)/(b^2-ac),
(20)

 

the semi-axes lengths are

 

a^' = sqrt((2(af^2+cd^2+gb^2-2bdf-acg))/((b^2-ac)[sqrt((a-c)^2+4b^2)-(a+c)]))
(21)
b^' = sqrt((2(af^2+cd^2+gb^2-2bdf-acg))/((b^2-ac)[-sqrt((a-c)^2+4b^2)-(a+c)])).
(22)

 

and the counterclockwise angle of rotation from the x-axis to the major axis of the ellipse is

 

 phi={0 for b=0 and a<c; 1/2pi for b=0 and a>c; 1/2cot^(-1)((a-c)/(2b)) for b!=0 and a<c; pi/2+1/2cot^(-1)((a-c)/(2b)) for b!=0 and a>c.
(23)

EllipseDirectrix

The ellipse can also be defined as the locus of points whose distance from the focus is proportional to the horizontal distance from a vertical line known as the conic section directrix, where the ratio is <1. Letting r be the ratio and d the distance from the center at which the directrix lies, then in order for this to be true, it must hold at the extremes of the major and minor axes, so

 

 r=(a-c)/(d-a)=(sqrt(b^2+c^2))/d.
(24)

 

Solving gives

 

d = (a^2)/(sqrt(a^2-b^2))=(a^2)/c
(25)
r = (sqrt(a^2-b^2))/a=c/a.
(26)

 

The focal parameter of the ellipse is

 

p = (b^2)/(sqrt(a^2-b^2))
(27)
= (a^2-c^2)/c
(28)
= (a(1-e^2))/e,
(29)

 

where e is a characteristic of the ellipse known as the eccentricity, to be defined shortly.

Ellipse4Points

An ellipse whose axes are parallel to the coordinate axes is uniquely determined by any four non-concyclic points on it, and the ellipse passing through the four points (x_1,y_1), (x_2,y_2), (x_3,y_3), and (x_4,y_4) has equation

 

 |x^2 y^2 x y 1; x_1^2 y_1^2 x_1 y_1 1; x_2^2 y_2^2 x_2 y_2 1; x_3^2 y_3^2 x_3 y_3 1; x_4^2 y_4^2 x_4 y_4 1|=0.
(30)

 

Let four points on an ellipse with axes parallel to the coordinate axes have angular coordinates t_i for i=1, 2, 3, and 4. Such points are concyclic when

 

 s_1s_2s_3+s_1s_2s_4+s_1s_3s_4+s_2s_3s_4-(s_1+s_2+s_3+s_4)=0,
(31)

 

where the intermediate variable s_i=tan(t_i/2) has been defined (Berger et al. 1984; Trott 2006, pp. 39-40). Rather surprisingly, this same relationship results after simplification of the above where s_i is now interpreted as s_i=sin(t_i). An equivalent, but more complicated, condition is given by

 

 (s_1^4+s_2^4+s_3^4+s_4^4)+(4s_2^2s_3^2s_1^2+s_2^2s_4^2s_1^2+s_3^2s_4^2s_1^2+s_2^2s_3^2s_4^2)-4(s_2s_3s_4s_1^3+s_2s_3s_4^3s_1+s_2s_3^3s_4s_1+s_2^3s_3s_4s_1)-2(s_2^2s_1^2+s_3^2s_1^2+s_4^2s_1^2+s_2^2s_3^2+s_2^2s_4^2-2s_3^2s_4^2)+8s_2s_3s_4s_1. =0
(32)

 

Like hyperbolas, noncircular ellipses have two distinct foci and two associated directrices, each conic section directrix being perpendicular to the line joining the two foci (Eves 1965, p. 275).

 

Define a new constant 0<=e<1 called the eccentricity (where e=0 is the case of a circle) to replace b

 

 e=sqrt(1-(b^2)/(a^2)),
(33)

 

from which it follows that

 

b = asqrt(1-e^2)
(34)
= sqrt(a^2-c^2)
(35)
c = sqrt(a^2-b^2)
(36)
= ae
(37)
e = sqrt(1-(b^2)/(a^2))
(38)
= c/a.
(39)

 

The eccentricity can therefore be interpreted as the position of the focus as a fraction of the semimajor axis.

EllipseFocus

If r and theta are measured from a focus F instead of from the center C (as they commonly are in orbital mechanics) then the equation s of the ellipse are

 

x = c+rcostheta
(40)
y = rsintheta,
(41)

 

and (◇) becomes

 

 ((c+rcostheta)^2)/(a^2)+(r^2sin^2theta)/(b^2)=1.
(42)

 

Clearing the denominators gives

 

 b^2(c^2+2crcostheta+r^2cos^2theta)+a^2r^2sin^2theta=a^2b^2.
(43)

 

Substituting in sin^2theta=1-cos^2theta gives

 

 b^2c^2+2rcb^2costheta+b^2r^2cos^2theta+a^2r^2-a^2r^2cos^2theta=a^2b^2.
(44)

 

Plugging in to re-express b and c in terms of a and e,

 

 a^2(1-e^2)a^2e^2+2aea^2(1-e^2)rcostheta+a^2(1-e^2)r^2cos^2theta +a^2r^2-a^2r^2cos^2theta=a^2[a^2(1-e^2)].
(45)

 

Dividing by -a^2 and simplifying gives

 

 -r^2+[ercostheta-a(1-e^2)]^2=0,
(46)

 

which can be solved for r to obtain

 

 r=+/-[ercostheta-a(1-e^2)].
(47)

 

The sign can be determined by requiring that r must be positive. When e=0, (47) becomes r=+/-(-a), but since a is always positive, we must take the negative sign, so (47) becomes

 

 r=a(1-e^2)-ercostheta
(48)

 

 r(1+ecostheta)=a(1-e^2)
(49)

 

 r=(a(1-e^2))/(1+ecostheta).
(50)

Ellipsex

The distance from a focus to a point with horizontal coordinate x (where the origin is taken to lie at the center of the ellipse) is found from

 

 costheta=(x-c)/r.
(51)

 

Plugging this into (50) yields

 

 r+e(x-c)=a(1-e^2)
(52)

 

 r=a(1-e^2)-e(x-c).
(53)

 

In pedal coordinates with the pedal point at the focus, the equation of the ellipse is

 

 (b^2)/(p^2)=(2a)/r-1.
(54)

 

The arc length of the ellipse is

 

s(t) = aE(t,e)
(55)
= aE(t,sqrt(1-(b^2)/(a^2)))
(56)
= bE(t,sqrt(1-(a^2)/(b^2))),
(57)

 

where E(t,e) is an incomplete elliptic integral of the second kind with elliptic modulus e (the eccentricity).

 

The relationship between the polar angle from the ellipse center theta and the parameter t follows from

 

 theta=tan^(-1)(y/x)=tan^(-1)(b/atant).
(58)

EllipseFunction

This function is illustrated above with theta shown as the solid curve and t as the dashed, with b/a=0.6. Care must be taken to make sure that the correct branch of the inverse tangent function is used. As can be seen, theta weaves back and forth around t, with crossings occurring at multiples of pi/2. The curvature and tangential angle of the ellipse are given by

 

kappa(t) = (ab)/((b^2cos^2t+a^2sin^2t)^(3/2))
(59)
phi(t) = tan^(-1)(a/btant).
(60)

 

The entire perimeter p of the ellipse is given by setting t=2pi (corresponding to theta=2pi), which is equivalent to four times the length of one of the ellipse's quadrants,

 

p = aE(2pi,e)
(61)
= 4aE(1/2pi,e)
(62)
= 4aE(e),
(63)

 

where E(e) is a complete elliptic integral of the second kind with elliptic modulus e (the eccentricity). The perimeter can be computed using the rapidly converging Gauss-Kummer series as

 

p = pi(a+b)sum_(n=0)^(infty)(1/2; n)^2h^(2n)
(64)
= pi(a+b)(1+1/4h^2+1/(64)h^4+1/(256)h^6+...)
(65)

 

(Sloane's A056981 and A056982), where (n; k) is a binomial coefficient and

 

 h=(a-b)/(a+b).
(66)

 

This can also be written analytically as

 

p = pi(a+b)_2F_1(-1/2,-1/2;1;h^2)
(67)
= 2(a+b)[2E(h)+2(h^2-1)K(h)],
(68)

 

where _2F_1(a,b;c;z) is a hypergeometric function, K(k) is a complete elliptic integral of the first kind.


ارسال نظر

سوال: کلمه"ریاضی" را در پایین وارد کنید
پررنگ کج خط دار خط دار در وسط | سمت چپ وسط سمت راست | قرار دادن شکلک قراردادن لینکقرار دادن لینک حفاظت شده انتخاب رنگ | پنهان کردن متن قراردادن نقل قول تبدیل نوشته ها به زبان روسی قراردادن Spoiler

پروژه دانلود مقاله