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) {
78 if (value != AltosRecord.MISSING && prev_value != AltosRecord.MISSING)
79 return value - prev_value;
80 return AltosRecord.MISSING;
85 double t = set_time - prev_set_time;
87 if (c != AltosRecord.MISSING && t != 0)
89 return AltosRecord.MISSING;
93 if (value == AltosRecord.MISSING)
94 return AltosRecord.MISSING;
95 if (prev_value == AltosRecord.MISSING)
96 return AltosRecord.MISSING;
98 return (value + prev_value) / 2 * (set_time - prev_set_time);
105 void set_derivative(AltosValue in) {
106 double n = in.rate();
108 if (n == AltosRecord.MISSING)
111 double p = prev_value;
112 double pt = prev_set_time;
114 if (p == AltosRecord.MISSING) {
116 pt = in.time() - 0.01;
119 /* Clip changes to reduce noise */
120 double ddt = in.time() - pt;
121 double ddv = (n - p) / ddt;
124 if (Math.abs(ddv) > 1000) {
134 filter_len = ascent_filter_len;
136 filter_len = descent_filter_len;
138 double f = 1/Math.exp(ddt/ filter_len);
139 n = p * f + n * (1-f);
144 void set_integral(AltosValue in) {
145 double change = in.integrate();
147 if (change != AltosRecord.MISSING)
148 set(prev_value + change, in.time());
151 void copy(AltosValue old) {
153 set_time = old.set_time;
154 prev_value = old.value;
155 prev_set_time = old.set_time;
156 max_value = old.max_value;
160 value = AltosRecord.MISSING;
161 prev_value = AltosRecord.MISSING;
162 max_value = AltosRecord.MISSING;
171 double v = measured.value();
172 if (v != AltosRecord.MISSING)
174 return computed.value();
177 boolean is_measured() {
178 return measured.value() != AltosRecord.MISSING;
182 double m = measured.max();
184 if (m != AltosRecord.MISSING)
186 return computed.max();
189 double prev_value() {
190 if (measured.value != AltosRecord.MISSING && measured.prev_value != AltosRecord.MISSING)
191 return measured.prev_value;
192 return computed.prev_value;
195 AltosValue altos_value() {
196 if (measured.value() != AltosRecord.MISSING)
202 double c = measured.change();
203 if (c == AltosRecord.MISSING)
204 c = computed.change();
209 double r = measured.rate();
210 if (r == AltosRecord.MISSING)
215 void set_measured(double new_value, double time) {
216 measured.set(new_value, time);
219 void set_computed(double new_value, double time) {
220 computed.set(new_value, time);
223 void set_derivative(AltosValue in) {
224 computed.set_derivative(in);
227 void set_derivative(AltosCValue in) {
228 set_derivative(in.altos_value());
231 void set_integral(AltosValue in) {
232 computed.set_integral(in);
235 void set_integral(AltosCValue in) {
236 set_integral(in.altos_value());
239 void copy(AltosCValue old) {
240 measured.copy(old.measured);
241 computed.copy(old.computed);
245 measured = new AltosValue();
246 computed = new AltosValue();
253 public boolean landed;
254 public boolean ascent; /* going up? */
255 public boolean boost; /* under power */
258 public int device_type;
259 public int config_major;
260 public int config_minor;
261 public int apogee_delay;
262 public int main_deploy;
263 public int flight_log_max;
265 private double pressure_to_altitude(double p) {
266 if (p == AltosRecord.MISSING)
267 return AltosRecord.MISSING;
268 return AltosConvert.pressure_to_altitude(p);
271 private AltosCValue ground_altitude;
273 public double ground_altitude() {
274 return ground_altitude.value();
277 public void set_ground_altitude(double a) {
278 ground_altitude.set_measured(a, time);
281 class AltosGroundPressure extends AltosValue {
282 void set(double p, double time) {
284 ground_altitude.set_computed(pressure_to_altitude(p), time);
288 private AltosGroundPressure ground_pressure;
290 public double ground_pressure() {
291 return ground_pressure.value();
294 public void set_ground_pressure (double pressure) {
295 ground_pressure.set(pressure, time);
298 class AltosAltitude extends AltosCValue {
300 private void set_speed(AltosValue v) {
301 if (!acceleration.is_measured() || !ascent)
302 speed.set_derivative(this);
305 void set_computed(double a, double time) {
306 super.set_computed(a,time);
311 void set_measured(double a, double time) {
312 super.set_measured(a,time);
318 private AltosAltitude altitude;
320 public double altitude() {
321 double a = altitude.value();
322 if (a != AltosRecord.MISSING)
326 return AltosRecord.MISSING;
329 public double max_altitude() {
330 double a = altitude.max();
331 if (a != AltosRecord.MISSING)
333 return AltosRecord.MISSING;
336 public void set_altitude(double new_altitude) {
337 altitude.set_measured(new_altitude, time);
340 class AltosPressure extends AltosValue {
341 void set(double p, double time) {
343 altitude.set_computed(pressure_to_altitude(p), time);
347 private AltosPressure pressure;
349 public double pressure() {
350 return pressure.value();
353 public void set_pressure(double p) {
354 pressure.set(p, time);
357 public double height() {
358 double k = kalman_height.value();
359 if (k != AltosRecord.MISSING)
362 double a = altitude();
363 double g = ground_altitude();
364 if (a != AltosRecord.MISSING && g != AltosRecord.MISSING)
366 return AltosRecord.MISSING;
369 public double max_height() {
370 double k = kalman_height.max();
371 if (k != AltosRecord.MISSING)
374 double a = altitude.max();
375 double g = ground_altitude();
376 if (a != AltosRecord.MISSING && g != AltosRecord.MISSING)
378 return AltosRecord.MISSING;
381 class AltosSpeed extends AltosCValue {
384 acceleration.set_derivative(this);
387 void set_derivative(AltosCValue in) {
388 super.set_derivative(in);
392 void set_computed(double new_value, double time) {
393 super.set_computed(new_value, time);
397 void set_measured(double new_value, double time) {
398 super.set_measured(new_value, time);
403 private AltosSpeed speed;
405 public double speed() {
406 return speed.value();
409 public double max_speed() {
413 class AltosAccel extends AltosCValue {
414 void set_measured(double a, double time) {
415 super.set_measured(a, time);
417 speed.set_integral(this.measured);
421 AltosAccel acceleration;
423 public double acceleration() {
424 return acceleration.value();
427 public double max_acceleration() {
428 return acceleration.max();
431 public AltosValue kalman_height, kalman_speed, kalman_acceleration;
433 public void set_kalman(double height, double speed, double acceleration) {
434 kalman_height.set(height, time);
435 kalman_speed.set(speed, time);
436 kalman_acceleration.set(acceleration, time);
439 public double battery_voltage;
440 public double pyro_voltage;
441 public double temperature;
442 public double apogee_voltage;
443 public double main_voltage;
445 public double ignitor_voltage[];
448 public AltosGPS temp_gps;
449 public boolean temp_gps_clear_sats_pending;
450 public boolean gps_pending;
451 public int gps_sequence;
456 public static final int MIN_PAD_SAMPLES = 10;
459 public int gps_waiting;
460 public boolean gps_ready;
464 public AltosGreatCircle from_pad;
465 public double elevation; /* from pad */
466 public double range; /* total distance */
468 public double gps_height;
470 public double pad_lat, pad_lon, pad_alt;
472 public int speak_tick;
473 public double speak_altitude;
475 public String callsign;
476 public String firmware_version;
478 public double accel_plus_g;
479 public double accel_minus_g;
481 public double ground_accel;
482 public double ground_accel_avg;
484 public int log_format;
486 public AltosMs5607 baro;
488 public AltosRecordCompanion companion;
490 public void set_npad(int npad) {
492 gps_waiting = MIN_PAD_SAMPLES - npad;
493 if (this.gps_waiting < 0)
495 gps_ready = gps_waiting == 0;
503 received_time = System.currentTimeMillis();
504 time = AltosRecord.MISSING;
505 time_change = AltosRecord.MISSING;
506 prev_time = AltosRecord.MISSING;
507 tick = AltosRecord.MISSING;
508 prev_tick = AltosRecord.MISSING;
509 boost_tick = AltosRecord.MISSING;
510 state = AltosLib.ao_flight_invalid;
511 flight = AltosRecord.MISSING;
514 rssi = AltosRecord.MISSING;
516 device_type = AltosRecord.MISSING;
517 config_major = AltosRecord.MISSING;
518 config_minor = AltosRecord.MISSING;
519 apogee_delay = AltosRecord.MISSING;
520 main_deploy = AltosRecord.MISSING;
521 flight_log_max = AltosRecord.MISSING;
523 ground_altitude = new AltosCValue();
524 ground_pressure = new AltosGroundPressure();
525 altitude = new AltosAltitude();
526 pressure = new AltosPressure();
527 speed = new AltosSpeed();
528 acceleration = new AltosAccel();
530 temperature = AltosRecord.MISSING;
531 battery_voltage = AltosRecord.MISSING;
532 pyro_voltage = AltosRecord.MISSING;
533 apogee_voltage = AltosRecord.MISSING;
534 main_voltage = AltosRecord.MISSING;
535 ignitor_voltage = null;
537 kalman_height = new AltosValue();
538 kalman_speed = new AltosValue();
539 kalman_acceleration = new AltosValue();
543 temp_gps_clear_sats_pending = false;
554 elevation = AltosRecord.MISSING;
555 range = AltosRecord.MISSING;
556 gps_height = AltosRecord.MISSING;
558 pad_lat = AltosRecord.MISSING;
559 pad_lon = AltosRecord.MISSING;
560 pad_alt = AltosRecord.MISSING;
562 speak_tick = AltosRecord.MISSING;
563 speak_altitude = AltosRecord.MISSING;
567 accel_plus_g = AltosRecord.MISSING;
568 accel_minus_g = AltosRecord.MISSING;
569 accel = AltosRecord.MISSING;
571 ground_accel = AltosRecord.MISSING;
572 ground_accel_avg = AltosRecord.MISSING;
574 log_format = AltosRecord.MISSING;
575 serial = AltosRecord.MISSING;
581 void copy(AltosState old) {
590 received_time = old.received_time;
594 prev_tick = old.tick;
595 boost_tick = old.boost_tick;
604 device_type = old.device_type;
605 config_major = old.config_major;
606 config_minor = old.config_minor;
607 apogee_delay = old.apogee_delay;
608 main_deploy = old.main_deploy;
609 flight_log_max = old.flight_log_max;
613 ground_altitude.copy(old.ground_altitude);
614 altitude.copy(old.altitude);
615 pressure.copy(old.pressure);
616 speed.copy(old.speed);
617 acceleration.copy(old.acceleration);
619 battery_voltage = old.battery_voltage;
620 pyro_voltage = old.pyro_voltage;
621 temperature = old.temperature;
622 apogee_voltage = old.apogee_voltage;
623 main_voltage = old.main_voltage;
624 ignitor_voltage = old.ignitor_voltage;
626 kalman_height.copy(old.kalman_height);
627 kalman_speed.copy(old.kalman_speed);
628 kalman_acceleration.copy(old.kalman_acceleration);
631 gps = old.gps.clone();
634 if (old.temp_gps != null)
635 temp_gps = old.temp_gps.clone();
638 temp_gps_clear_sats_pending = old.temp_gps_clear_sats_pending;
639 gps_sequence = old.gps_sequence;
640 gps_pending = old.gps_pending;
643 imu = old.imu.clone();
648 mag = old.mag.clone();
653 gps_waiting = old.gps_waiting;
654 gps_ready = old.gps_ready;
657 if (old.from_pad != null)
658 from_pad = old.from_pad.clone();
662 elevation = old.elevation;
665 gps_height = old.gps_height;
666 pad_lat = old.pad_lat;
667 pad_lon = old.pad_lon;
668 pad_alt = old.pad_alt;
670 speak_tick = old.speak_tick;
671 speak_altitude = old.speak_altitude;
673 callsign = old.callsign;
675 accel_plus_g = old.accel_plus_g;
676 accel_minus_g = old.accel_minus_g;
678 ground_accel = old.ground_accel;
679 ground_accel_avg = old.ground_accel_avg;
681 log_format = old.log_format;
685 companion = old.companion;
699 if (gps.locked && gps.nsat >= 4) {
700 /* Track consecutive 'good' gps reports, waiting for 10 of them */
701 if (state == AltosLib.ao_flight_pad) {
703 if (pad_lat != AltosRecord.MISSING) {
704 pad_lat = (pad_lat * 31 + gps.lat) / 32;
705 pad_lon = (pad_lon * 31 + gps.lon) / 32;
706 pad_alt = (pad_alt * 31 + gps.alt) / 32;
709 if (pad_lat == AltosRecord.MISSING) {
715 if (gps.lat != 0 && gps.lon != 0 &&
716 pad_lat != AltosRecord.MISSING &&
717 pad_lon != AltosRecord.MISSING)
721 if (h == AltosRecord.MISSING)
723 from_pad = new AltosGreatCircle(pad_lat, pad_lon, 0, gps.lat, gps.lon, h);
724 elevation = from_pad.elevation;
725 range = from_pad.range;
726 gps_height = gps.alt - pad_alt;
730 public void set_tick(int new_tick) {
731 if (new_tick != AltosRecord.MISSING) {
732 if (prev_tick != AltosRecord.MISSING) {
733 while (new_tick < prev_tick - 32767) {
742 public void set_boost_tick(int boost_tick) {
743 if (boost_tick != AltosRecord.MISSING)
744 this.boost_tick = boost_tick;
747 public String state_name() {
748 return AltosLib.state_name(state);
751 public void set_state(int state) {
752 if (state != AltosLib.ao_flight_invalid) {
754 ascent = (AltosLib.ao_flight_boost <= state &&
755 state <= AltosLib.ao_flight_coast);
756 boost = (AltosLib.ao_flight_boost == state);
761 public void set_device_type(int device_type) {
762 this.device_type = device_type;
765 public void set_config(int major, int minor, int apogee_delay, int main_deploy, int flight_log_max) {
766 config_major = major;
767 config_minor = minor;
768 this.apogee_delay = apogee_delay;
769 this.main_deploy = main_deploy;
770 this.flight_log_max = flight_log_max;
773 public void set_callsign(String callsign) {
774 this.callsign = callsign;
777 public void set_firmware_version(String version) {
778 firmware_version = version;
781 public void set_flight(int flight) {
783 /* When the flight changes, reset the state */
784 if (flight != AltosRecord.MISSING) {
785 if (this.flight != AltosRecord.MISSING &&
786 this.flight != flight) {
789 this.flight = flight;
793 public void set_serial(int serial) {
794 /* When the serial changes, reset the state */
795 if (serial != AltosRecord.MISSING) {
796 if (this.serial != AltosRecord.MISSING &&
797 this.serial != serial) {
800 this.serial = serial;
805 if (rssi == AltosRecord.MISSING)
810 public void set_rssi(int rssi, int status) {
811 if (rssi != AltosRecord.MISSING) {
813 this.status = status;
817 public void set_received_time(long ms) {
821 public void set_gps(AltosGPS gps, int sequence) {
823 this.gps = gps.clone();
824 gps_sequence = sequence;
830 public void make_baro() {
832 baro = new AltosMs5607();
835 public void set_ms5607(int pres, int temp) {
837 baro.set(pres, temp);
839 set_pressure(baro.pa);
840 set_temperature(baro.cc / 100.0);
844 public void make_companion (int nchannels) {
845 if (companion == null)
846 companion = new AltosRecordCompanion(nchannels);
849 public void set_companion(AltosRecordCompanion companion) {
850 this.companion = companion;
853 void update_accel() {
854 double ground = ground_accel;
856 if (ground == AltosRecord.MISSING)
857 ground = ground_accel_avg;
858 if (accel == AltosRecord.MISSING)
860 if (ground == AltosRecord.MISSING)
862 if (accel_plus_g == AltosRecord.MISSING)
864 if (accel_minus_g == AltosRecord.MISSING)
867 double counts_per_g = (accel_minus_g - accel_plus_g) / 2.0;
868 double counts_per_mss = counts_per_g / 9.80665;
869 acceleration.set_computed((ground - accel) / counts_per_mss, time);
872 public void set_accel_g(double accel_plus_g, double accel_minus_g) {
873 if (accel_plus_g != AltosRecord.MISSING) {
874 this.accel_plus_g = accel_plus_g;
875 this.accel_minus_g = accel_minus_g;
880 public void set_ground_accel(double ground_accel) {
881 if (ground_accel != AltosRecord.MISSING) {
882 this.ground_accel = ground_accel;
887 public void set_accel(double accel) {
888 if (accel != AltosRecord.MISSING) {
890 if (state == AltosLib.ao_flight_pad) {
891 if (ground_accel_avg == AltosRecord.MISSING)
892 ground_accel_avg = accel;
894 ground_accel_avg = (ground_accel_avg * 7 + accel) / 8;
900 public void set_temperature(double temperature) {
901 if (temperature != AltosRecord.MISSING) {
902 this.temperature = temperature;
907 public void set_battery_voltage(double battery_voltage) {
908 if (battery_voltage != AltosRecord.MISSING) {
909 this.battery_voltage = battery_voltage;
914 public void set_pyro_voltage(double pyro_voltage) {
915 if (pyro_voltage != AltosRecord.MISSING) {
916 this.pyro_voltage = pyro_voltage;
921 public void set_apogee_voltage(double apogee_voltage) {
922 if (apogee_voltage != AltosRecord.MISSING) {
923 this.apogee_voltage = apogee_voltage;
928 public void set_main_voltage(double main_voltage) {
929 if (main_voltage != AltosRecord.MISSING) {
930 this.main_voltage = main_voltage;
935 public void set_ignitor_voltage(double[] voltage) {
936 this.ignitor_voltage = voltage;
939 public double time_since_boost() {
940 if (tick == AltosRecord.MISSING)
943 if (boost_tick != AltosRecord.MISSING) {
944 return (tick - boost_tick) / 100.0;
949 public boolean valid() {
950 return tick != AltosRecord.MISSING && serial != AltosRecord.MISSING;
953 public AltosGPS make_temp_gps(boolean sats) {
954 if (temp_gps == null) {
955 temp_gps = new AltosGPS(gps);
959 temp_gps_clear_sats_pending = true;
960 else if (temp_gps_clear_sats_pending) {
961 temp_gps.cc_gps_sat = null;
962 temp_gps_clear_sats_pending = false;
967 public void set_temp_gps() {
968 set_gps(temp_gps, gps_sequence + 1);
973 public AltosState clone() {
974 AltosState s = new AltosState();
980 public void init (AltosRecord cur, AltosState prev_state) {
982 System.out.printf ("init\n");
984 cur = new AltosRecord();
988 /* Discard previous state if it was for a different board */
989 if (prev_state != null && prev_state.serial != cur.serial)
994 set_ground_altitude(cur.ground_altitude());
995 set_altitude(cur.altitude());
997 set_kalman(cur.kalman_height, cur.kalman_speed, cur.kalman_acceleration);
999 received_time = System.currentTimeMillis();
1001 set_temperature(cur.temperature());
1002 set_apogee_voltage(cur.drogue_voltage());
1003 set_main_voltage(cur.main_voltage());
1004 set_battery_voltage(cur.battery_voltage());
1006 set_pressure(cur.pressure());
1009 set_state(cur.state);
1011 set_accel_g (cur.accel_minus_g, cur.accel_plus_g);
1012 set_ground_accel(cur.ground_accel);
1013 set_accel (cur.accel);
1015 if (cur.gps_sequence != gps_sequence)
1016 set_gps(cur.gps, cur.gps_sequence);
1020 public AltosState(AltosRecord cur) {
1024 public AltosState (AltosRecord cur, AltosState prev) {
1029 public AltosState () {