Файл: Плотность воздуха против температуры.svg

import  csv ,  numpy , scipy , scipy.optimize import  matplotlib.pyplot  as  plt x  =  [] y  =  [] с  open ( "data.csv" ,  "r" )  as  f :  для  строки  в  csv . DictReader ( f ,  fieldnames = [ 'x' , 'y' ],  delimiter = ';' ):  x . добавить (int ( строка [ 'x' ]))  y . append ( float ( строка [ 'y' ])) x  =  numpy . массив ( x ) y  =  numpy . массив ( y )intx  =  numpy . linspace ( min ( x ),  max ( x ),  100 ) P1  =  numpy . poly1d ( NumPy . polyfit ( х ,  у ,  2 ))P2_data  =  scipy . оптимизировать . curve_fit ( lambda  t , a , b :  a * numpy . exp ( b * t ),  x ,  y ,  ) [ 0 ] P2  =  lambda  x  :  P2_data [ 0 ] * numpy . exp ( P2_data [ 1 ] * x )plt . scatter ( x , y ,  label = "Плотность воздуха при давлении 1 атм" ) plt . plot ( intx ,  P1 ( intx ),  label  =  "Полиномиальный тренд \ n $ x \ mapsto 10 ^ {{- 5}} x ^ 2 {} x + {} $" . format ( * numpy . round ( P1 , 4 ) [ 1 :])) plt . сюжет ( intx , P2 ( intx ),  label  =  "Экспоненциальный тренд \ n $ x \ mapsto \ exp ( {} x) {} $" . format ( *  numpy . round ( P2_data , 4 ))) plt . легенда () plt . xlabel ( "Температура [° C]" ) plt . ylabel ( «Плотность [кг / м³]» ) plt . title ( «Плотность воздуха в зависимости от температуры» ) plt .savefig ( "output.svg" , format = "svg" , transparent = True )