Étape 4: Écriture du code
Comme nous l’avons mentionné plus tôt, nous avons opté pour écrire une interface utilisateur graphique personnalisée pour le spiromètre. Parce que tous ceux qui lisent que ceci n’aient pas accès à MATLAB, nous simplement fournir les équations générales et laisser de côté les détails désordonnés de l’interface utilisateur dans ce guide. Nous énoncerons les bases du code en quelques étapes.
Dans une boucle qui parcourt le temps d’échantillon, où ii est l’itération :
- Lire tension A0 et appliquez un filtre moyens mobile pour réduire le bruit
- Voltage(II) = * lire la tension de la broche (application dépendra que vous utilisiez MATLAB ou Arduino *
- À l’aide de l’équation de la fiche technique du transducteur de pression, de convertir la tension en une différence de pression (psi).
- dP_psi(ii) = (14.6959 / 2) * (voltage(ii) - V_reference)
- Convertir la pression en Pascals pour des calculs plus faciles
- dP_Pa(ii) = dP_psi (ii) * 6894.76
- À l’aide de l’équation pour le débit volumétrique à travers un tube de Venturi, (IDs est le diamètre intérieur de la petite portion du tube, et bêta est le rapport entre les deux diamètres. Rho est la masse volumique de l’air à une température donnée)
Lucratives = A2*sqrt((2*(P1-P2))/(ρ*(1-(A2/A1)^2)))
Reportez-vous à l’étape 2 pour plus d’informations sur cette équation
- Si vous voulez calculer que le volume total d’air expiré dans votre gamme de l’échantillon, nous pouvons utiliser une simple intégration numérique :
- VolTotal += lucratives * dt
- Où dt est la fréquence d’échantillonnage