Producción del día 20 de Julio de 2015
Los datos del experimento:
Se desea comprobar si el filamento que podemos llegar a extruir con el sistema de la tractora puede llegar a ser bueno como para regularlo.
In [2]:
%pylab inline
#Importamos las librerías utilizadas
import numpy as np
import pandas as pd
import seaborn as sns
In [3]:
#Mostramos las versiones usadas de cada librerías
print ("Numpy v{}".format(np.__version__))
print ("Pandas v{}".format(pd.__version__))
print ("Seaborn v{}".format(sns.__version__))
In [4]:
#Abrimos el fichero csv con los datos de la muestra
datos = pd.read_csv('datos.csv')
In [5]:
#Almacenamos en una lista las columnas del fichero con las que vamos a trabajar
columns = ['Diametro X','VELOCIDAD']
In [6]:
#Mostramos un resumen de los datos obtenidoss
datos[columns].describe()
#datos.describe().loc['mean',['Diametro X [mm]', 'Diametro Y [mm]']]
Out[6]:
Representamos ambos diámetro y la velocidad de la tractora en la misma gráfica
In [26]:
#datos.ix[:, "Diametro X":"Diametro Y"].plot(secondary_y=['VELOCIDAD'],figsize=(16,10),ylim=(0.5,3)).hlines([1.85,1.65],0,3500,colors='r')
#graf = datos[columns].plot(secondary_y=['VELOCIDAD'],ylim=(1,2.5),figsize=(10,5),title='Relación entre la velocidad de tracción y el diámetro')
#graf.axhspan(1.65,1.85, alpha=0.2)
#datos['RPM TRAC'].plot(secondary_y='RPM TRAC')
fig, ax1 = plt.subplots()
ax1.plot(datos['VELOCIDAD'], 'b-')
ax1.set_xlabel('Tiempo (s)')
ax1.set_ylabel('Velocidad de tracción (RPM)', color='b')
ax1.set_ylim(0.5,3.1)
ax1.set_xlim(-1,200)
ax2 = ax1.twinx()
ax2.plot(datos['Diametro X'], 'r-')
ax2.set_ylabel('Diámetro (mm)', color='r')
ax2.set_ylim(0.5,3.1)
ax2.set_xlim(-1,200)
plt.axhspan(1.65,1.85, alpha=0.2)
plt.figure(figsize=(10,5))
plt.show()
In [27]:
graf = datos.ix[:, "Diametro X":"Diametro Y"].boxplot(return_type='axes')
graf.axhspan(1.65,1.85, alpha=0.2)
Out[27]:
Con esta segunda aproximación se ha conseguido estabilizar los datos. Se va a tratar de bajar ese porcentaje. Como cuarta aproximación, vamos a modificar las velocidades de tracción. El rango de velocidades propuesto es de 1.5 a 5.3, manteniendo los incrementos del sistema experto como en el actual ensayo.
Comparativa de Diametro X frente a Diametro Y para ver el ratio del filamento
In [9]:
plt.scatter(x=datos['Diametro X'], y=datos['Diametro Y'], marker='.')
Out[9]:
In [10]:
datos_filtrados = datos[(datos['Diametro X'] >= 0.9) & (datos['Diametro Y'] >= 0.9)]
In [11]:
#datos_filtrados.ix[:, "Diametro X":"Diametro Y"].boxplot(return_type='axes')
In [12]:
plt.scatter(x=datos_filtrados['Diametro X'], y=datos_filtrados['Diametro Y'], marker='.')
Out[12]:
In [13]:
ratio = datos_filtrados['Diametro X']/datos_filtrados['Diametro Y']
ratio.describe()
Out[13]:
In [14]:
rolling_mean = pd.rolling_mean(ratio, 50)
rolling_std = pd.rolling_std(ratio, 50)
rolling_mean.plot(figsize=(12,6))
# plt.fill_between(ratio, y1=rolling_mean+rolling_std, y2=rolling_mean-rolling_std, alpha=0.5)
ratio.plot(figsize=(12,6), alpha=0.6, ylim=(0.5,1.5))
Out[14]: