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=optimset('Display','iter','MaxIter',500,'LargeScale',...
    'off','GradObj','off','GradConstr','off','TolFun',1e-3,'MaxFunEvals',20000','Algorithm','active-set');
    %Optimization
x = fmincon(@error_medio_cuadratico_mcc,x0,[],[],[],[],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')