Как построить три графика одновременно на одной оси?

Question

Я новичок в этой области кодирования Chaos & MATLAB. У меня проблема, как описано ниже.

Ожидаемый результат: я пытаюсь получить точные графики фазовой плоскости и график временных рядов Aceng Sambas et.al. система, как показано на рисунке внизу его бумаги. Первый рисунок показывает фазовый простой портрет, а второй рисунок показывает отклик во временной области. Ссылка на этот документ приведена здесь [ 1 ]. {Этот материал был взят из сети}

Постановка задачи: я не могу получить несколько графиков на одной оси, и я не могу понять логику c за получением нескольких графиков по тем же осям.

Как я поступил: я начал записывать код для этой системы, так как есть три различных значения параметров K, я пытался решить их отдельно и позже я попытался получить график для всех трех различных значений K.

Результат, который я получил: я не могу получить желаемый график фазы Phase на графике 2-го и временного ряда (или временной области - это график для состояние v / s времени). Я основал определенный комментарий MATLAB, который был полезен при построении нескольких графиков вместе, и я включил их, но не повезло.

.

Я предоставил свои коды ниже, извините меня, если я предоставил слишком много кодов.

function  Aceng_Sambas_Parameter_Boosting
function  Aceng_Sambas_Parameter_Boosting_2
function  Aceng_Sambas_Parameter_Boosting_3
t0 = 0; % Initial point
tn = 600; % Terminal point
Y0 = [0.2, 0.2, 0.2]; % Initial values

a = 5;b = 2;k1=0;% system parameters
a = 5;b = 2;k2=2;% system parameters
a = 5;b = 2;k3=5;% system parameters

[t, Y] = ode45(@ Aceng_Sambas_Parameter_Boosting, [t0, tn], Y0); %Solving the equations 
[t2, Y2] = ode45(@ Aceng_Sambas_Parameter_Boosting_2, [t0, tn], Y0);%Solving the equations 
[t3, Y3] = ode45(@ Aceng_Sambas_Parameter_Boosting_3, [t0, tn], Y0);%Solving the equations 

% Create a(2-D)Plot 
plot(Y(:,2),Y(:,3),Y2(:,2),Y2(:,3),Y3(:,2),Y3(:,3),'Color', [0.72,0.27,1.00],'-r','-g') 
xlabel('y')
ylabel('z')
legend
title(' Aceng_Sambas System')
hold on
fig()

% Time series plot

plot(t,Y(:,3),t2,Y2(:,3),t3,Y3(:,3),'Color', [0.72,0.27,1.00],'-r','-g')
xlabel('t')
ylabel('z')
legend
title('Aceng_Sambas System')
hold on
fig()

  function dYdt =  Aceng_Sambas_Parameter_Boosting(~,X) %System equation for different values of k

    dxdt =  X(2)*(X(3)+k1);
    dydt = X(1)-X(2);
    dzdt = a*abs(X(1))- b*X(1)^2;
    dYdt = [dxdt; dydt; dzdt];

    function dYdt =  Aceng_Sambas_Parameter_Boosting_2(~,H) %System equation for different values of k

    dxdt =  H(2)*(H(3)+k2);
    dydt = H(1)-H(2);
    dzdt = a*abs(H(1))- b*H(1)^2;
    dYdt = [dxdt; dydt; dzdt];

    function dYdt =  Aceng_Sambas_Parameter_Boosting_3(~,P) %System equation for different values of k

    dxdt =  P(2)*(P(3)+k3);
    dydt = P(1)-P(2);
    dzdt = a*abs(P(1))- b*P(1)^2;
    dYdt = [dxdt; dydt; dzdt];

    end
    end

end
end
    end
    end