CN110765638A - Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance - Google Patents

Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance Download PDF

Info

Publication number
CN110765638A
CN110765638A CN201911070465.3A CN201911070465A CN110765638A CN 110765638 A CN110765638 A CN 110765638A CN 201911070465 A CN201911070465 A CN 201911070465A CN 110765638 A CN110765638 A CN 110765638A
Authority
CN
China
Prior art keywords
calculating
deviation
influence
platform
starlight
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.)
Granted
Application number
CN201911070465.3A
Other languages
Chinese (zh)
Other versions
CN110765638B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201911070465.3A priority Critical patent/CN110765638B/en
Publication of CN110765638A publication Critical patent/CN110765638A/en
Application granted granted Critical
Publication of CN110765638B publication Critical patent/CN110765638B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The invention discloses a method for calculating the influence of a disturbance gravitational field on starlight/inertia composite guidance, which comprises the steps of calculating total platform reference deviation according to a vertical deviation resolving error, a geodetic longitude and latitude error and an initial alignment (orientation) error, further obtaining a drop point deviation caused by the platform reference deviation, and obtaining the influence of the disturbance gravitational field on the starlight/inertia composite guidance of a single-star leveling platform according to the influence of the disturbance gravitational force of an active section and a passive section on the drop point deviation. The invention analyzes the influence of the disturbance gravitational field on the single-star leveling platform starlight/inertia composite guidance, and lays a foundation for improving the hit precision of the ballistic missile.

