METHOD FOR DETERMINING TRANSITION HEIGHT ELEMENTS IN FLIGHT CLIMBING STAGE BASED ON CONSTANT VALUE SEGMENT IDENTIFICATION
20220358846 · 2022-11-10
Inventors
Cpc classification
International classification
Abstract
A method for determining transition height elements in a flight climbing stage based on constant value segment identification comprises the steps of splitting a speed component and a Mach component from a flight track, and performing linear interpolation on the two respectively; discretizing the interpolated speed component, and setting a threshold for filtering to obtain a speed discrete value set; identifying a constant-speed segment, and acquiring a maximum constant-speed value and a maximum moment of the constant-speed segment; keeping the Mach component of the track with a time no less than the constant-speed maximum moment; discretizing the kept Mach components, and filtering to obtain a Mach discrete value set; identifying a constant-Mach segment, and acquiring a constant-Mach value corresponding to a minimum moment of the constant-Mach segment; and calculating a transition height in the flight climbing stage according to the constant-speed value and the constant-Mach value obtained.
Claims
1. A method for determining transition height elements in a flight climbing stage based on constant value segment identification, comprising the following steps of: step 1: for track data TR={tp.sub.ii=1, . . . ,n} of one flight, wherein an i.sup.th track point tp.sub.i is denoted by one vector, tp.sub.i=[ts.sub.i, sp.sub.i, ma.sub.i], ts.sub.i, sp.sub.i and ma.sub.i respectively denote a time, a speed and a Mach number of the current track point, respectively extracting a speed component and a Mach component from the track TR and recording the two components as a first speed component TR.sub.s_raw and a first Mach component TR.sub.m_raw; wherein, n denotes a total number of track points in the track data TR, and n is a positive integer; step 2: expanding the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw by adopting a linear interpolation method to obtain a second speed component TR.sub.s and a second Mach component TR.sub.m; step 3: discretizing the second speed component TR.sub.s of the track to obtain a discrete speed component TR.sub.sd; step 4: filtering each discrete value in the discrete speed component TR.sub.sd according to a threshold thr, and acquiring a speed discrete value set SP; step 5: identifying a constant-speed segment of the flight according to the speed discrete value set SP, and acquiring a maximum constant-speed value sp.sub.c and a maximum moment ts.sub.cs of the constant-speed segment; Step 6: keeping the Mach component of the track with a time no less than the ts.sub.cs, in the second Mach component TR.sub.m to obtain a third Mach component TR.sub.m_cut; step 7: discretizing the third Mach component TR.sub.m_cut of the track to obtain a discrete Mach component TR.sub.md; step 8: filtering each discrete value in the discrete Mach component TR.sub.md according to the threshold thr to acquire a Mach discrete value set MA; step 9: identifying a constant-Mach segment of the flight according to the Mach discrete value set MA, and acquiring a constant-Mach value ma.sub.c corresponding to a minimum moment; step 10: calculating a transition height H.sub.trans of the flight according to the maximum constant-speed value sp.sub.c and the constant-Mach value ma.sub.c corresponding to the minimum moment; and step 11: obtaining a real situation of flight track and adjusting flying parameters of the flight.
2. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 1, wherein in the step 1, the process of extracting the speed component and the Mach component is: the first speed component is that TR.sub.s_raw={s.sub.i,i=1, . . . , n}, wherein s.sub.i=[ts.sub.i,sp.sub.i]; and the first Mach component is that TR.sub.m_raw, ={m.sub.i,i=1, . . . ,n}, wherein m.sub.i=[ts.sub.i, ma.sub.i].
3. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 2, wherein in the step 2, the linear interpolation process is: step 2.1: arranging the track points in the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw in an ascending order according to the time ts.sub.i of the track points, wherein the time ts.sub.i of the track points is in a unit of second; step 2.2: when ts.sub.i+1−ts.sub.i<1, respectively interpolating ts.sub.i+1−ts.sub.i−1 speed values and Mach values respectively, wherein the p.sup.th interpolated speed value is that s.sub.interp_p=[ts.sub.i+p,sp.sub.i+p(sp.sub.i+1−sp.sub.i)/(ts.sub.i+1−ts.sub.i)], and the p.sup.th interpolated Mach value is that m.sub.interp_p=[ts.sub.i+p,ma.sub.i+p(ma.sub.i+1−ma.sub.i)/(ts.sub.i+1−ts.sub.i)], wherein p=1,2, . . . ,ts.sub.i+1−ts.sub.i−1; and step 2.3: when ts.sub.i+1−ts.sub.i≤1, no interpolation is needed; after interpolating the track points in the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw, acquiring a second speed component TR.sub.s={s.sub.idx,idx=1, . . . ,N}, wherein s.sub.idx=[ts.sub.idx,sp.sub.idx], and a second Mach component TR.sub.m={m.sub.idx,idx=1, . . . ,N}, wherein m.sub.idx=[ts.sub.idx,ma.sub.idx], and N denotes a sum of a total number of track points and a total number of interpolation points in the track data TR.
4. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 3, wherein in the step 3, the process of discretizing the second speed component TR.sub.s of the track is: for any speed sp.sub.idx in the second speed component TR.sub.s, in a unit of knot, when satisfying that qj−0.5q sp.sub.idx<qj+0.5q, then a discrete value of the speed is that sp.sub.idx.sup.d=qj, wherein q is a speed discrete precision, and q belongs to R.sup.+, j is an index variable, and j=0,1,2, . . . ; and a speed component discrete value is that TR.sub.sd={s.sub.idx.sup.d,idx=1, . . . ,N}, wherein s.sub.idx.sup.d=ts.sub.idx,sp.sub.idx.sup.d].
5. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 4, wherein in the step 4, the process of acquiring the speed discrete value set SP is: setting the threshold to be that thr=0.01N, then the speed discrete value set is that SP={sp.sub.idx.sup.d||TR.sub.sd(sp.sub.idx.sup.d|≥thr,idx=1,2, . . . ,N}, wherein |TR.sub.sd(sp.sub.idx.sup.d)| denotes a number of sp.sub.idx contained in the TR.sub.sd.
6. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 5, wherein in the step 5, the process of acquiring the maximum constant-speed value sp.sub.c and the maximum moment ts.sub.e, of the constant-speed segment is: Step 5.1: arranging the elements in the speed discrete value set SP in a descending order, and acquiring that SP=[sp.sub.1.sup.c,sp.sub.2.sup.c, . . . , sp.sub.k.sup.c, . . . , sp.sub.|SP|.sup.c], wherein |SP| denotes a number of elements in the speed discrete value set SP, 1≤k≤|SP|, and letting that k=1; step 5.2: acquiring a first track point set TR.sub.s.sup.k={s.sub.idx.sup.d|s.sub.idx.sup.d∈TR.sub.sd,sp.sub.idx.sup.d=sp.sub.k.sup.c,idx=1, . . . , N}; step 5.3: arranging the track points s.sub.idx.sup.d in the first track point set TR.sub.s.sup.k according to an ascending ordering of ts.sub.idx, and when a time difference of two continuous track points is less than or equal to 4 seconds, dividing the two track points into one track point set; if the time difference is greater than 4 seconds, dividing the previous track point into a current track point set and dividing the latter track point into next track point set, thus dividing the track points into g.sub.k track point sets, wherein TR.sub.s.sup.k={TR.sub.s1.sup.k,TR.sub.s2.sup.k, . . . TR.sub.sg1.sup.k}; step 5.4: detecting each track point set, and discarding a track point set if a total duration of the track point group is less than 30 seconds or a standard deviation of a speed value of the track point set is greater than 0.3q; otherwise, keeping the track point set; and step 5.5: when a number of the kept track point sets is greater than or equal to 1, then sp.sub.c=sp.sub.k.sup.c and the maximum moment in the track point set is ts.sub.cs, executing step 6; when the number of the kept track point sets is 0, and k+1|SP|, letting k.fwdarw.k+1, and skipping to step 5.2; and when k+1>|SP|, letting that sp.sub.c=−1, and ts.sub.cs=0, and then executing step 6.
7. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 6, wherein in the step 6, the process of acquiring the third Mach component TR.sub.m_cut is: recording a number of elements in the third Mach component TR.sub.m_cut as N.sub.cut wherein TR.sub.m_cut={m.sub.idx|m.sub.idx∈TR.sub.m,ts.sub.idx≥ts.sub.cs,idx=1,2, . . . ,N}; and in step 7: the process of discretizing the third Mach component TR.sub.m_cut, of the track is: for any Mach number ma.sub.index, the Mach number ma.sub.index is dimensionless; when the Mach number satisfies that uj−0.5u≤ma.sub.index<uj+0.5u, a discrete value of the Mach number is that ma.sub.index.sup.d=uj, wherein u denotes a Mach number discrete precision, u belongs to R.sup.+, j is an index variable, and j=0,1,2, . . . ; and a Mach component discrete value is that TR.sub.md={m.sub.index.sup.d,index=1, . . . , N.sub.cut}, wherein m.sub.index.sup.d=[ts.sub.index,ma.sub.index.sup.d].
8. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 7, wherein in the step 8, the process of acquiring the Mach discrete value set MA is: setting the threshold to be that thr=0.01N, then the Mach discrete value set is that MA={ma.sub.index.sup.d||TR.sub.md(ma.sub.index.sup.d)|≥thr,index=1,2, . . . N.sub.cut}, wherein |TR.sub.md(ma.sub.index.sup.d)| denotes a number of ma.sub.index.sup.d contained in the discrete Mach component TR.sub.md.
9. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 8, wherein in the step 9, the process of acquiring the constant-Mach value ma.sub.i corresponding to the minimum moment is: step 9.1: recording that MA=[ma.sub.1.sup.c,ma.sub.2.sup.c, . . . ,ma.sub.k1.sup.c, . . . ,ma|MA|c], wherein |MA| denotes a number of elements in the Mach discrete value set MA, and 1≤k1≤|MA|; letting that ts.sub.cm=+∞, and ma.sub.c=−1; and letting that k1=1; step 9.2: acquiring a second track point set TR.sub.m.sup.k1={m.sub.index.sup.d|m.sub.index.sup.d∈TR.sub.md,ma.sub.index.sup.d=ma.sub.k1.sup.c,index=1, . . . , N.sub.cut}; step 9.3: arranging the track points m.sub.index.sup.d in the second track point set TR.sub.m.sup.k1 according to an ascending ordering of ts.sub.index, and when a time difference of two continuous track points is less than or equal to 4 seconds, dividing the two track points into one track point set; if the time difference is greater than 4 seconds, dividing the previous track point into a current track point set and dividing the latter track point into next track point set, thus dividing the track points into g.sub.k1 track point sets, wherein TR.sub.m.sup.k1={TR.sub.m1.sup.k1,TR.sub.m2.sup.k1, . . . TR.sub.mgk1.sup.k1}; step 9.4: detecting each track point set, and discarding a track point set if a total duration of the track point group is less than 100 seconds or a standard deviation of a Mach value of the track point set is greater than 0.3u; otherwise, keeping the track point set; and step 9.5: when a number of the kept track point sets is greater than or equal to 1, and the minimum moment ts.sub.min in the track point set is less than ts.sub.cm, setting that ts.sub.cm=ts.sub.min, and ma.sub.c=ma.sub.k1.sup.c and then executing step 10; when the number of the kept track point sets is 0, and k1+1≤|MA|, letting k1.fwdarw.k1+1, and skipping to step 9.2; and when k1+1>|MA|, letting that ma.sub.c=−1, and ts.sub.cm=+∞, and then executing step 10.
10. The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to claim 9, wherein in the step 10, the process of calculating the transition height of the flight is: when sp.sub.c=−1 or ma.sub.c=−1, the transition height elements are not obtained, and the transition height cannot be calculated; otherwise, values of sp.sub.c and ma.sub.c are substituted according to a transition height calculation function provided by Base of Aircraft Data BADA to calculate the transition height.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0045] The advantages of the above and/or other aspects of the present invention will become more apparent by further explaining the present invention with reference to the following drawings and detailed description.
[0046]
[0047]
[0048]
DETAILED DESCRIPTION
[0049] The embodiments of the present invention will be described hereinafter with reference to the drawings.
[0050] With reference to
[0051] step 1: for track data TR={tp.sub.i i=1, . . . , n} of one flight, wherein an it track point tp.sub.i is denoted by one vector, tp.sub.i=[ts.sub.i,sp.sub.i, ma.sub.i], ts.sub.i, sp.sub.i and ma.sub.i respectively denote a time, a speed and a Mach number of the current track point, respectively extracting a speed component and a Mach component from the track TR and recording the two components as a first speed component TR.sub.s_raw and a first Mach component TR.sub.m_raw;
[0052] the process of extracting the speed component and the Mach component being: the first speed component is that TR.sub.s_raw={s.sub.i,i=1, . . . ,n}, wherein s.sub.i=[ts.sub.i, sp.sub.i]; the first Mach component is that TR.sub.m_raw={m.sub.i,i=1, . . . ,n}, wherein m.sub.i=[ts.sub.i, ma.sub.i]; and n denotes a total number of track points in the track data TR.sub.i and n is a positive integer.
[0053] Step 2: expanding the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw by adopting a linear interpolation method to obtain a second speed component TR.sub.s and a second Mach component TR.sub.m;
[0054] the linear interpolation process is:
[0055] step 2.1: arranging the track points in the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw in an ascending order according to the time ts.sub.i of the track points, wherein the time ts.sub.i of the track points is in a unit of second;
[0056] step 2.2: when ts.sub.i+1−ts.sub.i>1, respectively interpolating ts.sub.i+1−ts.sub.i−1 speed values and Mach values respectively, wherein the p.sup.th interpolated speed value is that s.sub.interp_p=[ts.sub.i+p,sp.sub.i+p(sp.sub.i+1−sp.sub.i)/(ts.sub.i+1−ts.sub.i)], and the p.sup.th interpolated Mach value is that m.sub.interp_p=[ts.sub.i+p,ma.sub.i+p(ma.sub.i+1−ma.sub.i)/(ts.sub.i+1−ts.sub.i)], wherein p=1,2, . . . , ts.sub.i+1−ts.sub.i−1; and
[0057] step 2.3: when ts.sub.i+1−ts.sub.i≤1, no interpolation is needed; after interpolating the track points in the first speed component TR.sub.s_raw and the first Mach component TR.sub.m_raw, acquiring a second speed component TR.sub.s={s.sub.idx,idx=1, . . . , N}, wherein s.sub.idx=[ts.sub.idx,sp.sub.idx], and a second Mach component TR.sub.m={m.sub.idx, idx=1, . . . , N}, wherein m.sub.idx=[ts.sub.idx, ma.sub.idx], and N denotes a sum of a total number of track points and a total number of interpolation points in the track data TR.
[0058] Step 3: discretizing the second speed component TR.sub.i of the track to obtain a discrete speed component TR.sub.sd;
[0059] for any speed sp.sub.idx in the second speed component TR.sub.s, in a unit of knot, when satisfying that qj−0.5q sp.sub.idx<qj+0.5q, then a discrete value of the speed is that sp.sub.idx.sup.d=qj, wherein q is a speed discrete precision, and q belongs to R.sup.+, j is an index variable, and j=0,1,2, . . . ; and a speed component discrete value is that TR.sub.sd={s.sub.idx.sup.d=1, . . . , N}, wherein s.sub.idx.sup.d=[ts.sub.idx.sup.d].
[0060] Step 4: filtering each discrete value in the discrete speed component TR.sub.sd according to a threshold thr, and acquiring a speed discrete value set SP;
[0061] setting the threshold to be that thr=0.01N, then the speed discrete value set is that SP={sp.sub.idx.sup.d||TR.sub.sd(sp.sub.idx.sup.d|≥thr,idx=1,2, . . . ,N}, wherein |TR.sub.sd (sp.sub.idx.sup.d)| denotes a number of sp.sub.idx.sup.d contained in the TR.sub.sd.
[0062] Step 5: identifying a constant-speed segment of the flight according to the speed discrete value set SP, and acquiring a maximum constant-speed value sp.sub.c and a maximum moment ts.sub.cs, of the constant-speed segment;
[0063] step 5.1: arranging the elements in the speed discrete value set SP in a descending order, and acquiring that SP=[sp.sub.1.sup.c,sp.sub.2.sup.c, . . . , sp.sub.k.sup.c, . . . sp.sub.|sp|.sup.c], wherein |SP| denotes a number of elements in the speed discrete value set SP, 1≤k≤|SP|, and letting that k=1;
[0064] step 5.2: acquiring a first track point set TR.sub.s.sup.k={s.sub.idx.sup.d|s.sub.idx.sup.d∈TR.sub.sd,sp.sub.idx.sup.d=sp.sub.k.sup.c,idx=1, . . . ,N};
[0065] step 5.3: arranging the track points s.sub.idx.sup.d in the first track point set TR.sub.s.sup.k according to an ascending ordering of ts.sub.idx, and when a time difference of two continuous track points is less than or equal to 4 seconds, dividing the two track points into one track point set; if the time difference x is greater than 4 seconds, dividing the previous track point into a current track point set and dividing the latter track point into next track point set, thus dividing the track points into g.sub.k track point sets, wherein TR.sub.s.sup.k={TR.sub.s1.sup.k,TR.sub.s2.sup.k, . . . TR.sub.sg1.sup.k};
[0066] step 5.4: detecting each track point set, and discarding a track point set if a total duration of the track point group is less than 30 seconds or a standard deviation of a speed value of the track point set is greater than 0.3 q, and q=6; otherwise, keeping the track point set; and
[0067] step 5.5: when a number of the kept track point sets is greater than or equal to 1, then sp.sub.c=sp.sub.k.sup.c and the maximum moment in the track point set is ts.sub.cs, executing step 6; when the number of the kept track point sets is 0, and k+1≤|SP|, letting k.fwdarw.k+1 and skipping to step 5.2; and when k+1>|SP|, letting that sp.sub.c=−1, and ts.sub.cs, =0, and then executing step 6.
[0068] Step 6: keeping the Mach component of the track with a time no less than the ts.sub.cs in the second Mach component TR.sub.m to obtain a third Mach component TR.sub.m_cut;
[0069] recording a number of elements in the third Mach component TR.sub.m_cut as N.sub.cute, wherein TR.sub.m_cut={m.sub.idx|m.sub.idx∈TR.sub.m,ts.sub.idx≥ts.sub.cs,idx=1,2, . . . ,N}.
[0070] Step 7: discretizing the third Mach component TR.sub.m_cut of the track to obtain a discrete Mach component TR.sub.md;
[0071] for any Mach number ma.sub.index, the Mach number ma.sub.index is dimensionless; when the Mach number satisfies that uj−0.5u≤ma.sub.index<uj+0.5u, a discrete value of the Mach number is that ma.sub.index.sup.d=uj, wherein u denotes a Mach number discrete precision, u belongs to R.sup.+, in this embodiment, u=0.01, j is an index variable, and j=0,1,2, . . . ; and a Mach component discrete value is that TR.sub.md={m.sub.index.sup.d,index=1, . . . , N}, wherein m.sub.index.sup.d=[ts.sub.index,ma.sub.index.sup.d].
[0072] Step 8: filtering each discrete value in the discrete Mach component TR.sub.md according to the threshold thr to acquire a Mach discrete value set MA;
[0073] setting the threshold to be that thr=0.01N, then the Mach discrete value set is that MA={ma.sub.index.sup.d||TR.sub.md(ma.sub.index.sup.d)|≥thr,index=1, 2, . . . ,N.sub.cut}, wherein TR.sub.md(ma.sub.index.sup.d)| denotes a number of ma.sub.index.sup.d contained in the discrete Mach component TR.sub.md.
[0074] Step 9: identifying a constant-Mach segment of the flight according to the Mach discrete value set MA, and acquiring a constant-Mach value ma.sub.c corresponding to a minimum moment;
[0075] step 9.1: recording that MA=[ma.sub.1.sup.c,ma.sub.2.sup.c, . . . , ma.sub.k1.sup.c, . . . ,ma.sub.|MA|.sup.c], wherein ≡MA| denotes a number of elements in the Mach discrete value set MA, and 1≤k1≤|MA|; letting that ts.sub.cm=+∞, and ma.sub.c=−1; and letting that k1=1;
[0076] step 9.2: acquiring a second track point set
TR.sub.m.sup.k1={m.sub.index.sup.d|m.sub.index.sup.d∈TR.sub.md,ma.sub.index.sup.d=ma.sub.k1.sup.c,index=1, . . . ,N.sub.cut};
[0077] step 9.3: arranging the track points m.sub.index.sup.d in the second track point set TR. according to an ascending ordering of ts.sub.index, and when a time difference of two continuous track points is less than or equal to 4 seconds, dividing the two track points into one track point set; if the time difference is greater than 4 seconds, dividing the previous track point into a current track point set and dividing the latter track point into next track point set, thus dividing the track points into g.sub.k1 track point sets, wherein TR.sub.m.sup.k1={TR.sub.m1.sup.k1,TR.sub.m2.sup.k1, . . . TR.sub.mgk1.sup.k1};
[0078] step 9.4: detecting each track point set, and discarding a track point set if a total duration of the track point group is less than 100 seconds or a standard deviation of a Mach value of the track point set is greater than 0.3 u (u=0.01); otherwise, keeping the track point set; and
[0079] step 9.5: when a number of the kept track point sets is greater than or equal to 1, and the minimum moment ts.sub.min in the track point set is less than ts.sub.cm, setting that ts.sub.cm=ts.sub.min, and ma.sub.c=ma.sub.k1.sup.c, and then executing step 10; when the number of the kept track point sets is 0, and k1+1≤|MA|, letting k1.fwdarw.k1+1, and skipping to step 9.2; and when k1+1>|MA|, letting that ma.sub.c=−1, and ts.sub.cm=+∞, and then executing step 10.
[0080] Step 10: calculating a transition height H.sub.trans of the flight according to the maximum constant-speed value sp.sub.c and the constant-Mach value ma.sub.c corresponding to the minimum moment.
[0081] When sp.sub.c=−1 or ma.sub.c=−1, the transition height elements are not obtained, and the transition height cannot be calculated; otherwise, values of sp.sub.c and ma.sub.c are substituted according to a transition height calculation function provided by Base of Aircraft Data BADA to calculate the transition height.
[0082] With reference to
[0083] In this embodiment, as shown in
[0084] According to the result of step 10, the real situation of the flight track is obtained, which is used for adjusting flying parameters of the flight.
[0085] The method for determining the transition height elements in the flight climbing stage based on constant value segment identification according to this embodiment is loaded and operated in a processing server this embodiment an air traffic flow management system (ATFM system) or a corresponding computer of an air traffic control system (ATC system).
[0086] In a specific implementation, the present application provides a computer storage medium and a corresponding data processing unit, wherein the computer storage medium is capable of storing a computer program, and the computer program, when executed by the data processing unit, can run the inventive contents of the method for determining the transition height elements in the flight climbing stage based on constant value segment identification provided by the present invention and some or all steps in various embodiments. The storage medium may be a magnetic disk, an optical disk, a Read Only Storage (ROM) or a Random Access Storage (RAM), and the like.
[0087] Those skilled in the art can clearly understand that the technical solutions in the embodiments of the present invention can be realized by means of a computer program and a corresponding general hardware platform thereof. Based on such understanding, the essence of the technical solutions in the embodiments of the present invention or the part contributing to the prior art, may be embodied in the form of a computer program, i.e., a software product. The computer program, i.e., the software product is stored in a storage medium comprising a number of instructions such that a device (which may be a personal computer, a server, a singlechip, a MUU or a network device, and the like) comprising the data processing unit executes the methods described in various embodiments or some parts of the embodiments of the present invention.
[0088] The present invention provides the method for determining the transition height elements in the flight climbing stage based on constant value segment identification. There are many methods and ways to realize the technical solutions. The above is only the specific embodiments of the present invention. It should be pointed out that those of ordinary skills in the art can make some improvements and embellishments without departing from the principle of the present invention, and these improvements and embellishments should also be regarded as falling with the scope of protection of the present invention. All the unspecified components in the embodiments can be realized by the prior art.