METHOD AND SYSTEM FOR DETERMINING A POSITION ERROR OF AN AIRCRAFT

Abstract

A method and system for determining a position error of an aircraft. The system comprises a set of collection modules implemented in an iterative manner, including an offset collection module, a position vector collection module, a set of sequentially implemented computation modules, including an angular alignment bias computation module, an observable position bias computation module and a position error computation module. The system thus allows a position error to be determined that is used to perform a landing in good conditions.

Claims

1. A method for determining a position error of an aircraft during an approach procedure of the aircraft with a view to landing on a runway, comprising at least the following steps implemented in an iterative manner: a first set of collection steps implemented in each iteration, the first set of steps comprising: an offset collection step, implemented by an offset collection module, involving collecting a current offset measured by an offset measurement module based on a reference guidance signal transmitted by a transmitting station of a landing assistance device, the landing assistance device comprising at least the ground-based transmitting station configured to transmit the reference guidance signal defining a reference guidance axis, the runway having a longitudinal axis and a threshold; a position vector collection step, implemented by a position vector collection module, involving collecting a current position vector of the aircraft measured by a position vector measurement module based on signals transmitted by a geopositioning system, the current position vector comprising a current longitudinal position and at least one from among a current lateral position and a current vertical position; a second set of computation steps comprising: a first computation step, implemented by a first computation module in a first iteration, involving computing an angular alignment bias of the reference guidance axis in relation to the approach axis of the runway based on the current offset, the current longitudinal position and said at least one from among the current lateral position and the current vertical position; a second computation step, implemented by a second computation module in a second iteration subsequent to the first iteration, involving computing an observable position bias based on the angular alignment bias computed in the first iteration, the current offset, the current longitudinal position and said at least one from among the current lateral position and the current vertical position; a third computation step, implemented by a third computation module in a third iteration subsequent to the second iteration, involving computing the position error based on the observable position bias computed in the second iteration, the current longitudinal position and the longitudinal position measured in the first iteration; a first transmission step, implemented by a first transmission module, involving transmitting the position error to a first user device.

2. The method according to claim 1, further comprising a filtering step, implemented by a filtering module, involving filtering at least one of the following two biases: the angular alignment bias computed in the first computation step; and the observable position bias computed in the second computation step.

3. The method according to claim 1, further comprising: a fourth computation step, implemented by a fourth computation module, involving computing a corrected angular alignment bias based on the current offset, the current longitudinal position, said at least one from among the current lateral position and the current vertical position and the position error computed in the third computation step; a second transmission step, implemented by a second transmission module, involving transmitting the corrected angular alignment bias to a second user device.

4. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current lateral position, the angular alignment bias is computed in the first computation step using the following relation: $b_{{\mathrm{LOC}}_{\mathrm{align}}}=\arcsin \left(\frac{b_{{\mathrm{LOC}}_{\mathrm{offset}}}-\frac{{\eta }_{{\mathrm{LOC}}_{\mu A}}\times 0.7\left(L-x_{\mathrm{RWY}}\right)}{L}-y_{\mathrm{RWY}}}{L-x_{\mathrm{RWY}}}\right),$ wherein: b.sub.LOC.sub.offset corresponds to the lateral position of the transmitting station in relation to the longitudinal axis of the runway; η.sub.LOC.sub.μA corresponds to a current lateral offset measured in microamperes by the offset measurement module; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; y.sub.RWY corresponds to the current lateral position collected by the position vector collection module; L corresponds to the longitudinal position of the transmitting station in relation to the threshold of the runway.

5. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current lateral position, an observable lateral position bias is computed in the second computation step using the following relation: $b_{y_{\mathrm{obs}}}=y_{\mathrm{RWY}}-b_{{\mathrm{LOC}}_{\mathrm{offset}}}-\left(L-x_{\mathrm{RWY}}\right)\sin \left(b_{{\mathrm{LOC}}_{\mathrm{align}}}\right)-\frac{{\eta }_{{\mathrm{LOC}}_{\mu A}}\times 0.7\times \left(L-x_{\mathrm{RWY}}\right)}{L},$ wherein: η.sub.GS.sub.rad corresponds to a current lateral offset measured in radians by the offset measurement module; b.sub.LOC.sub.offset corresponds to the lateral position of the transmitting station in relation to the longitudinal axis of the runway; b.sub.LOC.sub.align corresponds to the angular alignment bias computed in the first computation step; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; y.sub.RWY corresponds to the current lateral position collected by the position vector collection module; L corresponds to the longitudinal position of the transmitting station in relation to the threshold of the runway.

6. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current lateral position, a lateral position error is computed in the third computation step in one of the following two manners: by computing a temporal gradient of the observable lateral position bias, then extrapolating the value of a bias at a transmitting beacon as a function of said temporal gradient; or using a spatial extrapolation, in particular defined by the following relation: $b_y=b_{y_{\mathrm{obs}}}\frac{L+x_{b_{{\mathrm{LOC}}_{\mathrm{align}}}}}{x_{\mathrm{RWY}}-x_{b_{{\mathrm{LOC}}_{\mathrm{align}}}}},$ wherein: $x_{b_{{\mathrm{LOC}}_{\mathrm{align}}}}$ corresponds to the longitudinal position measured by a position vector module during the first iteration; b.sub.y.sub.obs corresponds to the observable lateral position bias computed in the second step; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; L corresponds to the longitudinal position of the transmitting station in relation to the threshold of the runway.

