close all;
clear ;
%%
    set(0,'DefaultAxesFontSize',12, 'DefaultAxesFontName','Times New Roman');
    set(0,'DefaultTextFontSize',12, 'DefaultTextFontName','Times New Roman');
    set(0,'DefaultAxesXGrid','on', 'DefaultAxesYGrid','on');
    set(0,'DefaultLineLineWidth',2);
    set(0,'DefaultAxesColorOrder',[0 0 0]);
%%
% Xticks = [0:0.2:1];
%%

% model = 'vc_direct_ad.mdl';
model = 'vc_direct_ad_robast.slx';

load_system(model)
tpx=0.48;
tpy=0.999;

variationMc = [0.25,0.5,0.75,1]*0.484;
variationRangeR2 = [0.5:0.2:2];
variationColour = ['k', 'r', 'g', 'y'];
variationLegend = ['0,25*M_n', '0,5*M_n', '0,75*M_n', 'M_n']
    vR1 = 1;
    vR2 = 1;

%% Figure 1
A=figure;
set(A, 'Units','centimeters', 'Position',[0 0 17 11.5]);
movegui(A, 'center');
q = 0;
for Mcr = variationMc
    q=q+1;

    arrayiq = [];
    arraypa = [];
    i = 0;

    for vR2 = variationRangeR2
        i = i + 1;
        sim(model)
        arrayiq(i) = iq(end);
        arraypa(i) = Pa(end);

    end


    %%
    sub1 =subplot(1,2,1), plot(variationRangeR2,arrayiq,variationColour(q))
    hold on;
    grid on;
    %xlim([t(1) t(end)]); ylim([min(iq) max(iq)]);
    title('Струм ротора i_q, А','Units','normalized','Position',[tpx tpy]);


    %%
    sub2 = subplot(1,2,2), plot(variationRangeR2,arraypa,variationColour(q))
    hold on;
    grid on;
    %xlim([t(1) t(end)]); ylim([min(pa)-10 max(pa)+10]);
    title('Активна потужність, Вт','Units','normalized','Position',[tpx tpy]);

end