如何在同一轴上绘制三个同时图形?

时间:2020-02-08 12:27:10

标签: matlab matrix plot

我是Chaos&MATLAB编码领域的新手。我有一个如下所述的问题。

预期的输出:我正在尝试获取Aceng Sambas等人的确切相平面图和时间序列图。系统如下图所示。 第一个图显示相位平面肖像,第二个图显示时域响应。本文的链接位于[1]。 {此资料是从网上取得的}

问题陈述:我无法在同一轴上获得多个图,并且我也无法理解在同一轴上获得多个图的逻辑。

我如何进行:由于存在三个不同的参数K值,因此我开始写下该系统的代码,我尝试分别求解它们,后来又尝试获取K的所有三个不同值的图

我获得的结果:我无法在2-d和时间序列图(或时域-这是状态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


Phase plain Portrait as shown in their research paper- the desired result that I need Time domain plot as shown in their research paper- the desired result that I need

0 个答案:

没有答案