7. The method according to claim 3, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current lateral position, the corrected angular alignment bias is computed in the fourth computation step using the following relation: $b_{{\mathrm{LOC}}_{{\mathrm{align}}_{\mathrm{TRUE}}}}=\arcsin \left(\frac{b_{{\mathrm{LOC}}_{\mathrm{offset}}}-\frac{{\eta }_{{\mathrm{LOC}}_{\mu A}}\times 0.7\left(L-x_{\mathrm{RWY}}\right)}{L}-y_{\mathrm{RWY}}+b_y}{L-x_{\mathrm{RWY}}}\right).$

8. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current vertical position, the angular alignment bias is computed in the first computation step using the following relation: $b_{\mathrm{GPA}}=\frac{\begin{array}{c}\mathrm{GPA}-\mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\frac{0.12\times \mathrm{GPA}\times {\eta }_{{\mathrm{GC}}_{\mathrm{ddm}}}}{0.0875}-\\ \arctan \left(\frac{z_{\mathrm{RWY}}-b_z}{x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan \mathrm{GPA}}}\right)\end{array}}{1+\frac{0.12\times \mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\times {\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}}{0.0875}},$ wherein: GPA corresponds to an angle of the approach axis with the longitudinal axis of the runway; η.sub.GS.sub.ddm corresponds to the current vertical offset measured as DDM by the offset measurement module; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; z.sub.RWY corresponds to the current vertical position collected by the position vector collection module; TCH corresponds to a height at which the vertical offset is zero at a position in line with the threshold of the runway.

9. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current vertical position, the observable vertical position bias is computed in the second computation step using the following relation: $b_{z_{\mathrm{obs}}}=z_{\mathrm{RWY}}-\left(x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan \left(\mathrm{GPA}\right)}\right)\tan \left(\mathrm{GPA}-b_{\mathrm{GPA}}-{\eta }_{{\mathrm{GS}}_{\mathrm{rad}}}\right),$ wherein: η.sub.GS.sub.rad corresponds to the current vertical offset measured in radians by the offset measurement module; b.sub.GPA corresponds to the angular alignment bias computed in the first computation step; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; z.sub.RWY corresponds to the current vertical position collected by the position vector collection module.

10. The method according to claim 1, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current vertical position, the vertical position error is computed in the third computation step in one of the following two manners: by computing a temporal gradient of the observable vertical position bias, then extrapolating the value of a bias at a transmitting beacon as a function of said temporal gradient; or using a spatial extrapolation, in particular defined by the following relation: $b_z=b_{z_{\mathrm{obs}}}\frac{L_{\mathrm{GS}}+x_{b_{\mathrm{GPA}}}}{x_{\mathrm{RWY}}-x_{b_{\mathrm{GPA}}}},$ wherein: x.sub.b.sub.GPA corresponds to the longitudinal position measured by a position vector module during the first iteration; b.sub.GPA corresponds to the angular alignment bias computed in the first computation step; x.sub.RWY corresponds to the current longitudinal position collected by the position vector collection module; L.sub.GS corresponds to a longitudinal distance between the transmitting station and the threshold of the runway.

11. The method according to claim 3, wherein, when said at least one from among the current lateral position and the current vertical position corresponds to the current vertical position, the corrected angular alignment bias is computed in the fourth step using the following relation: $b_{{\mathrm{GPA}}_{\mathrm{TRUE}}}=\frac{\begin{array}{c}\mathrm{GPA}-\mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\frac{0.12\times \mathrm{GPA}\times {\eta }_{{\mathrm{GC}}_{\mathrm{ddm}}}}{0.0875}-\\ \arctan \left(\frac{z_{\mathrm{RWY}}-b_z}{x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan \mathrm{GPA}}}\right)\end{array}}{1+\frac{0.12\times \mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\times {\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}}{0.0875}}.$

12. A system for determining a position error of an aircraft during an approach procedure of the aircraft with a view to landing on a runway, comprises at least the following modules implemented in an iterative manner: a first set of collection modules implemented in each iteration, the first set of modules comprising: an offset collection module configured to collect a current offset measured by an offset measurement module based on a reference guidance signal transmitted by a transmitting station of a landing assistance device, the landing assistance device comprising at least the ground-based transmitting station configured to transmit the reference guidance signal defining a reference guidance axis, the runway having a longitudinal axis and a threshold; a position vector collection module configured to collect a current position vector of the aircraft measured by a position vector measurement module based on signals transmitted by a geopositioning system, the current position vector comprising a current longitudinal position and at least one from among a current lateral position and a current vertical position; a second set of computation modules comprising: a first computation module, implemented in a first iteration, configured to compute an angular alignment bias of the reference guidance axis in relation to the approach axis of the runway based on the current offset, the current longitudinal position and said at least one from among the current lateral position and the current vertical position; a second computation module, implemented in a second iteration subsequent to the first iteration, configured to compute an observable position bias based on the angular alignment bias computed in the first iteration, the current offset, the current longitudinal position and said at least one from among the current lateral position and the current vertical position; a third computation module, implemented in a third iteration subsequent to the second iteration, configured to compute the position error based on the observable position bias computed in the second iteration, the current longitudinal position and the longitudinal position measured in the first iteration; a first transmission module configured to transmit the position error to a first user device.

