数学建模案例-大学物理-阻尼振动.pdfVIP

数学模型案例-阻尼振动的数学模型 问题背景与描述： 简谐振动是一种无阻尼振动，而实际上，任何振动物体都要受到阻力的作用，这种振动 叫阻尼振动。阻尼振动有摩擦阻尼和辐射阻尼两种。 根据实验证实，当物体以不太大的速率在粘性介质中运动时，介质对物体的阻力与物体 的运动速率成正比，方向与运动方向相反。 数学模型 阻力与速率的关系： 阻尼系数，它与物体的形状、大小及介质有关。 对弹簧振子，在弹性力及阻力的作用下，物体的运动方程为 （此处自行思考：得到该模型的基本原理是什么？） 令 ， 无阻尼时振子的固有角频率， 为阻尼因子， 代入上式后可得等价形式： 数学模型的解 在阻尼作用较小时，即 ，此微分方程的解为 t x  A e cos(t  ) 0 0 其中 此时阻尼振动的周期为 当 时，这种阻尼振动称为欠阻尼； 当 时，称为过阻尼； 当 时，称为临界阻尼。 *****进一步的学习资料*******： 解上述二阶微分方程，Matlab 只要用一个函数ode23 就能解决问题。那我 们现在就开始学吧！ode23 函数可以用来解微分方程组，但只能是一阶的。所以对于二阶 的微分方程，我们可以将其分解成两个一阶的方程来解。而且使用此函数，必须提前编写一 个ode 文件 （ode 文件）。 ode 文件是一种函数M 文件，即：在m-file 编辑器中，文件必须以function 开头，后 面跟的函数名必须与将来保存的M 文件名一致。函数文件第一句的具体形式为： function y=name (x) 其中为y 输出变量名，x 为变量，name 为函数名，保存时文件名必须是name.m。当函 数有一个以上的输出函数时，输出参数包含在方括号里，变参包含在圆括号里，即函数文件 的第一行的具体形式为： function [out1,out2,…] name (in1,in2,…) 例：计算数组元素的平均值 解：functiony=average(x) y=sum(x)/length(x) 解微分方程的函数ode23 调用格式为： ode23(FUN,TSPAN,YO) 其中， FUN 为字符串，表示微分方程的ode 文件名， TSPAN=[TO,TFINAL]表示积分区间， YO 为初始条件。此函数表示在初始条件下Y0 从T0 到TFINAL 对微分方程y=F(t,y)进 行积分。 函数在图形窗口返回具体的数值解，则可以采用以下的格式： [T,Y]=ode23(FUN,TSPAN,Y0) 这样在T 向量返回积分点，Y 返回对应积分点的积分函数值。 例：编写m-file 描述物体作阻尼振动的振动曲线，分欠阻尼、过阻尼和临界阻尼三种情况。 解：首先将阻尼振动的微分方程 分解为一阶微分方程：令 则原方程化为 打开m-file 编辑器， 在其中编写ode 文件为 functionf=zuni(t,x,flag,beita,wo); f=[x(2); -2*beita.*x(2)-wo.^2*x(1)]; 保存文件名为zuni.m 另外打开一个新m-file，在其中编写程序代码： beita=input(欠阻尼系数p=); wo=input(固有频率wo(pwo)=); [t,x]=ode23(zuni,[0:0.001:5*pi],[6,8],[],beita,wo);%调用ode23,[6,8]为初 始条件， %即x=6,x=8,[]说明以后是输入参数，beita,wo 是参数名 subplot(3,1,1),plot(t,x(:,1),r),title(欠阻尼)%取x(1),绘制t-x(1)位移曲 线，红色 beita=input(过阻尼系数p=)