倒立摆实验程序

巧了我就是萌 提交于 2019-12-29 20:22:02

这里写自定义目录标题

倒立摆实验程序

clear;
t0=0; tf=20;
x0=[0.5;-0.2;0;0];
[t,x]=ode45(‘dlfun’,[t0,tf],x0);
figure;
plot(t,x(:,1),‘r’,t,x(:,2),‘b’);
title(‘输出曲线’);
legend(‘旋臂角度’,‘摆杆角度’);

function xdot=dlfun(t,x);
m1=0.2; m2=0.052; l1=0.10; l2=0.12; r=0.2; km=0.0236; ke=0.2865;
g=9.8; j1=0.004; j2=0.001; f1=0.01; f2=0.001;
a=j1+m2rr; b=m2rl2; c=j2; d=f1+kmke;
e=(m1
l1+m2r)g; f=f2; h=m2l2g;
u=1;
xdot=zeros(4,1);
xdot(1)=x(3);
xdot(2)=x(4);
xdot(3)=((-dc).x(3)+(fbcos(x(2)-x(1))).x(4)+bbsin(x(2)-x(1)).cos(x(2)-x(1)).x(3).x(3)-bcsin(x(1)-x(2)).x(4).x(4)+ecsin(x(1))-hbsin(x(2)).cos(x(2)-x(1))+kmcu)/(ac-bb.cos(x(1)-x(2)).cos(x(2)-x(1)));
xdot(4)=((d
b
cos(x(1)-x(2))).x(3)-(af).x(4)-ab
sin(x(2)-x(1)).x(3).x(3)+bbsin(x(1)-x(2)).cos(x(1)-x(2)).x(4).x(4)-ebsin(x(1)).cos(x(1)-x(2))+ahsin(x(2))-bcos(x(1)-x(2))kmu)/(ac-b*b.*cos(x(2)-x(1)).*cos(x(2)-x(1)));
end

导出

如果你想尝试使用此编辑器, 你可以在此篇文章任意编辑。当你完成了一篇文章的写作, 在上方工具栏找到 文章导出 ,生成一个.md文件或者.html文件进行本地保存。

导入

如果你想加载一篇你写过的.md文件,在上方工具栏可以选择导入功能进行对应扩展名的文件导入,
继续你的创作。

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!