13. The system according to claim 12, further comprising at least one of the following filtering modules: a first filtering module configured to filter the angular alignment bias computed by the first computation module; and a second filtering module configured to filter the observable position bias computed by the second computation module.

14. The system according to claim 12, further comprising: a fourth computation module configured to compute a corrected angular alignment bias based on the current offset, the current longitudinal position, said at least one from among the current lateral position and the current vertical position and the position error computed by the third computation module; a second transmission module configured to transmit the corrected angular alignment bias to a second user device.

15. An aircraft comprising a system for determining a position error, according to claim 12.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0097] The invention, with its features and advantages, will become more clearly apparent from reading the description, which is provided with reference to the accompanying drawings, in which:

[0098] FIG. 1 schematically shows the determination system;

[0099] FIG. 2 shows a profile view of an aircraft with the determination system on board;

[0100] FIG. 3 schematically shows the determination method;

[0101] FIG. 4A shows a top view of a runway and of an aircraft in an approach phase with a view to landing on the runway;

[0102] FIG. 4B shows a profile view of a runway and of an aircraft in an approach phase with a view to landing on the runway;

[0103] FIG. 5A shows a top view of a runway and of an aircraft at positions measured by the position vector measurement module and actual positions of the aircraft;

[0104] FIG. 5B shows a profile view of a runway and of an aircraft at positions measured by the position vector measurement module and actual positions of the aircraft.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0105] FIG. 1 shows the system 1 for determining a position error of an aircraft AC, for example, a transport plane, in particular, during an approach phase with a view to landing on a runway RWY. The position error is determined using a landing assistance device.

[0106] The runway RWY has a longitudinal axis B1 and a threshold T (FIG. 4A).

[0107] The determination system 1 comprises at least the following modules implemented in an iterative manner:

[0108] a first set of collection modules implemented in each iteration of a set of iterations i, i+N1, i+N1+N2 (with N1 and N2 being strictly positive whole numbers); and

[0109] a second set of computation modules.

[0110] The first set of modules comprises:

[0111] an offset collection module 2 COLL1; and

[0112] a position vector collection module 3 COLL2.

[0113] Each of the modules 2 and 3 of the first set of modules is implemented in each iteration. An offset or a position is considered to be current when it is collected in the current iteration.

[0114] The offset collection module 2 is configured to collect a current offset measured by an offset measurement module 21 MEAS1 based on a reference guidance signal transmitted by the ground-based transmitting station 22 of a landing assistance device.

[0115] The landing assistance device can correspond to a device of the ILS system type. The landing assistance device comprises at least one ILS transmitting station 22 (ground-based) of a runway RWY configured to transmit a reference guidance signal. This guidance signal defines a reference guidance axis A.

[0116] In a first embodiment, the considered reference guidance axis A is a lateral guidance axis. The transmitting station 22 is then a station of the LOC (“Localizer beacon”) type configured to transmit a reference lateral guidance signal. This lateral guidance signal defines the reference lateral guidance axis A. The offset measured by the offset measurement module 21 is then a current lateral offset η.sub.LOC.sub.μA.

[0117] The current lateral offset η.sub.LOC.sub.μA can be measured in microamperes. The lateral offset is zero when the aircraft follows a path aligned with the reference lateral guidance axis A.

[0118] In a second embodiment, the considered reference guidance axis A is a vertical guidance axis. The transmitting station 22 is then a station of the Glide (“Glide slope beacon”) type configured to transmit a reference vertical guidance signal, corresponding to a signal of the “Glide Slope signal” type. This vertical guidance signal defines the reference vertical guidance axis A. The offset measured by the offset measurement module 21 is then a current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad.

[0119] The current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad can be measured in radians in relation to the reference vertical guidance axis A. It also can be measured as DDM (“Difference in Depth of Modulation”). In practice, it is expressed as a percentage or in microamperes. This vertical guidance signal has two parts. Each of the parts is transmitted at a different frequency to one another. The DDM corresponds to a difference between the depth of modulation of a part transmitted at one frequency and the depth of modulation of the other part transmitted at another frequency. The vertical offset in radians or as DDM is zero when the aircraft follows a path aligned with the reference vertical guidance axis A.

[0120] The position vector collection module 3 is configured to collect a current position vector of the aircraft AC measured by a position vector measurement module 31 MEAS2 based on signals transmitted by a geopositioning system 32 GNSS. Advantageously, the position vector is expressed in a coordinate system associated with the runway RWY. The current position vector comprises at least one current longitudinal position x.sub.RWY and a current lateral position y.sub.RWY in the first embodiment (respectively a current vertical position z.sub.RWY in the second embodiment). The geopositioning system can correspond to a geolocation and navigation system using a GNSS (“Global Navigation Satellite System”), such as a GPS (“Global Positioning System”).

[0121] The determination system 1 also comprises a second set of computation modules comprising:

