In [1]:
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
#Импортируем функцию calc_dryair_dens
from calcdryairdens import calc_dryair_dens
In [2]:
ts = np.linspace(-100, 1000, 221)
ps = np.array([0.101325, 0.5, 1, 2, 3, 4, 5, 10, 15, 20]) * 1e6
denss = np.zeros((ts.shape[0], ps.shape[0]), dtype=float)
zs = np.zeros((ts.shape[0], ps.shape[0]), dtype=float)
In [3]:
for i in range(ts.shape[0]):
for j in range(ps.shape[0]):
denss[i,j], zs[i,j] = calc_dryair_dens(ts[i], ps[j])
In [4]:
plt.rcParams['figure.figsize'] = [15, 7]
plt.rcParams.update({'font.size': 16})
plt.title('Зависимость коэффициента сжимаемости сухого воздуха $z$ \nот температуры (C) \
при различных давлениях (МПа)')
plt.xlabel('$t, C$'); plt.ylabel('$z$')
plt.plot(ts, zs)
plt.legend(ps/1e6)
plt.xlim (-100,1000);
#plt.ylim (0.95,1.05)
plt.grid()
In [5]:
plt.title('Зависимость коэффициента сжимаемости сухого воздуха $z$ \nот температуры (C) \
при различных давлениях (МПа)')
plt.xlabel('$t, C$'); plt.ylabel('$z$')
plt.plot(ts, zs)
plt.legend(ps/1e6)
plt.xlim (-100,200);
plt.ylim (0.95,1.05)
plt.grid()
In [6]:
table = pd.DataFrame(denss, columns = ps/1e6, index = ts)
In [7]:
table
Out[7]:
In [8]:
#Сохранение результатов расчёта в csv-файл:
#table.to_csv('Плотность_сухого_воздуха.csv', sep = ';')
In [9]:
#Сохранение результатов расчёта в html-файл:
#table.to_html('Плотность_сухого_воздуха.html')