2 * Copyright © 2010 Keith Packard <keithp@keithp.com>
4 * This program is free software; you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation; version 2 of the License.
8 * This program is distributed in the hope that it will be useful, but
9 * WITHOUT ANY WARRANTY; without even the implied warranty of
10 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 * General Public License for more details.
13 * You should have received a copy of the GNU General Public License along
14 * with this program; if not, write to the Free Software Foundation, Inc.,
15 * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
19 * Track flight state from telemetry or eeprom data stream
24 import altosui.AltosRecord;
25 import altosui.AltosGPS;
27 public class AltosState {
39 boolean ascent; /* going up? */
41 double ground_altitude;
52 double max_acceleration;
61 static final int MIN_PAD_SAMPLES = 10;
68 AltosGreatCircle from_pad;
69 double elevation; /* from pad */
70 double range; /* total distance */
75 double speak_altitude;
78 void init (AltosRecord cur, AltosState prev_state) {
84 ground_altitude = data.ground_altitude();
85 height = data.filtered_altitude() - ground_altitude;
87 report_time = System.currentTimeMillis();
89 acceleration = data.acceleration();
90 speed = data.accel_speed();
91 temperature = data.temperature();
92 drogue_sense = data.drogue_voltage();
93 main_sense = data.main_voltage();
94 battery = data.battery_voltage();
98 if (prev_state != null) {
100 /* Preserve any existing gps data */
101 npad = prev_state.npad;
102 ngps = prev_state.ngps;
103 gps = prev_state.gps;
104 pad_lat = prev_state.pad_lat;
105 pad_lon = prev_state.pad_lon;
106 pad_alt = prev_state.pad_alt;
107 max_height = prev_state.max_height;
108 max_acceleration = prev_state.max_acceleration;
109 max_speed = prev_state.max_speed;
111 /* make sure the clock is monotonic */
112 while (tick < prev_state.tick)
115 time_change = (tick - prev_state.tick) / 100.0;
117 /* compute barometric speed */
119 double height_change = height - prev_state.height;
121 baro_speed = (prev_state.baro_speed * 3 + (height_change / time_change)) / 4.0;
123 baro_speed = prev_state.baro_speed;
132 if (state == Altos.ao_flight_pad) {
134 /* Track consecutive 'good' gps reports, waiting for 10 of them */
135 if (data.gps != null && data.gps.locked && data.gps.nsat >= 4)
140 /* Average GPS data while on the pad */
141 if (data.gps != null && data.gps.locked && data.gps.nsat >= 4) {
143 /* filter pad position */
144 pad_lat = (pad_lat * 31.0 + data.gps.lat) / 32.0;
145 pad_lon = (pad_lon * 31.0 + data.gps.lon) / 32.0;
146 pad_alt = (pad_alt * 31.0 + data.gps.alt) / 32.0;
148 pad_lat = data.gps.lat;
149 pad_lon = data.gps.lon;
150 pad_alt = data.gps.alt;
156 gps_waiting = MIN_PAD_SAMPLES - npad;
160 gps_ready = gps_waiting == 0;
162 ascent = (Altos.ao_flight_boost <= state &&
163 state <= Altos.ao_flight_coast);
165 /* Only look at accelerometer data on the way up */
166 if (ascent && acceleration > max_acceleration)
167 max_acceleration = acceleration;
168 if (ascent && speed > max_speed)
171 if (height > max_height)
173 if (data.gps != null) {
174 if (gps == null || !gps.locked || data.gps.locked)
176 if (ngps > 0 && gps.locked) {
177 from_pad = new AltosGreatCircle(pad_lat, pad_lon, gps.lat, gps.lon);
183 gps_height = gps.alt - pad_alt;
184 if (from_pad != null) {
185 elevation = Math.atan2(height, from_pad.distance) * 180 / Math.PI;
186 range = Math.sqrt(height * height + from_pad.distance * from_pad.distance);
193 public AltosState(AltosRecord cur) {
197 public AltosState (AltosRecord cur, AltosState prev) {