ITD模态识别怎么编程

66奋斗江南

急求ITD模态识别怎么编程？请高手请教，谢谢

勤奋的wwf

高手，你搞定了没？等着你的回信

chybeyond

勤奋的wwf 发表于 2015-4-6 15:44
高手，你搞定了没？等着你的回信

1. 网上程序，可作参考
   
2. %用ITD方法对随机减量后数据进行整体模态参数识别
   
3. clear
   
4. clc
   
5. close all hidden
   
6. for i=1:17   %共17组数据
   
7. fni=input('ITD法模态参数识别－输入数据文件名：','s');
   
8. fid1=fopen(fni,'r');
   
9. x=fscanf(fid1,'%f',[8,inf]);%每组数据8个通道（有两个参考点）
   
10. y(1+8*(i-1):8+8*(i-1),:)=x;%将各组数据逐一存入同一矩阵
    
11. status=fclose(fid1);
    
12. end
    
13. y=y';
    
14. %根据实际情况对数据矩阵重新排列
    
15. z(:,1)=y(:,11);
    
16. z(:,2)=y(:,13);
    
17. z(:,3)=y(:,27);
    
18. z(:,4)=y(:,45);
    
19. z(:,5)=y(:,49);
    
20. z(:,6)=y(:,61);
    
21. z(:,7)=y(:,69);
    
22. z(:,8)=y(:,77);
    
23. z(:,9)=y(:,93);
    
24. z(:,10)=y(:,101);
    
25. z(:,11)=y(:,113);
    
26. z(:,12)=y(:,117);
    
27. z(:,13)=y(:,129);
    
28. z(:,14)=y(:,3);
    
29. z(:,15)=y(:,21);
    
30. z(:,16)=y(:,29);
    
31. z(:,17)=y(:,33);
    
32. z(:,18)=y(:,41);
    
33. z(:,19)=y(:,57);
    
34. z(:,20)=y(:,65);
    
35. z(:,21)=y(:,81);
    
36. z(:,22)=y(:,89);
    
37. z(:,23)=y(:,97);
    
38. z(:,24)=y(:,105);
    
39. z(:,25)=y(:,121);
    
40. z(:,26)=y(:,119);
    
41. z(:,27)=y(:,1);
    
42. z(:,28)=y(:,9);
    
43. z(:,29)=y(:,19);
    
44. z(:,30)=y(:,37);
    
45. z(:,31)=y(:,51);
    
46. z(:,32)=y(:,53);
    
47. z(:,33)=y(:,67);
    
48. z(:,34)=y(:,75);
    
49. z(:,35)=y(:,85);
    
50. z(:,36)=y(:,109);
    
51. z(:,37)=y(:,115);
    
52. z(:,38)=y(:,123);
    
53. z(:,39)=y(:,133);
    
54. z(:,40)=y(:,5);
    
55. z(:,41)=y(:,17);
    
56. z(:,42)=y(:,25);
    
57. z(:,43)=y(:,35);
    
58. z(:,44)=y(:,43);
    
59. z(:,45)=y(:,59);
    
60. z(:,46)=y(:,73);
    
61. z(:,47)=y(:,83);
    
62. z(:,48)=y(:,91);
    
63. z(:,49)=y(:,99);
    
64. z(:,50)=y(:,107);
    
65. z(:,51)=y(:,125);
    
66. z(:,52)=y(:,131);
    
67. %ITD法
    
68. h=z';
    
69. n=length(h(1,:));
    
70. n1=n-2;
    
71. mn=26; %模态阶数
    
72. sf=100; %采样频率
    
73. nm=2*mn; %建立特征方程阶数
    
74. dt=1/sf;%计算时间间隔
    
75. t=0:dt:(n1-1)*dt;%建立离散时间向量
    
76. %计算自由振动响应矩阵
    
77. for k=1:nm
    
78.     x0(k,:)=h(k,1:n1);
    
79.     x2(k,:)=h(k,2:n1+1);
    
80. end
    
81. b2=x2/x0;
    
82. [A,V]=eig(b2);%计算特征值和特征向量
    
83. %变换特征值对角阵为一向量
    
84. for k=1:nm
    
85.     u(k)=V(k,k);
    
86. end
    
87. f1=abs(log(u'))./(2*pi*dt);%计算模态频率向量
    
88. d1=sqrt(1./(((imag(log(u'))./real(log(u'))).^2)+1));%计算阻尼比向量
    
89. %计算振型系数向量
    
90. for i=1:nm
    
91.     A1(:,i)=A(:,i)/log(u(i));
    
92. end
    
93. s=imag(A1)*10;
    
94. u=u';
    
95. for k=0:n1-1
    
96. va(:,k+1)=u.^k;
    
97. end
    
98. h1=real(A*va);
    
99. for k0=1:13
    
100.     figure(k0)
     
101.     for i0=4*k0-3:4*k0
     
102.         for n0=1:4
     
103.     subplot(2,2,n0)
     
104. plot(t,x0((k0-1)*4+n0,:),':',t,h1((k0-1)*4+n0,:));
     
105.         grid on
     
106.         end
     
107.     end
     
108. end

复制代码

勤奋的wwf

chybeyond 发表于 2015-4-6 15:52

大神，我现在有频率-加速度的频响曲线数据，问能否识别出其模态质量，刚度，阻尼比

chybeyond

勤奋的wwf 发表于 2015-4-6 16:23
大神，我现在有频率-加速度的频响曲线数据，问能否识别出其模态质量，刚度，阻尼比

未研究过相关领域，你可以百度一下相关文献资料

勤奋的wwf

chybeyond 发表于 2015-4-6 16:47
未研究过相关领域，你可以百度一下相关文献资料

我查过文献，里面都是一带而过，具体怎么实现的都没有