[0122] a computation module 4 CALC1 implemented in a first iteration i;

[0123] a computation module 5 CALC2 implemented in a second iteration i+N1 subsequent to the first iteration i;

[0124] a computation module 7 CALC3 implemented in a third iteration i+N1+N2 subsequent to the second iteration i+N1.

[0125] Therefore, the modules 4, 5 and 7 are implemented sequentially.

[0126] Numbers N1 and N2 are selected so that the distances covered by the aircraft between iterations i and i+N1, on the one hand, and between iterations i+N1 and i+N1+N2, on the other hand, are sufficient to allow sufficient precision for the various computations. For example, if iteration i corresponds to a distance from the aircraft to the runway threshold that is equal to 5 nautical miles (approximately 9.3 km), the number N1 can be selected so that iteration i+N1 corresponds to a distance from the aircraft to the runway threshold that is equal to 4 nautical miles (approximately 7.4 km) and the number N2 can be selected so that iteration i+N1+N2 corresponds to a distance from the aircraft to the runway threshold that is equal to 3 nautical miles (approximately 5.6 km). Numbers N1 and N2 are particularly determined as a function of a computation cycle time (corresponding to the transition from an iteration i to an iteration i+1) of the one or more avionics computers implementing modules 4, 5 and 7, so as to correspond to the aforementioned distances.

[0127] In the first embodiment, the computation module 4 is configured to compute an angular alignment bias b.sub.LOC.sub.align of the reference lateral guidance axis A in relation to the longitudinal axis B1 of the runway RWY based on the current lateral offset η.sub.LOC.sub.μA, the current longitudinal position x.sub.RWY and the current lateral position y.sub.RWY. The longitudinal axis B1 of the runway corresponds to the lateral component of the approach axis of the runway.

[0128] The longitudinal axis B1 of the runway corresponds to an axis separating the runway RWY into two substantially identical parts along the length of the runway RWY. The threshold T of the runway RWY corresponds to the origin of an orthonormal frame of reference (T, .sub.RWY, .sub.RWY, .sub.RWY) covered by the longitudinal axis B1 of the runway RWY.

[0129] In practice, the actual reference guidance axis is not perfectly straight, in particular due to multiple paths of the guidance signals transmitted by the station 22 due to the environment. The reference guidance axis A considered in the remainder of the invention is straight, as shown in FIGS. 4A and 4B. It corresponds to an approximation of the actual reference guidance axis that is obtained, for example, by an average, by smoothing or by a linear interpolation of the actual reference guidance axis.

[0130] According to the preferred embodiment, the angular alignment bias b.sub.LOC.sub.align is computed by the computation module 4 using the following relation:

[00011]$b_{LO{C}_{\mathrm{align}}}=\arcsin \left(\frac{b_{{\mathrm{LOC}}_{\mathrm{offse}t}}-\frac{{\eta }_{{\mathrm{LOC}}_{\mu A}\times 0.7\left(L-x_{\mathrm{RWY}}\right)}}{L}-y_{\mathrm{RWY}}}{L-x_{\mathrm{RWY}}}\right),$

[0131] wherein:

[0132] b.sub.LOC.sub.offset corresponds to the lateral position (in meters) of the transmitting station 22 in relation to the longitudinal axis B1 of the runway RWY;

[0133] η.sub.LOC.sub.μA corresponds to the current lateral offset measured in microamperes by the offset measurement module 21;

[0134] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0135] y.sub.RWY corresponds to the current lateral position collected (in meters) by the position vector collection module 3;

[0136] L corresponds to the longitudinal position (in meters) of the transmitting station 22 in relation to the threshold T of the runway RWY.

[0137] The angular alignment bias b.sub.LOC.sub.align is tainted by the lateral position error by generated by the position vector measurement module 31. The value of the angular alignment bias b.sub.LOC.sub.align that is computed thus is set for the remainder of the computations.

[0138] The computation module 5 is configured to compute an observable lateral position bias b.sub.y.sub.obs based on the angular alignment bias b.sub.LOC.sub.align computed in the first iteration i, the current lateral offset η.sub.LOC.sub.align the current longitudinal position x.sub.RWY and the current lateral position y.sub.RWY.

[0139] According to the preferred embodiment, the observable lateral position bias b.sub.y.sub.obs is computed by the computation module 5 using the following relation:

[00012]$b_{y_{obs}}=y_{RWY}-b_{LO{C}_{\mathrm{offse}t}}-\left(L-x_{RWY}\right)\sin \left(b_{LO{C}_{\mathrm{align}}}\right)-\frac{{{\eta }_{\mathrm{LOC}}}_{\mu A}\times 0.7\times \left(L-x_{RWY}\right)}{L},$

[0140] wherein:

[0141] η.sub.LOC.sub.μA corresponds to the current lateral offset measured in microamperes by the offset measurement module 21;

[0142] b.sub.LOC.sub.offset corresponds to the lateral position (in meters) of the transmitting station 22 in relation to the longitudinal axis B of the runway RWY;

[0143] b.sub.LOC.sub.align corresponds to the angular alignment bias (in radians) computed by the computation module 4;

[0144] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0145] y.sub.RWY corresponds to the current lateral position collected (in meters) by the position vector collection module 3;

