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
22 package org.altusmetrum.AltosLib;
24 public class AltosState {
25 public AltosRecord data;
29 public long report_time;
32 public double time_change;
36 public boolean landed;
37 public boolean ascent; /* going up? */
38 public boolean boost; /* under power */
40 public double ground_altitude;
41 public double altitude;
43 public double acceleration;
44 public double battery;
45 public double temperature;
46 public double main_sense;
47 public double drogue_sense;
48 public double accel_speed;
49 public double baro_speed;
51 public double max_height;
52 public double max_acceleration;
53 public double max_accel_speed;
54 public double max_baro_speed;
61 public static final int MIN_PAD_SAMPLES = 10;
65 public int gps_waiting;
66 public boolean gps_ready;
68 public AltosGreatCircle from_pad;
69 public double elevation; /* from pad */
70 public double range; /* total distance */
72 public double gps_height;
74 public double pad_lat, pad_lon, pad_alt;
76 public int speak_tick;
77 public double speak_altitude;
79 public double speed() {
86 public double max_speed() {
87 if (max_accel_speed != 0)
88 return max_accel_speed;
89 return max_baro_speed;
92 public void init (AltosRecord cur, AltosState prev_state) {
95 ground_altitude = data.ground_altitude();
97 altitude = data.altitude();
99 if (data.kalman_height != AltosRecord.MISSING)
100 height = data.kalman_height;
102 if (prev_state != null)
103 height = (prev_state.height * 15 + altitude - ground_altitude) / 16.0;
106 report_time = System.currentTimeMillis();
108 if (data.kalman_acceleration != AltosRecord.MISSING)
109 acceleration = data.kalman_acceleration;
111 acceleration = data.acceleration();
112 temperature = data.temperature();
113 drogue_sense = data.drogue_voltage();
114 main_sense = data.main_voltage();
115 battery = data.battery_voltage();
119 if (prev_state != null) {
121 /* Preserve any existing gps data */
122 npad = prev_state.npad;
123 ngps = prev_state.ngps;
124 gps = prev_state.gps;
125 pad_lat = prev_state.pad_lat;
126 pad_lon = prev_state.pad_lon;
127 pad_alt = prev_state.pad_alt;
128 max_height = prev_state.max_height;
129 max_acceleration = prev_state.max_acceleration;
130 max_accel_speed = prev_state.max_accel_speed;
131 max_baro_speed = prev_state.max_baro_speed;
132 imu = prev_state.imu;
133 mag = prev_state.mag;
135 /* make sure the clock is monotonic */
136 while (tick < prev_state.tick)
139 time_change = (tick - prev_state.tick) / 100.0;
141 if (data.kalman_speed != AltosRecord.MISSING) {
142 baro_speed = accel_speed = data.kalman_speed;
144 /* compute barometric speed */
146 double height_change = height - prev_state.height;
148 baro_speed = (prev_state.baro_speed * 3 + (height_change / time_change)) / 4.0;
150 baro_speed = prev_state.baro_speed;
152 if (acceleration == AltosRecord.MISSING) {
153 /* Fill in mising acceleration value */
154 accel_speed = baro_speed;
156 acceleration = (accel_speed - prev_state.accel_speed) / time_change;
158 acceleration = prev_state.acceleration;
160 /* compute accelerometer speed */
161 accel_speed = prev_state.accel_speed + acceleration * time_change;
172 if (acceleration == AltosRecord.MISSING)
178 if (cur.new_gps && (state == AltosLib.ao_flight_pad || state == AltosLib.ao_flight_idle)) {
180 /* Track consecutive 'good' gps reports, waiting for 10 of them */
181 if (data.gps != null && data.gps.locked && data.gps.nsat >= 4)
186 /* Average GPS data while on the pad */
187 if (data.gps != null && data.gps.locked && data.gps.nsat >= 4) {
189 /* filter pad position */
190 pad_lat = (pad_lat * 31.0 + data.gps.lat) / 32.0;
191 pad_lon = (pad_lon * 31.0 + data.gps.lon) / 32.0;
192 pad_alt = (pad_alt * 31.0 + data.gps.alt) / 32.0;
194 pad_lat = data.gps.lat;
195 pad_lon = data.gps.lon;
196 pad_alt = data.gps.alt;
201 pad_alt = ground_altitude;
203 data.new_gps = false;
205 gps_waiting = MIN_PAD_SAMPLES - npad;
209 gps_ready = gps_waiting == 0;
211 ascent = (AltosLib.ao_flight_boost <= state &&
212 state <= AltosLib.ao_flight_coast);
213 boost = (AltosLib.ao_flight_boost == state);
215 /* Only look at accelerometer data under boost */
216 if (boost && acceleration > max_acceleration && acceleration != AltosRecord.MISSING)
217 max_acceleration = acceleration;
218 if (boost && accel_speed > max_accel_speed && accel_speed != AltosRecord.MISSING)
219 max_accel_speed = accel_speed;
220 if (boost && baro_speed > max_baro_speed && baro_speed != AltosRecord.MISSING)
221 max_baro_speed = baro_speed;
223 if (height > max_height && height != AltosRecord.MISSING)
225 if (data.gps != null) {
226 if (gps == null || !gps.locked || data.gps.locked)
228 if (ngps > 0 && gps.locked) {
229 from_pad = new AltosGreatCircle(pad_lat, pad_lon, gps.lat, gps.lon);
235 gps_height = gps.alt - pad_alt;
236 if (from_pad != null) {
237 elevation = Math.atan2(height, from_pad.distance) * 180 / Math.PI;
238 range = Math.sqrt(height * height + from_pad.distance * from_pad.distance);
245 public AltosState(AltosRecord cur) {
249 public AltosState (AltosRecord cur, AltosState prev) {