matlab双摆程序,matlab-双摆仿真

在物理学和数学中，在动力系统领域，双摆是一个摆锤，另一个摆锤连接在其末端，是一个简单的物理系统，具有丰富的动态特性，对初始条件双摆的运动由一组耦合的常微分方程控制并且是混沌的。

由于双摆运动涉及到二阶微分方程组，在matlab中对双摆的仿真需要用到ode45求其数值解。

clear all;

%控制r1 r2两个角位移 就可改变双摆初状态

r1=0.8;

r2=0.8;

m1=1;

m2=1;

L1=1;

L2=1;

g=9.8;

Da=inline(['[x(3);x(4);',...

'inv([(m1+m2)*L1,m2*L2*cos(x(1)-x(2));',...

'm1*L1*cos(x(1)-x(2)),m1*L2])*'...

'[m2*L2*x(4)^2*sin(x(2)-x(1))-(m1+m2)*g*sin(x(1));',...

'm2*L1*x(3)^2*sin(x(1)-x(2))-m2*g*sin(x(2))]]'],'t','x',...

'flag','m1','m2','L1','L2','g');

set(gcf,'DoubleBuffer','on');

[t,x]=ode