为什么要用matlab做振动波动仿真
倍频程的计算？
倍频程的计算？
将全频域按几何等比级数的间隔划分，使得中心频率fc取做带宽上、下限f1、f2的几何平均值，且带宽hf2-f1 总是和中心频率fc保持一常数关系，hv×fc。如果v等于根号二的倒数(0.707)，那么f22*f1，则定义这样的频率带宽叫倍频程带宽；如果v等于三倍根号二的倒数(0.236)，那么h0.236fc，则定义这样的频率带宽为1/3倍频程带宽。
1/3倍频程作用主要是分析噪声能量的频率分布。另外做分析的时候加了计权网络可起到滤波功能。
每个倍频程或者1/3倍频程的获得是通过带通滤波实现的。但是作为总的倍频程或者1/3倍频程分析来看，主要是为了研究信号能量在不同频带的分布。
使用1/3倍频程主要是因为人耳对声音的感觉，其频率分辨能力不是单一频率，而是频带，而1/3倍频程曾经被认为是比较符合人耳特性的频带划分方法，不过现在心理声学里提出了Critical Band这么个频带划分方法，听说更符合人耳特性。
先要知道1/3倍频程的划分方法，相关的书和国标都有公式和现成的数据表格，然后，你将时间域的声信号fft变换到频率域，对定义的每个1/3倍频带的声压计算等效连续声压级。这就是1/3倍频程声压级。
FFT后再进行1/3倍频程分析，在王济和胡晓编“MATLAB在振动信号处理中的应用”(中国水利水电出版社)一书中有一节用介绍1/3倍频程分析，它是在FFT之后用1/3倍频程滤波器对信号进行分析处理，求出1/3倍频程滤波器输出的均方根值，并提供了MATLAB程序。
质量矩阵和刚度矩阵怎么求频率？
可以使用MATLAB求解这个问题，
用Matlab 数值工具箱中库函数命令eig 来实现。
[ x ,d] eig(ik ,im) ; % ik 和im 分别为结构的刚度矩阵和质量矩阵
d sqrt (d)
for i 1 :cn % cn 为结构的层数,即质点数
[ dl (i) ,j ] min (d) ;
Xgd( : ,i) x ( : ,j)
d(j) max (d) 1
end %以此循环对所求频率和振型进行排序
w dl ; %所求自振频率
x xgd ; %所求结构主振型
举个例子，可能会更清楚。
某三层钢筋混凝土结构,结构的各层特性参数为:第一层到第
三层质量m 分别为2 762 kg ,2 760 kg ,2 300 kg ,第一层到第三层
刚度k 分别为2. 485 ×104 N/ m , 1. 921 ×104 N/ m , 1. 522 ×
104 N/ m。地震波采用200gal El Centro 波,采样周期为0. 02 s。