TwoSensorsSkeleton.c is functional

Lab_1/TwoSensorsSkeleton.c

@@ -0,0 +1,179 @@
+#include <stdio.h>
+#include <math.h>
+
+// Forward TC function
+float NISTdegCtoMilliVoltsKtype(float tempDegC);  // returns EMF in millivolts
+
+// Inverse TC function
+float NISTmilliVoltsToDegCKtype(float tcEMFmV);  // returns temp in degC assuming 0 degC cold jcn
+
+float adc_to_voltage(float v_ref, int n_adc);
+float kelvin_to_c(float k);
+float resistance_to_temperature(float r);
+float voltage_to_resistance(float v);
+float voltage_to_erc(float v);
+
+int main()
+{
+	// Define VRef
+	float v_ref = 5, e_rc, e_comp;
+	float thermistor_temp, thermocouple_temp;
+
+	// User input for pins A0 and A1
+	int thermistor_pin, thermocouple_pin;
+	printf("Enter values for thermistor pin, thermocouple pin: ");
+	scanf("%d %d", &thermistor_pin, &thermocouple_pin);
+
+	// Calculate thermistor temperature in degrees C ( Part b, i,ii,iii & v)
+	thermistor_temp = kelvin_to_c(resistance_to_temperature(voltage_to_resistance(adc_to_voltage(v_ref, thermistor_pin))));
+
+	// Calculate thermocouple temperature in degrees C ( Part c, i - iv)
+	e_rc = 1000*voltage_to_erc(adc_to_voltage(v_ref, thermistor_pin));
+	e_comp = NISTdegCtoMilliVoltsKtype(thermistor_temp);
+	thermocouple_temp = NISTmilliVoltsToDegCKtype(e_rc + e_comp);
+
+
+	// Output results
+	printf("Thermistor temperature (deg C): %f \n", thermistor_temp);
+	printf("Thermocouple temperature with CJC (deg C): %f \n", thermocouple_temp);
+
+	return 0;
+}
+
+
+/* Write a function here to convert ADC value to voltages. (Part a, equation 1)
+Call it from the main() function above */
+float adc_to_voltage(float v_ref, int n_adc) {
+	return (float)n_adc*v_ref/1024.0;
+}
+
+
+/* Write a function to convert degrees K to degrees C  (Part b, (iv))
+Call it from the main() function above */
+
+float kelvin_to_c(float k) {
+	return k-273.15;
+}
+
+
+float resistance_to_temperature(float r) {
+	// Define Thermistor constants
+	float t_0 = 298.15;
+	float r_0 = 10000;
+	float b = 3975;
+
+	return 1 / ( (1/t_0) + (1/b)*log(r/r_0));
+}
+
+
+float voltage_to_resistance(float v) {
+	return 1000*((10.0*3.3/v)-10.0);
+}
+
+
+float voltage_to_erc(float v) {
+	return (v-0.35)/54.4;
+}
+
+
+/* returns EMF in millivolts */
+float NISTdegCtoMilliVoltsKtype(float tempDegC)
+{
+    int i;
+    float milliVolts = 0;
+    if(tempDegC >= -170 && tempDegC < 0)
+    {
+        const float coeffs[11] =
+        {
+            0.000000000000E+00,
+            0.394501280250E-01,
+            0.236223735980E-04,
+            -0.328589067840E-06,
+            -0.499048287770E-08,
+            -0.675090591730E-10,
+            -0.574103274280E-12,
+            -0.310888728940E-14,
+            -0.104516093650E-16,
+            -0.198892668780E-19,
+            -0.163226974860E-22
+        };
+        for (i=0; i<=10; i++)
+        {
+            milliVolts += coeffs[i] * pow(tempDegC,i);
+        }
+    }
+    else if(tempDegC >= 0 && tempDegC <= 1372)
+    {
+        const float coeffs[10] =
+        {
+            -0.176004136860E-01,
+            0.389212049750E-01,
+            0.185587700320E-04,
+            -0.994575928740E-07,
+            0.318409457190E-09,
+            -0.560728448890E-12,
+            0.560750590590E-15,
+            -0.320207200030E-18,
+            0.971511471520E-22,
+            -0.121047212750E-25
+        };
+        const float a0 =  0.118597600000E+00;
+        const float a1 = -0.118343200000E-03;
+        const float a2 =  0.126968600000E+03;
+
+        for (i=0; i<=9; i++)
+        {
+            milliVolts += coeffs[i] * pow(tempDegC,i);
+        }
+
+        milliVolts += a0*exp(a1*(tempDegC - a2)*(tempDegC - a2));
+    }
+    else
+    {
+        milliVolts = 99E99;
+    }
+    return milliVolts;
+}
+
+// returns temperature in deg C.
+float NISTmilliVoltsToDegCKtype(float tcEMFmV)
+{
+
+        int i, j;
+        float tempDegC = 0;
+        const float coeffs[11][3] =
+        {
+          {0.0000000E+00,  0.000000E+00, -1.318058E+02},
+         {2.5173462E+01,  2.508355E+01,  4.830222E+01},
+         {-1.1662878E+00,  7.860106E-02, -1.646031E+00},
+         {-1.0833638E+00, -2.503131E-01,  5.464731E-02},
+         {-8.9773540E-01,  8.315270E-02, -9.650715E-04},
+         {-3.7342377E-01, -1.228034E-02,  8.802193E-06},
+         {-8.6632643E-02,  9.804036E-04, -3.110810E-08},
+         {-1.0450598E-02, -4.413030E-05,  0.000000E+00},
+         {-5.1920577E-04,  1.057734E-06,  0.000000E+00},
+         {0.0000000E+00, -1.052755E-08,  0.000000E+00}
+       };
+       if(tcEMFmV >=-5.891 && tcEMFmV <=0 )
+       {
+           j=0;
+       }
+       else if (tcEMFmV > 0 && tcEMFmV <=20.644  )
+       {
+           j=1;
+       }
+       else if (tcEMFmV > 20.644 && tcEMFmV <=54.886  )
+       {
+           j=2;
+       }
+       else
+       {
+           return 99E9;
+       }
+
+       for (i=0; i<=9; i++)
+        {
+            tempDegC += coeffs[i][j] * pow(tcEMFmV,i);
+        }
+    return tempDegC;
+}