2012-10-15 79 views
1

有2個信號A(t),B(t)表示輸出功率。 A是時間延遲輸出。如何實現這個公式計算時移信號的互相關

Cross_cor(delay)=[A(t + delay) − <A(t)>]*[B(t) − <B(t)>]/{[A(t)-<A(t)]^2 * [B(t)-<B(t)]^2} ^2 
  1. 我不明白分子意味着什麼?
  2. 如何繪製彩色地圖或類似的衝浪地塊,從而可以用不同的顏色看到不同的地區。
  3. 如何獲得作爲延遲偏移時間(延遲X軸,Y軸上的Cross_corr)的函數的互相關圖,使峯可以可視化?
+0

那你嘗試做?你爲什麼認爲這是一個3-D信號?它只是一個需要繪製的1-D信號(延遲是唯一改變的)。你有沒有嘗試過'doc xcorr'? – bla

+0

我知道xcorr內置函數。但是這個函數只會給出A,B之間的相關性,而我想觀察(a)作爲延遲(b)的函數的相關性,使得峯值指示存在強相關性。如何將延遲參數合併到xcorr中?關於3D事物,有幾個影響相關性的參數。說,平均相關維度vs延遲vs噪聲的圖。 – Chaitali

+0

查看本維基百科文章中的「時間序列分析」條目:http://en.wikipedia.org/wiki/Cross-correlation。互相關通常用於從統計角度來看靜止的數據。如果將互相關應用於平穩信號(具有相關噪聲),則互相關僅取決於方程的延遲項。在這種情況下,您可以獲得二維圖,延遲與互相關。 我從來沒有處理過不穩定的信號。但是,如果你不得不做一個三維圖,我想這將是:(x)t,(y)延遲(z)互相關。 – jespestana

回答

0

我會說[X]是X的,則平均:

R(DT)= [(A(T + DT) - [A(T)])*( B(t) - [B(t)])] /(std(A)* std(B))

分子取兩個函數A(t)和B(t)這是什麼減去平均值。我可以創建兩個新的功能,使這一切看起來有點容易:

C(T)= A(T) - [A(T)] d(T)= B(T) - [ B(t)的]

然後:

R(DT)= [C(T + DT)* d(T)] /(標準(C)* STD(d))

現在我正在研究這兩個函數,並採取意思。如果這兩個函數完全相互移動,則這將是1,這意味着正相關。

你可以做這樣做在Matlab如下:

a = randn(1,50); 
b = randn(1,50); 

x1 = mean((a - mean(a)) .* (b - mean(b)))/(std(a) * std(b)) 

這給我在DT = 0的crosscor,另一個DT做:

x2 = mean((a(2:50) - mean(a(2:50))) .* (b(1:49) - mean(b(1:49))))/(std(a) * std(b))