AlGO DONE FOR GAUSSIAN

This commit is contained in:
pavanvpatil 2022-10-28 19:32:32 +05:30
parent 25dfbc2444
commit b10a48516c
2 changed files with 148 additions and 21 deletions

View File

@ -1,8 +1,8 @@
/* pavan changes */ /* pavan changes */
#include <time.h>
#include <gtk/gtk.h> #include <gtk/gtk.h>
/* variables declared of attack_detect function */ /* variables declared of attack_detect function */
/* function declared */ /* function declared */
gboolean attack_detect(struct data_frame *df); gboolean attack_detect_freq(struct data_frame *df);
gboolean attack_detect_vol(struct data_frame *df)
/* pavan changes */ /* pavan changes */

View File

@ -6,14 +6,51 @@
#include <math.h> #include <math.h>
#include "parser.h" #include "parser.h"
long double AVERAGE_OF_FREQUENCY = 50; struct freqlist
unsigned long long int COUNT = 500;
gboolean attack_detect(struct data_frame *df)
{ {
float CURR_FREQ = 50 + to_intconvertor(df->dpmu[0]->freq) * 1e-3; int idcode;
long double AVERAGE_OF_FREQUENCY;
unsigned long long int COUNT;
struct freqlist* next;
};
struct vollist
{
int idcode;
long double AVERAGE_OF_VOLTAGE;
unsigned long long int COUNT;
};
struct freqlist *head = NULL;
struct vollist *headvol = NULL
gboolean attack_detect_freq(struct data_frame *df)
{
if (head == NULL)
{
head = (struct freqlist *)malloc(sizeof(struct freqlist));
head->AVERAGE_OF_FREQUENCY = 50;
head->COUNT = 500;
return TRUE;
}
else
{
struct freqlist *temp = head;
struct freqlist *previous = NULL;
while (temp != NULL)
{
if (to_intconvertor(df->idcode) == temp->idcode)
{
float CURR_FREQ;
if (df->dpmu[0]->fmt->freq == '0'){
CURR_FREQ = 50 + to_intconvertor(df->dpmu[0]->freq) * 1e-3;
}
else{
CURR_FREQ = decode_ieee_single(df->dpmu[0]->freq);
}
printf("Current freq: %f\n", CURR_FREQ); printf("Current freq: %f\n", CURR_FREQ);
float DETECT_PERCENT = (fabs(AVERAGE_OF_FREQUENCY - CURR_FREQ) / (AVERAGE_OF_FREQUENCY * 1.0f)) * 100; float DETECT_PERCENT = (fabs(temp->AVERAGE_OF_FREQUENCY - CURR_FREQ) / (temp->AVERAGE_OF_FREQUENCY * 1.0f)) * 100;
/* detecting based on thershold frequency can vary only around 0.2 hz*/ /* detecting based on thershold frequency can vary only around 0.2 hz*/
if (DETECT_PERCENT > 0.9) if (DETECT_PERCENT > 0.9)
@ -30,10 +67,100 @@ gboolean attack_detect(struct data_frame *df)
printf("\033[0m"); printf("\033[0m");
/* updating mean of frequency */ /* updating mean of frequency */
AVERAGE_OF_FREQUENCY = ((AVERAGE_OF_FREQUENCY * COUNT) + CURR_FREQ) / ++COUNT; temp->AVERAGE_OF_FREQUENCY = ((temp->AVERAGE_OF_FREQUENCY * temp->COUNT) + CURR_FREQ) / ++temp->COUNT;
printf("avg freq: %Lf\n", AVERAGE_OF_FREQUENCY); printf("avg freq: %Lf\n", temp->AVERAGE_OF_FREQUENCY);
return TRUE; return TRUE;
} }
} }
previous=temp;
temp=temp->next;
}
if(temp==NULL)
{
struct freqlist* bring = (struct freqlist *)malloc(sizeof(struct freqlist));
bring = (struct freqlist *)malloc(sizeof(struct freqlist));
bring->AVERAGE_OF_FREQUENCY = 50;
bring->COUNT = 500;
previous->next=bring;
return TRUE;
}
}
}
gboolean attack_detect_vol(struct data_frame *df)
{
float CURR_vol;
if (df->dpmu[0]->fmt->phasor == '0')
{
unsigned char s1[2];
unsigned char s2[2];
strncpy(s1, df->dpmu[0]->phasors[0], 2);
strncpy(s2, df->dpmu[0]->phasors[0] + 2, 2);
long double v1 = to_intconvertor(s1);
long double v2 = to_intconvertor(s2);
CURR_vol = sqrt((v1 * v1) + (v2 * v2));
}
else
{
unsigned char s1[2];
unsigned char s2[2];
strncpy(s1, df->dpmu[0]->phasors[0], 4);
strncpy(s2, df->dpmu[0]->phasors[0] + 2, 4);
long double v1 = decode_ieee_single(s1);
long double v2 = decode_ieee_single(s2);
CURR_vol = sqrt((v1 * v1) + (v2 * v2));
}
if (head == NULL)
{
head = (struct vollist *)malloc(sizeof(struct vollist));
head->AVERAGE_OF_VOLTAGE = CURR_vol;
head->COUNT = 500;
return TRUE;
}
else
{
struct vollist *temp = head;
struct vollist *previous = NULL;
while (temp != NULL)
{
if (to_intconvertor(df->idcode) == temp->idcode)
{
printf("Current vol: %f\n", CURR_vol);
float DETECT_PERCENT = (fabs(temp->AVERAGE_OF_VOLTAGE - CURR_vol) / (temp->AVERAGE_OF_VOLTAGE * 1.0f)) * 100;
/* detecting based on thershold VOLTAGE can vary only around 0.2 hz*/
if (DETECT_PERCENT > 5)
{
printf("\033[0;31m");
printf("ATTACK DETECTED!");
printf("\033[0m");
return FALSE;
}
else
{
printf("\033[0;32m");
printf("NO PROBLEM :)\n");
printf("\033[0m");
/* updating mean of VOLTAGE */
temp->AVERAGE_OF_VOLTAGE = ((temp->AVERAGE_OF_VOLTAGE * temp->COUNT) + CURR_vol) / ++temp->COUNT;
printf("avg vol: %Lf\n", temp->AVERAGE_OF_VOLTAGE);
return TRUE;
}
}
previous = temp;
temp = temp->next;
}
if (temp == NULL)
{
struct vollist *bring = (struct vollist *)malloc(sizeof(struct vollist));
bring = (struct vollist *)malloc(sizeof(struct vollist));
bring->AVERAGE_OF_VOLTAGE = CURR_vol;
bring->COUNT = 500;
previous->next = bring;
return TRUE;
}
}
}
/* pavan changes */ /* pavan changes */