关于广义互相关求时延的问题

野菊

TestGuru 发表于 2016-10-18 12:24
用互相关做时延测量做过几种，比如音响要调节各方位距离远近不同的喇叭到听者处的时间延迟等，还有几个相 ...

感谢提供资料

lebronze

TestGuru 发表于 2016-10-18 12:24
用互相关做时延测量做过几种，比如音响要调节各方位距离远近不同的喇叭到听者处的时间延迟等，还有几个相 ...

好吧，我再找找问题

Posion

楼主问题研究咋样了？

lebronze

Posion 发表于 2016-10-20 09:01
楼主问题研究咋样了？

没解决，反而发现更多问题了。
比如：在求出互相关函数后，求峰值是求数据的最大值还是数据绝对值得最大值？
请问你有什么经验没

zghitboy

楼主问题解决了吗

xiujuan507

楼主，能交流一下吗？我最近也在做这个

eastar

xiujuan507 发表于 2016-12-27 09:24
楼主，能交流一下吗？我最近也在做这个

你哪里出现问题了

xiujuan507

eastar 发表于 2016-12-27 09:28
你哪里出现问题了

看了帖子，有用xcorr,FFT的，有的人用小波变换的，不知道小波变换该从何入手？

怪咖先生

TestGuru 发表于 2016-10-18 12:24
用互相关做时延测量做过几种，比如音响要调节各方位距离远近不同的喇叭到听者处的时间延迟等，还有几个相 ...

挺好的视频

林夕sq123

lebronze 发表于 2016-10-17 17:44
多谢前辈回复。
1.我目前已经确定用广义互相关（GCC）做时延估计了，时域就不考虑了
2.频域互相关要先 ...

您好 能给一些广义互相关的程序吗、最近在做这个

sd

1. 
   
2. %t=0.001:0.001:1;
   
3. t=1:1000;
   
4. s1=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));
   
5. 
   
6. for t=1:500;
   
7. s2(t)=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));
   
8. end
   
9. 
   
10. for t=501:1000;
    
11. s2(t)=sin(2*pi*200*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));
    
12. end
    
13. 
    
14. 
    
15. subplot(9,2,1)
    
16. plot(s1)
    
17. title('原始信号')
    
18. ylabel('S1')
    
19. 
    
20. subplot(9,2,2)
    
21. plot(s2)
    
22. title('故障信号')
    
23. ylabel('S2')
    
24. 
    
25. wpt=wpdec(s1,3,'db1','shannon');
    
26. 
    
27. %plot(wpt);
    
28. s130=wprcoef(wpt,[3,0]);
    
29. s131=wprcoef(wpt,[3,1]);
    
30. s132=wprcoef(wpt,[3,2]);
    
31. s133=wprcoef(wpt,[3,3]);
    
32. s134=wprcoef(wpt,[3,4]);
    
33. s135=wprcoef(wpt,[3,5]);
    
34. s136=wprcoef(wpt,[3,6]);
    
35. s137=wprcoef(wpt,[3,7]);
    
36. s10=norm(s130);
    
37. s11=norm(s131);
    
38. s12=norm(s132);
    
39. s13=norm(s133);
    
40. s14=norm(s134);
    
41. s15=norm(s135);
    
42. s16=norm(s136);
    
43. s17=norm(s137);
    
44. 
    
45. st10=std(s130);
    
46. st11=std(s131);
    
47. st12=std(s132);
    
48. st13=std(s133);
    
49. st14=std(s134);
    
50. st15=std(s135);
    
51. st16=std(s136);
    
52. st17=std(s137);
    
53. 
    
54. disp('正常信号的特征向量');
    
55. snorm1=[s10,s11,s12,s13,s14,s15,s16,s17]
    
56. std1=[st10,st11,st12,st13,st14,st15,st16,st17]
    
57. 
    
58. subplot(9,2,3);plot(s130);
    
59. ylabel('S130');
    
60. subplot(9,2,5);plot(s131);
    
61. ylabel('S131');
    
62. subplot(9,2,7);plot(s132);
    
63. ylabel('S132');
    
64. subplot(9,2,9);plot(s133);
    
65. ylabel('S133');
    
66. subplot(9,2,11);plot(s134);
    
67. ylabel('S134');
    
68. subplot(9,2,13);plot(s135);
    
69. ylabel('S135');
    
70. subplot(9,2,15);plot(s136);
    
71. ylabel('S136');
    
72. subplot(9,2,17);plot(s137);
    
73. ylabel('S137');
    
74. 
    
75. 
    
76. wpt=wpdec(s2,3,'db1','shannon');
    
77. 
    
78. %plot(wpt);
    
79. s230=wprcoef(wpt,[3,0]);
    
80. s231=wprcoef(wpt,[3,1]);
    
81. s232=wprcoef(wpt,[3,2]);
    
82. s233=wprcoef(wpt,[3,3]);
    
83. s234=wprcoef(wpt,[3,4]);
    
84. s235=wprcoef(wpt,[3,5]);
    
85. s236=wprcoef(wpt,[3,6]);
    
86. s237=wprcoef(wpt,[3,7]);
    
87. 
    
88. s20=norm(s230);
    
89. s21=norm(s231);
    
90. s22=norm(s232);
    
91. s23=norm(s233);
    
92. s24=norm(s234);
    
93. s25=norm(s235);
    
94. s26=norm(s236);
    
95. s27=norm(s237);
    
96. 
    
97. st20=std(s230);
    
98. st21=std(s231);
    
