clc
clear
close all
t0=0;tf=250;
step_qw=0.01;
p1=20;
q0=[(rand(1)-0.5)*p1;(rand(1)-0.5)*p1;0.2;0.2;
(rand(1)-0.5)*p1;(rand(1)-0.5)*p1;0.1;0.1;
(rand(1)-0.5)*p1;(rand(1)-0.5)*p1;0.1;0.1;
(rand(1)-0.5)*p1;(rand(1)-0.5)*p1;0.1;0.1;
(rand(1)-0.5)*p1;(rand(1)-0.5)*p1;0.1;0.1];
%这一组是赋初值,里面有些给的是随机值,也可以改成定值
%save('q0.mat','q0');
load('q0.mat','q0');
[t,q]=ode45('Whole_Twodimension_total_state',t0:step_qw:tf,q0);
figure(1);
hold on
xlabel('x position','FontName','Times New Roman','FontSize',14);ylabel('y position','FontName','Times New Roman','FontSize',14);
plot(q(:,1),q(:,2),'g');
hold on
plot(q(:,5),q(:,6),'r');
hold on
plot(q(:,9),q(:,10),'b');
hold on
plot(q(:,13),q(:,14),'k');
hold on
plot(q(:,17),q(:,18),'c');
hold on
legend('leader ','follower 1','follower 2','follower 3','follower 4',4);
hold on
plot(q(1,1),q(1,2),'ro');
hold on
plot(q(1,5),q(1,6),'bo');
hold on
plot(q(1,9),q(1,10),'bo');
hold on
plot(q(1,13),q(1,14),'bo');
hold on
plot(q(1,17),q(1,18),'bo');
hold on
plot(q(end,1),q(end,2),'rp');
hold on
box on
figure(2)
hold on
xlabel('t/s','FontName','Times New Roman','FontSize',14);ylabel('$\bar{v}_i$','interpreter','latex','fontsize',14);
hold on
plot(t,q(:,3),'r:',t,q(:,4),'r');
hold on
plot(t,q(:,7),'b:',t,q(:,8),'b');
hold on;
plot(t,q(:,11),'k:',t,q(:,12),'k');
hold on;
plot(t,q(:,15),'m:',t,q(:,16),'m');
hold on
plot(t,q(:,19),'m:',t,q(:,20),'m');
box on
legend('leader ','follower 1','follower 2','follower 3','follower 4',4);
figure(3)
hold on
xlabel('t/s','FontName','Times New Roman','FontSize',14);ylabel('xi-x0','interpreter','latex','fontsize',14);
hold on
plot(t,q(:,5)-q(:,1),'r:',t,q(:,6)-q(:,2),'r');
hold on
plot(t,q(:,9)-q(:,1),'b:',t,q(:,10)-q(:,2),'b');
hold on;
plot(t,q(:,13)-q(:,1),'k:',t,q(:,14)-q(:,2),'k');
hold on;
plot(t,q(:,17)-q(:,1),'m:',t,q(:,18)-q(:,2),'m');
hold on
legend('leader ','follower 1','follower 2','follower 3','follower 4',4);
box on
figure(4)
hold on
xlabel('t/s','FontName','Times New Roman','FontSize',14);ylabel('vi-v0','interpreter','latex','fontsize',14);
hold on
plot(t,q(:,3)-q(:,7),'r:',t,q(:,4)-q(:,8),'r');
hold on
plot(t,q(:,3)-q(:,11),'b:',t,q(:,4)-q(:,12),'b');
hold on;
plot(t,q(:,3)-q(:,15),'k:',t,q(:,4)-q(:,16),'k');
hold on;
plot(t,q(:,3)-q(:,19),'m:',t,q(:,4)-q(:,20),'m');
hold on
box on
legend('leader ','follower 1','follower 2','follower 3','follower 4',4);
grid off

心梓
- 粉丝: 865
- 资源: 8040
最新资源
- 基于51单片机protues仿真的的智能窗帘有三个模式设计(仿真图、源代码)
- 多配送中心选址与车辆路径优化的集成策略:遗传算法在MDVRPTW中的应用及其Matlab代码解析,多配送中心选址与车辆路径优化问题的遗传算法研究:Matlab完整代码实现及数据可修改,多配送中心车辆路
- 基于51单片机protues仿真的电子钟仿真程序(7个功能)
- 两级式单相光伏并网仿真研究:MATLAB 2021a版本下的DC-DC变换与桥式逆变技术实现功率跟踪与并网效果优化,基于Matlab 2021a的两级式单相光伏并网仿真研究:实现最大功率跟踪与稳定的直
- WinCC高级报表功能强大:自定义模版、打印预览与PDF、Excel导出、自定义打印区域及数据统计(最大值、最小值、平均值、求和),适用于WinCC 7.4/7.5/8.0及博图WinCC ,WinC
- 基于51单片机protues仿真的多功能音乐播放设计(仿真图、AD原理图、源代码)
- 基于BP神经网络对一阶线性系统的拟合分析与实验结果,基于BP神经网络拟合一阶线性系统的算法实践与测试结果,BP神经网络拟合一阶线性系统 (A)Train.mlx、Test.mlx、Weight.mat
- GVIM,WINDOWS版本的VIM
- 光伏并网储能逆变器的技术原理及应用前景探索:高效、智能、可持续能源解决方案,光伏并网储能逆变器的技术原理及应用前景探索:高效、智能、可持续能源解决方案,光伏并网储能逆变器 ,光伏并网; 储能; 逆变器
- 基于51单片机protues仿真的多模式带音乐跑马灯(仿真图、源代码)
- Simulink下二极管钳位型三电平SVPWM闭环并网系统:三相波形对称,五电平输出稳定,高效功率跟踪与高质量并网性能,Simulink下的三电平二极管钳位型SVPWM闭环并网系统:800V直流输入
- 基于51单片机protues仿真的多种信号发生器、波形发生器设计(仿真图、源代码)
- 三相桥式两电平逆变器并网仿真:LCL滤波器下SPWM调制与电流内环PI解耦策略,功率输出稳定且并网性能优异,三相桥式两电平逆变器并网仿真:LCL滤波器下SPWM调制与PI解耦控制,稳定功率输出,低TH
- 基于51单片机protues仿真的楼道人数计数器设计(仿真图、源代码、论文)
- 《三菱PLC高效编程利器:自编ST语言库,通用性强且提高工作效率》,《三菱PLC高效编程利器:自编ST语言库,通用性强且提高工作效率》,三菱PLC结构化编程,自己为了方便编程编写的ST语言库,结构清晰
- yolov8图像分割苹果数据集,包括标注json文件
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈



- 1
- 2
前往页