Description

Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance
Technical Field
The invention relates to the field of guidance control, can be applied to single-star leveling station starlight/inertia composite guidance of ballistic missiles/carrier rockets, and is particularly suitable for analyzing the influence of a disturbance gravitational field on the single-star leveling station starlight/inertia composite guidance.
Background
The starlight/inertia composite guidance utilizes the measurement information of the fixed star to calibrate an error angle between a platform coordinate system and a launching inertia coordinate system, corrects the drop point deviation of the ballistic missile according to the error angle, and comprehensively utilizes the inertia guidance and the starlight information to improve the guidance precision of the ballistic missile.
According to the different star modes of the star sensor arranged on the platform, the starlight/inertia composite guidance can be divided into two schemes of adjusting the platform and not adjusting the platform, but on the basis of meeting the requirements of precision and quick maneuvering, the single-star leveling platform scheme has better reliability and feasibility. The single-star-adjustment platform starlight/inertia composite guidance method corrects the drift error of the inertia platform by using the measurement information of the optimal star vector, can greatly improve the guidance precision and the rapid launching capability of the ballistic missile, simultaneously reduces the cost of a weapon system, and has strong environment adaptability.
Factors of the disturbance gravitational field influencing the drop point precision of the ballistic missile mainly include the vertical deviation of a launching point and the disturbance gravitational force in the flight process, and the disturbance gravitational force can be divided into an active section disturbance gravitational force and a passive section disturbance gravitational force. In order to improve the hit precision of the ballistic missile, the disturbance gravitational field must be compensated, and the premise is to analyze the influence of the disturbance gravitational field on the single-star leveling platform starlight/inertia composite guidance.
Disclosure of Invention
The invention aims to solve the technical problem that aiming at the defects of the prior art, the invention provides a method for calculating the influence of a disturbance gravitational field on the single-star leveling platform star light/inertia composite guidance, and the hit precision of a ballistic missile is improved.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: a method for calculating the influence of a disturbance gravitational field on starlight/inertia composite guidance comprises the following steps:
1) astronomical longitude λ from ideal transmission pointTAstronomical latitude BTConversion matrix from platform coordinate system to transmit inertial coordinate system
Figure BDA0002260785790000011
The resolving errors delta ξ and delta η of the vertical deviation and the measuring errors delta lambda of the geodetic longitude and latitude0、ΔB0The platform reference deviation α caused by the initial positioning error is obtained0(ii) a According to the angle of rotation in pitching and yawing directions in the process of leveling the platform
Figure BDA0002260785790000012
Psi, calculating from the initial alignment error [ epsilon ]0xε0yε0z]TInduced platform reference deviation αd
2) Platform reference deviation α due to initial positioning error0And platform reference deviation α due to initial alignment errordCalculating a total platform reference deviation α;
3) calculate the landing offset Δ L due to the platform reference offset αC、ΔHC(ii) a Calculating the influence delta L of the disturbance gravitation of the active segment and the passive segment on the deviation of the drop pointK、ΔHK
4) And calculating the influences delta L and delta H of the disturbance gravitational field on the starlight/inertia composite guidance of the single-star leveling platform.
Platform reference deviation α due to initial positioning error0The calculation process of (2) includes: according to the astronomical longitude and latitude lambda of the ideal transmitting pointT、BTConversion matrix from platform coordinate system to transmit inertial coordinate system
Figure BDA0002260785790000021
Astronomical positioning error delta lambdaT、ΔBTObtaining platform reference deviation α caused by astronomical positioning error0(ii) a Based on the relationship between astronomical positioning error and geodetic positioning error and vertical deviation resolving error when measuring independently, the method obtains the angle from delta lambda0、ΔB0α by Δ ξ, Δ η0
Initial alignment error [ epsilon ]0xε0yε0z]TInduced platform reference deviation αdMeter (2)The calculation formula is as follows:
Figure BDA0002260785790000022
the platform reference deviation α is calculated as α ═ α0d
Landing deviation Δ L due to platform reference deviation αC、ΔHCThe calculation process of (2) is as follows: according to the inertial navigation principle, calculating the navigation error delta v at the shutdown point of the missiletC、ΔrtCCalculating the landing point deviation delta L caused by the platform reference deviation α according to the partial derivative of the landing point deviation to the shutdown point state and based on the navigation error of the missile at the shutdown pointC、ΔHC
Influence delta L of disturbance gravitation of active segment and passive segment on drop point deviationK、ΔHKThe calculation process of (2) includes: calculating the required velocity increment delta v of the disturbance gravitation of the active segment and the passive segmentK1、ΔvK2(ii) a According to the state (v)KK,tK) Calculating the deviation Delta L of the falling point caused by the disturbance gravitation of the active segment and the passive segment according to the partial derivatives of the longitudinal stroke and the transverse stroke to the speedK、ΔHK
The calculation formulas of Δ L and Δ H are:
Figure BDA0002260785790000031
compared with the prior art, the invention has the beneficial effects that: the method calculates the total platform reference deviation according to the vertical deviation calculation error, the geodetic longitude and latitude error and the initial alignment (orientation) error, further obtains the drop point deviation caused by the platform reference deviation, and obtains the influence of the disturbance gravitational field on the single-star leveling stage starlight inertial composite guidance according to the influence of the disturbance gravitational force of the active section and the passive section on the drop point deviation. The invention analyzes the influence of the disturbance gravitational field on the single-star leveling platform starlight/inertia composite guidance, and lays a foundation for improving the hit precision of the ballistic missile.
Drawings
FIG. 1 is a flow chart of the influence of a disturbance gravitational field on single-star modulation platform starlight/inertia composite guidance;
FIG. 2a is a graph of the disturbance gravitational field effect and compensation results for region 1;
FIG. 2b is a graph of the disturbance gravitational field effect and compensation results for region 2;
FIG. 2c is a graph of the disturbance gravitational field effect and compensation results for region 3;
fig. 2d is a graph of the effect of the disturbance gravitational field and the compensation results for region 4.
Detailed Description
For the spacecraft of the embodiment, the specific implementation steps of the invention are as follows:
s1, calculating the influence of the initial positioning error on the platform reference deviation according to the following steps:
(S1-a) according to the astronomical longitude and latitude lambda of the ideal transmitting pointT、BTTransformation matrix from platform coordinate system to emission inertial coordinate system
Figure BDA0002260785790000032
Astronomical positioning error delta lambdaT、ΔBTCalculating platform reference deviation caused by astronomical positioning error α0
Figure BDA0002260785790000033
Note the book
Figure BDA0002260785790000034
Figure BDA0002260785790000041
(S1-b) when the independent measurement is not performed, the relationship between the astronomical positioning error and the geodetic positioning error and the vertical deviation resolving error is
Figure BDA0002260785790000042
Thus can utilize Δ λ0、ΔB0、Δξ, Δ η calculate a platform reference deviation α due to initial positioning error0
Figure BDA0002260785790000043
S2, calculating the error [ epsilon ] from initial alignment (orientation)0xε0yε0z]TInduced platform reference deviation αd
The initial alignment (orientation) error can be represented by misalignment angles of the platform body coordinate system to three axes of the emission inertial coordinate system, and a conversion mode of yaw before pitch is adopted in the platform regulation, namely
Figure BDA0002260785790000044
Wherein psi is the angle of the yaw direction;
Figure BDA0002260785790000045
is the angle of the pitch direction rotation.
S3, calculating total platform reference deviation α
α=α0d(5)
S4, calculating the influence of the initial positioning error on the platform reference deviation according to the following steps:
(S4-a) according to the inertial navigation principle, calculating the navigation error delta v at the shutdown point of the missiletC、ΔrtC
Figure BDA0002260785790000046
Wherein, WPTo look at the acceleration.
(S4-b) calculating a landing deviation DeltaL caused by the platform reference deviation α based on a partial derivative of the landing deviation with respect to the shutdown stateC、ΔHC
Figure BDA0002260785790000051
Wherein,
Figure BDA0002260785790000052
the deviation of the falling point is the partial derivative of the position and the speed of the shutdown point.
S5, calculating the drop point deviation caused by the disturbance gravity of the active section and the passive section according to the following steps:
(S5-a) calculating the required speed increment delta v of the active segment and the passive segmentK1、ΔvK2
Figure BDA0002260785790000053
Wherein,
Figure BDA0002260785790000054
to disregard the required velocity of the passive section disturbing gravity, Freq(vKK,tK) To account for the required speed of the passive segment perturbing gravity.
(S5-b) calculating the deviation of the falling point caused by the disturbance of the gravity by the active segment and the passive segment:
s6, calculating the falling point deviations delta L and delta H caused by the disturbance gravitational field
Figure BDA0002260785790000056
To further illustrate the operation of the present invention, a simulation example is given. The emission points in the simulation are shown in table 1, the emission azimuth angle is traversed to 360 degrees from 0 degrees, and the step length is 30 degrees; the vertical deviation resolving error of the launching point is 3 ' and 3 ', the earth longitude and latitude error is 1 ' and 1 ', the initial alignment (orientation) error is [01000] ' (2.7 sigma), and the disturbance gravity in the calculation considers the X order.
TABLE 1 emission points with large topographic disparity
Figure BDA0002260785790000057
The influence results of the disturbance gravitational field on the single-star leveling platform starlight/inertia composite guidance are shown in the figures 2 a-2 d. For the convenience of comparison, table 2 shows the influence of the disturbance gravitational field on the star light inertial composite guidance of the single-star leveling station when the azimuth angle is 0 °. The influence of the disturbance gravitational field on the single-star leveling station starlight/inertia composite guidance precision is in the order of hundred meters, and the influence of the disturbance gravitational field on the ballistic missile is not only related to the position of a launching point, but also related to a launching azimuth angle.
TABLE 2 influence of disturbance gravitational field on the Star light inertia composite guidance of Single Star levelling station at 0 degree azimuth
Figure BDA0002260785790000061
Reference to the literature
[1] Zhengwei, national soup construction, disturbance gravitational field zhongzhao missile aeromechanics [ M ]. beijing: national defense industry press, 2009.
[2] Wanglie, a strategic missile trajectory rapid forecasting and guidance method based on a state space perturbation method is researched by [ D ]. Changsha: national defense science and technology university, 2018.
[3] An analysis and compensation method for the influence of an earth disturbance gravitational field on the guidance precision of a ballistic missile researches [ D ]. Changsha: national defense science and technology university, 2017.

Claims (7)

1. A method for calculating the influence of a disturbance gravitational field on starlight/inertia composite guidance is characterized by comprising the following steps:
1) astronomical longitude λ from ideal transmission pointTAstronomical latitude BTConversion matrix from platform coordinate system to transmit inertial coordinate system
Figure FDA0002260785780000011
The resolving errors delta ξ and delta η of the vertical deviation and the measuring errors delta lambda of the geodetic longitude and latitude0、ΔB0To obtainTo platform reference deviation α caused by initial positioning error0(ii) a According to the angle of rotation in pitching and yawing directions in the process of leveling the platform
Figure FDA0002260785780000012
Psi, calculating from the initial alignment error [ epsilon ]0xε0yε0z]TInduced platform reference deviation αd
2) Platform reference deviation α due to initial positioning error0And platform reference deviation α due to initial alignment errordCalculating a total platform reference deviation α;
3) calculate the landing offset Δ L due to the platform reference offset αC、ΔHC(ii) a Calculating the influence delta L of the disturbance gravitation of the active segment and the passive segment on the deviation of the drop pointK、ΔHK
4) And calculating the influences delta L and delta H of the disturbance gravitational field on the starlight/inertia composite guidance of the single-star leveling platform.
2. The method for calculating the influence of the disturbance gravitational field on the starlight/inertia composite guidance according to claim 1, wherein the platform reference deviation α caused by the initial positioning error0The calculation process of (2) includes: according to the astronomical longitude and latitude lambda of the ideal transmitting pointT、BTConversion matrix from platform coordinate system to transmit inertial coordinate system
Figure FDA0002260785780000013
Astronomical positioning error delta lambdaT、ΔBTObtaining platform reference deviation α caused by astronomical positioning error0(ii) a Based on the relationship between astronomical positioning error and geodetic positioning error and vertical deviation resolving error when measuring independently, the method obtains the angle from delta lambda0、ΔB0α by Δ ξ, Δ η0
3. The method for calculating the influence of a perturbed gravitational field on a starlight/inertial composite guidance according to claim 1, characterized in that the initial alignment error isDifference [ epsilon ]0xε0yε0z]TInduced platform reference deviation αdThe calculation formula of (2) is as follows:
4. the method for calculating the influence of the disturbance gravitational field on the starlight/inertia composite guidance according to claim 1, wherein the calculation formula of the platform reference deviation α is α - α0d
5. The method for calculating the influence of the disturbance gravitational field on the starlight/inertia composite guidance according to claim 1, wherein the landing point deviation Δ L caused by the platform reference deviation αC、ΔHCThe calculation process of (2) is as follows: according to the inertial navigation principle, calculating the navigation error delta v at the shutdown point of the missiletC、ΔrtCCalculating the landing point deviation delta L caused by the platform reference deviation α according to the partial derivative of the landing point deviation to the shutdown point state and based on the navigation error of the missile at the shutdown pointC、ΔHC
6. The method for calculating the influence of the disturbance gravitational field on the starlight/inertia composite guidance according to claim 1, wherein the influence of the disturbance gravitational field of the active segment and the passive segment on the landing point deviation is Δ LK、ΔHKThe calculation process of (2) includes: calculating the required velocity increment delta v of the disturbance gravitation of the active segment and the passive segmentK1、ΔvK2(ii) a According to the state (v)KK,tK) Calculating the deviation Delta L of the falling point caused by the disturbance gravitation of the active segment and the passive segment according to the partial derivatives of the longitudinal stroke and the transverse stroke to the speedK、ΔHK
7. The method for calculating the influence of the disturbance gravitational field on the starlight/inertia composite guidance according to claim 1, wherein the calculation formulas of Δ L and Δ H are as follows:
Figure FDA0002260785780000022
CN201911070465.3A 2019-11-05 2019-11-05 Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance Active CN110765638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911070465.3A CN110765638B (en) 2019-11-05 2019-11-05 Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911070465.3A CN110765638B (en) 2019-11-05 2019-11-05 Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance

