MATLAB计算方法和技巧10_1阻尼频谱.docVIP

连续阻尼函数的频谱 [问题]物体作初速度为零的阻尼运动的位移为 ，(t ≥ 0) 其中，A0为t = 0时的初位移，β为阻尼因子，ω为准圆频率，ω0是无阻尼简谐振动的圆频率。试求物体作阻尼运动的频谱，并画出频谱曲线。 [数学模型] 方法一：用解析公式。非周期性的振动可以分解为一系列简谐振动，它可以当作周期T→∞或基频ω→0的周期性运动，因此分解的一系列简谐振动的频率是连续分布的。根据傅立叶积分，可将振动x(t)按频率ω分解为 ， (10.1.1) 可以证明：频率为ω的振幅为 。 (10.1.2) 其模就是频谱。 阻尼振动在t 0时可认为x = 0，因此频谱为 。 (10.1.3) 利用复数积分 ， 和欧拉公式eix = cos x + isin x，令上式两边实部和虚部分别相等可得 ， (10.1.4a) 。 (10.1.4b) 因此频谱为 ， 当t→∞时|e-iωt| ≤ 1，e-βt→0，因此 ， 利用可得振幅为 。 (10.1.5) 阻尼运动的频谱为 。 (10.1.6) 如果β = 0，当ω→ω0时，|F(ω)|→∞。 [算法]由(10.1.5)式可得简化振幅 。 (10.1.5*) 其中，ω* = ω/ω0，β* = β/ω0。物体作阻尼运动的简化频谱为 。 (10.1.6*) [程序].m如下。 %阻尼振动(非周期性振动)的傅立叶变换(用解析式) clear %清除变量 b=input(请输入相对阻尼因子:); %键盘输入阻尼因子与固有角频率的倍数 if b==1,b=1-eps;end %1则改小一点 t=0:0.1:20; %固有角频率与时间的乘积w0t向量 ww=sqrt(1-b^2); %准角频率向量 x=exp(-b*t).*(cos(ww*t)+b/ww*sin(ww*t));%位移函数 f1=figure; %创建图形窗口 plot(t,x) %画曲线 grid on %加网格 fs=16; %字体大小 xlabel(\it\omega\rm_0\itt,FontSize,fs)%标记横坐标 ylabel(\itx/A\rm_0,FontSize,fs) %标记纵坐标 title(质点阻尼运动的位移曲线,FontSize,fs)%标题 text(0,0,[\it\beta/\omega\rm_0=,num2str(b)],FontSize,fs)%参数 w=0:0.001:3; %简化角频率 f=sqrt((4*b^2+w.^2)./((1-w.^2).^2+4*b^2*w.^2));%频谱 f2=figure; %创建图形窗口 plot(w,f,LineWidth,2) %画频谱曲线 grid on %加网格 xlabel(\it\omega/\omega\rm_0,FontSize,fs)%标记横坐标 ylabel(\itF\omega\rm_0/\itA\rm_0,FontSize,fs)%标记纵坐标 title(质点阻尼运动的频谱曲线,FontSize,fs)%标题 text(0,0,[\it\beta/\omega\rm_0=,num2str(b)],FontSize,fs)%参数 P10_1_1a图 P10_1_1b图 P10_1_1c图 P10_1_1d图 方法二：用快速傅立叶变换指令。任何一个非周期函数都可以看成是由某个周期函数当周期趋于无穷大转化而来的。为了便于计算，往往取足够长的时间当作周期T，并且当t 0时，非周期函数当零处理，即 。 (10.1.17) 在(0