PROGRAMA PARA LEER ARCHIVOS
PROGRAMA DE RESISTENCIA
este programa usa el método de archivo para y así poder sacar los valores de las medidas de varia resistencias de un bloc de notas para poder aplicarlas en un problema estadístico. y también da la opción de que el usuario introduzca sus propias medidas en el bloc de notas si el así lo requiere.
from pylab import *
import numpy as np
respuesta=input(str("si quiere usar los datos guardados introduzca OKEY si quiere sus datos introduzca MIO\n"))
if(respuesta == "MIO"):
archivo=open("PROYECTO.txt","w")
muestra=0
final=[]
pequeña=[]
otra=[]
muestra=0
c=int(input("Introduzca el numero de muestras\n"))
b=int(input("introduzca el numero de lecturas\n"))
for j in range(c):
muestra=muestra+1
print("Muestra ",muestra)
suma=0
for i in range(b):
a=float(input(""))
a=str(a)
archivo.write(a)
archivo.write(" ")
archivo.write("\n")
archivo.close()
archivo=open("PROYECTO.txt","r+")
for linea in archivo.readlines():
elen=linea.split()
otra.append(elen)
for i in range (c):
for j in range(b):
pequeña.append(float(otra[i][j]))
final.append(pequeña)
pequeña=[]
promedio2=0
rango2=0
for i in range(len(final)):
for j in range(b):
suma=suma+final[i][j]
promedio=suma/b
promedio2=promedio2+promedio
final[i].append(promedio)
rango=max(final[i])-min(final[i])
final[i].append(rango)
rango2=rango2+rango
suma=0
rangof=rango2/c
promediof=promedio2/c
print()
print(" Muestras Lecturas(OHMS) Promedio Rango" )
for i in range(len(final)):
print(" \t",i+1," \t",end="")
for j in range((b+2)):
print(" \t%.1f"%(final[i][j]),end="")
print()
print(" PROM \t%.1f"%(promediof)," \t%.1f"%(rangof))
LCI=promediof-(.729*rangof)
LCS=promediof+(.729*rangof)
rango3=rangof*0
rango4=rangof*2.282
print(" LCS \t%.1f"%(LCS)," \t%.1f"%(rango4))
print(" LCI \t%.1f"%(LCI)," \t",rango3)
x=[]
PROM=[]
graficalCS=[]
graficalCI=[]
for i in range (len(final)):
graficalCI.append(LCI)
graficalCS.append(LCS)
x.append(promediof)
PROM.append(final[i][b])
plt.title("RESISTENCIAS")
plt.ylabel("OHMS")
plt.xlabel("N° Muestra ")
plt.plot(PROM,label="PROM", marker="s", linestyle="-", color="b")
plt.plot(x,label="x",marker="D", linestyle="-",color="g")
plt.plot(graficalCS,label="LCS", marker="^", linestyle="-", color="r")
plt.plot(graficalCI,label="LCI", marker="x", linestyle="-", color="c")
plt.legend(loc="lower left")
plt.show()
graficarango=[]
PROMrango=[]
graficalCS2=[]
graficalCI2=[]
for i in range (len(final)):
graficalCI2.append(rango4)
graficalCS2.append(rango3)
graficarango.append(rangof)
PROMrango.append(final[i][b+1])
plt.title("RANGOS")
plt.ylabel("OHMS")
plt.xlabel("N° Muestra ")
plt.plot(PROMrango,label="PROM", marker="s", linestyle="-", color="b")
plt.plot(graficarango,label="x",marker="D", linestyle="-",color="g")
plt.plot(graficalCS2,label="LCS", marker="^", linestyle="-", color="r")
plt.plot(graficalCI2,label="LCI", marker="x", linestyle="-", color="c")
plt.legend(loc="lower left")
plt.show()
archivo.close()
if(respuesta=="OKEY"):
archivo=open("PROYECTOLIMPIO.txt")
muestra=0
final=[]
pequeña=[]
otra=[]
muestra=0
c=12
b=4
archivo=open("PROYECTOLIMPIO.txt","r+")
for linea in archivo.readlines():
elen=linea.split()
otra.append(elen)
for i in range (c):
for j in range(b):
pequeña.append(float(otra[i][j]))
final.append(pequeña)
pequeña=[]
promedio2=0
rango2=0
suma=0
for i in range(len(final)):
for j in range(b):
suma=suma+final[i][j]
promedio=suma/b
promedio2=promedio2+promedio
final[i].append(promedio)
rango=max(final[i])-min(final[i])
final[i].append(rango)
rango2=rango2+rango
suma=0
rangof=rango2/c
promediof=promedio2/c
print()
print(" Muestras Lecturas(OHMS) Promedio Rango" )
for i in range(len(final)):
print(" \t",i+1," \t",end="")
for j in range((b+2)):
print(" \t%.1f"%(final[i][j]),end="")
print()
print(" PROM \t%.1f"%(promediof)," \t%.1f"%(rangof))
LCI=promediof-(.729*rangof)
LCS=promediof+(.729*rangof)
rango3=rangof*0
rango4=rangof*2.282
print(" LCS \t%.1f"%(LCS)," \t%.1f"%(rango4))
print(" LCI \t%.1f"%(LCI)," \t",rango3)
x=[]
PROM=[]
graficalCS=[]
graficalCI=[]
for i in range (len(final)):
graficalCI.append(LCI)
graficalCS.append(LCS)
x.append(promediof)
PROM.append(final[i][b])
plt.title("RESISTENCIAS")
plt.ylabel("OHMS")
plt.xlabel("N° Muestra ")
plt.plot(PROM,label="PROM", marker="s", linestyle="-", color="b")
plt.plot(x,label="x",marker="D", linestyle="-",color="g")
plt.plot(graficalCS,label="LCS", marker="^", linestyle="-", color="r")
plt.plot(graficalCI,label="LCI", marker="x", linestyle="-", color="c")
plt.legend(loc="lower left")
plt.show()
graficarango=[]
PROMrango=[]
graficalCS2=[]
graficalCI2=[]
for i in range (len(final)):
graficalCI2.append(rango4)
graficalCS2.append(rango3)
graficarango.append(rangof)
PROMrango.append(final[i][b+1])
plt.title("RANGOS")
plt.ylabel("OHMS")
plt.xlabel("N° Muestra ")
plt.plot(PROMrango,label="PROM", marker="s", linestyle="-", color="b")
plt.plot(graficarango,label="x",marker="D", linestyle="-",color="g")
plt.plot(graficalCS2,label="LCS", marker="^", linestyle="-", color="r")
plt.plot(graficalCI2,label="LCI", marker="x", linestyle="-", color="c")
plt.legend(loc="lower left")
plt.show()
Comentarios
Publicar un comentario