##// END OF EJS Templates
Monitoreo-de-Potencia
RSS
Se agrego el threshold...
JesusTapia-dev -
r20:96b1842fbedf schottkyExp
parent child
Show More
Show file before | Show file after | Hide comments
@@ -1,2 +1,2
1 IPP(km) DutyCicle(%) Wu Wd Nd Nu Eu Ed Sd Su
2 1500 0.1 250 0 0 0 0 0 200 200 No newline at end of file
1 IPP(km) DutyCicle(%) threshold(%) Wu Wd Nd Nu Eu Ed Sd Su
2 1500 0.1 10 250 0 0 0 0 0 200 200 No newline at end of file
Show file before | Show file after | Hide comments
@@ -1,125 +1,126
1 1 import paho.mqtt.client as mqtt
2 2 import paho.mqtt.publish as publish
3 3 from datetime import datetime, timezone, timedelta
4 4 import numpy as np
5 5 import json
6 6 import time
7 7 subscribe_topic="tesis/analogRaw"
8 8 publish_topic="atrad/test3"#"tesis/potencia"
9 9 broker_address = "10.10.10.102"#10.10.10.102"#"192.168.43.149"
10 10 analogRawMatriz=[[],[],[],[],[],[],[],[]]#guardamos los valores del
11 11 numMaxdatos=1000 #definimos el numero máximo de datos a promediar
12 12 contador = 0
13 13 interval_start_time = time.time()
14 14 interval_duration=30 #duracion del intervalo (en segundos) para el envio de datos
15 15 def potenciaAncho10us(AnalogRaw):
16 16 potencia=[0,0,0,0,0,0,0,0]
17 17 m=1.5476
18 18 b=-91.898
19 19 for i in range(8):
20 20 potencia[i]=m*AnalogRaw[i]+b
21 21 return potencia
22 22
23 23 def potenciaAncho20us(AnalogRaw):
24 24 potencia=[0,0,0,0,0,0,0,0]
25 25 m=0.6233
26 26 b=62.891
27 27 for i in range(8):
28 28 potencia[i]=m*AnalogRaw[i]+b
29 29 return potencia
30 30 #la potencia default considera en general valores mayores a 100 us
31 31 def potenciaDefault(AnalogRaw):
32 32 potencia=[0,0,0,0,0,0,0,0]
33 33 m=0.6233
34 34 b=62.891
35 35 for i in range(8):
36 36 potencia[i]=m*AnalogRaw[i]+b
37 37 return potencia
38 38
39 39 def leer_datos_desde_txt(archivo):
40 40 try:
41 41 with open(archivo, 'r') as file:
42 42 lineas = file.readlines()
43 43 # Obtiene la segunda línea (índice 1)
44 44 segunda_linea = lineas[1]
45 45 datos = segunda_linea.strip().split()
46 46 datos_numericos = [float(dato) for dato in datos]
47 47 return datos_numericos
48 48 except FileNotFoundError:
49 49 print(f"El archivo '{archivo}' no fue encontrado.")
50 50 except Exception as e:
51 51 print(f"Error al leer el archivo: {e}")
52 52
53 53 def procesamiento_data(analogRawMatriz):
54 54 if analogRawMatriz[0]:
55 55 average_raw=[0,0,0,0,0,0,0,0]
56 56 for i in range(8):
57 57 average_raw[i] = sum(analogRawMatriz[i]) / len(analogRawMatriz[i])
58 58 if(ancho<15):
59 59 potencia=potenciaAncho10us(average_raw)
60 60 if ancho>15 and ancho<25:
61 61 potencia=potenciaAncho20us(average_raw)
62 62 timestamp_actual = int(time.time())
63 63 fecha_utc=datetime.utcfromtimestamp(timestamp_actual)
64 64 zona_horaria_utc_menos_5 = timezone(timedelta(hours=-5))
65 65 fecha_utc_menos_5 = fecha_utc.replace(tzinfo=timezone.utc).astimezone(zona_horaria_utc_menos_5)
66 66 fecha_legible=fecha_utc_menos_5.strftime('%d-%m-%Y %H:%M:%S')
67 67 estado=[1 if x > 0 else 0 for x in potenciaNominal]
68 processed_data = {"Ancho_us":IPP,"pow": potencia, "time": fecha_legible,"potenciaNominal":potenciaNominal,"status":estado}
68 processed_data = {"Ancho_us":IPP,"pow": potencia, "time": fecha_legible,"potenciaNominal":potenciaNominal,"status":estado,"threshold":threshold}
69 69 client.publish(publish_topic, json.dumps(processed_data))
70 70 print("---------------------------------------")
71 71 print(average_raw)
72 72 print(processed_data)
73 73 print(type(processed_data["pow"]))
74 74
75 75 def on_connect(client, userdata, flags, rc):
76 76 if rc == 0:
77 77 print("Conexión exitosa con el broker")
78 78 client.subscribe(subscribe_topic)
79 79 else:
80 80 print("Error de conexión. Código de retorno =", rc)
81 81 # Callback cuando se recibe un mensaje en el tópico suscrito
82 82 def on_message(client, userdata, msg):
83 83 global analogRawMatriz
84 84 global contador
85 85 mensaje = msg.payload.decode()
86 86 lista=json.loads(mensaje)
87 87 if contador<numMaxdatos:
88 88 for i in range(8):
89 89 analogRawMatriz[i].append(lista[i])
90 90 contador=contador+1
91 91 #Leo los datos de la configuracion y halla el ancho
92 92 archivo_txt = 'commandPotencia.txt'
93 93 datos_numericos = leer_datos_desde_txt(archivo_txt)
94 if len(datos_numericos)!=10:
94 if len(datos_numericos)!=11:
95 95 print("Hay más(o menos) valores de los que debería")
96 96 else:
97 97 IPP_km=datos_numericos[0]
98 98 IPP=IPP_km*1/150#IPP en ms
99 99 Dutty=datos_numericos[1]
100 threshold=datos_numericos[2]
100 101 ancho=IPP*Dutty*pow(10,3)/100
101 potenciaNominal=datos_numericos[2:]
102 potenciaNominal=datos_numericos[3:]
102 103 # Configurar el cliente MQTT
103 104 client = mqtt.Client()
104 105 # Configurar los callbacks
105 106 client.on_connect = on_connect
106 107 client.on_message = on_message
107 108 client.connect(broker_address, port=1883, keepalive=60)
108 109
109 110 client.loop_start()
110 111 try:
111 112 while True:
112 113 # Verifica si ha pasado el intervalo de tiempo
113 114 current_time = time.time()
114 115 elapsed_time = current_time - interval_start_time
115 116 if elapsed_time >= interval_duration:
116 117 procesamiento_data(analogRawMatriz) # Calcula la potencia y envía los datos
117 118 interval_start_time = current_time # Reinicia el temporizador del intervalo
118 119 analogRawMatriz=[[],[],[],[],[],[],[],[]]
119 120 time.sleep(1)
120 121
121 122 except KeyboardInterrupt:
122 123 print("Programa detenido por el usuario.")
123 124 client.disconnect()
124 125 client.loop_stop()
125 126 #client.loop_forever() No newline at end of file
General Comments 0
You need to be logged in to leave comments. Login now