Publications (2)

Publication Number Publication Date
CN110765638A true CN110765638A (en) 2020-02-07
CN110765638B CN110765638B (en) 2024-07-09

Family

ID=69336365

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911070465.3A Active CN110765638B (en) 2019-11-05 2019-11-05 Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance

Country Status (1)

Country Link
CN (1) CN110765638B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111572817A (en) * 2020-06-08 2020-08-25 北京航天自动控制研究所 Optimization calculation method for platform starlight correction coefficient

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5396326A (en) * 1989-04-03 1995-03-07 Northrop Grumman Corporation Two gimbal error averaging astro-inertial navigator
CN107036626A (en) * 2016-11-08 2017-08-11 哈尔滨工业大学 A kind of long-range rocket initial alignment orientation error impact analysis method
CN108507569A (en) * 2017-11-10 2018-09-07 中国人民解放军国防科技大学 Missile-borne fixed star library rapid generation method for starlight/inertia composite guidance

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5396326A (en) * 1989-04-03 1995-03-07 Northrop Grumman Corporation Two gimbal error averaging astro-inertial navigator
CN107036626A (en) * 2016-11-08 2017-08-11 哈尔滨工业大学 A kind of long-range rocket initial alignment orientation error impact analysis method
CN108507569A (en) * 2017-11-10 2018-09-07 中国人民解放军国防科技大学 Missile-borne fixed star library rapid generation method for starlight/inertia composite guidance

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张艳玲;辛万青;王增寿;: "平台单星方案原理及仿真研究", 导弹与航天运载技术, no. 05, 10 October 2007 (2007-10-10), pages 26 - 29 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111572817A (en) * 2020-06-08 2020-08-25 北京航天自动控制研究所 Optimization calculation method for platform starlight correction coefficient

