bonzer 发表于 2006-12-2 20:56

求 ：模态识别ITD法程序，急！感激！！

恳请高手帮忙，小弟感激不尽！

pengweicai 发表于 2006-12-3 16:08

这个是我以前做作业时留下的。MATLAB程序。
    信号是一组正弦，最后一项 1e-5*ones(4,2048)....为噪声。
   以前不怎么会用，所以程序不是很好。
         仅供参考。


%ITD(DLS)
%M--测点个数   %R--采样点个数   %ts<0.004   delt<0.008ks=4
SSg=; %阻尼理论值
WWd=         %频率理论值
M=5;R=512;ts=0.001;ks=4;delt=ks*ts;R2=2*R;M2=2*M;pi=3.1415926;
%t=;
%t=;
t=;
y1(:,:)=3.5e-4*exp(-0.2*t).*sin(10*t+0.13*pi)+2.0e-4*exp(-0.15*t).*sin(26*t+0.20*pi)+1.00e-4*exp(-0.26*t).*sin(32*t+0.40*pi)+1.2e-4*exp(-0.13*t).*sin(45*t+0.50*pi)+1.4e-4*exp(-0.14*t).*sin(62*t+0.15*pi)+1e-5*ones(4,2048);
y2(:,:)=1.5e-4*exp(-0.2*t).*sin(10*t+0.24*pi)+2.5e-4*exp(-0.15*t).*sin(26*t+0.40*pi)+1.57e-4*exp(-0.26*t).*sin(32*t+0.40*pi)+0.9e-4*exp(-0.13*t).*sin(45*t+0.75*pi)+1.7e-4*exp(-0.14*t).*sin(62*t+0.35*pi)+3e-6*ones(4,2048);
y3(:,:)=1.1e-4*exp(-0.2*t).*sin(10*t+0.32*pi)+1.6e-4*exp(-0.15*t).*sin(26*t+0.20*pi)+2.20e-4*exp(-0.26*t).*sin(32*t+0.60*pi)+1.4e-4*exp(-0.13*t).*sin(45*t+0.35*pi)+1.0e-4*exp(-0.14*t).*sin(62*t+0.50*pi)+1e-6*ones(4,2048);
y4(:,:)=0.9e-4*exp(-0.2*t).*sin(10*t+0.72*pi)+1.4e-4*exp(-0.15*t).*sin(26*t+0.82*pi)+1.80e-4*exp(-0.26*t).*sin(32*t+0.06*pi)+2.4e-4*exp(-0.13*t).*sin(45*t+0.35*pi)+1.3e-4*exp(-0.14*t).*sin(62*t+0.15*pi)+5e-6*ones(4,2048);
y5(:,:)=1.4e-4*exp(-0.2*t).*sin(10*t+0.32*pi)+1.0e-4*exp(-0.15*t).*sin(26*t+0.20*pi)+0.90e-4*exp(-0.26*t).*sin(32*t+0.45*pi)+1.6e-4*exp(-0.13*t).*sin(45*t+0.23*pi)+1.7e-4*exp(-0.14*t).*sin(62*t+0.80*pi)+7e-6*ones(4,2048);
D(1,:)=y1(1,:);D(2,:)=y2(1,:);D(3,:)=y3(1,:);D(4,:)=y4(1,:);D(5,:)=y5(1,:);
D(6,:)=y1(2,:);D(7,:)=y2(2,:);D(8,:)=y3(2,:);D(9,:)=y4(2,:);D(10,:)=y5(2,:);
DD(1,:)=y1(3,:);DD(2,:)=y2(3,:);DD(3,:)=y3(3,:);DD(4,:)=y4(3,:);DD(5,:)=y5(3,:);
DD(6,:)=y1(4,:);DD(7,:)=y2(4,:);DD(8,:)=y3(4,:);DD(9,:)=y4(4,:);DD(10,:)=y5(4,:);
A=((DD*D')*inv(D*D')+(DD*DD')*inv(D*DD'))/2;
=eig(A);
for j=1:M
U(j)=real(Y(2*j-1,2*j-1));
V(j)=imag(Y(2*j-1,2*j-1));
Sg(j)=-log(U(j)*U(j)+V(j)*V(j))/(2*delt);   
Wd(j)=atan(V(j)/U(j))/delt;
end
MCF=abs(F(1,1:2:9)./F(6,1:2:9))
Sg            % 阻尼识别值 0.1583    0.1566    0.2442    0.1543    0.2226
Wd            % 频率识别值 61.9866   44.9597    9.9330   26.0067   31.9731
(Sg-SSg)./SSg*100   % 误差(%) 13.0567   20.4311   22.0982    2.8976-14.4016   
(Wd-WWd)./WWd*100   % 误差(%) -0.0216   -0.0895   -0.6698    0.0257   -0.0842

bonzer 发表于 2006-12-3 16:23

谢谢主任！ 我自己编了一个频率和振型识别的还不错，但是阻尼比识别误差很大！（没加噪声）不知道为什么

pengweicai 发表于 2006-12-3 16:26

这是很正常的。
   阻尼比的识别目前都还没有非常好的方法，加之ITD法本身也不是很完善。
         想要得到更好的结果，就需要用别的识别方法。

bonzer 发表于 2006-12-3 16:36

目前什么识别方法比较好呢？ERA？

pengweicai 发表于 2006-12-3 17:06

不好意思 已经很久没有关注了。

bonzer 发表于 2006-12-3 17:06

要用ITD识别，需要自由衰减信号。
这就需要一些技术（随机减量？。。。。）对得到的随机信号进行分析

66奋斗江南 发表于 2012-8-7 10:53

高手，看到希望指教，我现在也遇到这样的问题，现在急求，谢谢
页: [1]
查看完整版本: 求 ：模态识别ITD法程序，急！感激！！