我是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