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) の関数として相互相関のプロットを取得する方法は?
4

1 に答える 1

0

[X]はXの平均であると言います。

R（dt）= [（A（t + dt）-[A（t）]）*（B（t）-[B（t）]）] /（std（A）* std（B））

分子は2つの関数A（t）とB（t）を取り、それらを中央に配置します。これは、平均を減算することです。これを少し簡単に見せるための2つの新しい関数を作成できます。

C（t）= A（t）-[A（t）] D（t）= B（t）-[B（t）]

それから：

R（dt）= [C（t + dt）* D（t）] /（std（C）* std（D））

今、私はこれら2つの機能を組み合わせて、平均を取っています。これらの2つの関数が互いに正確に移動している場合、これは1になり、正の相関を意味します。

次のようにして、Matlabでこれを行うことができます。

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

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

これにより、dt = 0でのクロスコアが得られます。別のdtでは、次のようになります。

x2 = mean((a(2:50) - mean(a(2:50))) .* (b(1:49) - mean(b(1:49))))/(std(a) * std(b))
于 2012-10-15T23:15:21.167 に答える