[0146] L corresponds to the longitudinal position (in meters) of the transmitting station 22 in relation to the threshold T of the runway RWY.

[0147] According to an alternative embodiment, the determination system 1 further comprises a filtering module 6A FILT1 configured to filter the angular alignment bias b.sub.LOC.sub.align computed by the computation module 4.

[0148] The filtering is computed, for example, on k1 successive iterations, with 0<k1<N2. This filtering can correspond to a computation of the average of the angular alignment bias b.sub.LOC.sub.align based on angular alignment bias b.sub.LOC.sub.align computed in previous iterations.

[0149] According to another alternative embodiment, the determination system 1 further comprises a filtering module 6B FILT2 configured to filter the observable lateral position bias b.sub.y.sub.obs computed by the computation module 5.

[0150] The filtering is computed, for example, on m1 successive iterations, with 0<m1. This filtering can correspond to a computation of the average of the observable lateral position bias b.sub.y.sub.obs based on observable lateral position bias b.sub.y.sub.obs computed in previous iterations.

[0151] The computation module 7 is configured to compute the lateral position error by based on the observable lateral position bias b.sub.y.sub.obs computed in the second iteration i+N1, the current longitudinal position x.sub.RWY and the longitudinal position measured in the first iteration i.

[0152] According to the preferred embodiment, the lateral position error b.sub.y is computed by the computation module 7 using a spatial extrapolation that can be defined, in a non-limiting manner, by the following relation:

[00013]$b_y=b_{y_{obs}}\frac{L+x_{b_{LO{C}_{\mathrm{align}}}}}{x_{RWY}-x_{b_{LO{C}_{\mathrm{align}}}}},$

[0153] wherein:

[00014]$X_{b_{LO{C}_{\mathrm{align}}}}$

corresponds to the longitudinal position measured (in meters) by the position vector measurement module 31 during the first iteration i;

[0154] b.sub.y.sub.obs corresponds to the observable lateral position bias (in meters) computed by the computation module 5;

[0155] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0156] L corresponds to the longitudinal position (in meters) of the transmitting station 22 in relation to the threshold T of the runway RWY.

[0157] In an alternative embodiment, the lateral position error b.sub.y is determined by computing a gradient of the observable lateral position b.sub.y.sub.obs. This gradient corresponds to a temporal gradient. The value of the position bias is then extrapolated in relation to the time remaining until landing.

[0158] The determination system 1 further comprises a transmission module 8 TRANS1 configured to transmit the lateral position error b.sub.y to a user device 81 US1.

[0159] The determination system 1 can also comprise:

[0160] a computation module 9 CALC4; and

[0161] a transmission module 10 TRANS2.

[0162] The computation module 9 is configured to compute a corrected angular alignment bias

[00015]$b_{LO{C}_{{\mathrm{align}}_{TRUE}}}$

based on the current lateral offset η.sub.LOC.sub.μA, the current longitudinal position x.sub.RWY, the current lateral position y.sub.RWY and the lateral position error b.sub.y computed by the computation module 7.

[0163] The corrected angular alignment bias b

[00016]$b_{LO{C}_{{\mathrm{align}}_{TRUE}}}$

can be computed by the computation module 9 using the following relation:

[00017]$b_{LO{C}_{{\mathrm{align}}_{TRUE}}}=\arcsin \left(\frac{b_{{\mathrm{LOC}}_{\mathrm{offse}t}}-\frac{{\eta }_{{\mathrm{LOC}}_{\mu A}\times 0.7\left(L-x_{\mathrm{RWY}}\right)}}{L}-y_{\mathrm{RWY}}+b_y}{L-x_{\mathrm{RWY}}}\right).$

[0164] The transmission module 10 is configured to transmit the corrected angular alignment bias

[00018]$\left(b_{LO{C}_{{\mathrm{align}}_{TRUE}}}\right)$

to a second user device 101 US2.

[0165] Determining the lateral position error b.sub.y makes it possible to know whether a realignment of the aircraft AC is necessary. If the lateral position error b.sub.y is non-zero, a realignment is necessary, either by the pilot of the aircraft AC, or by the automatic pilot of the aircraft AC. This can also indicate a failure in the reference lateral guidance signal or erroneous data from information from the database that is used that relates to the runway RWY. In this case, a “go round’ may need to be performed if the value of the lateral position error b.sub.y is too high.

[0166] FIG. 5A shows several positions P1, P2, P3, P4 and P5 of an aircraft AC. Positions P1 and P3 represent positions determined in two different iterations (with t.sub.0 corresponding to the first iteration i and t.sub.1 corresponding to the second iteration i+N1) and by the position vector measurement module 31. Positions P2 and P4 represent actual positions of the aircraft AC in the same iterations. The determination system 1 allows the lateral position error b.sub.y to be determined between positions P1 and P2 and between positions P3 and P4.

[0167] Position P5 in the second iteration i+N1 corresponds to a position determined with the observable lateral position bias b.sub.y.sub.obs by considering the angular alignment bias b.sub.LOC.sub.align computed in the first iteration i.

