《常微分方程》实验报告四.doc

《常微分方程》实验报告四 专业 信息与计算科学 班级 　 　 姓名 　 学号 　 实验地点 实验室　 实验时间2015.12.11　 实验名称：机械振动中的单摆问题 实验目的：1．与计算机软件Maple结合起来，学会求常微分方程组解析解的方法； 2．会用Maple软件做简单的计算机仿真； 3．能够理论联系实际，从理论结论中客观分析实践中的现象。 实验内容：（给出实验程序与运行结果） 问题：如图所示为一单摆示意图，物块的质量为，摆线长为，如图建立坐标系， （1）建立系统的运动函数，即物块坐标随时间的函数、； （2）设=10kg,L =1m，系统作微摆动，存在空气阻力，空气阻力与速度成正比，比例系数为u，u=2， 求在初始条件x[0]=0 , [0]=0 下的系统运动函数、； （3）当系统作微摆动，不存在空气阻力时系统的运动函数、； （4）合理运用Matlab分别画出（2）、（3）中确定的各个函数的图象； （5）对（2）、（3）所得结果进行分析。根据图象说明物块是做阻尼运动、无阻尼运动，还是周期运动；并分析给定参量具体数值后对反映系统的运动规律有无影响。 LM L M x 可以发现物体运动遵循牛顿第二运动定律。 2．对内容（2）的分析：与（1）相比，物块做圆周运动是一样的；不同的是，物块作微摆动时，受力分析后我们作一个近似，，从而对在x方向上平衡建立一个方程。 利用Maltab软件求方程组的解析解：一般情况下，对高阶方程首先引进新参量对方程进行降阶，之后直接使用dsolve命令直接求出方程组的解。但是这里求出的解可能并不简单明了，故而必须求它的数值解。 3．对内容（3），与（2）类似，而且相对简单。由于这里不再受空气阻力。 4.对内容（4），用Matlab软件中简单的作图命令plot分别对（2）、（3）中所得函数进行做图。 5.对内容（5）通过分析（2）、（3）所得函数的图象，判断物块分别做什么样的运动。是阻尼、无阻尼运动？还是周期运动？运动有无衰减现象？我们给定具体参数值对反映系统的运动规律有无影响？ 实验分析： 1．（1）的模型的建立 使物块运动的力F是物块的重力在切线方向的分力： . 而物块的加速度为 . 根据牛顿第二定律，得 . 上式右端的负号是由于力F与位移s的正方向相反的缘故。上式可简化为 . 初始条件为： 2．的模型的建立 由于系统作微摆动，故。又存在空气阻力，空气阻力与速度成正比，比例系数为u，所以摆动的方程为 . 将=10kg,L =1m，u=2代入得： 3．的模型的建立 当系统作微摆动，不存在空气阻力时系统的运动函数为 用Matlab分别画出（2）、（3）中确定的各个函数的图象， 对（2）先看看有没有解析解，运行Maltlab程序（1）得s =0，ans =0 知原方程没有解析解，下面求解数值解。令可将原方程化为如下方程组 建立M文件fun9.m如下即%M文件 fun9.m function f=fun9(t,y) f=[y(2), -9.8*y(1)-0.2*y(2)]; 运行MATLAB代码（2）结果见图1： 图1 数值解 由图1可见，随时间有衰减现象。 对（3）用同样的方法可得结论。 令可将原方程化为如下方程组 建立M文件fun9.m如下： %M文件 fun9.m function f=fun9(t,y) f=[y(2), -9.8*y(1)]; %f向量必须为一列向量 运行MATLAB代码（3）结果见图2 图2 数值解 由图2可见，随时间周期变化。 5．根据图1和图2可知，（2）的物块做阻尼运动，（3）的物块做周期运动； 并分析给定参量具体数值后对反映系统的运动规律有无影响。对于（2），摆线越长频率越小；空气阻力比例系数越大，衰减现象越明显；质量越大，衰减速度变缓，对于（3），初始值影响振幅。 源代码 1.有空气阻力 clear; s=dsolve(D2y+0.2*Dy+9.8*y=0,y(0)=0,Dy(0)=0,t) simplify(s) function?f=fun9(t,y)f=[y(2),?-9.8*y(1)-0.2*y(2)];? clear;close;? [t,y]=ode45(fun9,[-10,10],[0.5,0]);plot(t,y(:,1));?xlabel(t),ylabel(y1) 2.无空气阻力 clear; s=dsolve(D2y=9.8*sin(y),y(0)=15,Dy(0)=0,t) simplify(s)