Summary of the invention
The present invention aims to provide the solution that a kind of implementation is flexible, reliability is high, precision is high, compact to design, can solve the problem in the above-mentioned correlation technique.
According to an aspect of the present invention, provide a kind of multipath searching device that is used for cdma system, having comprised: the slip correlation module, it is relevant to obtain the time delay expansion to be used for the data that receive are slided; Time delay expansion interpolation and mean value calculation module are used for that interpolation is carried out in the relevant time delay expansion that obtains of sliding and handle, and expansion is asked on average to obtain time delay expansion mean value to the time delay after the interpolation; The peak value searching module is used to search for the peak value that each activates the cell time delay expansion, and sorts by the peak value size; The sequential adjusting module is used for the sequential of down physical channel and uplink physical channel is adjusted, and guarantees that largest peaks is positioned at search window central authorities; And the peak allocation module, be used to cell allocation to refer to the peak.
In above-mentioned multipath searching device, time delay expansion interpolation and mean value calculation module comprise: time delay expansion interpolation unit is used for interpolation is carried out in the relevant time delay expansion that obtains of sliding; And average calculation unit, be used for the time delay after interpolation expansion is averaged to obtain time delay expansion mean value.
In said apparatus, time delay expansion interpolation unit comprises: interpolation unit, be used between the sampling of input, inserting 0, and obtain 2 times sampling rate; And filter unit, FIR filtering is carried out in the sampling after the interpolation.
In said apparatus, interpolation unit is set a sampling and adjusted sign: be masked as sampling adjustment under 0 the situation, odd positions is deposited the position from the relevant back sampling of sliding, and the position of interpolation sampling is deposited in the even number position; And adjust under 1 the situation of being masked as in sampling, the position from the relevant back sampling of sliding is deposited in the even number position, and odd positions is deposited the position of interpolation sampling.
In said apparatus, described sequential adjusting module comprises: the down physical channel sequential is adjusted computing module, is used to calculate the exploratory sequential adjusted value of down physical channel; Adjust computing module with the uplink physical channel sequential, be used to calculate the sequential adjusted value of uplink physical channel; And the level and smooth module of down physical channel sequential adjustment, be used for according to the level and smooth down physical channel sequential of the exploratory sequential adjusted value of the down physical channel adjusted value that calculates.
In said apparatus, refer to that the peak distribution module refers to the peak distribution module for fixing peak distribution module or the self adaptation of referring to.
In said apparatus, fixing finger peak distribution module comprises: peak value meter and finger table normalization unit, peak are used for peak value meter and the index unification that refers to the peak table, and peak value are done normalization to carry out peakedness ratio; First refers to the peak protected location, is used to write down the peak-peak of current area, and checks in peak value meter and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; First refers to the peak refresh unit, is used to seek the candidate who is in synchronous regime and refers to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; The primary importance amending unit is used for referring to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And first refer to the peak map unit, is used for obtaining referring to the peak allocation result according to distributing pattern to extract the bigger finger peak of peak value from referring to the peak table.
In said apparatus, self adaptation refers to that the peak distribution module comprises: the peak value meter integral unit, be used for the peak value of all activated sub-district is sorted, and the peak value of getting the front some forms the integration table; Integrate table and peak value meter normalization unit, be used for index unification, and peak value is carried out normalization to carry out peakedness ratio integration table and peak value meter; Second refers to the peak protected location, is used for the peak-peak of recording integrating table, and checks in described integration table and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Second refers to the peak refresh unit, is used to seek the candidate who is in synchronous regime and refers to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Second place amending unit is used for referring to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And second refer to the peak map unit, is used for from referring to that the peak table extracts the bigger finger peak of peak value, obtains referring to the peak allocation result.
According to another aspect of the present invention, provide a kind of multipath searching method that is used for cdma system, this method may further comprise the steps: step S202, and it is relevant to obtain the time delay expansion that the data that receive are slided; Step S204 carries out interpolation to the relevant time delay expansion that obtains of sliding and handles, and expansion is asked on average to obtain time delay expansion mean value to the time delay after the interpolation; Step S206 searches for the peak value that each activates the cell time delay expansion, and sorts by the peak value size; Step S208 adjusts the sequential of down physical channel and uplink physical channel, guarantees that largest peaks is positioned at search window central authorities; And step S210, refer to that with fixing peak apportion design or the self adaptation of referring to the peak apportion design is that cell allocation refers to the peak.
In above-mentioned multipath searching method, step S204 comprises: step a, and interpolation is carried out in the relevant time delay expansion that obtains to sliding; And step b, expansion averages to obtain time delay expansion mean value to the time delay after the interpolation.
In above-mentioned multipath searching method, step a comprises: insert 0 between the sampling of input, obtain 2 times sampling rate; And the sampling after the interpolation carried out FIR filtering.
In above-mentioned multipath searching method, in the step of carrying out interpolation, set a sampling and adjust sign:
Be masked as in described sampling adjustment under 0 the situation, odd positions is deposited the position from the relevant back sampling of sliding, and the position of interpolation sampling is deposited in the even number position; And
Adjust under 1 the situation of being masked as in sampling, the position from the relevant back sampling of sliding is deposited in the even number position, and odd positions is deposited the position of interpolation sampling.
In above-mentioned multipath searching method, step S208 comprises: the sequential adjusted value that calculates exploratory sequential adjusted value of down physical channel and uplink physical channel; And according to the level and smooth down physical channel sequential of the exploratory sequential adjusted value of the down physical channel that calculates adjusted value.
In above-mentioned multipath searching method, refer to the peak by the following steps fixed allocation: step a, unify the index of peak value meter and finger peak table, and peak value done normalization to carry out peakedness ratio; Step b, the peak-peak of record current area, and in described peak value meter, check and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Step c seeks the candidate who is in synchronous regime and refers to the peak, refers to described candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Steps d refers to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And step e, according to distributing pattern from refer to the peak table, to extract the bigger finger peak of peak value, obtain referring to the peak allocation result.
In above-mentioned multipath searching method, distribute by the following steps self adaptation to refer to peak: step a, the peak value of all activated sub-district is sorted, the peak value of getting the front some forms the integration table; Step b will integrate the index unification of table and peak value meter, and peak value is carried out normalization to carry out peakedness ratio; Step c, the peak-peak of recording integrating table, and in the integration table, check and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Steps d is sought the candidate who is in synchronous regime and is referred to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Step e refers to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And step f, from refer to the peak table, extract the bigger finger peak of peak value, obtain referring to the peak allocation result.
By technique scheme, the present invention has realized following technique effect:
1) improved system accuracy by increasing time delay expansion interpolation step, in this step, realized the function of time delay expansion mean value calculation in the conventional method simultaneously, made program structure compact more, reasonable.
2) time delay expansion interpolation and peak value searching can be selected implementation by the order of host computer, and like this, not only implementation is flexible, has increased reliability, and implementation can backup each other, and complements each other, and selects according to actual conditions.
3) support two kinds to refer to the peak distribution method: fixing peak apportion design and the self adaptation of referring to refers to the peak apportion design, and the method for salary distribution is flexible, and adaptability is strong.
4) utilize the result of Multipath searching, finished the sequential adjustment simultaneously.Not only make program structure compact more, reasonable, and reduced the scheduling between task, thereby saved overhead.
Other features and advantages of the present invention will be set forth in the following description, and, partly from specification, become apparent, perhaps understand by implementing the present invention.Purpose of the present invention and other advantages can realize and obtain by specifically noted structure in the specification of being write, claims and accompanying drawing.
Embodiment
Below with reference to the accompanying drawings and in conjunction with the embodiments, describe the present invention in detail.
Figure 1A is the block diagram according to multipath searching device of the present invention.With reference to Figure 1A, multipath searching device according to the present invention comprises: slip correlation module 102, and it is relevant to obtain the time delay expansion to be used for the data that receive are slided; Time delay expansion interpolation and mean value calculation module 104 are used for that interpolation is carried out in the relevant time delay expansion that obtains of sliding and handle, and expansion is asked on average to obtain time delay expansion mean value to the time delay after the interpolation; Peak value searching module 106 is used to search for the peak value that each activates the cell time delay expansion, and sorts by the peak value size; Sequential adjusting module 108 is used for the sequential of down physical channel and uplink physical channel is adjusted, and guarantees that largest peaks is positioned at search window central authorities; And peak allocation module 110, be used to cell allocation to refer to the peak.
In above-mentioned multipath searching device, time delay expansion interpolation and mean value calculation module 104 comprise: time delay expansion interpolation unit 104-2 is used for interpolation is carried out in the relevant time delay expansion that obtains of sliding; And average calculation unit 104-4, be used for the time delay after interpolation expansion is averaged to obtain time delay expansion mean value.
In said apparatus, time delay expansion interpolation unit 104-2 comprises: interpolation unit 104-2-2, be used between the sampling of input, inserting 0, and obtain 2 times sampling rate; And filter unit 104-2-4, FIR filtering is carried out in the sampling after the interpolation.
In said apparatus, interpolation unit is set a sampling and adjusted sign: be masked as sampling adjustment under 0 the situation, odd positions is deposited the position from the relevant back sampling of sliding, and the position of interpolation sampling is deposited in the even number position; And adjust under 1 the situation of being masked as in sampling, the position from the relevant back sampling of sliding is deposited in the even number position, and odd positions is deposited the position of interpolation sampling.
In said apparatus, described sequential adjusting module 108 comprises: the down physical channel sequential is adjusted computing module 108-2, is used to calculate the exploratory sequential adjusted value of down physical channel; Adjust computing module 108-3 with the uplink physical channel sequential, be used to calculate the sequential adjusted value of uplink physical channel; And the level and smooth module 108-4 of down physical channel sequential adjustment, be used for sequential adjusted value according to the level and smooth down physical channel of the exploratory sequential adjusted value of down physical channel that calculates.
In said apparatus, refer to that peak distribution module 110 refers to that for fixing peak distribution module 110B or self adaptation refer to peak distribution module 110C.
Figure 1B shows the fixing block diagram that refers to the peak distribution module.Shown in Figure 1B, fixing finger peak distribution module 110B comprises: peak value meter and finger peak table normalization unit 110B-2 are used for the index unification with peak value meter, and peak value are done normalization to carry out peakedness ratio; First refers to peak protected location 110B-4, is used to write down the peak-peak of current area, and checks in peak value meter and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; First refers to peak refresh unit 110B-6, is used to seek the candidate who is in synchronous regime and refers to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Primary importance amending unit 110B-8 is used for referring to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And first refer to peak map unit 110B-10, is used for obtaining referring to the peak allocation result according to distributing pattern to extract the bigger finger peak of peak value from referring to the peak table.
Fig. 1 C shows the block diagram that self adaptation refers to the peak distribution module.Shown in Fig. 1 C, in said apparatus, self adaptation refers to that peak distribution module 110C comprises: peak value meter integral unit 110C-2, be used for the peak value of all activated sub-district is sorted, and the peak value of getting the front some forms the integration table; Integrate table and peak value meter normalization unit 110C-4, be used for index unification, and peak value is carried out normalization to carry out peakedness ratio integration table and peak value meter; Second refers to peak protected location 110C-6, is used for the peak-peak of recording integrating table, and checks in described integration table and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Second refers to peak refresh unit 110C-8, is used to seek the candidate who is in synchronous regime and refers to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Second place amending unit 110C-10 is used for referring to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And second refer to peak map unit 110C-12, is used for from referring to that the peak table extracts the bigger finger peak of peak value, obtains referring to the peak allocation result.
Fig. 2 is the flow chart according to multipath searching method of the present invention.With reference to Fig. 2, multipath searching method according to the present invention may further comprise the steps: step S202, and it is relevant to obtain the time delay expansion that the data that receive are slided; Step S204 carries out interpolation to the relevant time delay expansion that obtains of sliding and handles, and expansion is asked on average to obtain time delay expansion mean value to the time delay after the interpolation; Step S206 searches for the peak value that each activates the cell time delay expansion, and sorts by the peak value size; Step S208 adjusts the sequential of down physical channel and uplink physical channel, guarantees that largest peaks is positioned at search window central authorities; And step S210, for cell allocation refers to the peak.
In above-mentioned multipath searching method, step S204 comprises: step a, and interpolation is carried out in the relevant time delay expansion that obtains to sliding; And step b, expansion averages to obtain time delay expansion mean value to the time delay after the interpolation.
In above-mentioned multipath searching method, step a comprises: insert 0 between the sampling of input, obtain 2 times sampling rate; And the sampling after the interpolation carried out FIR filtering.
In above-mentioned multipath searching method, set a sampling and adjust sign in the step of carrying out interpolation: be masked as in described sampling adjustment under 0 the situation, odd positions is deposited the position from the relevant back sampling of sliding, and the position of interpolation sampling is deposited in the even number position; And
Adjust under 1 the situation of being masked as in sampling, the position from the relevant back sampling of sliding is deposited in the even number position, and odd positions is deposited the position of interpolation sampling.
In above-mentioned multipath searching method, step S208 comprises: the sequential adjusted value that calculates exploratory sequential adjusted value of down physical channel and uplink physical channel; And according to the sequential adjusted value of the level and smooth down physical channel of the exploratory sequential adjusted value of the down physical channel that calculates.
In above-mentioned multipath searching method, refer to the peak by the following steps fixed allocation: step a, unify the index of peak value meter and finger peak table, and peak value done normalization to carry out peakedness ratio; Step b, the peak-peak of record current area, and in described peak value meter, check and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Step c seeks the candidate who is in synchronous regime and refers to the peak, refers to described candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Steps d refers to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And step e, according to distributing pattern from refer to the peak table, to extract the bigger finger peak of peak value, obtain referring to the peak allocation result.
In above-mentioned multipath searching method, distribute by the following steps self adaptation to refer to peak: step a, the peak value of all activated sub-district is sorted, the peak value of getting the front some forms the integration table; Step b will integrate the index unification of table and peak value meter, and peak value is carried out normalization to carry out peakedness ratio; Step c, the peak-peak of recording integrating table, and in the integration table, check and refer to that whether the path, peak exists and checks the variation that refers to the peak position, refreshes the state at finger peak according to the path detection result; Steps d is sought the candidate who is in synchronous regime and is referred to the peak, refers to the candidate refer to that with the candidate peak replaces a certain finger peak under the situation of peak value greater than the peak value at other a certain finger peak at peak; Step e refers to that according to the sequential adjusted value correction of down physical channel peak and candidate refer to the value at peak; And step f, from refer to the peak table, extract the bigger finger peak of peak value, obtain referring to the peak allocation result.
Below in conjunction with Fig. 3 to Fig. 8 embodiments of the invention are described.Fig. 3 is time delay expansion interpolation and the mean value calculation module flow chart according to the embodiment of the invention.
The implementation of time delay expansion interpolation and mean value calculation and peak value searching can be selected, and during specific implementation, selects by the host computer configures operation parameters.
The resolution of the time delay growth data that sliding is correlated with obtains for example is every sample sampling, and this time delay expansion will reach the resolution of every 2*sample sampling after the interpolation.
Interpolation process comprises two steps: the first step, between the input sampling, insert 0, and obtain 2 times sampling rate; Second step is with the filtering of FIR interpolation filter.
Here established a sampling and adjusted sign: if sampling adjustment is masked as 0: odd positions is deposited the position from the relevant back sampling of sliding; The position of interpolation sampling is deposited in the even number position.If sampling adjustment is masked as 1: the position from the relevant back sampling of sliding is deposited in the even number position; Odd positions is deposited the position of interpolation sampling.
In the peak value searching process, the average noise of calculation delay expansion, the average of the just necessary calculation delay expansion of calculating average noise, the interpolation process of this process and front can be finished together, and realization flow is as shown in Figure 3.With reference to Fig. 3, this flow process may further comprise the steps:
Step S302 at first judges whether to be hardware interrupts, if words execution in step S304 then, if not, then return and rejudge;
Step S304 judges whether the time delay growth data is finished, if words execution in step S306 then, if not, then return and rejudge;
Step S306 obtains the time delay growth data and mark is adjusted in sampling;
Step S308 judge to adjust to mark whether being zero, if words execution in step S310 then, if not, execution in step S312 then;
Step S310, odd number count point and are directly output, and even number point is first filter output;
Step S312, even number count point and are directly output, and odd point is first filter output;
Step S314 exports average even data;
Step S316 exports average odd data; And
Step S318 averages all data from the sample survey.
Fig. 4 adjusts calculation flow chart according to the down physical channel sequential of a plurality of sub-districts of the embodiment of the invention.
The down physical channel sequential changes according to the fluctuation in path delay.The uplink physical channel sequential is corresponding with the down physical channel sequential.Phase place with the path of peak value maximum refreshes down physical channel with reference to sequential.Refresh down physical channel with reference to sequential, the uplink physical channel sequential also adjusts accordingly at every turn, yet its adjustment has certain limitation, and the scope that its sequential moves in our monitor elapsed 1 second is just no longer revised sequential if the sequential fluctuation surpasses certain value.Identical sequential restriction also is applicable to the down physical channel sequential.
This step is divided into three sub-steps and realizes: calculate the adjustment of down physical channel sequential; Calculate the adjustment of uplink physical channel sequential; Level and smooth down physical channel sequential is adjusted.
The purpose that the down physical channel sequential is adjusted is the central authorities that the path of peak value maximum placed search window.All calculate a tentative down physical channel sequential respectively for each sub-district of opening and adjust numerical value.The adjustment calculation flow chart of a plurality of sub-districts as shown in Figure 4, the specific implementation of step S406 as shown in Figure 5, this is the exploratory adjustment calculation flow process of each sub-district.With reference to Fig. 4, calculation process may further comprise the steps:
Step S402 judges whether that first sub-district opened, if, execution in step S404 then, otherwise execution in step S408;
Step S404 points to corresponding peak value meter;
Step S406 calls " down physical channel sequential adjustment exploration value calculation procedure ";
Step S408 switches to another sub-district;
Step S410 judges whether that all sub-districts all calculate, if then return, otherwise returns execution in step S402.
Fig. 5 is the down physical channel sequential adjustment exploration value calculation flow chart according to a sub-district of the embodiment of the invention.With reference to Fig. 5, calculate according to the down physical channel sequential adjustment exploration value of a sub-district of the embodiment of the invention and may further comprise the steps:
Step S502 is provided with the initial value of exploratory adjusted value;
Step S504, peak and symbol middle position are relatively; And
Step S506 sets exploratory adjusted value according to comparative result.Because the peak value in the peak value meter is arranged from big to small by peak value, so the calculating that can adjust with the position of first peak value in the peak value meter.According to highest peak value position and symbol middle position result relatively, determine exploratory adjusted value by certain rule.
Adjust for the uplink physical channel sequential, scope with uplink physical channel sequential adjustment in the monitor elapsed certain hour, the interior all adjusted values of the certain hour of promptly adding up over, this accumulated value and adjustment limited field are compared, according to result relatively, determine the sequential adjusted value by certain rule.Flow chart as shown in Figure 6.
Fig. 6 adjusts calculation flow chart according to the uplink physical channel sequential of the embodiment of the invention.With reference to Fig. 6, calculate according to the uplink physical channel sequential adjustment of the embodiment of the invention and may further comprise the steps:
Step S602 sets the exploratory adjusted value of uplink physical channel according to the exploratory adjusted value of the down physical channel of the strongest sub-district of peak value;
Step S604, the adjusted value in the certain hour that adds up;
Step S606, accumulated value compares with the adjustment limited field; And
Step S608 is according to comparative result decision adjusted value.
For the variation of level and smooth down physical channel sequential, adjusted value is restricted.The method of down physical channel sequential restriction and the method for above-mentioned uplink physical channel are similar, and only the restriction of down physical channel sequential relates to a plurality of sub-districts.Here repeat no more.
Fig. 7 is the fixedly finger peak apportion design flow chart according to the embodiment of the invention.The peak distribution of fixing finger is meant that for the sub-district of opening some, each cell allocation is distributed to the finger peak numbering of each sub-district and also fixed to the finger peak of fixed number.To each sub-district of opening, the fixing peak allocation order that refers to is called following steps, finishes to refer to refreshing and safeguarding of peak information.With reference to Fig. 7, the fixing peak apportion design that refers to may further comprise the steps:
Step S702 judges whether current area is opened, if execution in step S704 then, otherwise execution in step S714;
Step S704, peak value meter and finger peak table normalization are about to peak value meter and refer to that the index of peak table is unified, and peak value is done normalization, are convenient to peakedness ratio, and write down the peak-peak of current area;
Step S706 refers to the peak protection, promptly in current peak value meter, checks to refer to whether the path, peak also exists, and the variation that refers to the peak position.Result and finger peak counter (record refers to the number of times or the undiscovered number of times of peak path discovery) according to path detection refresh the state that refers to the peak;
Step S708 refers to that the peak refreshes, and promptly seeks the candidate who is in synchronous regime and refers to the peak, if its peak value refers to the peak value at peak greater than certain, refers to that then the peak replacement refreshes;
Step S710, position correction is promptly revised according to the sequential adjustment of down physical channel and is referred to that peak and candidate refer to the position at peak;
Step S712 refers to the peak mapping, promptly according to distributing pattern, extracts the bigger finger peak of peak value from refer to the peak table, obtains referring to the peak allocation result; And
Step S714 has judged whether last sub-district, if then return, otherwise execution in step S716; And
Step S716 carries out the flow process to next sub-district.
Fig. 8 is that the self adaptation according to the embodiment of the invention refers to peak allocation flow figure.With reference to Fig. 8, self adaptation refers to that the peak distribution method may further comprise the steps:
Step S802, peak value meter integrate, and promptly to integrate be that peak value with all activated sub-district sorts by certain rule to peak value meter, gets the peak value of the some of its front, forms and integrate table;
Step S804, integration table and refer to peak table normalization is about to the integration table and refers to that the index of peak table is unified, and peak value is done normalization, is convenient to peakedness ratio, and the peak-peak of recording integrating table;
Step S806 refers to the peak protection, promptly in current integration table, checks to refer to whether the path, peak also exists, and the variation that refers to the peak position.Result and finger peak counter (record refers to the number of times or the undiscovered number of times of peak path discovery) according to path detection refresh the state that refers to the peak;
Step S808 refers to that the peak refreshes, and promptly seeks the candidate who is in synchronous regime and refers to the peak, if its peak value refers to the peak value at peak greater than certain, refers to that then the peak replacement refreshes;
Step S810, position correction is promptly revised according to the sequential adjustment of down physical channel and is referred to that peak and candidate refer to the position at peak; And
Step S812 refers to the peak mapping, promptly extracts the bigger finger peak of peak value from refer to the peak table, obtains referring to the peak allocation result.
The above is the preferred embodiments of the present invention only, is not limited to the present invention, and for a person skilled in the art, the present invention can have various changes and variation.Within the spirit and principles in the present invention all, any modification of being done, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.