[0168] In the second embodiment, the computation module 4 is configured to compute an angular alignment bias b.sub.GPA of the reference vertical guidance axis A in relation to a nominal approach descent axis B2 based on the current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad, the current longitudinal position x.sub.RWY and the current vertical position z.sub.RWY. The axis of the nominal approach descent path B2, as shown in FIG. 4B, corresponds to the vertical component of the approach axis of the runway. It generally forms an angle GPA of 3° with the longitudinal axis B1 of the runway RWY.

[0169] According to a preferred embodiment, the angular alignment bias b.sub.GPA is computed by the computation module 4 using the following relation:

[00019]$b_{\mathrm{GPA}}=\frac{\begin{array}{c}\mathrm{GPA}-\mathrm{sign}\left({\eta }_{G{S}_{ddm}}\right)\frac{0.12\times GPA\times {\eta }_{G{C}_{ddm}}}{0.0875}-\\ \arctan \left(\frac{z_{\mathrm{RWY}}}{x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan GPA}}\right)\end{array}}{1+\frac{0.12\times \mathrm{sign}\left({\eta }_{G{S}_{ddm}}\right)\times {\eta }_{G{S}_{ddm}}}{0.0875}},$

[0170] wherein:

[0171] GPA corresponds to an angle (in radians) of the nominal descent approach axis B2 with the longitudinal axis B1 of the runway RWY;

[0172] η.sub.GS.sub.ddm corresponds to the current vertical offset measured as DDM by the offset measurement module 21;

[0173] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0174] z.sub.RWY corresponds to the current vertical position collected (in meters) by the position vector collection module 3;

[0175] TCH corresponds to a height at which the vertical offset is zero at a position in line with the threshold T of the runway RWY.

[0176] The angular alignment bias b.sub.GPA is tainted by the vertical position error b.sub.z generated by the position vector measurement module 31. The value of the angular alignment bias b.sub.GPA that is thus computed is set for the remainder of the computations.

[0177] The computation module 5 is configured to compute an observable vertical position bias b.sub.z.sub.obs based on the angular alignment bias b.sub.GPA computed in the first iteration i, the current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad, the current longitudinal position x.sub.RWY and the current vertical position z.sub.RWY.

[0178] According to the preferred embodiment, the observable vertical position bias b.sub.z.sub.obs is computed by the computation module 5 using the following relation:

[00020]$b_{z_{\mathrm{obs}}}=z_{\mathrm{RWY}}-\left(x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan \left(\mathrm{GPA}\right)}\right)\tan \left(\mathrm{GPA}-b_{\mathrm{GPA}}-{\eta }_{{\mathrm{GS}}_{\mathrm{rad}}}\right),$

[0179] wherein:

[0180] η.sub.GS.sub.rad corresponds to the current vertical offset measured in radians by the offset measurement module 21;

[0181] b.sub.GPA corresponds to the angular alignment bias (in radians) computed by the computation module 4;

[0182] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0183] z.sub.RWY corresponds to the current vertical position collected (in meters) by the position vector collection module 3.

[0184] According to an alternative embodiment, the determination system 1 further comprises a filtering module 6A FILT1 configured to filter the angular alignment bias b.sub.GPA computed by the computation module 4.

[0185] The filtering is computed, for example, on k2 successive iterations, with 0<k2<N2. This filtering can correspond to a computation of the average of the angular alignment bias b.sub.GPA based on angular alignment bias b.sub.GPA computed in previous iterations.

[0186] According to an alternative embodiment, the determination system 1 further comprises a filtering module 6B FILT2 configured to filter the observable vertical position bias b.sub.z.sub.obs computed by the computation module 5.

[0187] The filtering is computed, for example, on m2 successive iterations, with 0<m2. This filtering can correspond to a computation of the average of the observable vertical position bias b.sub.z.sub.obs based on observable vertical position bias b.sub.z.sub.obs computed in previous iterations.

[0188] The computation module 7 is configured to compute the vertical position error b.sub.z based on the observable vertical position bias b.sub.z.sub.obs computed in the second iteration i+N1, the current longitudinal position x.sub.RWY and the longitudinal position measured in the first iteration i.

[0189] According to the preferred embodiment, the vertical position error b.sub.z is computed by the computation module 7 using a spatial extrapolation that can be defined, in a non-limiting manner, by the following relation:

[00021]$b_z=b_{z_{\mathrm{obs}}}\frac{L_{\mathrm{GS}}+x_{b_{\mathrm{GPA}}}}{x_{\mathrm{RWY}}-x_{b_{\mathrm{GPA}}}},$

[0190] wherein:

[0191] x.sub.b.sub.GPA corresponds to the longitudinal position measured (in meters) by the position vector measurement module 31 during the first iteration i;

[0192] b.sub.GPA corresponds to the angular alignment bias computed (in radians) by the computation module 4 in the first iteration i;

[0193] x.sub.RWY corresponds to the current longitudinal position collected (in meters) by the position vector collection module 3;

[0194] L.sub.GS corresponds to a longitudinal distance (in meters) between the transmitting station and the threshold of the runway.

[0195] The longitudinal distance L.sub.GS between the transmitting station 22 and the threshold T of the runway RWY can be computed based on the following relation:

[00022]$L_{\mathrm{GS}}=\frac{\mathrm{TCH}}{\tan \left(\mathrm{GPA}-{\hat{b}}_{\mathrm{GPA}}\right)},$

