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_1;
24 public class AltosState implements Cloneable {
25 public AltosRecord record;
27 public static final int set_position = 1;
28 public static final int set_gps = 2;
29 public static final int set_data = 4;
33 static final double ascent_filter_len = 0.1;
34 static final double descent_filter_len = 2.0;
38 public long received_time;
41 public double prev_time;
42 public double time_change;
44 private int prev_tick;
45 public int boost_tick;
49 private double prev_value;
50 private double max_value;
51 private double set_time;
52 private double prev_set_time;
53 private double max_rate = 1000.0;
55 void set(double new_value, double time) {
56 if (new_value != AltosRecord.MISSING) {
58 if (max_value == AltosRecord.MISSING || value > max_value) {
65 void set_filtered(double new_value, double time) {
66 if (prev_value != AltosRecord.MISSING)
67 new_value = (prev_value * 15.0 + new_value) / 16.0;
84 if (value != AltosRecord.MISSING && prev_value != AltosRecord.MISSING)
85 return value - prev_value;
86 return AltosRecord.MISSING;
91 double t = set_time - prev_set_time;
93 if (c != AltosRecord.MISSING && t != 0)
95 return AltosRecord.MISSING;
99 if (value == AltosRecord.MISSING)
100 return AltosRecord.MISSING;
101 if (prev_value == AltosRecord.MISSING)
102 return AltosRecord.MISSING;
104 return (value + prev_value) / 2 * (set_time - prev_set_time);
111 void set_derivative(AltosValue in) {
112 double n = in.rate();
114 if (n == AltosRecord.MISSING)
117 double p = prev_value;
118 double pt = prev_set_time;
120 if (p == AltosRecord.MISSING) {
122 pt = in.time() - 0.01;
125 /* Clip changes to reduce noise */
126 double ddt = in.time() - pt;
127 double ddv = (n - p) / ddt;
130 if (Math.abs(ddv) > 1000) {
140 filter_len = ascent_filter_len;
142 filter_len = descent_filter_len;
144 double f = 1/Math.exp(ddt/ filter_len);
145 n = p * f + n * (1-f);
150 void set_integral(AltosValue in) {
151 double change = in.integrate();
153 if (change != AltosRecord.MISSING) {
154 double prev = prev_value;
155 if (prev == AltosRecord.MISSING)
157 set(prev + change, in.time());
161 void copy(AltosValue old) {
163 set_time = old.set_time;
164 prev_value = old.value;
165 prev_set_time = old.set_time;
166 max_value = old.max_value;
170 value = AltosRecord.MISSING;
171 prev_value = AltosRecord.MISSING;
172 max_value = AltosRecord.MISSING;
181 double v = measured.value();
182 if (v != AltosRecord.MISSING)
184 return computed.value();
187 boolean is_measured() {
188 return measured.value() != AltosRecord.MISSING;
192 double m = measured.max();
194 if (m != AltosRecord.MISSING)
196 return computed.max();
199 double prev_value() {
200 if (measured.value != AltosRecord.MISSING && measured.prev_value != AltosRecord.MISSING)
201 return measured.prev_value;
202 return computed.prev_value;
205 AltosValue altos_value() {
206 if (measured.value() != AltosRecord.MISSING)
212 double c = measured.change();
213 if (c == AltosRecord.MISSING)
214 c = computed.change();
219 double r = measured.rate();
220 if (r == AltosRecord.MISSING)
225 void set_measured(double new_value, double time) {
226 measured.set(new_value, time);
229 void set_computed(double new_value, double time) {
230 computed.set(new_value, time);
233 void set_derivative(AltosValue in) {
234 computed.set_derivative(in);
237 void set_derivative(AltosCValue in) {
238 set_derivative(in.altos_value());
241 void set_integral(AltosValue in) {
242 computed.set_integral(in);
245 void set_integral(AltosCValue in) {
246 set_integral(in.altos_value());
249 void copy(AltosCValue old) {
250 measured.copy(old.measured);
251 computed.copy(old.computed);
255 measured = new AltosValue();
256 computed = new AltosValue();
263 public boolean landed;
264 public boolean ascent; /* going up? */
265 public boolean boost; /* under power */
268 public int device_type;
269 public int config_major;
270 public int config_minor;
271 public int apogee_delay;
272 public int main_deploy;
273 public int flight_log_max;
275 private double pressure_to_altitude(double p) {
276 if (p == AltosRecord.MISSING)
277 return AltosRecord.MISSING;
278 return AltosConvert.pressure_to_altitude(p);
281 private AltosCValue ground_altitude;
283 public double ground_altitude() {
284 return ground_altitude.value();
287 public void set_ground_altitude(double a) {
288 ground_altitude.set_measured(a, time);
291 class AltosGroundPressure extends AltosCValue {
292 void set_filtered(double p, double time) {
293 computed.set_filtered(p, time);
294 ground_altitude.set_computed(pressure_to_altitude(computed.value()), time);
297 void set_measured(double p, double time) {
298 super.set_measured(p, time);
299 ground_altitude.set_computed(pressure_to_altitude(p), time);
303 private AltosGroundPressure ground_pressure;
305 public double ground_pressure() {
306 return ground_pressure.value();
309 public void set_ground_pressure (double pressure) {
310 ground_pressure.set_measured(pressure, time);
313 class AltosAltitude extends AltosCValue {
315 private void set_speed(AltosValue v) {
316 if (!acceleration.is_measured() || !ascent)
317 speed.set_derivative(this);
320 void set_computed(double a, double time) {
321 super.set_computed(a,time);
326 void set_measured(double a, double time) {
327 super.set_measured(a,time);
333 private AltosAltitude altitude;
335 public double altitude() {
336 double a = altitude.value();
337 if (a != AltosRecord.MISSING)
341 return AltosRecord.MISSING;
344 public double max_altitude() {
345 double a = altitude.max();
346 if (a != AltosRecord.MISSING)
348 return AltosRecord.MISSING;
351 public void set_altitude(double new_altitude) {
352 altitude.set_measured(new_altitude, time);
355 class AltosPressure extends AltosValue {
356 void set(double p, double time) {
358 if (state == AltosLib.ao_flight_pad)
359 ground_pressure.set_filtered(p, time);
360 double a = pressure_to_altitude(p);
361 altitude.set_computed(a, time);
365 private AltosPressure pressure;
367 public double pressure() {
368 return pressure.value();
371 public void set_pressure(double p) {
372 pressure.set(p, time);
375 public double height() {
376 double k = kalman_height.value();
377 if (k != AltosRecord.MISSING)
380 double a = altitude();
381 double g = ground_altitude();
382 if (a != AltosRecord.MISSING && g != AltosRecord.MISSING)
384 return AltosRecord.MISSING;
387 public double max_height() {
388 double k = kalman_height.max();
389 if (k != AltosRecord.MISSING)
392 double a = altitude.max();
393 double g = ground_altitude();
394 if (a != AltosRecord.MISSING && g != AltosRecord.MISSING)
396 return AltosRecord.MISSING;
399 class AltosSpeed extends AltosCValue {
402 acceleration.set_derivative(this);
405 void set_derivative(AltosCValue in) {
406 super.set_derivative(in);
410 void set_computed(double new_value, double time) {
411 super.set_computed(new_value, time);
415 void set_measured(double new_value, double time) {
416 super.set_measured(new_value, time);
421 private AltosSpeed speed;
423 public double speed() {
424 return speed.value();
427 public double max_speed() {
431 class AltosAccel extends AltosCValue {
432 void set_measured(double a, double time) {
433 super.set_measured(a, time);
435 speed.set_integral(this.measured);
439 AltosAccel acceleration;
441 public double acceleration() {
442 return acceleration.value();
445 public double max_acceleration() {
446 return acceleration.max();
449 public AltosValue kalman_height, kalman_speed, kalman_acceleration;
451 public void set_kalman(double height, double speed, double acceleration) {
452 kalman_height.set(height, time);
453 kalman_speed.set(speed, time);
454 kalman_acceleration.set(acceleration, time);
457 public double battery_voltage;
458 public double pyro_voltage;
459 public double temperature;
460 public double apogee_voltage;
461 public double main_voltage;
463 public double ignitor_voltage[];
466 public AltosGPS temp_gps;
467 public int temp_gps_sat_tick;
468 public boolean gps_pending;
469 public int gps_sequence;
474 public static final int MIN_PAD_SAMPLES = 10;
477 public int gps_waiting;
478 public boolean gps_ready;
482 public AltosGreatCircle from_pad;
483 public double elevation; /* from pad */
484 public double range; /* total distance */
486 public double gps_height;
488 public double pad_lat, pad_lon, pad_alt;
490 public int speak_tick;
491 public double speak_altitude;
493 public String callsign;
494 public String firmware_version;
496 public double accel_plus_g;
497 public double accel_minus_g;
499 public double ground_accel;
500 public double ground_accel_avg;
502 public int log_format;
504 public AltosMs5607 baro;
506 public AltosRecordCompanion companion;
508 public void set_npad(int npad) {
510 gps_waiting = MIN_PAD_SAMPLES - npad;
511 if (this.gps_waiting < 0)
513 gps_ready = gps_waiting == 0;
521 received_time = System.currentTimeMillis();
522 time = AltosRecord.MISSING;
523 time_change = AltosRecord.MISSING;
524 prev_time = AltosRecord.MISSING;
525 tick = AltosRecord.MISSING;
526 prev_tick = AltosRecord.MISSING;
527 boost_tick = AltosRecord.MISSING;
528 state = AltosLib.ao_flight_invalid;
529 flight = AltosRecord.MISSING;
532 rssi = AltosRecord.MISSING;
534 device_type = AltosRecord.MISSING;
535 config_major = AltosRecord.MISSING;
536 config_minor = AltosRecord.MISSING;
537 apogee_delay = AltosRecord.MISSING;
538 main_deploy = AltosRecord.MISSING;
539 flight_log_max = AltosRecord.MISSING;
541 ground_altitude = new AltosCValue();
542 ground_pressure = new AltosGroundPressure();
543 altitude = new AltosAltitude();
544 pressure = new AltosPressure();
545 speed = new AltosSpeed();
546 acceleration = new AltosAccel();
548 temperature = AltosRecord.MISSING;
549 battery_voltage = AltosRecord.MISSING;
550 pyro_voltage = AltosRecord.MISSING;
551 apogee_voltage = AltosRecord.MISSING;
552 main_voltage = AltosRecord.MISSING;
553 ignitor_voltage = null;
555 kalman_height = new AltosValue();
556 kalman_speed = new AltosValue();
557 kalman_acceleration = new AltosValue();
561 temp_gps_sat_tick = 0;
572 elevation = AltosRecord.MISSING;
573 range = AltosRecord.MISSING;
574 gps_height = AltosRecord.MISSING;
576 pad_lat = AltosRecord.MISSING;
577 pad_lon = AltosRecord.MISSING;
578 pad_alt = AltosRecord.MISSING;
580 speak_tick = AltosRecord.MISSING;
581 speak_altitude = AltosRecord.MISSING;
585 accel_plus_g = AltosRecord.MISSING;
586 accel_minus_g = AltosRecord.MISSING;
587 accel = AltosRecord.MISSING;
589 ground_accel = AltosRecord.MISSING;
590 ground_accel_avg = AltosRecord.MISSING;
592 log_format = AltosRecord.MISSING;
593 serial = AltosRecord.MISSING;
599 void copy(AltosState old) {
608 received_time = old.received_time;
612 prev_tick = old.tick;
613 boost_tick = old.boost_tick;
622 device_type = old.device_type;
623 config_major = old.config_major;
624 config_minor = old.config_minor;
625 apogee_delay = old.apogee_delay;
626 main_deploy = old.main_deploy;
627 flight_log_max = old.flight_log_max;
631 ground_altitude.copy(old.ground_altitude);
632 altitude.copy(old.altitude);
633 pressure.copy(old.pressure);
634 speed.copy(old.speed);
635 acceleration.copy(old.acceleration);
637 battery_voltage = old.battery_voltage;
638 pyro_voltage = old.pyro_voltage;
639 temperature = old.temperature;
640 apogee_voltage = old.apogee_voltage;
641 main_voltage = old.main_voltage;
642 ignitor_voltage = old.ignitor_voltage;
644 kalman_height.copy(old.kalman_height);
645 kalman_speed.copy(old.kalman_speed);
646 kalman_acceleration.copy(old.kalman_acceleration);
649 gps = old.gps.clone();
652 if (old.temp_gps != null)
653 temp_gps = old.temp_gps.clone();
656 temp_gps_sat_tick = old.temp_gps_sat_tick;
657 gps_sequence = old.gps_sequence;
658 gps_pending = old.gps_pending;
661 imu = old.imu.clone();
666 mag = old.mag.clone();
671 gps_waiting = old.gps_waiting;
672 gps_ready = old.gps_ready;
675 if (old.from_pad != null)
676 from_pad = old.from_pad.clone();
680 elevation = old.elevation;
683 gps_height = old.gps_height;
684 pad_lat = old.pad_lat;
685 pad_lon = old.pad_lon;
686 pad_alt = old.pad_alt;
688 speak_tick = old.speak_tick;
689 speak_altitude = old.speak_altitude;
691 callsign = old.callsign;
693 accel_plus_g = old.accel_plus_g;
694 accel_minus_g = old.accel_minus_g;
696 ground_accel = old.ground_accel;
697 ground_accel_avg = old.ground_accel_avg;
699 log_format = old.log_format;
703 companion = old.companion;
717 if (gps.locked && gps.nsat >= 4) {
718 /* Track consecutive 'good' gps reports, waiting for 10 of them */
719 if (state == AltosLib.ao_flight_pad) {
721 if (pad_lat != AltosRecord.MISSING) {
722 pad_lat = (pad_lat * 31 + gps.lat) / 32;
723 pad_lon = (pad_lon * 31 + gps.lon) / 32;
724 pad_alt = (pad_alt * 31 + gps.alt) / 32;
727 if (pad_lat == AltosRecord.MISSING) {
733 if (gps.lat != 0 && gps.lon != 0 &&
734 pad_lat != AltosRecord.MISSING &&
735 pad_lon != AltosRecord.MISSING)
739 if (h == AltosRecord.MISSING)
741 from_pad = new AltosGreatCircle(pad_lat, pad_lon, 0, gps.lat, gps.lon, h);
742 elevation = from_pad.elevation;
743 range = from_pad.range;
744 gps_height = gps.alt - pad_alt;
748 public void set_tick(int new_tick) {
749 if (new_tick != AltosRecord.MISSING) {
750 if (prev_tick != AltosRecord.MISSING) {
751 while (new_tick < prev_tick - 32767) {
760 public void set_boost_tick(int boost_tick) {
761 if (boost_tick != AltosRecord.MISSING)
762 this.boost_tick = boost_tick;
765 public String state_name() {
766 return AltosLib.state_name(state);
769 public void set_state(int state) {
770 if (state != AltosLib.ao_flight_invalid) {
772 ascent = (AltosLib.ao_flight_boost <= state &&
773 state <= AltosLib.ao_flight_coast);
774 boost = (AltosLib.ao_flight_boost == state);
779 public void set_device_type(int device_type) {
780 this.device_type = device_type;
783 public void set_config(int major, int minor, int apogee_delay, int main_deploy, int flight_log_max) {
784 config_major = major;
785 config_minor = minor;
786 this.apogee_delay = apogee_delay;
787 this.main_deploy = main_deploy;
788 this.flight_log_max = flight_log_max;
791 public void set_callsign(String callsign) {
792 this.callsign = callsign;
795 public void set_firmware_version(String version) {
796 firmware_version = version;
799 public void set_flight(int flight) {
801 /* When the flight changes, reset the state */
802 if (flight != AltosRecord.MISSING) {
803 if (this.flight != AltosRecord.MISSING &&
804 this.flight != flight) {
807 this.flight = flight;
811 public void set_serial(int serial) {
812 /* When the serial changes, reset the state */
813 if (serial != AltosRecord.MISSING) {
814 if (this.serial != AltosRecord.MISSING &&
815 this.serial != serial) {
818 this.serial = serial;
823 if (rssi == AltosRecord.MISSING)
828 public void set_rssi(int rssi, int status) {
829 if (rssi != AltosRecord.MISSING) {
831 this.status = status;
835 public void set_received_time(long ms) {
839 public void set_gps(AltosGPS gps, int sequence) {
841 this.gps = gps.clone();
842 gps_sequence = sequence;
848 public void make_baro() {
850 baro = new AltosMs5607();
853 public void set_ms5607(int pres, int temp) {
855 baro.set(pres, temp);
857 set_pressure(baro.pa);
858 set_temperature(baro.cc / 100.0);
862 public void make_companion (int nchannels) {
863 if (companion == null)
864 companion = new AltosRecordCompanion(nchannels);
867 public void set_companion(AltosRecordCompanion companion) {
868 this.companion = companion;
871 void update_accel() {
872 double ground = ground_accel;
874 if (ground == AltosRecord.MISSING)
875 ground = ground_accel_avg;
876 if (accel == AltosRecord.MISSING)
878 if (ground == AltosRecord.MISSING)
880 if (accel_plus_g == AltosRecord.MISSING)
882 if (accel_minus_g == AltosRecord.MISSING)
885 double counts_per_g = (accel_minus_g - accel_plus_g) / 2.0;
886 double counts_per_mss = counts_per_g / 9.80665;
887 acceleration.set_computed((ground - accel) / counts_per_mss, time);
890 public void set_accel_g(double accel_plus_g, double accel_minus_g) {
891 if (accel_plus_g != AltosRecord.MISSING) {
892 this.accel_plus_g = accel_plus_g;
893 this.accel_minus_g = accel_minus_g;
898 public void set_ground_accel(double ground_accel) {
899 if (ground_accel != AltosRecord.MISSING) {
900 this.ground_accel = ground_accel;
905 public void set_accel(double accel) {
906 if (accel != AltosRecord.MISSING) {
908 if (state == AltosLib.ao_flight_pad) {
909 if (ground_accel_avg == AltosRecord.MISSING)
910 ground_accel_avg = accel;
912 ground_accel_avg = (ground_accel_avg * 7 + accel) / 8;
918 public void set_temperature(double temperature) {
919 if (temperature != AltosRecord.MISSING) {
920 this.temperature = temperature;
925 public void set_battery_voltage(double battery_voltage) {
926 if (battery_voltage != AltosRecord.MISSING) {
927 this.battery_voltage = battery_voltage;
932 public void set_pyro_voltage(double pyro_voltage) {
933 if (pyro_voltage != AltosRecord.MISSING) {
934 this.pyro_voltage = pyro_voltage;
939 public void set_apogee_voltage(double apogee_voltage) {
940 if (apogee_voltage != AltosRecord.MISSING) {
941 this.apogee_voltage = apogee_voltage;
946 public void set_main_voltage(double main_voltage) {
947 if (main_voltage != AltosRecord.MISSING) {
948 this.main_voltage = main_voltage;
953 public void set_ignitor_voltage(double[] voltage) {
954 this.ignitor_voltage = voltage;
957 public double time_since_boost() {
958 if (tick == AltosRecord.MISSING)
961 if (boost_tick != AltosRecord.MISSING) {
962 return (tick - boost_tick) / 100.0;
967 public boolean valid() {
968 return tick != AltosRecord.MISSING && serial != AltosRecord.MISSING;
971 public AltosGPS make_temp_gps(boolean sats) {
972 if (temp_gps == null) {
973 temp_gps = new AltosGPS(gps);
977 if (tick != temp_gps_sat_tick)
978 temp_gps.cc_gps_sat = null;
979 temp_gps_sat_tick = tick;
984 public void set_temp_gps() {
985 set_gps(temp_gps, gps_sequence + 1);
990 public AltosState clone() {
991 AltosState s = new AltosState();
997 public void init (AltosRecord cur, AltosState prev_state) {
999 System.out.printf ("init\n");
1001 cur = new AltosRecord();
1005 /* Discard previous state if it was for a different board */
1006 if (prev_state != null && prev_state.serial != cur.serial)
1011 set_ground_altitude(cur.ground_altitude());
1012 set_altitude(cur.altitude());
1014 set_kalman(cur.kalman_height, cur.kalman_speed, cur.kalman_acceleration);
1016 received_time = System.currentTimeMillis();
1018 set_temperature(cur.temperature());
1019 set_apogee_voltage(cur.drogue_voltage());
1020 set_main_voltage(cur.main_voltage());
1021 set_battery_voltage(cur.battery_voltage());
1023 set_pressure(cur.pressure());
1026 set_state(cur.state);
1028 set_accel_g (cur.accel_minus_g, cur.accel_plus_g);
1029 set_ground_accel(cur.ground_accel);
1030 set_accel (cur.accel);
1032 if (cur.gps_sequence != gps_sequence)
1033 set_gps(cur.gps, cur.gps_sequence);
1037 public AltosState(AltosRecord cur) {
1041 public AltosState (AltosRecord cur, AltosState prev) {
1046 public AltosState () {