特殊弦振动问题的MATLAB仿真
特殊弦振动问题的MATLAB仿真
【摘要】本文针对一个特殊弦振动问题进行模型求解MATLAB进行仿真。
【关键词】弦振动;MATLAB;仿真
引言
本文主要讨论了如下弦振动问题:一根均匀弦两端分别在处固定,设初始速度为零,初始时刻弦的形状为一抛物线,抛物线的顶点为,讨论其弦振动的位移的情况。
问题求解
由所给出的问题可以得到下面的定界问题:
可以设想该泛定方程具有分离变量形式的特解:
将该式代入原泛定方程得:
因为等式两边的变量不同所以只有当两边同时等于一个常数时才能成立,设此常数为,则有:
由此可得到两个常微分方程:
现在将代入边界条件得:
由于不是所需要的解,所以不恒为零,则有
下面,我们来求解特征值和特征函数的问题:
为了求解这个特征值问题,需分一下三种情形讨论:
当时,方程的通解为:
其中,A,B为两任意常数,代入边界条件得:,即,舍去舍去
所以可得弦振动位移函数
问题的仿真
第二节中,得到的函数表达式
0项,观察仿真结果。
3.1仿真代码
function?z?=?My_Function?(?x,t,L,a,h,k?)z?=?0;for?n?=?0:k????z?=?z?+?32*h/pi^3*(1/(2*n+1)^3)*cos((2*n+1)*pi*a*t/L)*sin((2*n+1)*pi*x/L);endend%%clear;clc;X?=?0:0.2:10;T?=?0:0.2:5;LX?=?length(X);LT?=?length(T);Z?=?zeros(LX,LT);L?=?10;a?=?4;h?=?8;k?=?20;for?i?=?1:LX????for?j?=?1:LT?????????Z(i,j)?=?My_Function(X(i),T(j),L,a,h,k);????endendfiguresurf(T,X,Z);
3.2仿真结果
图1 一个周期()图像
周期(10)图像
从图1中我们可以看出一根长为10单位,初始时刻为顶点在(5,8)的抛物线的弦为