[0196] wherein {circumflex over (b)}.sub.GPA corresponds to an estimate of the angular alignment bias b.sub.GPA.

[0197] In an alternative embodiment, the vertical position error b.sub.z is determined by computing an observable vertical position gradient b.sub.z.sub.obs. This gradient corresponds to a temporal gradient in relation to the time remaining until landing. The value of the position bias is subsequently extrapolated.

[0198] The determination system 1 further comprises a transmission module 8 TRANS1 configured to transmit the vertical position error b.sub.z to a user device 81 US1.

[0199] The determination system 1 can also comprise:

[0200] a computation module 9 CALC4; and

[0201] a transmission module 10 TRANS2.

[0202] The computation module 9 is configured to compute a corrected angular alignment bias b.sub.GPA.sub.TRUE based on the current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad, the current longitudinal position x.sub.RWY, the current vertical position z.sub.RWY and the vertical position error b.sub.z computed by the computation module 7.

[0203] The corrected angular alignment bias b.sub.GPA.sub.TRUE can be computed by the computation module 9 using the following relation:

[00023]$b_{{\mathrm{GPA}}_{\mathrm{TRUE}}}=\frac{\begin{array}{c}\mathrm{GPA}-\mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\frac{0.12\times \mathrm{GPA}\times {\eta }_{{\mathrm{GC}}_{\mathrm{ddm}}}}{0.0875}-\\ \arctan \left(\frac{z_{\mathrm{RWY}}-b_z}{x_{\mathrm{RWY}}+\frac{\mathrm{TCH}}{\tan \mathrm{GPA}}}\right)\end{array}}{1+\frac{0.12\times \mathrm{sign}\left({\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}\right)\times {\eta }_{{\mathrm{GS}}_{\mathrm{ddm}}}}{0.0875}}$

[0204] The transmission module 10 is configured to transmit the corrected angular alignment bias b.sub.GPA.sub.TRUE to a user device 101 US2.

[0205] Determining the vertical position error b.sub.z makes it possible to know whether a realignment of the aircraft AC is necessary. If the vertical position error b.sub.z is non-zero, a realignment is necessary, either by the pilot of the aircraft AC, or by the automatic pilot of the aircraft AC. This can also indicate a failure in the reference vertical guidance signal or erroneous data from information from the database that is used that relates to the runway RWY. In this case, a “go round” may need to be performed if the value of the vertical position error b.sub.z is too high.

[0206] FIG. 5B shows several positions P1′, P2′, P3′, P4′ and P5′ of an aircraft AC. Positions P1′ and P3′ represent positions determined in two different iterations (with t.sub.0 corresponding to the first iteration i and t.sub.1 corresponding to the second iteration i+N1) and by the position vector measurement module 31. Positions P2′ and P4′ represent actual positions of the aircraft AC in the same iterations. The determination system 1 allows the vertical position error b.sub.z to be determined between positions P1′ and P2′ and between positions P3′ and P4′.

[0207] Position P5′ in the second iteration i+N1 corresponds to a position determined with the observable vertical position bias b.sub.z.sub.obs by considering the angular alignment bias b.sub.GPA computed in the first iteration i.

[0208] The invention also relates to a method for determining a position error of an aircraft AC (FIG. 3).

[0209] The determination method comprises at least the following steps implemented in an iterative manner:

[0210] a first set of collection steps implemented in each iteration (i, i+N1, i+N1+N2), the first set of steps comprising:

[0211] an offset collection step E1, implemented by the offset collection module 2, involving collecting a current lateral offset η.sub.LOC.sub.μA in the first embodiment (respectively a current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad in the second embodiment) measured by the offset measurement module 21 based on the reference guidance signal transmitted by the transmitting station 22 of the landing assistance device;

[0212] a position vector collection step E2, implemented by the position vector collection module 3, involving collecting a current position vector of the aircraft AC measured by the position vector measurement module 31 based on signals transmitted by a geopositioning system 32, the current position vector comprising at least one current longitudinal position x.sub.RWY and a current lateral position y.sub.RWY in the first embodiment (respectively a current vertical position z.sub.RWY in the second embodiment);

[0213] a second set of computation steps comprising:

[0214] a first computation step E3, implemented by the computation module 4 in a first iteration i, involving computing an angular alignment bias b.sub.LOC.sub.align (respectively b.sub.GPA) of the reference guidance axis A in relation to the approach axis B1 (respectively B2) of the runway RWY based on the current offset η.sub.LOC.sub.μA (respectively η.sub.GS.sub.ddm, η.sub.GS.sub.rad), the current longitudinal position x.sub.RWY and the current lateral position y.sub.RWY (respectively the current vertical position z.sub.RWY);

[0215] a computation step E4, implemented by the computation module 5 in a second iteration i+N1 subsequent to the first iteration i, involving computing an observable position bias b.sub.y.sub.obs (respectively b.sub.z.sub.obs) based on the angular alignment bias b.sub.LOC.sub.align (respectively b.sub.GPA) computed in the first iteration i, the current offset η.sub.LOC.sub.μA (respectively η.sub.GS.sub.ddm, η.sub.GS.sub.rad), the current longitudinal position x.sub.RWY and the current lateral position y.sub.RWY (respectively the current vertical position z.sub.RWY);