99. st22=std(s232);
    
100. st23=std(s233);
     
101. st24=std(s234);
     
102. st25=std(s235);
     
103. st26=std(s236);
     
104. st27=std(s237);
     
105. 
     
106. disp('故障信号的特征向量');
     
107. snorm2=[s20,s21,s22,s23,s24,s25,s26,s27]
     
108. std2=[st20,st21,st22,st23,st24,st25,st26,st27]
     
109. 
     
110. 
     
111. subplot(9,2,4);plot(s230);
     
112. ylabel('S230');
     
113. subplot(9,2,6);plot(s231);
     
114. ylabel('S231');
     
115. subplot(9,2,8);plot(s232);
     
116. ylabel('S232');
     
117. subplot(9,2,10);plot(s233);
     
118. ylabel('S233');
     
119. subplot(9,2,12);plot(s234);
     
120. ylabel('S234');
     
121. subplot(9,2,14);plot(s235);
     
122. ylabel('S235');
     
123. subplot(9,2,16);plot(s236);
     
124. ylabel('S236');
     
125. subplot(9,2,18);plot(s237);
     
126. ylabel('S237');
     
127. 
     
128. %fft
     
129. 
     
130. figure
     
131. y1=fft(s1,1024);
     
132. py1=y1.*conj(y1)/1024;
     
133. y2=fft(s2,1024);
     
134. py2=y2.*conj(y2)/1024;
     
135. 
     
136. y130=fft(s130,1024);
     
137. py130=y130.*conj(y130)/1024;
     
138. y131=fft(s131,1024);
     
139. py131=y131.*conj(y131)/1024;
     
140. y132=fft(s132,1024);
     
141. py132=y132.*conj(y132)/1024;
     
142. y133=fft(s133,1024);
     
143. py133=y133.*conj(y133)/1024;
     
144. y134=fft(s134,1024);
     
145. py134=y134.*conj(y134)/1024;
     
146. y135=fft(s135,1024);
     
147. py135=y135.*conj(y135)/1024;
     
148. y136=fft(s136,1024);
     
149. py136=y136.*conj(y136)/1024;
     
150. y137=fft(s137,1024);
     
151. py137=y137.*conj(y137)/1024;
     
152. 
     
153. y230=fft(s230,1024);
     
154. py230=y230.*conj(y230)/1024;
     
155. y231=fft(s231,1024);
     
156. py231=y231.*conj(y231)/1024;
     
157. y232=fft(s232,1024);
     
158. py232=y232.*conj(y232)/1024;
     
159. y233=fft(s233,1024);
     
160. py233=y233.*conj(y233)/1024;
     
161. y234=fft(s234,1024);
     
162. py234=y234.*conj(y234)/1024;
     
163. y235=fft(s235,1024);
     
164. py235=y235.*conj(y235)/1024;
     
165. y236=fft(s236,1024);
     
166. py236=y236.*conj(y236)/1024;
     
167. y237=fft(s237,1024);
     
168. py237=y237.*conj(y237)/1024;
     
169. 
     
170. f=1000*(0:511)/1024;
     
171. subplot(1,2,1);
     
172. plot(f,py1(1:512));
     
173. ylabel('P1');
     
174. title('原始信号的功率谱')
     
175. subplot(1,2,2);
     
176. plot(f,py2(1:512));
     
177. ylabel('P2');
     
178. title('故障信号的功率谱')
     
179. 
     
180. figure
     
181. 
     
182. subplot(4,2,1);
     
183. plot(f,py130(1:512));
     
184. ylabel('P130');
     
185. title('S130的功率谱')
     
186. subplot(4,2,2);
     
187. plot(f,py131(1:512));
     
188. ylabel('P131');
     
189. title('S131的功率谱')
     
190. subplot(4,2,3);
     
191. plot(f,py132(1:512));
     
192. ylabel('P132');
     
193. subplot(4,2,4);
     
194. plot(f,py133(1:512));
     
195. ylabel('P133');
     
196. subplot(4,2,5);
     
197. plot(f,py134(1:512));
     
198. ylabel('P134');
     
199. subplot(4,2,6);
     
200. plot(f,py135(1:512));
     
201. ylabel('P135');
     
202. subplot(4,2,7);
     
203. plot(f,py136(1:512));
     
204. ylabel('P136');
     
205. subplot(4,2,8);
     
206. plot(f,py137(1:512));
     
207. ylabel('P137');
     
208. figure
     
209. 
     
210. subplot(4,2,1);
     
211. plot(f,py230(1:512));
     
212. ylabel('P230');
     
213. title('S230的功率谱')
     
214. subplot(4,2,2);
     
215. plot(f,py231(1:512));
     
216. ylabel('P231');
     
217. title('S231的功率谱')
     
218. subplot(4,2,3);
     
219. plot(f,py232(1:512));
     
220. ylabel('P232');
     
221. subplot(4,2,4);
     
222. plot(f,py233(1:512));
     
223. ylabel('P233');
     
224. subplot(4,2,5);
     
225. plot(f,py234(1:512));
     
226. ylabel('P234');
     
227. subplot(4,2,6);
     
228. plot(f,py235(1:512));
     
229. ylabel('P235');
     
230. subplot(4,2,7);
     
231. plot(f,py236(1:512));
     
232. ylabel('P236');
     
233. subplot(4,2,8);
     
234. plot(f,py237(1:512));
     
235. ylabel('P237');
     
236. figure
     
237. plottree(wpt)
     
238. 
     
239. 
     

复制代码