En una entrada anterior de este blog, publique el proceso de control estadístico en una maquina que fabricaba resistencias de 1000, con una tolerancias de +/- 100 , y en el programa se mostraba la captura de datos, como estos datos se metían a una matriz para trabajarlos y hacer cálculos con ellos, para finalmente imprimir la gráfica del control estadístico la cual es la parte visual de este programa
Pero en el anterior programa había un inconveniente, este era que al momento de ingresar los datos no se guardaban en ningún sitio por ende los datos que se ingresaban ya después de cerrar la corrida se perdían para siempre y si se quería volver a revisar la misma corrida después de haberla finalizado tenias que volver a ingresar los datos.
Ahora les mostrare el programa ya corregido, donde te da la opción de usar datos nuevos, ingresandolos y guardándolos en un bloc de notas, o el poder extraer datos de un bloc de notas, siempre y cuando este bloc de notas este guardado en la misma carpeta donde este guardado el programa del control estadistico
Explicación del programa.
El programa en esencia es el mismo que el anterior mostrado, solo que se le agrego el método "open-close" que sirve para abrir archivos de texto, donde existen varias notaciones como "r" donde da permiso de lectura , "w" el cual da permiso de sobre escritura en un determinado archivo.
Al inicio del programa te pregunta si deseas utilizar datos nuevos y guardarlos o utilizar datos de la anterior corrida, cuando utilizas datos nuevos estos se sobrescriben en un archivo de texto predefinido y si son datos anteriores se abre el archivo predestinado como el default (ejemplo de la anterior entrada ) o si ya se guardaron datos se abren los datos de la ultima corrida.
PROGRAMA
Cada linea del programa esta comentada para su entendimiento
import numpy as np
from pylab import *
print("Control Estadistico")
print("Puedes usar datos nuevos o usar los datos de la ultima entrada")
resp=input("¿Quieres poner datos nuevos? ")
respuesta="si"
respuesta1="Si"
if (resp in respuesta) or (resp in respuesta1 ): #en caso de usar datos nuevos
print("LLENADO DE DATOS")
archivo=open("matriz.txt","w")
renglon=int(input("Dame el numero de muestras: "))
columna=int(input("Dame el numero de lecturas: "))
#guarda los datos en un archivo
for i in range(renglon):
print("MUESTRA", i+1)
for j in range(columna):
val=int(input("Lectura "+str(j+1)+" "))
archivo.write(str(val))
archivo.write(" ")
archivo.write("\n")
archivo.close()
#guarda los datos en la matriz
archivo=open("matriz.txt","r")
print()
listass=[]
for i in range(renglon):
lisp=archivo.readline() #en esta variable guardo las lineas del documento
llena=lisp.split() #con esta funcion separo los datos de cada renglon y los guardo en una lista se podria decir
listass.append(llena)#por ultimo ya separados los datos y en una lista, guardo la lista en listass
archivo.close()
#llena matriz
A=np.zeros((renglon,columna))
for i in range(renglon):
for j in range(columna):
A[i][j]= int(listass[i][j])
else: #en caso de datos anteriores
print("USANDO LOS DATOS ANTERIORES")
archivo=open("matriz.txt","r")
#dice numero de renglones para la matriz
renglon=0
for n in archivo: #con este for cuento los renglones del archivo
renglon=renglon+1
#print(renglon)
#print()
archivo.close()
archivo=open("matriz.txt","r")
listass=[] #declaro lista vacia
for i in range(renglon):
lisp=archivo.readline() #en esta variable guardo las lineas del documento
llena=lisp.split() #con esta funcion separo los datos de cada renglon y los guardo en una lista se podria decir
listass.append(llena)#por ultimo lla separados los datos y en una lista, guardo la lista en listass
#print(listass)
#print()
#print(listass[0])
#print()
#print(listass[0][0])
#print()
#dice la cantidad de columnas para la matriz
uno=listass[0]
columna=len(uno)
#print(columna)
archivo.close()
#llena matriz
A=np.zeros((renglon,columna))
for i in range(renglon):
for j in range(columna):
A[i][j]= int(listass[i][j])
#imprime la matriz
print("MUESTRAS")
for i in range(renglon):
print("|",i+1,"\t|",end=" ")
for j in range(columna):
print(A[i][j], end=" ")
print()
print()
#Calis
#CALCULA LOS PROMEDIOS
prom=[]
for i in range(renglon):
p=np.sum(A[i])/columna
prom.append(p)
print("Promedios",prom)
print()
#Calcula el promedio total de promedios
suma=0
i=0
for elemento in prom:
suma+=elemento
i+=1
promedio=suma/i
print("Promedio de promedios",promedio)
print()
#CALCULO DE LOS RANGOS
lisrang=[]
for i in range(renglon):
M=np.max(A[i])
m=np.min(A[i])
R=M-m
lisrang.append(R)
print("Rangos", lisrang)
print()
#PROMEDIO DE RANGOS
suma1=0
z=0
for elemento in lisrang:
suma1+=elemento
z+=1
promedio1=suma1/z
print("Promedio de rangos",promedio1)
print()
print()
#CALCULA LINEA CENTRAL EN PROMEDIOS +
lineaCentralP=[]
Lc=0
while (Lc<renglon):
lineaCentralP.append(promedio)
Lc=Lc+1
print("linea central en promedios",lineaCentralP[0])
print()
#Tabla para sacar los A2 D3 D4
A2=[0,0,1.880,1.023,0.729,0.577,0.486,0.429,0.373,0.377,0.308,0.285,0.266,0.249,0.235,0.223]
D3=[0,0,0,0,0,0,0,0.076,0.136,0.184,0.223,0.256,0.284,0.308,0.329,0.348]
D4=[0,0,3.268,2.574,2.282,2.114,2.004,1.924,1.864,1.816,1.777,1.744,1.717,1.692,1.671,1.652]
#LIMITE CONTROL INFERIOR EN PROMEDIOS
lci=promedio-(A2[columna]*promedio1)
LCI=0
lineainferiorP=[]
while(LCI<renglon):
lineainferiorP.append(lci)
LCI=LCI+1
print("linea inferior en promedios",lineainferiorP[0])
print()
#LIMITE CONTROL SUPERIOR EN PROMEDIOS
lcs=promedio+(A2[columna]*promedio1)
LCS=0
LSP=[]
while(LCS<renglon):
LSP.append(lcs)
LCS=LCS+1
print("linea superior en promedios",LSP[0])
print()
print()
print()
#print("D3",D3[lec])
#print("D4",D4[lec])
#LINEA CENTRAL EN RANGOS
LCR=[]
y=0
while(y<renglon):
LCR.append(promedio1)
y=y+1
print("linea central en rangos",LCR[0])
print()
#LIMITE CONTROL INFERIOR EN RANGOS
lcir=0*promedio1
v=0
LCIR=[]
while(v<renglon):
LCIR.append(lcir)
v=v+1
print("linea inferior en rangos",LCIR[0])
print()
#LIMITE CONTROL SUPERIOR EN RANGOS
lcsr=2.282*promedio1
w=0
LSR=[]
while(w<renglon):
LSR.append(lcsr)
w=w+1
print("linea superior en rangos",LSR[0])
print()
#GRAFICA PARA PROMEDIOS
plt.plot(lineaCentralP, label="PROM", marker='s') #linea centra
plt.plot(lineainferiorP, label="LCI", marker='s') #linea inferior
plt.plot(LSP, label="LCS", marker='D') #linea superior
plt.plot(prom, label="Promedios",marker='s')#promedios
plt.grid(True)
plt.legend()
plt.title("Control en promedios")
plt.xlabel("Muestras")
plt.ylabel("Promedios")
plt.show()
#GRAFICA PARA RANGOS
plt.plot(LCR, label="prom", marker='s')
plt.plot(LCIR, label="LCI", marker='s')
plt.plot(LSR, label="LCS", marker='s')
plt.plot(lisrang, label="Rangos", marker='D')
plt.grid(True)
plt.legend()
plt.title("Control en rangos")
plt.xlabel("Muestras")
plt.ylabel("Rangos")
plt.show()
EJEMPLO USANDO DATOS DE LA ANTERIOR CORRIDA
En la siguiente imagen se muestra los datos guardados en un bloc de notas
Como podrás ver son los datos del ejemplo de la entrada anterior
CORRIDA
En la anterior imagen se muestra la corrida del programa en el caso de querer usar datos de un archivo, entonces el programa toma los datos del archivo y los utiliza para hacer el control estadístico, y en la corrida ya no pide los datos porque los esta sacando del archivo de texto
GRAFICAS
EJEMPLO USANDO DATOS NUEVOS
CORRIDA
En esta corrida de ve que pregunta si deseaba usar datos nuevos y al contestar si, se empiezan a llenar los datos en este caso fueron datos aleatorios.
Después de haber llenado los datos despliega la tabla y los cálculos
GRÁFICAS
DATOS GUARDADOS EN EL ARCHIVO
Estos datos ya quedaron guardados, listos para volverse a usar o para imprimirse y revisarlos en el control estadístico.
NOTA: LOS ARCHIVOS DE TEXTO DEBEN DE ESTAR EN LA MISMA CARPETA DONDE SE GUARDO EL PROGRAMA ADEMAS, SI SE DESEA USAR OTRO NOMBRE PARA EL ARCHIVO EL CÓDIGO FUENTE DEBE DE MODIFICARSE PONIENDO EL NOMBRE NUEVO.
Por ultimo me gustaría agregar que la programación es de estar depurando y quitando errores constantemente, ya que sobre la marcha saldrán errores o surgirán nuevas necesidades y por ende debemos de estar modificando los programas constantemente para que no sean obsoletos.
Comentarios
Publicar un comentario