CN106126834B - A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement - Google Patents

A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement Download PDF

Info

Publication number
CN106126834B
CN106126834B CN201610497850.6A CN201610497850A CN106126834B CN 106126834 B CN106126834 B CN 106126834B CN 201610497850 A CN201610497850 A CN 201610497850A CN 106126834 B CN106126834 B CN 106126834B
Authority
CN
China
Prior art keywords
curvature
curve
linear
point
mileage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610497850.6A
Other languages
Chinese (zh)
Other versions
CN106126834A (en
Inventor
元向辉
李晓雯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610497850.6A priority Critical patent/CN106126834B/en
Publication of CN106126834A publication Critical patent/CN106126834A/en
Application granted granted Critical
Publication of CN106126834B publication Critical patent/CN106126834B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention discloses a kind of best linear determining methods of orbit plane based on inertial navigation angle measurement, purpose is, meet that railroad track demand of technical standard, result precision are high, treatment process calculation amount is small, data processing time is short, used technical solution are as follows: on the basis of being based on inertial navigation angle measurement, in conjunction with the best wiresizing optimization model of orbit plane, outline segmentation is carried out to track according to fixed length curvature curve least square fitting slope variation first, then carry out accurate segmentation using combined iteration algorithm and is fitted best linear parameter.

Description

A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement
Technical field
The invention belongs to the data processing fields in orbit geometry parameter detection, and in particular to one kind is based on inertial navigation angular amount The best linear determining method of the orbit plane of survey.
Background technique
Rail track is to geometric parameter, such as track is linear and track alignment irregularity etc., with climatic environment and fortune The change of battalion's load is constantly changing, and has an impact to train operating safety.The track rail difference linear and linear ideal to reality As track alignment irregularity, track alignment irregularity will cause vehicle pitching, rolling stock wheel increases the horizontal force of rail, seriously When may cause that widening of gauge, the section of track be traversing or even train derailing.Therefore, reinforce track rail to detect to geometric parameter to guarantee Traffic safety is one of the emphasis of railway maintenance work.
Orbit parameter detection method widely used at present has inertial reference method, chord measurement and method of absolute coordinate system.Inertia base Quasi- method be according in car body mass block displacement and its relative displacement calculates track irregularity between axle box method, be chiefly used in big Type net rail detects vehicle, and the examined speed of precision is affected, and is not suitable for low speed line detection;Chord measurement is chiefly used in rail inspection instrument, Have the advantages that the not examined speed of measured value influences, but this method " pushing away greatly with small ", all can not really reflect track State, therefore error is larger;Method of absolute coordinate system cooperates the absolute position of observation point on CPIII control point survey track using total station Coordinate is set, is then fitted orbit geometry parameter using least-squares algorithm, this method detection speed is slow, needs to establish CPIII control Making point and carrying out long term maintenance, Data processing also needs Reference Design data, adds somewhat to railway maintenance work Cost and complexity.
Different from the above traditional detection method using position measurement and position versus as means, researcher proposes a kind of base In the track detection method of inertia measurement, odometer is cooperated to detect track rail to geometric parameter by measurement direction angle, This method makes full use of the high precision of gyroscope measurement angle, has high detection efficiency, precision height, testing result not examined The advantages that rate limitation, is with a wide range of applications.Currently, the track irregularity parameter based on the detection method calculates mostly Rely on original design data, and for some design datas missing or track is linear need re-optimization whole positive in the case of, then require The best linear parameter of orbit plane is determined from the orbital direction angle comprising track irregularity and measurement noise, mileage.
Summary of the invention
In order to solve the problems in the prior art, the present invention propose it is a kind of have meet railroad track demand of technical standard, As a result the orbit plane based on inertial navigation angle measurement for the advantage that precision is high, treatment process calculation amount is small, data processing time is short Best linear determining method.
In order to achieve the goal above, the technical scheme adopted by the invention is as follows: the following steps are included:
1) inertial navigation angle measurement is carried out to orbit plane and obtains inertial navigation deflection data sequence, by inertial navigation direction angular data sequence Column obtain the curvature Variation at each test point along route mileage difference, and form orderly song with the mileage at test point Rate-mileage sequence, to obtain curvature-mileage curve;
2) to two segment data of front and back of the test point i in the curvature of step 1)-mileage curveWithIt is quasi- using least-squares line respectively Hop algorithm is fitted, miIndicate the mileage at measuring point i, kiIt indicates the curvature at measuring point i, calculates the last period matched curve with after The slope differences of one section of matched curve are based on least square fitting slope variation according to the linear curvature feature of track and curvature curve Rule is segmented to obtain multiple groups waypoint coordinate, so that orbit plane is divided into multistage fundamental line tuple to linear;
3) linear according to track linear curvature feature fitting first to the waypoint of each section of fundamental line tuple, then group It is optimal to the fundamental line tuple to close iteration, next fundamental line tuple is reprocessed, after finally integrating to multistage fundamental line tuple The best linear parameter of whole track is obtained, so that it is determined that best linear.
Step 1) the mean curvature variation characteristic should be similar to the linear curvature feature of ideal track, and wherein ideal track is linear Curvature feature is linear for the symmetrical basic model track as made of straight line-easement curve-circular curve-easement curve-Straight Combination, The inertial navigation deflection of straightway remains unchanged, and straightway curvature is 0;Circular curve segment inertial navigation deflection changes linearly, circular curve It is circular curve radius that section curvature, which is steady state value K ∝ a 1/R, R,;Easement curve section inertial navigation deflection changes in second-degree parabola, bent Rate linear change between 0~K with mileage.
Curvature obtained in the step 1)-mileage curve need to carry out abandoning rough error and median filter process.
The slope differences of the case where being positive in the step 2) for circular curve curvature, two sections of matched curves are gentle in point of tangent to spiral There is minimum at straight point, has maximum at point of spiral to curve and point of curve to spiral;The case where being negative for circular curve curvature, two sections of fittings are bent The slope differences of line have minimum at point of spiral to curve and point of curve to spiral, have maximum at point of tangent to spiral and point of spiral to tangent.
Being fitted linear step according to the linear curvature feature of track in the step 3) includes: starting point determining first, one group Waypoint and terminal abscissa (M0,M1,M2,M3,M4,Me), then in conjunction with the feature of the linear curvature of orbit plane, according to following Principle determines its curvature value:Wherein, straightway curvature is 0, circle Curved section curvature is proportional to 1/R, MjThe abscissa of measuring point j is represented, the value of ω is that abscissa is M2To M3Measuring point curvature value it is flat Mean, E (k (M2),k(M2+1),…,k(M3)) represent to k (M2),k(M2+1),…,k(M3) it is a series of number be averaging, finally will Curvature-mileage curve that each point that coordinate determines is fitted after being gradually connected, the deflection-being fitted after mileage integral Mileage curve, and calculate the target function value under each waypoint combination.
The specific steps of combined iteration in the step 3) are as follows: determine the region of search range that iteration carries out first, substitute into general It is initial target solution (dS, X) that slightly waypoint track is linear, calculates initial target functional value f (dS, X) using objective function, described Objective function are as follows: min (f (dS, X))=∑ Δi 2, dS is independent variable, and S can be obtained by integrating to dS along mileage, and wherein S is for rail The ordered set of deflection is fitted at each measuring point mileage in road and the mileage, X is the ordered set of metric data mileage, deflection It closes, Δiii' indicate measuring point direction angle alphaiWith fitting direction angle alphai' deviation;
Then change waypoint combination in region of search, redefine the linear (dS of trackn, X), calculating target function value f (dSn, X), if f (dSn, X) and < f (dS, X), then receive new explanation, otherwise receives former solution;
Finally judge whether all waypoints combination in traversal search domain, be then algorithm termination, otherwise repeatedly step 2).
Constraint condition is introduced in the objective function constitutes the best wiresizing optimization model of orbit plane, constraint condition are as follows:
In formula: kZ=0 indicates that the curvature of straightway is 0, kciiThe curvature for indicating i-th section of circular curve is certain value εi, kh(i-1)=a* εi+b(i-1)Indicate that each measuring point curvature is 0 to ε on the preceding easement curve of i-th section of circular curveiBetween linear change, kh(i-1)=-a* εi+b(i+1)Indicate that each measuring point curvature is in ε on the rear easement curve of i-th section of circular curveiTo 0 linear change;Respectively smallest curve radius as defined in design specification, smallest circle length of curve and minimum length of transition curve; lformer,llatterRespectively indicate the length of the two sections of easement curves in same circular curve front and back.
Compared with prior art, track curvature curve of the invention is based on least square fitting slope variation rule, in conjunction with Symmetrical basic model track curvature changing rule can carry out outline segmentation to track, use group according to previous step outline segmentation result Close the best linear parameter that iterative algorithm gradually handles multistage fundamental line tuple available whole track after integration.This hair It is bright that track curvature changing rule is based in linear fitting, therefore meet Railway Design code requirement.The present invention is changed by combination Optimal trajectory parameter is determined for algorithm, and the region of search traversed is limited, therefore result precision is high, treatment process calculation amount is small, number According to short processing time.
It further, include measurement error and track irregularity in curvature-mileage curve that metric data obtains, by abandoning Rough error and median filtering can obtain and the more similar curvature-mileage curve of the linear curvature feature of ideal track.
Detailed description of the invention
Fig. 1 is that the ideal orientation angle that circular curve curvature is the symmetrical basic model railroad track horizontal alignment of timing changes with mileage Curve graph;
Fig. 2 is that the desired curvature that circular curve curvature is the symmetrical basic model railroad track horizontal alignment of timing changes song with mileage Line chart;
Fig. 3 is that circular curve curvature is timing track curvature curve least square fitting slope variation rule figure;
Fig. 4 is that the ideal orientation angle of symmetrical basic model railroad track horizontal alignment when circular curve curvature is negative changes song with mileage Line chart;
Fig. 5 is the desired curvature of symmetrical basic model railroad track horizontal alignment when circular curve curvature is negative with mileage change curve Figure;
Fig. 6 is track curvature curve least square fitting slope variation rule figure when circular curve curvature is negative;
Fig. 7 is the comparison of the linear and best linear definitive result of the linear measurement of orbit plane, to measure curvature, filtering curvature And optimization curvature indicates;
Fig. 8 is the comparison of the linear and best linear definitive result of the linear measurement of orbit plane, with measurement direction angle and optimization Deflection indicates;
Fig. 9 is flow chart of the method for the present invention.
Specific embodiment
Below with reference to specific embodiment and Figure of description the present invention will be further explained explanation.
The present invention specifically includes the following steps:
The first step, data prediction: carrying out inertial navigation angle measurement to orbit plane and obtain inertial navigation deflection data sequence, will Inertial navigation deflection data sequence obtained in detection process obtains the Curvature varying at each test point: K=d along route mileage difference α/dm, positive and negative expression circular curve bending direction form orderly curvature-mileage sequence with mileage at the point, it is inner to obtain curvature- Journey curve.It include measurement error and track irregularity in the curvature that metric data obtains-mileage curve, by abandoning rough error in Value filtering can obtain curvature-mileage curve more similar with the linear curvature feature of ideal track.
Wherein, the linear curvature feature of the ideal track is directed to straight by straight line-easement curve-circular curve-easement curve- The symmetrical basic model track that line is composed is linear, and the analysis linear direction corner characteristics along mileage of track have: straightway deflection It remains unchanged, circular curve segment deflection linear change, easement curve section deflection changes in second-degree parabola, as shown in Figure 1.Its It is 0 that curvature, which changes with mileage in straightway curvature, and circular curve segment curvature is a steady state value K ∝ 1/R (R is circular curve radius), is delayed And curved section curvature with mileage the linear change between 0 and K, as shown in Figure 2.
Second step, the outline segmentation based on track curvature curve least square fitting slope variation track: one section of definition is complete Whole straight line-easement curve-circular curve-easement curve-straight line linear combination is a fundamental line tuple.The song that the first step is obtained Two segment datas before and after measuring point i in rate-mileage curveWithIt is fitted respectively using least squares line fitting algorithm, miIt indicates at measuring point i Mileage, kiIt indicates the curvature at measuring point i, the slope differences of the last period matched curve and latter section of matched curve is calculated, according to rail The linear curvature feature in road and curvature curve are based on least square fitting slope variation rule can be to linear carry out outline segmentation Multiple groups waypoint coordinate is obtained, track is divided into multistage fundamental line tuple.
Wherein, it is to curvature-mileage curve that the track curvature curve, which is based on least square fitting slope variation rule, In measuring point i before and after two segment datasWith It is fitted respectively using least squares line fitting algorithm, calculates the slope of the last period matched curve and latter section of matched curve Difference, slope differences have extreme value at each waypoint.The case where being positive for circular curve curvature, it is gentle in point of tangent to spiral by taking Fig. 1 as an example There is minimum at straight point, has maximum at point of spiral to curve and point of curve to spiral, as shown in Figure 3;The case where being negative for circular curve curvature, There is minimum at point of spiral to curve and point of curve to spiral, has maximum at point of tangent to spiral and point of spiral to tangent.Track turn to it is opposite with Fig. 1 when deflection with Mileage change curve as shown in figure 4, curvature with mileage change curve as shown in figure 5, track curvature curve at this time is based on minimum Two to multiply fit slope changing rule as shown in Figure 6.
Third step, the linear optimization algorithm of track based on combined iteration: for the waypoint of each fundamental line tuple, first root Linear according to the linear curvature feature fitting of ideal track, then combined iteration is optimal to the fundamental line tuple, reprocesses next basic Line element group, algorithm flow chart are as shown in Figure 9.The linear fitting of track should carry out as follows: determine starting point, a set of segmentation point and terminal Abscissa (M0,M1,M2,M3,M4,Me) after, it, should be true according to following principle in conjunction with the feature of Fig. 1-2 middle orbit horizontal alignment curvature Its fixed curvature value:Wherein, straightway curvature is 0, circular curve segment Curvature is proportional to 1/R, MjThe abscissa of measuring point j is represented, the value of ω is that abscissa is M2To M3Measuring point curvature value average, E (k(M2),k(M2+1),…,k(M3)) represent to k (M2),k(M2+1),…,k(M3) it is a series of number be averaging, by coordinate determine Curvature-mileage curve that each point is fitted after being gradually connected, the deflection-mileage curve being fitted after mileage integral, It is possible thereby to calculate the target function value under each waypoint combination.Use the multistage fundamental line of combined iteration algorithm process whole Tuple, the best linear parameter of available whole track after integration, wherein optimal curvatures and measuring value compare as shown in fig. 7, Optimum orientation angle and the comparison of original measuring value are as shown in Figure 8.
Wherein, the combined iteration algorithm is to be used as the waypoint combination in the certain threshold range of outline waypoint to search Rope domain, all combinations determine the optimal linear algorithms of track in traversal search domain, and steps are as follows: step 1: determining that iteration carries out Region of search range, substitute into outline waypoint track it is linear be initial target solution (dS, X), calculate initial target functional value f (dS, X).Step 2: change waypoint combination in region of search, redefine the linear (dS of trackn, X), calculating target function value f (dSn, X), if f (dSn, X) and < f (dS, X), then receive new explanation, otherwise receives former solution.Step 3: judge whether traversal search domain Interior all waypoints combination is then algorithm termination, otherwise repeatedly step 2.
Wherein, the objective function are as follows: min (f (dS, X))=∑ Δi 2, wherein dS is independent variable, and meaning is track The ordered set { (m of regressive curvature value at each measuring point mileage and the mileage1,k1),(m2,k2),…,(mn,kn), to the edge dS Mileage, which carries out integral, can obtain S, and meaning is that the ordered set of deflection is fitted at each measuring point mileage of track and the mileage {(m11),(m22),…,(mnn), it can differentiate that track is linear by S and dS, calculate track linear parameter (linear length, circle Sweep etc.), the mileage at linear start, end, therefore available all fronts Track desigh parameter;X be metric data mileage, Ordered set the X={ (m of deflection11),(m22),…,(mnn)};Use Δiii' indicate measuring point deflection and intend Close the deviation of deflection.Objective function using minimize each measuring point deflection and be fitted deflection sum of square of deviations as target, The best wiresizing optimization model of orbit plane is constituted together with constraint condition, solution is each for the track for meeting railroad track design specification Optimum orientation angle sequence at measuring point can calculate the best linear parameter of track by it.
Wherein, the constraint condition should meet the linear curvature variation requirement of track and railroad track design rule All kinds of parameter minimum values requirement of model defined, and should to meet front and back length of transition curve equal for symmetrical basic alignment.Therefore Prescribed Properties:
In formula: kZ=0 indicates that the curvature of straightway is 0, kciiThe curvature for indicating i-th section of circular curve is certain value εi, kh(i-1)=a* εi+b(i-1)Indicate that each measuring point curvature is 0 to ε on the preceding easement curve of i-th section of circular curveiBetween linear change, kh(i-1)=-a* εi+b(i+1)Indicate that each measuring point curvature is in ε on the rear easement curve of i-th section of circular curveiTo 0 linear change;Respectively smallest curve radius as defined in design specification, smallest circle length of curve and minimum length of transition curve; lformer,llatterRespectively indicate the length of the two sections of easement curves in same circular curve front and back.

Claims (5)

1. a kind of best linear determining method of orbit plane based on inertial navigation angle measurement, which comprises the following steps:
1) inertial navigation angle measurement is carried out to orbit plane and obtains inertial navigation deflection data sequence, by inertial navigation deflection data sequence edge Route mileage difference obtains the curvature Variation at each test point, and it is inner with the mileage at test point to form orderly curvature- Program column, to obtain curvature-mileage curve;
2) to two segment data of front and back of the test point i in the curvature of step 1)-mileage curve WithIt is fitted respectively using least squares line fitting algorithm, miIndicate measuring point i The mileage at place, kiIt indicates the curvature at measuring point i, calculates the slope differences of the last period matched curve and latter section of matched curve, according to The linear curvature feature of track and curvature curve are based on least square fitting slope variation rule, to it is linear be segmented to obtain it is more Set of segmentation point coordinate, so that orbit plane is divided into multistage fundamental line tuple;
3) to the waypoint of each section of fundamental line tuple, linear according to the linear curvature feature fitting of track first, then combination changes In generation, is optimal to the fundamental line tuple, reprocesses next fundamental line tuple, after finally integrating to multistage fundamental line tuple to obtain the final product To the best linear parameter of whole track, so that it is determined that best linear;Linear step is fitted according to the linear curvature feature of track It include: starting point determining first, a set of segmentation point and terminal abscissa (M0,M1,M2,M3,M4,Me), then in conjunction with orbit plane The feature of linear curvature determines its curvature value according to following principle: Wherein, straightway curvature is 0, and circular curve segment curvature is proportional to 1/R, MjThe abscissa of measuring point j is represented, the value of ω is that abscissa is M2To M3Measuring point curvature value average, E (k (M2),k(M2+1),…,k(M3)) represent to k (M2),k(M2+1),…,k(M3) A series of numbers are averaging, the curvature-mileage curve being fitted after finally each point that coordinate determines gradually is connected, along mileage product Deflection-mileage the curve being fitted after point, and calculate the target function value under each waypoint combination;
The specific steps of combined iteration are as follows: determine the region of search range that iteration carries out first, it is linear to substitute into outline waypoint track For initial target solution (dS, X), initial target functional value f (dS, X) is calculated using objective function, the objective function are as follows: min (f (dS, X))=∑ Δi 2, dS is independent variable, and S can be obtained by integrating to dS along mileage, and wherein S is for each measuring point mileage of track And the ordered set of deflection is fitted at the mileage, X is the ordered set of metric data mileage, deflection, Δiii' table Show measuring point direction angle alphaiWith fitting direction angle alphai' deviation;
Then change waypoint combination in region of search, redefine the linear (dS of trackn, X), calculating target function value f (dSn, X), if f (dSn, X) and < f (dS, X), then receive new explanation, otherwise receives former solution;
Finally judge whether all waypoints combination in traversal search domain, be then algorithm termination, otherwise repeatedly step 2).
2. a kind of best linear determining method of orbit plane based on inertial navigation angle measurement according to claim 1, special Sign is that step 1) the mean curvature variation characteristic should be similar to the linear curvature feature of ideal track, and wherein ideal track is linear Curvature feature is linear for the symmetrical basic model track as made of straight line-easement curve-circular curve-easement curve-Straight Combination, The inertial navigation deflection of straightway remains unchanged, and straightway curvature is 0;Circular curve segment inertial navigation deflection changes linearly, circular curve It is circular curve radius that section curvature, which is steady state value K ∝ a 1/R, R,;Easement curve section inertial navigation deflection changes in second-degree parabola, bent Rate linear change between 0~K with mileage.
3. a kind of best linear determining method of orbit plane based on inertial navigation angle measurement according to claim 2, special Sign is that curvature obtained in the step 1)-mileage curve need to carry out abandoning rough error and median filter process.
4. a kind of best linear determining method of orbit plane based on inertial navigation angle measurement according to claim 1, special The slope differences of the case where sign is, is positive in the step 2) for circular curve curvature, two sections of matched curves are gentle in point of tangent to spiral There is minimum at straight point, has maximum at point of spiral to curve and point of curve to spiral;The case where being negative for circular curve curvature, two sections of fittings are bent The slope differences of line have minimum at point of spiral to curve and point of curve to spiral, have maximum at point of tangent to spiral and point of spiral to tangent.
5. a kind of best linear determining method of orbit plane based on inertial navigation angle measurement according to claim 1, special Sign is, constraint condition is introduced in the objective function and constitutes the best wiresizing optimization model of orbit plane, constraint condition are as follows:
In formula: kZ=0 indicates that the curvature of straightway is 0, kciiThe curvature for indicating i-th section of circular curve is certain value εi, kh(i-1) =a* εi+b(i-1)Indicate that each measuring point curvature is 0 to ε on the preceding easement curve of i-th section of circular curveiBetween linear change, kh(i-1)=- a*εi+b(i+1)Indicate that each measuring point curvature is in ε on the rear easement curve of i-th section of circular curveiTo 0 linear change;Rmin,lminPoint It Wei not smallest curve radius, smallest circle length of curve as defined in design specification and minimum length of transition curve;lformer,llatterPoint The length of the two sections of easement curves in same circular curve front and back is not indicated.
CN201610497850.6A 2016-06-28 2016-06-28 A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement Active CN106126834B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610497850.6A CN106126834B (en) 2016-06-28 2016-06-28 A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610497850.6A CN106126834B (en) 2016-06-28 2016-06-28 A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement

Publications (2)

Publication Number Publication Date
CN106126834A CN106126834A (en) 2016-11-16
CN106126834B true CN106126834B (en) 2019-07-19

Family

ID=57285426

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610497850.6A Active CN106126834B (en) 2016-06-28 2016-06-28 A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement

Country Status (1)

Country Link
CN (1) CN106126834B (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609300B (en) * 2017-09-27 2018-12-14 中南大学 A kind of existing railway planar line position integrative reconstruction design method
CN109857978B (en) * 2019-01-17 2024-03-26 上海上塔软件开发有限公司 Online self-adaptive segmentation algorithm suitable for power curve of household electrical appliance
CN110015445B (en) * 2019-02-15 2020-12-11 北京空间飞行器总体设计部 Earth-moon L2 point Halo track maintaining method
CN110880202B (en) * 2019-12-02 2023-03-21 中电科特种飞机系统工程有限公司 Three-dimensional terrain model creating method, device, equipment and storage medium
CN111340908B (en) * 2020-02-10 2023-07-25 兰州交通大学 Rail electronic map generation method
CN112733223B (en) * 2020-12-28 2021-09-14 中铁二局集团有限公司 Method for back-calculating corresponding line mileage and offset distance by using known coordinate points
CN113722812B (en) * 2021-07-16 2023-12-08 华东交通大学 Calculation method for widened subway circle slow section
CN114169040B (en) * 2021-10-28 2024-04-02 北京自动化控制设备研究所 Existing line design line type fitting method based on inertia and odometer information
CN114111767B (en) * 2021-10-28 2023-11-03 北京自动化控制设备研究所 Method for optimizing line design line type based on multi-information fusion
CN114201850B (en) * 2021-10-29 2024-05-03 北京自动化控制设备研究所 Method for determining line type parameters of existing line type correction design
CN114546330B (en) * 2022-04-26 2022-07-12 成都登临科技有限公司 Function implementation method, approximation interval segmentation method, chip, device and medium
CN114912159B (en) * 2022-07-18 2022-09-13 中国铁路设计集团有限公司 Method for fitting geometric line shape of rail transit line plane
CN115310029B (en) * 2022-10-12 2023-02-03 中国铁建高新装备股份有限公司 Method and device for resolving geometrical parameters of railway theoretical line track

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103399849A (en) * 2013-06-24 2013-11-20 中南大学 Road three-dimensional linear automatic optimization method based on improved particle swarm optimization
CN104636562A (en) * 2015-02-12 2015-05-20 铁道第三勘察设计院集团有限公司 High-speed railway line design method based on train-track system dynamics

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103399849A (en) * 2013-06-24 2013-11-20 中南大学 Road three-dimensional linear automatic optimization method based on improved particle swarm optimization
CN104636562A (en) * 2015-02-12 2015-05-20 铁道第三勘察设计院集团有限公司 High-speed railway line design method based on train-track system dynamics

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Track Irregularities Estimation Based on Acceleration measurements;Wang Y, Qin Y, Wei X;《International Conference on Measurement, Information and Control. IEEE, 2012》;20120520;83-87

Also Published As

Publication number Publication date
CN106126834A (en) 2016-11-16

Similar Documents

Publication Publication Date Title
CN106126834B (en) A kind of best linear determining method of the orbit plane based on inertial navigation angle measurement
CN109798844B (en) Rail profile automatic registration method in dynamic measurement of rail abrasion
CN102251451B (en) Track geometric state measurement system and method based on multisource information fusion technology
CN109977599B (en) Integral intelligent reconstruction method for existing railway longitudinal plane line position
CN109823362A (en) A kind of orbital forcing determination method
CN108820039A (en) A kind of automatic driving vehicle bend crosswise joint system and method
CN105510632B (en) Obtain the method and apparatus of pickup data
CN109655052B (en) Target motion trajectory smoothness evaluation method
CN109594439A (en) A kind of track geometric parameter measurement method being applied in combination based on GNSS/INS and total station
CN110736999B (en) Railway turnout detection method based on laser radar
CN114201850B (en) Method for determining line type parameters of existing line type correction design
CN114912159B (en) Method for fitting geometric line shape of rail transit line plane
CN108229803B (en) Feature road section extraction method and dangerous road section judgment method and system
CN112233233B (en) Three-dimensional linear recovery method for tunnel relative point cloud data
Yang et al. An accurate and fast method to inspect rail wear based on revised global registration
CN107339982A (en) High ferro wire plotting method
CN110530300A (en) A kind of vertical Irregularity detection method of track
CN110017823A (en) Exempt from the tracks positioned Coordinate calculation method and device that horizontalization sets station based on total station
CN114169040B (en) Existing line design line type fitting method based on inertia and odometer information
Zhou et al. Onboard train localization based on railway track irregularity matching
RU2242392C2 (en) Method of and device for correcting errors in location of rail vehicle
CN113609565B (en) Acceptance method and device for static long wave irregularity of railway bridge track
Zhihua et al. A new approach to map-matching and parameter correcting for vehicle navigation system in the area of shadow of GPS signal
CN106123866A (en) Existing Railway Line mileage setting out method based on Aerial Surveying Technology
KR101240578B1 (en) Method and device of measuring the angle of inclination of a plan where a vehicle travels

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant