itd ITD模态参数识别的matlab源程序

ITD模态参数识别的matlab源程序，比较好用

function [FDS]=ITD（xmnfs）
%建立特征方程矩阵的阶数（为模态阶数的2倍）
b=x;
nm=2*mn;
%取时域数据1/2的长度
n=fix（length（b）/2）;
%将输入时域数据赋值给列向量h
h=b（1:2*n）;
%计算时间间隔
dt=1/fs;
%建立离散时间向量
%t=0:dt:（2*n-1）*dt;
%计算自由振动响应矩阵
L=length（h）;
for k=1:nm
    x1（k:）=h（k:L-（nm-k+1））‘;
    x2（k:）=h（k+1:L-（nm-k））‘;
end
%用最小二乘法求解特征方程矩阵
B=x1/x2;
%计算特征值及特征向量
[AV]=eig（B）;
%变换特征值对角阵为一向量
for k=1:nm
    U（k）=V（kk）;
end
%计算模态频率向量
F1=abs（log（U‘））./（2*pi*dt）;
%w=atan（imag（U‘）./real（U‘））/dt;
%d=-log（imag（U‘）.^2+real（U‘）.^2）./（2*dt）;
%F1=（w.^2+d.^2）.^0.5;
%计算阻尼比向量
D1=sqrt（1./（（imag（log（U‘））./real（（log（U‘）））.^2）+1））;
%D1=d./w;
%计算振型系数向量
for k=0:（2*n-1）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1494  2010-05-23 13:49  itd\ITD.m

     目录          0  2010-12-01 17:01  itd

----------- ---------  ---------- -----  ----

                 1494                    2