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