case118.py

import numpy as np
import matplotlib.pyplot as plt
from pandapower.networks import case118
import tqdm

import functional as F


def main():
    sep118 = case118()

    init = F.initialize_pes(sep118)

    optimization_iters = 10
    iterations = 20
    num_individuals = 50
    fitness = np.zeros((iterations, optimization_iters))
    perdas = np.zeros((iterations, optimization_iters))
    penaltygen = np.zeros((iterations, optimization_iters))
    penaltyv = np.zeros((iterations, optimization_iters))
    tempos = np.zeros((iterations, optimization_iters))
    solutions = np.zeros((iterations, 83))

    j, perdas_, pen_v, pen_gq, pen_tap, pen_bsh, global_best, tempo = F.optimize_ga(
        sep118, 
        pt=0.7, 
        rgap=0, 
        zeta=0.01, 
        psi=0.05, 
        sigma=0.01, 
        omega=0.01, 
        max_iter=optimization_iters, 
        num_individuals=num_individuals, 
        c1=1, 
        c2=0, 
        v_amp=0.1, 
        valor_inicial=0, 
        step=0.01, 
        wmax=0.9, 
        relatorio=True, 
        inicial=False
    )

    breakpoint()


if __name__ == "__main__":
    main()