AND
Exluding ZMOD
Dependencies: iAQ_Core Adafruit_SGP30_mbed mbed BME680
Diff: main.cpp
- Revision:
- 5:44f7dbc5c75d
- Parent:
- 3:3d51f8870e91
- Child:
- 6:96de1304b9e7
diff -r 54dc2a95c130 -r 44f7dbc5c75d main.cpp
--- a/main.cpp Fri Jul 26 09:31:58 2019 +0000
+++ b/main.cpp Fri May 01 14:31:45 2020 +0000
@@ -1,222 +1,61 @@
#include "mbed.h"
-#include "flow.h"
+#include "iAQ_Core.h"
+
+iAQ_Core myiAQ_Core ( PB_11,PB_10, iAQ_Core::iAQ_Core_ADDRESS );
+Serial ble(PC_12,PD_2);
+
+Ticker newReading;
-/////////////////////////
-// In this version of the program developed for the Breath project, flow and CO2, as well as 8 channel sensors,
-// are measured in a separate .h file called: "flow.h" which is included in the
-// main code. So 10 sets of data is streamed using a serial connection (TTL cable or Bluetooth)
-// without any intruption.
-// This version is especially suitable to be used for KST.
-// Also, a solenoid would be turned on and off based on calculating the standard deviation in CO2 profile.
-//
-// START POINT: calculates SD for 9 samples of CO2, if it's grater than 0.02 it enables the solenoid.
-// END POINT: calculates SD for 9 samples of CO2, if it's grater thatn 0.05 it disables the solenoid.
-//
-// You can easily change the threshold of Standard deviation to detect plateau
-//
-// Generated by: Mehrnaz Javadipour
-//////////////////////////
+iAQ_Core::iAQ_Core_status_t aux;
+iAQ_Core::iAQ_Core_data_t myiAQ_Core_data;
+uint32_t myState = 0;
+
-Serial ttl(PC_12,PD_2); //TTL cable TX,RX
-DigitalOut sol(PC_5); //Solenoid: Digital Output
-PwmOut led(PB_6);
-
+void changeDATA ( void )
+{
+ myState = 1;
+}
int main()
{
- ttl.baud(9600); //baudrate for the serial connection
- flow(); //calling flow from flow.h
- carbon(); //calling CO2 from flow.h
- s1(); //calling 8 channels from flow.h
- s2();
- s3();
- s4();
- s5();
- s6();
- s7();
- s8();
- getTemp(); //calling Temperature from flow.h
-
- //////////////////////////////
- // I defined a flag for each section of specific functions, so by enabling the
- // flag the section starts and by disabling the flag it finishes the section.
- // at the end of the program, I reset the flags so it would be ready for the next loop.
- /////////////////////////////
-
- int bf=0; //FLAG for detecting base flow
- int i=0;
- float bfArray[4]; //sampling flow for finding the average base flow
- float sf=0; //sum of flow samples for calculating base flow
- float fv=0; //final value of base flow
+ ble.baud(9600);
- int measurement_started=0; //FLAG for starting calculations after detecting breath
-
- int solstart=0; //FLAG for starting calculations for detecting plateau
- int m=0;
- int myArray[9]; //sampling 9 values of CO2
- unsigned int sum=0; //sum of 9 samples of CO2
- int avg=0; //average of 9 samples of CO2
- int difSum=0; //used for the Standard deviation algorithm
- long double var=0.0; //used for the Standard deviation algorithm
- float sigma=0.0; //final value for standar deviation
- int flags=0; //FLAG for keep taking samples from CO2 profile when it's too early to detect plateau
-
- int solend=0; //FLAG for ending calculations for detecting plateau
- unsigned int sum2=0; //same as before; used for finding standard deviation
- long double var2=0.0;
- float sigma2=0.0;
- int difSum2=0;
- int avg2=0;
- int flage=0; //FLAG for keep taking samples from CO2 profile when it's too early to finish plateau
-
- int fin=0;
-
+ uint32_t myWarmUpCounter = 0;
+
+
-
+ // iAQ-Core warm up is at least 5 minutes ( 300 * 1s ) or when the sensor is ready
+ do
+ {
+ aux = myiAQ_Core.iAQ_Core_GetNewReading ( &myiAQ_Core_data );
+ wait(1);
+ myWarmUpCounter++;
+ } while( ( myWarmUpCounter < 300 ) && ( myiAQ_Core_data.status == iAQ_Core::iAQ_Core_STATUS_RUNIN ) );
+
+
+ newReading.attach( &changeDATA, 0.1); // the address of the function to be attached ( changeDATA ) and the interval ( 1s )
+
+ // Let the callbacks take care of everything
while(1)
{
- led=1.00f; //an LED is fully turned on at the beginning, the brightness will be reduced when the plateau is detected.
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
+ sleep();
+
+
- if (bf==0) //finding base flow before breath
- {
- for(i=0; i<4; i++)
- {
- bfArray[i]=flow();
- sf+=bfArray[i];
- }
- fv=sf/4;
- fv=fv+0.2;
- //ttl.printf("set\n");
- bf=1;
- }
-
- //Starts calculations when it detects breathing into the device:
-
- if ((flow()>fv) and (measurement_started ==0))
- {
- measurement_started = 1;
- }
-
- //Starts detecting plateau:
-
- if ((measurement_started == 1) and (solstart==0))
- {
- //Takes 9 samples of CO2:
- //I have also included printing the values inside the loops so we don't loose any data during calculatins.
- for(m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- myArray[m]=carbon();
- }
- while(flags==0)
- {
- //While "flags" is enabled, keeps calculating the standard deviation.
- for(int m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- sum+=myArray[m];
- }
- avg=sum/9;
- for(int m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- difSum+=(myArray[m]-avg)*(myArray[m]-avg); //Find sum of difference between value X and mean
- }
- var=difSum/9;
- sigma=sqrt(var);
- if (sigma<0.02)
- {
-
- //if SD is less than 0.02 it means that it is too early to start the plateau
- //So we shift all but the first sample and define the new set of arrays:
-
- for(int m=0;m<8;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- myArray[m]=myArray[m+1]; //Shift all CO2 values to the left by 1 value
- }
- myArray[8]=carbon(); //assign a new value for the 9th sample
- }
- //The new set of arrays are now generated and is sent back to be used for preveious SD calculations.
- //If sigma for the new set is still small, a newer set will be generated and replaced
- //Otherwise, it's accepted and will turn on the solenoid:
- else
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- sol=1; //Solenoid is ON
- led=0.50f; //The brightness is reduced to half during the plateau
- flags=1; //breakes the while loop
- }
- }
- solend=1; //prepares the next section for finishing the plateau
- solstart =1;
- }
- if ((measurement_started == 1) and (solend==1))
- {
- // same process happens for finishing the plateau:
-
- for(m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- myArray[m]=carbon();
- }
- while(flage==0)
- {
- for(int m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- sum2+=myArray[m];
- }
- avg2=sum2/9;
- for(int m=0;m<9;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- difSum2+=(myArray[m]-avg2)*(myArray[m]-avg2);
- }
- var2=difSum2/9;
- sigma2=sqrt(var2);
- if (sigma2<0.05)
- {
- // here we defined the end threshold to be 0.05, it can be changed later based on experiment results
- for(int m=0;m<8;m++)
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- myArray[m]=myArray[m+1];
- }
- myArray[8]=carbon();
- }else
- {
- ttl.printf("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",getTemp(),flow(), carbon(),s1(),s2(),s3(),s4(),s5(),s6(),s7(),s8());
- sol=0; //Solenoid is OFF
- flage=1; //breakes the loop
- }
- }
- solend =0; //end of this section
- led=1.00f; //LED is back to full brightness
- bf=0; //reset the detecting base flow flag
- fin=1; //enables the next section flag
- }
- if((carbon()<10000) and (fin ==1))
- {
- //User has to wait for the CO2 level to drop less than 1% before testing again.
- //Once it is less than 1%, all the flags and parameters used in calculations are reset
- measurement_started =0;
- solstart=0;
- sum=0;
- var=0.0;
- sigma=0.0;
- difSum=0;
- sum2=0;
- var2=0.0;
- sigma2=0.0;
- difSum2=0;
- avg2=0;
- avg=0;
- flags=0;
- flage=0;
- fin=0;
- }
+ if ( myState == 1 ) {
+ // New reading
+ do {
+ aux = myiAQ_Core.iAQ_Core_GetNewReading ( &myiAQ_Core_data );
+ wait_ms(1);
+ } while( myiAQ_Core_data.status != iAQ_Core::iAQ_Core_STATUS_OK ); // [TODO] Dangerous!!! The uC may get stuck here if something goes wrong!
+ // [WORKAROUND] Insert a counter.
+
+ // Send data through the UART
+ ble.printf( "Pred: %d | Tvoc: %d | Resistance: %d\r\n", myiAQ_Core_data.pred, myiAQ_Core_data.Tvoc, myiAQ_Core_data.resistance );
+ wait(0.01);
+ myState = 0; // Reset the variable
+ }
+
}
}
\ No newline at end of file