Also Published As

Publication number Publication date
CN110765638B (en) 2024-07-09

Similar Documents

Publication Publication Date Title
CN107966156B (en) Guidance law design method suitable for carrier rocket vertical recovery section
CN111099045B (en) Full physical simulation method for double super satellite dynamics and control air floatation platform
CN101793523B (en) Combined navigation and photoelectric detection integrative system
CN107121929B (en) Robust reentry guidance method based on linear covariance model predictive control
CN105157705B (en) A kind of half strapdown radar seeker line of sight rate extracting method
CN102425980B (en) Control method for realizing overload autopilot by using accelerometer
CN108180910B (en) One kind being based on the uncertain aircraft quick high accuracy method of guidance of aerodynamic parameter
CN109633724B (en) Passive target positioning method based on single-satellite and multi-ground-station combined measurement
CN110672128B (en) Starlight/inertia combined navigation and error online calibration method
CN104089529B (en) Use the method and apparatus that fibre optic gyroscope is calibrated fighter plane armament systems
CN105069311B (en) A kind of long-range rocket emits initial state error propagation method of estimation
CN113900448B (en) Aircraft prediction correction composite guidance method based on sliding mode interference observer
CN105486307A (en) Line-of-sight angular rate estimating method of maneuvering target
CN111591470A (en) Aircraft accurate soft landing closed-loop guidance method adapting to thrust adjustable mode
CN108917764A (en) A kind of Double Satellite only ranging Relative Navigation
CN109724624A (en) A kind of airborne adaptive Transfer alignment algorithm suitable for wing flexure deformation
CN110160523A (en) Initial alignment orientation angle compensation method and vehicular weapons system based on orientation bookbinding
CN109141802A (en) Analogy method for the store Combinations control law in captive trajectory testing
CN110765638A (en) Method for calculating influence of disturbance gravitational field on starlight/inertia composite guidance
CN105354380B (en) The glide trajectories rapid correction method of compensation is influenced towards perturbation factors
CN110895418A (en) Low-speed rotating aircraft control method and system for compensating dynamic lag of steering engine
RU2423658C2 (en) Method of controlling and stabilising mobile carrier, integrated system, device for turning antenna reflector in two mutually-perpendicular planes and device for actuating differential aerodynamic controllers for realising said method
CN203928892U (en) The equipment that uses fibre optic gyroscope to calibrate fighter plane armament systems
CN114111805B (en) High-precision alignment method for carrier rocket multisource multiclass measurement data position references
CN115326008A (en) Dynamic estimation method for attitude error and time delay of shipborne inertial navigation system based on fixed star observation

Legal Events

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