2 * Copyright © 2009 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.
29 #define AO_ADC_RING 64
30 #define ao_adc_ring_next(n) (((n) + 1) & (AO_ADC_RING - 1))
31 #define ao_adc_ring_prev(n) (((n) - 1) & (AO_ADC_RING - 1))
33 #define AO_M_TO_HEIGHT(m) ((int16_t) (m))
34 #define AO_MS_TO_SPEED(ms) ((int16_t) ((ms) * 16))
35 #define AO_MSS_TO_ACCEL(mss) ((int16_t) ((mss) * 16))
38 * One set of samples read from the A/D converter
41 uint16_t tick; /* tick when the sample was read */
42 int16_t accel; /* accelerometer */
43 int16_t pres; /* pressure sensor */
44 int16_t temp; /* temperature sensor */
45 int16_t v_batt; /* battery voltage */
46 int16_t sense_d; /* drogue continuity sense */
47 int16_t sense_m; /* main continuity sense */
56 #define to_fix16(x) ((int16_t) ((x) * 65536.0 + 0.5))
57 #define to_fix32(x) ((int32_t) ((x) * 65536.0 + 0.5))
58 #define from_fix(x) ((x) >> 16)
61 * Above this height, the baro sensor doesn't work
63 #define AO_MAX_BARO_HEIGHT 12000
66 * Above this speed, baro measurements are unreliable
68 #define AO_MAX_BARO_SPEED 200
70 #define ACCEL_NOSE_UP (ao_accel_2g >> 2)
72 enum ao_flight_state {
73 ao_flight_startup = 0,
85 extern enum ao_flight_state ao_flight_state;
90 struct ao_adc ao_adc_ring[AO_ADC_RING];
96 #define ao_timer_set_adc_interval(i)
97 #define ao_wakeup(wchan) ao_dump_state()
98 #define ao_cmd_register(c)
99 #define ao_usb_disable()
100 #define ao_telemetry_set_interval(x)
101 #define ao_rdf_set(rdf)
102 #define ao_packet_slave_start()
103 #define ao_packet_slave_stop()
106 ao_igniter_drogue = 0,
110 struct ao_adc ao_adc_static;
119 static int32_t ao_k_height;
122 ao_ignite(enum ao_igniter igniter)
124 double time = (double) (ao_adc_static.tick + tick_offset) / 100;
126 if (igniter == ao_igniter_drogue) {
128 drogue_height = ao_k_height >> 16;
131 main_height = ao_k_height >> 16;
139 #define ao_add_task(t,f,n)
141 #define ao_log_start()
142 #define ao_log_stop()
144 #define AO_MS_TO_TICKS(ms) ((ms) / 10)
145 #define AO_SEC_TO_TICKS(s) ((s) * 100)
147 #define AO_FLIGHT_TEST
154 double emulator_error_max = 4;
160 ao_sleep(void *wchan);
162 const char const * const ao_state_names[] = {
163 "startup", "idle", "pad", "boost", "fast",
164 "coast", "drogue", "main", "landed", "invalid"
172 #include "ao_convert.c"
175 uint16_t main_deploy;
176 int16_t accel_plus_g;
177 int16_t accel_minus_g;
180 #define ao_config_get()
182 struct ao_config ao_config;
184 #define DATA_TO_XDATA(x) (x)
192 #define HAS_ACCEL_REF 0
195 #define GRAVITY 9.80665
196 extern int16_t ao_ground_accel, ao_flight_accel;
197 extern int16_t ao_accel_2g;
199 extern uint16_t ao_sample_tick;
201 extern int16_t ao_sample_height;
202 extern int16_t ao_sample_accel;
203 extern int32_t ao_accel_scale;
205 int ao_sample_prev_tick;
208 #include "ao_kalman.c"
209 #include "ao_sample.c"
210 #include "ao_flight.c"
212 #define to_double(f) ((f) / 65536.0)
214 static int ao_records_read = 0;
215 static int ao_eof_read = 0;
216 static int ao_flight_ground_accel;
217 static int ao_flight_started = 0;
218 static int ao_test_max_height;
219 static double ao_test_max_height_time;
220 static int ao_test_main_height;
221 static double ao_test_main_height_time;
229 if (!ao_test_main_height_time) {
230 ao_test_main_height_time = ao_test_max_height_time;
231 ao_test_main_height = ao_test_max_height;
233 drogue_error = fabs(ao_test_max_height_time - drogue_time);
234 main_error = fabs(ao_test_main_height_time - main_time);
235 if (drogue_error > emulator_error_max || main_error > emulator_error_max) {
237 emulator_app, emulator_name);
238 printf ("\tApogee error %g\n", drogue_error);
239 printf ("\tMain error %g\n", main_error);
240 printf ("\tActual: apogee: %d at %7.2f main: %d at %7.2f\n",
241 ao_test_max_height, ao_test_max_height_time,
242 ao_test_main_height, ao_test_main_height_time);
243 printf ("\tComputed: apogee: %d at %7.2f main: %d at %7.2f\n",
244 drogue_height, drogue_time, main_height, main_time);
255 ao_adc_ring[ao_adc_head] = ao_adc_static;
256 ao_adc_head = ao_adc_ring_next(ao_adc_head);
257 if (ao_flight_state != ao_flight_startup) {
258 double height = ao_pres_to_altitude(ao_sample_pres) - ao_ground_height;
259 double accel = ((ao_flight_ground_accel - ao_adc_static.accel) * GRAVITY * 2.0) /
260 (ao_config.accel_minus_g - ao_config.accel_plus_g);
263 tick_offset = -ao_adc_static.tick;
264 if ((prev_tick - ao_adc_static.tick) > 0x400)
265 tick_offset += 65536;
266 prev_tick = ao_adc_static.tick;
267 time = (double) (ao_adc_static.tick + tick_offset) / 100;
269 if (ao_test_max_height < height) {
270 ao_test_max_height = height;
271 ao_test_max_height_time = time;
273 if (height > ao_config.main_deploy) {
274 ao_test_main_height_time = time;
275 ao_test_main_height = height;
279 printf("%7.2f height %g accel %g state %s k_height %g k_speed %g k_accel %g drogue %d main %d error %d\n",
283 ao_state_names[ao_flight_state],
284 ao_k_height / 65536.0,
285 ao_k_speed / 65536.0 / 16.0,
286 ao_k_accel / 65536.0 / 16.0,
290 if (ao_flight_state == ao_flight_landed)
297 ao_sleep(void *wchan)
299 if (wchan == &ao_adc_head) {
311 if (ao_records_read > 2 && ao_flight_state == ao_flight_startup)
313 ao_adc_static.accel = ao_flight_ground_accel;
318 if (!fgets(line, sizeof (line), emulator_in)) {
319 if (++ao_eof_read >= 1000) {
321 printf ("no more data, exiting simulation\n");
324 ao_adc_static.tick += 10;
329 for (nword = 0; nword < 64; nword++) {
330 words[nword] = strtok_r(l, " \t\n", &saveptr);
332 if (words[nword] == NULL)
337 tick = strtoul(words[1], NULL, 16);
338 a = strtoul(words[2], NULL, 16);
339 b = strtoul(words[3], NULL, 16);
340 } else if (nword >= 6 && strcmp(words[0], "Accel") == 0) {
341 ao_config.accel_plus_g = atoi(words[3]);
342 ao_config.accel_minus_g = atoi(words[5]);
343 } else if (nword >= 4 && strcmp(words[0], "Main") == 0) {
344 ao_config.main_deploy = atoi(words[2]);
345 } else if (nword >= 36 && strcmp(words[0], "CALL") == 0) {
346 tick = atoi(words[10]);
347 if (!ao_flight_started) {
350 ao_flight_started = 1;
357 if (type != 'F' && !ao_flight_started)
362 ao_flight_ground_accel = a;
363 if (ao_config.accel_plus_g == 0) {
364 ao_config.accel_plus_g = a;
365 ao_config.accel_minus_g = a + 530;
367 if (ao_config.main_deploy == 0)
368 ao_config.main_deploy = 250;
369 ao_flight_started = 1;
374 ao_adc_static.tick = tick;
375 ao_adc_static.accel = a;
376 ao_adc_static.pres = b;
381 ao_adc_static.tick = tick;
382 ao_adc_static.temp = a;
383 ao_adc_static.v_batt = b;
396 #define COUNTS_PER_G 264.8
403 static const struct option options[] = {
404 { .name = "summary", .has_arg = 0, .val = 's' },
405 { .name = "debug", .has_arg = 0, .val = 'd' },
409 void run_flight_fixed(char *name, FILE *f, int summary)
411 emulator_name = name;
413 ao_summary = summary;
419 main (int argc, char **argv)
430 while ((c = getopt_long(argc, argv, "sd", options, NULL)) != -1) {
442 run_flight_fixed("<stdin>", stdin, summary);
444 for (i = optind; i < argc; i++) {
445 FILE *f = fopen(argv[i], "r");
450 run_flight_fixed(argv[i], f, summary);