clear;clc;close all
global R L km J D
Tres=0;
Vdc=240;
x0=[1 10.4e-3 1.8];
% Limites de las variables R=x(1)>=0, y C=x(2)>=0
lb = x0/2;
ub = 2*x0;

options = optimoptions('ga','MaxGenerations', 50, 'PopulationSize', 20);

x = ga(@error_medio_cuadratico_mcc,3,[],[],[],[],lb,ub,[],options);

R=x(1);L=x(2);km=x(3);
sim MyS2_ag_mdc_mdl;

plot(t,ia,t,ia1,'--r')
grid on
legend('modelo','experimental')

figure
plot(t,omega,t,omega1)
grid on
legend('modelo','experimental')