汕头公司做网站,山西网络网站建设销售公司,自己搭建服务器做网站,thinkphp官方网站ode45求解多自由度动力学方程实例Ode45函数调用形式如下#xff1a;[T,Y]ode45(odefun,tspan,y0)相关参数介绍如下#xff1a;参数名称参数说明odefun用于存放待求解的方程的m文件名#xff0c;方程必须用y’f(t,y)的形式存放tspan指定自变量范围的向量#xff0c;通常用[t…ode45求解多自由度动力学方程实例Ode45函数调用形式如下[T,Y]ode45(odefun,tspan,y0)相关参数介绍如下参数名称参数说明odefun用于存放待求解的方程的m文件名方程必须用y’f(t,y)的形式存放tspan指定自变量范围的向量通常用[t0tf]指定y0函数的边界条件即y0y(t0)对于方程组y0也可以是向量例若一三自由度多体动力学系统方程如下初始条件由于方程必须用y’f(t,y)的形式存放因此需要对方程组进行降阶处理。令则方程组可化为因此建立M函数文件来定义此方程组如下function dyfunc(t,y)dyzeros(6,1);dy(1)y(2);dy(2)0.5*(-3*y(2)1.5*y(4)-100*y(1)50*y(3)2.0*sin(3.754*t)); dy(3)y(4);dy(4)0.5*(1.5*y(2)-3*y(4)1.5*y(6)50*y(1)-100*y(3)50*y(5)-2.0*cos(2.2*t));dy(5)y(6); dy(6)0.5*(-1.5*y(4)-3*y(6)50*y(3)-100*y(5)1.0*sin(2.8*t));end在matlab命令窗口里输入一下命令y0[1 1 1 1 1 1];tspan[0 30];[t,y]ode45(func,tspan,y0);figure(1)plot(t,y(:,1),t,y(:,3),t,y(:,5));legend(x1,x2,x3);xlabel(时间(s),FontSize,10);ylabel(振动位移曲线,FontSize,10);figure(2)plot(t,y(:,2),t,y(:,4),t,y(:,6));legend(v1,v2,v3);xlabel(时间(s),FontSize,10);ylabel(‘振动速度曲线’,’FontSize’,10);运行结果