[0216] a computation step E6, implemented by the computation module 7 in a third iteration i+N1+N2 subsequent to the second iteration i+N1, involving computing the lateral position error b.sub.y (respectively the vertical position error b.sub.z) based on the observable lateral position bias b.sub.y.sub.obs (respectively the observable vertical position bias b.sub.z.sub.obs) computed in the second iteration i+N1, the current longitudinal position x.sub.RWY and the longitudinal position measured in the first iteration i;

[0217] a first transmission step E7, implemented by the transmission module 8, involving transmitting the lateral position error b.sub.y (respectively the vertical position error b.sub.z) to a user device 81.

[0218] Advantageously, the position vector of the aircraft is expressed in a coordinate system associated with the runway RWY.

[0219] The method can further comprise a filtering step Ef, implemented by the filtering module 6A, involving filtering the angular alignment position bias b.sub.LOC.sub.align (respectively the angular alignment bias b.sub.GPA) computed in the computation step E3.

[0220] The method can further comprise a filtering step E5, implemented by the filtering module 6B, involving filtering the observable lateral position bias b.sub.y.sub.obs(respectively the observable vertical position bias b.sub.z.sub.obs) computed in the computation step E4.

[0221] The method can also comprise:

[0222] a computation step E8, implemented by the computation module 9, involving computing a corrected angular alignment bias

[00024]$b_{{\mathrm{LOC}}_{{\mathrm{align}}_{\mathrm{TRUE}}}}$

(respectively b.sub.GPA.sub.TRUE) based on the current lateral offset η.sub.LOC.sub.μA (respectively based on the current vertical offset η.sub.GS.sub.ddm, η.sub.GS.sub.rad), the current longitudinal position x.sub.RWY, the current lateral position y.sub.RWY (respectively the current vertical position z.sub.RWY) and the lateral position error b.sub.y (respectively the vertical position error b.sub.z) computed in computation step E6;

[0223] a transmission step E9, implemented by a transmission module 10, involving transmitting the corrected angular alignment bias

[00025]$b_{{\mathrm{LOC}}_{{\mathrm{align}}_{\mathrm{TRUE}}}}$

(respectively b.sub.GPA.sub.TRUE) to a user device 101.

[0224] The systems and devices described herein may include a controller or a computing device comprising a processing unit and a memory which has stored therein computer-executable instructions for implementing the processes described herein. The processing unit may comprise any suitable devices configured to cause a series of steps to be performed so as to implement the method such that instructions, when executed by the computing device or other programmable apparatus, may cause the functions/acts/steps specified in the methods described herein to be executed. The processing unit may comprise, for example, any type of general-purpose microprocessor or microcontroller, a digital signal processing (DSP) processor, a central processing unit (CPU), an integrated circuit, a field programmable gate array (FPGA), a reconfigurable processor, other suitably programmed or programmable logic circuits, or any combination thereof.

[0225] The memory may be any suitable known or other machine-readable storage medium. The memory may comprise non-transitory computer readable storage medium such as, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. The memory may include a suitable combination of any type of computer memory that is located either internally or externally to the device such as, for example, random-access memory (RAM), read-only memory (ROM), compact disc read-only memory (CDROM), electro-optical memory, magneto-optical memory, erasable programmable read-only memory (EPROM), and electrically-erasable programmable read-only memory (EEPROM), Ferroelectric RAM (FRAM) or the like. The memory may comprise any storage means (e.g., devices) suitable for retrievably storing the computer-executable instructions executable by processing unit.

[0226] The methods and systems described herein may be implemented in a high-level procedural or object-oriented programming or scripting language, or a combination thereof, to communicate with or assist in the operation of the controller or computing device. Alternatively, the methods and systems described herein may be implemented in assembly or machine language. The language may be a compiled or interpreted language. Program code for implementing the methods and systems for detecting skew in a wing slat of an aircraft described herein may be stored on the storage media or the device, for example a ROM, a magnetic disk, an optical disc, a flash drive, or any other suitable storage media or device. The program code may be readable by a general or special-purpose programmable computer for configuring and operating the computer when the storage media or device is read by the computer to perform the procedures described herein.

[0227] Computer-executable instructions may be in many forms, including program modules, executed by one or more computers or other devices. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform particular tasks or implement particular abstract data types. Typically, the functionality of the program modules may be combined or distributed as desired in various embodiments.

[0228] While at least one exemplary embodiment of the present invention(s) is disclosed herein, it should be understood that modifications, substitutions and alternatives may be apparent to one of ordinary skill in the art and can be made without departing from the scope of this disclosure. This disclosure is intended to cover any adaptations or variations of the exemplary embodiment(s). In addition, in this disclosure, the terms “comprise” or “comprising” do not exclude other elements or steps, the terms “a” or “one” do not exclude a plural number, and the term “or” means either or both. Furthermore, characteristics or steps which have been described may also be used in combination with other characteristics or steps and in any order unless the disclosure or context suggests otherwise. This disclosure hereby incorporates by reference the complete disclosure of any patent or application from which it claims benefit or priority.