import numpy as np
import matplotlib.pyplot as plt

import soundFunctions as sound

plotsPath = '14-sound/plots/'


########################################
#   Data processing
########################################

temperature = 22.6 # C
absoluteHumidity = 0.006 # kg/m^3 at realtiveHumidity = 31.5 and temperature 22.6 C. Check at https://planetcalc.com/2167/

co2Max = 5 # %

co2x, soundSpeed = sound.speedOfSound(temperature, absoluteHumidity, co2Max)

linearFit = np.polyfit(soundSpeed, co2x, 1)

dtAir = 3.36 / 1000 # sec
dtBreath = 3.40 / 1000 # sec

L = 1.158 # m

uAir = L / dtAir
uBreath = L / dtBreath

co2xAir = np.polyval(linearFit, uAir)
co2xBreath = np.polyval(linearFit, uBreath)

# Plot

plt.figure('Скорость звука', figsize=(8, 5), dpi=200)

plt.title('Зависимость скорости звука\nот концентрации углекислого газа')
plt.xlabel('Концентрация CO$_{2}$ [%]')
plt.ylabel('Скорость звука [м/c]')

plt.plot(co2x, soundSpeed, label='Аналитическая зависимость')
plt.plot(co2xAir, uAir, '*', label='Значения в воздухе: {:.1f} [м/с], {:.1f} [%]'.format(uAir, co2xAir))
plt.plot(co2xBreath, uBreath, '*', label='Значения в выдохе: {:.1f} [м/с], {:.1f} [%]'.format(uBreath, co2xBreath))
plt.legend()

plt.grid(which='major')
plt.grid(which='minor', linestyle='--')
plt.minorticks_on()

plt.savefig(plotsPath + 'speed-of-sound.png')

# plt.show()
import numpy as np
import matplotlib.pyplot as plt

import soundFunctions as sound

plotsPath = '14-sound/plots/'


########################################
#   Data processing
########################################

temperature = 22.6 # C
absoluteHumidity = 0.006 # kg/m^3 at realtiveHumidity = 31.5 and temperature 22.6 C. Check at https://planetcalc.com/2167/

co2Max = 5 # %

co2x, soundSpeed = sound.speedOfSound(temperature, absoluteHumidity, co2Max)

linearFit = np.polyfit(soundSpeed, co2x, 1)

dtAir = 3.36 / 1000 # sec
dtBreath = 3.40 / 1000 # sec

L = 1.158 # m

uAir = L / dtAir
uBreath = L / dtBreath

co2xAir = np.polyval(linearFit, uAir)
co2xBreath = np.polyval(linearFit, uBreath)

# Plot

plt.figure('Скорость звука', figsize=(8, 5), dpi=200)

plt.title('Зависимость скорости звука\nот концентрации углекислого газа')
plt.xlabel('Концентрация CO$_{2}$ [%]')
plt.ylabel('Скорость звука [м/c]')

plt.plot(co2x, soundSpeed, label='Аналитическая зависимость')
plt.plot(co2xAir, uAir, '*', label='Значения в воздухе: {:.1f} [м/с], {:.1f} [%]'.format(uAir, co2xAir))
plt.plot(co2xBreath, uBreath, '*', label='Значения в выдохе: {:.1f} [м/с], {:.1f} [%]'.format(uBreath, co2xBreath))
plt.legend()

plt.grid(which='major')
plt.grid(which='minor', linestyle='--')
plt.minorticks_on()

plt.savefig(plotsPath + 'speed-of-sound.png')

# plt.show()
######################################## #   Data processing ######################################## # C # kg/m^3 at realtiveHumidity = 31.5 and temperature 22.6 C. Check at https://planetcalc.com/2167/ # % # sec # sec # m # Plot # plt.show()