diff options
Diffstat (limited to 'apps/systemlib')
-rw-r--r-- | apps/systemlib/airspeed.c | 79 | ||||
-rw-r--r-- | apps/systemlib/airspeed.h | 86 | ||||
-rw-r--r-- | apps/systemlib/conversions.c | 11 | ||||
-rw-r--r-- | apps/systemlib/conversions.h | 8 |
4 files changed, 183 insertions, 1 deletions
diff --git a/apps/systemlib/airspeed.c b/apps/systemlib/airspeed.c new file mode 100644 index 000000000..e213b66c2 --- /dev/null +++ b/apps/systemlib/airspeed.c @@ -0,0 +1,79 @@ +/**************************************************************************** + * + * Copyright (C) 2012-2013 PX4 Development Team. All rights reserved. + * Author: Lorenz Meier <lm@inf.ethz.ch> + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in + * the documentation and/or other materials provided with the + * distribution. + * 3. Neither the name PX4 nor the names of its contributors may be + * used to endorse or promote products derived from this software + * without specific prior written permission. + * + * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS + * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT + * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS + * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE + * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, + * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, + * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS + * OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED + * AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT + * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN + * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + * POSSIBILITY OF SUCH DAMAGE. + * + ****************************************************************************/ + +/** + * @file airspeed.c + * Airspeed estimation + * + * @author Lorenz Meier <lm@inf.ethz.ch> + * + */ + +#include "math.h" + + +float calc_indicated_airspeed(float pressure_front, float pressure_ambient, float temperature) +{ + return sqrtf((2.0f*(pressure_front - pressure_ambient)) / air_density_sea_level); +} + +/** + * Calculate true airspeed from indicated airspeed. + * + * Note that the true airspeed is NOT the groundspeed, because of the effects of wind + * + * @param speed current indicated airspeed + * @param pressure_ambient pressure at the side of the tube/airplane + * @param temperature air temperature in degrees celcius + * @return true airspeed in m/s + */ +float calc_true_airspeed_from_indicated(float speed, float pressure_ambient, float temperature) +{ + return speed * sqrtf(air_density_sea_level / get_air_density(pressure_ambient, temperature)); +} + +/** + * Directly calculate true airspeed + * + * Note that the true airspeed is NOT the groundspeed, because of the effects of wind + * + * @param pressure_front pressure inside the pitot/prandl tube + * @param pressure_ambient pressure at the side of the tube/airplane + * @param temperature air temperature in degrees celcius + * @return true airspeed in m/s + */ +float calc_true_airspeed(float pressure_front, float pressure_ambient, float temperature) +{ + return sqrtf((2.0f*(pressure_front - pressure_ambient)) / get_air_density(pressure_ambient, temperature)); +}
\ No newline at end of file diff --git a/apps/systemlib/airspeed.h b/apps/systemlib/airspeed.h new file mode 100644 index 000000000..62acfe2b0 --- /dev/null +++ b/apps/systemlib/airspeed.h @@ -0,0 +1,86 @@ +/**************************************************************************** + * + * Copyright (C) 2012-2013 PX4 Development Team. All rights reserved. + * Author: Lorenz Meier <lm@inf.ethz.ch> + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in + * the documentation and/or other materials provided with the + * distribution. + * 3. Neither the name PX4 nor the names of its contributors may be + * used to endorse or promote products derived from this software + * without specific prior written permission. + * + * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS + * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT + * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS + * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE + * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, + * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, + * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS + * OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED + * AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT + * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN + * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + * POSSIBILITY OF SUCH DAMAGE. + * + ****************************************************************************/ + +/** + * @file airspeed.c + * Airspeed estimation + * + * @author Lorenz Meier <lm@inf.ethz.ch> + * + */ + +#include "math.h" +#include "conversions.h" + + __BEGIN_DECLS + +/** + * Calculate indicated airspeed. + * + * Note that the indicated airspeed is not the true airspeed because it + * lacks the air density compensation. Use the calc_true_airspeed functions to get + * the true airspeed. + * + * @param pressure_front pressure inside the pitot/prandl tube + * @param pressure_ambient pressure at the side of the tube/airplane + * @param temperature air temperature in degrees celcius + * @return indicated airspeed in m/s + */ +__EXPORT float calc_indicated_airspeed(float pressure_front, float pressure_ambient, float temperature); + +/** + * Calculate true airspeed from indicated airspeed. + * + * Note that the true airspeed is NOT the groundspeed, because of the effects of wind + * + * @param speed current indicated airspeed + * @param pressure_ambient pressure at the side of the tube/airplane + * @param temperature air temperature in degrees celcius + * @return true airspeed in m/s + */ +__EXPORT float calc_true_airspeed_from_indicated(float speed, float pressure_ambient, float temperature); + +/** + * Directly calculate true airspeed + * + * Note that the true airspeed is NOT the groundspeed, because of the effects of wind + * + * @param pressure_front pressure inside the pitot/prandl tube + * @param pressure_ambient pressure at the side of the tube/airplane + * @param temperature air temperature in degrees celcius + * @return true airspeed in m/s + */ +__EXPORT float calc_true_airspeed(float pressure_front, float pressure_ambient, float temperature); + +__END_DECLS
\ No newline at end of file diff --git a/apps/systemlib/conversions.c b/apps/systemlib/conversions.c index 32df446d9..b08b9fc69 100644 --- a/apps/systemlib/conversions.c +++ b/apps/systemlib/conversions.c @@ -42,6 +42,10 @@ #include "conversions.h" +#define air_gas_constant 8.31432f +#define air_density_sea_level 1.225f +#define absolute_null_kelvin 273.15f + int16_t int16_t_from_bytes(uint8_t bytes[]) { @@ -143,4 +147,9 @@ void quat2rot(const float Q[4], float R[9]) R[2] = 2.0F * (Q[1] * Q[3] - Q[0] * Q[2]); R[5] = 2.0F * (Q[2] * Q[3] + Q[0] * Q[1]); R[8] = ((q0_2 + q3_2) - q1_2) - q2_2; -}
\ No newline at end of file +} + +float get_air_density(float static_pressure, float temperature_celsius) +{ + return static_pressure/(air_gas_constant * (temperature_celsius + absolute_null_kelvin)); +} diff --git a/apps/systemlib/conversions.h b/apps/systemlib/conversions.h index 3a24ca9b7..4db6de772 100644 --- a/apps/systemlib/conversions.h +++ b/apps/systemlib/conversions.h @@ -78,6 +78,14 @@ __EXPORT void rot2quat(const float R[9], float Q[4]); */ __EXPORT void quat2rot(const float Q[4], float R[9]); +/** + * Calculates air density. + * + * @param static_pressure ambient pressure in millibar + * @param temperature_celcius air / ambient temperature in celcius + */ +__EXPORT float get_air_density(float static_pressure, float temperature_celsius); + __END_DECLS #endif /* CONVERSIONS_H_ */ |