2019-01-26 17:02:19 +00:00
|
|
|
/*
|
|
|
|
* wx_handler.c
|
|
|
|
*
|
|
|
|
* Created on: 26.01.2019
|
|
|
|
* Author: mateusz
|
|
|
|
*/
|
|
|
|
|
2019-12-16 21:40:42 +00:00
|
|
|
#include "wx_handler.h"
|
|
|
|
|
2019-01-26 22:18:25 +00:00
|
|
|
#include <rte_wx.h>
|
2019-12-16 21:40:42 +00:00
|
|
|
#include <stm32f10x.h>
|
2020-01-04 13:15:33 +00:00
|
|
|
#include <math.h>
|
2019-01-26 17:02:19 +00:00
|
|
|
#include "drivers/_dht22.h"
|
|
|
|
#include "drivers/ms5611.h"
|
2020-01-04 13:15:33 +00:00
|
|
|
#include "drivers/analog_anemometer.h"
|
2019-01-26 22:18:25 +00:00
|
|
|
|
|
|
|
#include "station_config.h"
|
2019-01-26 17:02:19 +00:00
|
|
|
|
2019-12-16 21:40:42 +00:00
|
|
|
#include "delay.h"
|
2019-04-15 20:46:36 +00:00
|
|
|
#include "telemetry.h"
|
2019-12-16 21:40:42 +00:00
|
|
|
#include "main.h"
|
|
|
|
|
|
|
|
#define WX_WATCHDOG_PERIOD (SYSTICK_TICKS_PER_SECONDS * SYSTICK_TICKS_PERIOD * 90)
|
|
|
|
#define WX_WATCHDOG_RESET_DURATION (SYSTICK_TICKS_PER_SECONDS * SYSTICK_TICKS_PERIOD * 3)
|
|
|
|
|
2019-12-22 11:34:03 +00:00
|
|
|
#define WX_MAX_TEMPERATURE_SLEW_RATE 4.0f
|
|
|
|
|
2019-12-16 21:40:42 +00:00
|
|
|
uint32_t wx_last_good_temperature_time = 0;
|
|
|
|
uint32_t wx_last_good_wind_time = 0;
|
|
|
|
wx_pwr_state_t wx_pwr_state;
|
2019-12-22 11:34:03 +00:00
|
|
|
uint8_t wx_inhibit_slew_rate_check = 1;
|
2019-04-15 20:46:36 +00:00
|
|
|
|
2020-02-01 22:20:17 +00:00
|
|
|
static const float direction_constant = M_PI/180.0f;
|
|
|
|
|
2019-01-26 17:02:19 +00:00
|
|
|
void wx_get_all_measurements(void) {
|
|
|
|
|
|
|
|
int32_t return_value = 0;
|
|
|
|
|
2019-09-08 21:42:59 +00:00
|
|
|
#ifdef _METEO
|
|
|
|
// quering MS5611 sensor for temperature
|
2019-10-27 19:50:01 +00:00
|
|
|
return_value = ms5611_get_temperature(&rte_wx_temperature_ms, &rte_wx_ms5611_qf);
|
2019-09-08 21:42:59 +00:00
|
|
|
|
|
|
|
if (return_value == MS5611_OK) {
|
2019-10-27 19:50:01 +00:00
|
|
|
rte_wx_temperature_ms_valid = rte_wx_temperature_ms;
|
2019-09-08 21:42:59 +00:00
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
#endif
|
|
|
|
|
2019-01-26 22:18:25 +00:00
|
|
|
#if defined _METEO || defined _DALLAS_AS_TELEM
|
|
|
|
|
2019-01-26 17:02:19 +00:00
|
|
|
// quering dallas DS12B20 thermometer for current temperature
|
2019-04-15 20:46:36 +00:00
|
|
|
rte_wx_temperature_dallas = dallas_query(&rte_wx_current_dallas_qf);
|
2019-01-26 17:02:19 +00:00
|
|
|
|
|
|
|
// checking if communication was successfull
|
2019-04-15 20:46:36 +00:00
|
|
|
if (rte_wx_temperature_dallas != -128.0f) {
|
|
|
|
|
2019-12-22 11:34:03 +00:00
|
|
|
// calculate the slew rate
|
|
|
|
rte_wx_temperature_dalls_slew_rate = rte_wx_temperature_dallas - rte_wx_temperature_dallas_valid;
|
|
|
|
|
|
|
|
// chcecking the positive (ascending) slew rate of the temperature measuremenets
|
|
|
|
if (rte_wx_temperature_dalls_slew_rate > WX_MAX_TEMPERATURE_SLEW_RATE && wx_inhibit_slew_rate_check == 0) {
|
|
|
|
|
|
|
|
// if temeperature measuremenet has changed more than maximum allowed slew rate set degradadet QF
|
|
|
|
rte_wx_error_dallas_qf = DALLAS_QF_DEGRADATED;
|
|
|
|
|
|
|
|
// and increase the temperature only by 1.0f to decrease slew rate
|
|
|
|
rte_wx_temperature_dallas += 1.0f;
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
// chcecking the negaive (descending) slew rate of the temperature measuremenets
|
|
|
|
if (rte_wx_temperature_dalls_slew_rate < -WX_MAX_TEMPERATURE_SLEW_RATE && wx_inhibit_slew_rate_check == 0) {
|
|
|
|
|
|
|
|
// if temeperature measuremenet has changed more than maximum allowed slew rate set degradadet QF
|
|
|
|
rte_wx_error_dallas_qf = DALLAS_QF_DEGRADATED;
|
|
|
|
|
|
|
|
// and decrease the temperature only by 1.0f to decrease slew rate
|
|
|
|
rte_wx_temperature_dallas -= 1.0f;
|
|
|
|
|
|
|
|
}
|
|
|
|
|
2019-10-27 19:50:01 +00:00
|
|
|
// store current value
|
2019-01-26 17:02:19 +00:00
|
|
|
rte_wx_temperature_dallas_valid = rte_wx_temperature_dallas;
|
2019-04-15 20:46:36 +00:00
|
|
|
|
2019-10-27 19:50:01 +00:00
|
|
|
// include current temperature into the average
|
|
|
|
dallas_average(rte_wx_temperature_dallas, &rte_wx_dallas_average);
|
|
|
|
|
|
|
|
// update the current temperature with current average
|
|
|
|
rte_wx_temperature_average_dallas_valid = dallas_get_average(&rte_wx_dallas_average);
|
|
|
|
|
|
|
|
// update current minimal temperature
|
|
|
|
rte_wx_temperature_min_dallas_valid = dallas_get_min(&rte_wx_dallas_average);
|
|
|
|
|
|
|
|
// and update maximum also
|
|
|
|
rte_wx_temperature_max_dallas_valid = dallas_get_max(&rte_wx_dallas_average);
|
2019-12-16 21:40:42 +00:00
|
|
|
|
|
|
|
// updating last good measurement time
|
|
|
|
wx_last_good_temperature_time = master_time;
|
2019-04-15 20:46:36 +00:00
|
|
|
}
|
|
|
|
else {
|
2019-04-28 07:36:24 +00:00
|
|
|
// if there were a communication error set the error to unavaliable
|
2019-04-15 20:46:36 +00:00
|
|
|
rte_wx_error_dallas_qf = DALLAS_QF_NOT_AVALIABLE;
|
2019-04-28 07:36:24 +00:00
|
|
|
|
2019-04-15 20:46:36 +00:00
|
|
|
}
|
2019-12-22 11:34:03 +00:00
|
|
|
|
|
|
|
// enabling slew rate checking after first power up
|
|
|
|
wx_inhibit_slew_rate_check = 0;
|
2019-01-26 22:18:25 +00:00
|
|
|
#endif
|
2019-01-26 17:02:19 +00:00
|
|
|
|
2019-01-26 22:18:25 +00:00
|
|
|
#ifdef _METEO
|
2019-01-26 17:02:19 +00:00
|
|
|
// quering MS5611 sensor for pressure
|
|
|
|
return_value = ms5611_get_pressure(&rte_wx_pressure, &rte_wx_ms5611_qf);
|
|
|
|
|
|
|
|
if (return_value == MS5611_OK) {
|
|
|
|
rte_wx_pressure_valid = rte_wx_pressure;
|
|
|
|
}
|
|
|
|
|
|
|
|
// if humidity sensor is idle trigger the communiction & measuremenets
|
|
|
|
if (dht22State == DHT22_STATE_DONE || dht22State == DHT22_STATE_TIMEOUT)
|
|
|
|
dht22State = DHT22_STATE_IDLE;
|
|
|
|
|
2019-01-26 22:18:25 +00:00
|
|
|
#endif
|
2019-01-26 17:02:19 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
void wx_pool_dht22(void) {
|
|
|
|
|
|
|
|
dht22_timeout_keeper();
|
|
|
|
|
|
|
|
switch (dht22State) {
|
|
|
|
case DHT22_STATE_IDLE:
|
|
|
|
dht22_comm(&rte_wx_dht);
|
|
|
|
break;
|
|
|
|
case DHT22_STATE_DATA_RDY:
|
|
|
|
dht22_decode(&rte_wx_dht);
|
|
|
|
break;
|
|
|
|
case DHT22_STATE_DATA_DECD:
|
|
|
|
rte_wx_dht_valid = rte_wx_dht; // powrot do stanu DHT22_STATE_IDLE jest w TIM3_IRQHandler
|
2019-01-27 11:34:43 +00:00
|
|
|
//rte_wx_dht_valid.qf = DHT22_QF_FULL;
|
2019-01-26 17:02:19 +00:00
|
|
|
dht22State = DHT22_STATE_DONE;
|
|
|
|
|
|
|
|
#ifdef _DBG_TRACE
|
|
|
|
trace_printf("DHT22: temperature=%d,humi=%d\r\n", dht_valid.scaledTemperature, dht_valid.humidity);
|
|
|
|
#endif
|
|
|
|
break;
|
2019-01-27 11:34:43 +00:00
|
|
|
case DHT22_STATE_TIMEOUT:
|
|
|
|
rte_wx_dht_valid.qf = DHT22_QF_UNAVALIABLE;
|
|
|
|
break;
|
2019-01-26 17:02:19 +00:00
|
|
|
default: break;
|
|
|
|
}
|
|
|
|
|
|
|
|
}
|
2019-12-16 21:40:42 +00:00
|
|
|
|
2020-01-04 13:15:33 +00:00
|
|
|
void wx_pool_analog_anemometer(void) {
|
|
|
|
|
|
|
|
// locals
|
|
|
|
uint32_t average_windspeed = 0;
|
2020-02-01 22:20:17 +00:00
|
|
|
int32_t wind_direction_x_avg = 0;
|
|
|
|
int32_t wind_direction_y_avg = 0;
|
|
|
|
int16_t wind_direction_x = 0;
|
|
|
|
int16_t wind_direction_y = 0;
|
|
|
|
volatile float dir_temp = 0;
|
|
|
|
volatile float arctan_value = 0.0f;
|
2020-01-04 13:15:33 +00:00
|
|
|
short i = 0;
|
|
|
|
|
2020-02-23 06:50:22 +00:00
|
|
|
#ifdef _ANEMOMETER_ANALOGUE
|
2020-01-04 13:15:33 +00:00
|
|
|
// this windspeed is scaled * 10. Example: 0.2 meters per second is stored as 2
|
|
|
|
uint16_t scaled_windspeed = analog_anemometer_get_ms_from_pulse(rte_wx_windspeed_pulses);
|
2020-02-23 06:50:22 +00:00
|
|
|
#else
|
|
|
|
uint16_t scaled_windspeed = 0;
|
|
|
|
#endif
|
2020-01-04 13:15:33 +00:00
|
|
|
|
|
|
|
// putting the wind speed into circular buffer
|
|
|
|
rte_wx_windspeed[rte_wx_windspeed_it] = scaled_windspeed;
|
|
|
|
|
|
|
|
// increasing the iterator to the windspeed buffer
|
|
|
|
rte_wx_windspeed_it++;
|
|
|
|
|
|
|
|
// checking if iterator reached an end of the buffer
|
|
|
|
if (rte_wx_windspeed_it >= WIND_AVERAGE_LEN)
|
|
|
|
rte_wx_windspeed_it = 0;
|
|
|
|
|
|
|
|
// calculating the average windspeed
|
|
|
|
for (i = 0; i < WIND_AVERAGE_LEN; i++)
|
|
|
|
average_windspeed += rte_wx_windspeed[i];
|
|
|
|
|
|
|
|
average_windspeed /= WIND_AVERAGE_LEN;
|
|
|
|
|
|
|
|
// store the value in rte
|
|
|
|
rte_wx_average_windspeed = average_windspeed;
|
|
|
|
|
|
|
|
// reuse the local variable to find maximum value
|
|
|
|
average_windspeed = 0;
|
|
|
|
|
|
|
|
// looking for gusts
|
|
|
|
for (i = 0; i < WIND_AVERAGE_LEN; i++) {
|
|
|
|
if (average_windspeed < rte_wx_windspeed[i])
|
|
|
|
average_windspeed = rte_wx_windspeed[i];
|
|
|
|
}
|
|
|
|
|
|
|
|
// storing wind gusts value in rte
|
|
|
|
rte_wx_max_windspeed = average_windspeed;
|
|
|
|
|
2020-02-01 22:20:17 +00:00
|
|
|
// adding last wind direction to the buffers
|
|
|
|
if (rte_wx_winddirection_it >= WIND_AVERAGE_LEN)
|
|
|
|
rte_wx_winddirection_it = 0;
|
|
|
|
|
|
|
|
rte_wx_winddirection[rte_wx_winddirection_it++] = rte_wx_winddirection_last;
|
|
|
|
|
2020-01-04 13:15:33 +00:00
|
|
|
// calculating average wind direction
|
|
|
|
for (i = 0; i < WIND_AVERAGE_LEN; i++) {
|
2020-02-01 22:20:17 +00:00
|
|
|
|
|
|
|
dir_temp = (float)rte_wx_winddirection[i];
|
|
|
|
|
2020-01-04 13:15:33 +00:00
|
|
|
// split the wind direction into x and y component
|
2020-02-01 22:20:17 +00:00
|
|
|
wind_direction_x = (int16_t)(100.0f * cosf(dir_temp * direction_constant));
|
|
|
|
wind_direction_y = (int16_t)(100.0f * sinf(dir_temp * direction_constant));
|
2020-01-04 13:15:33 +00:00
|
|
|
|
|
|
|
// adding components to calculate average
|
|
|
|
wind_direction_x_avg += wind_direction_x;
|
|
|
|
wind_direction_y_avg += wind_direction_y;
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
// dividing to get average of x and y componen
|
|
|
|
wind_direction_x_avg /= WIND_AVERAGE_LEN;
|
|
|
|
wind_direction_y_avg /= WIND_AVERAGE_LEN;
|
|
|
|
|
|
|
|
// converting x & y component of wind direction back to an angle
|
2020-02-01 22:20:17 +00:00
|
|
|
arctan_value = atan2f(wind_direction_y_avg , wind_direction_x_avg);
|
|
|
|
|
|
|
|
rte_wx_average_winddirection = (int16_t)(arctan_value * (180.0f/M_PI));
|
2020-01-04 13:15:33 +00:00
|
|
|
|
|
|
|
if (rte_wx_average_winddirection < 0)
|
|
|
|
rte_wx_average_winddirection += 360;
|
|
|
|
}
|
|
|
|
|
2019-12-16 21:40:42 +00:00
|
|
|
void wx_pwr_init(void) {
|
|
|
|
GPIO_InitTypeDef GPIO_InitStructure;
|
|
|
|
GPIO_InitStructure.GPIO_Pin = GPIO_Pin_8;
|
|
|
|
GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;
|
|
|
|
GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_OD;
|
|
|
|
GPIO_Init(GPIOB, &GPIO_InitStructure);
|
|
|
|
|
|
|
|
wx_pwr_state = WX_PWR_OFF;
|
|
|
|
|
|
|
|
GPIO_ResetBits(GPIOB, GPIO_Pin_8);
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
void wx_pwr_periodic_handle(void) {
|
|
|
|
// check when last measuremenets was sent by wind or temperature sensor
|
|
|
|
if ( (master_time - wx_last_good_temperature_time >= WX_WATCHDOG_PERIOD ||
|
|
|
|
master_time - wx_last_good_wind_time >= WX_WATCHDOG_PERIOD) &&
|
|
|
|
wx_pwr_state == WX_PWR_ON) {
|
|
|
|
|
|
|
|
// if timeout watchod expired there is a time to reset the supply voltage
|
|
|
|
wx_pwr_state = WX_PWR_UNDER_RESET;
|
|
|
|
|
|
|
|
// pulling the output down to switch the relay
|
|
|
|
GPIO_ResetBits(GPIOB, GPIO_Pin_8);
|
|
|
|
|
|
|
|
// setting the last_good timers to current value to prevent reset loop
|
|
|
|
wx_last_good_temperature_time = master_time;
|
|
|
|
wx_last_good_wind_time = master_time;
|
|
|
|
|
|
|
|
return;
|
|
|
|
}
|
|
|
|
|
|
|
|
// service actual supply state
|
|
|
|
switch (wx_pwr_state) {
|
|
|
|
case WX_PWR_OFF:
|
|
|
|
|
|
|
|
// one second delay
|
|
|
|
delay_fixed(2000);
|
|
|
|
|
|
|
|
GPIO_SetBits(GPIOB, GPIO_Pin_8);
|
|
|
|
|
|
|
|
// power is off after power-up and needs to be powered on
|
|
|
|
wx_pwr_state = WX_PWR_ON;
|
|
|
|
break;
|
|
|
|
case WX_PWR_ON:
|
|
|
|
break;
|
|
|
|
case WX_PWR_UNDER_RESET:
|
|
|
|
|
|
|
|
GPIO_SetBits(GPIOB, GPIO_Pin_8);
|
|
|
|
|
|
|
|
wx_pwr_state = WX_PWR_ON;
|
|
|
|
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
}
|