• ベストアンサー

相互相関関数の問題について教えてください

一対の定常過程X(t),Y(t)を考える。 相互相関関数Rxy(τ)、Ryx(τ)は以下の性質を持つ。 (A)Rxy(τ)=Ryx(τ) (B)|Rxy(τ)|≦1/2[Rx(0)+Ry(0)] Rx(τ)とRy(τ)が、X(t),Y(t)の相互相関関数となるのはどのようなときか? よろしくお願いします。

質問者が選んだベストアンサー

  • ベストアンサー
noname#221368
noname#221368
回答No.1

「ここにすぐ書いて出せ!」と言われたら、自宅へノートを取りに逃げ帰りますが、たしか相互相関関数の定義にまで遡って地道に式を見ていけば、できるような気がします。

trot_hero
質問者

お礼

ありがとうございました。 相互相関関数の定義 Rxy=〈X(t),Y(t)〉         1 +k    =lim ―∫  X(t)・Y(t+τ)dt     k→∞ 2k -k を使って地道にといてみます。

関連するQ&A

  • 相互相関関数 使い方

    x(t)、y(t)の信号があるとき、互いの距離は22cmであり、行路差はLとする。行路差Lを求めたい。 相互相関関数は Φxy(t1,t2)=E[x(t1),y(t1)] 定常性、エルゴード性より、 Φxy(τ)=x(t)y(t+τ)の平均     =lim(T→∞)1/T*∫(0からT)x(t)y(t+τ)dt となると習いました。 結果は(0からτまで5.8msなので)L=13.9cm だそうです。 どのような計算を行ったのか”おそらく”でもいいので教えていただきたいです! よろしくお願いします。。

  • 【数学】偏相関係数の式の導出過程を教えてください

    偏相関係数の定義式の導出過程を教えてください。 rxy・z=(rxy-rxz*ryz)/(√1-r^2xz)*(√1-r^2yz) rxy=xとyの相関係数 ryz=yとzの相関係数 rxz=xとzの相関係数 rxy・z=zの影響を取り除いたxとyの偏相関係数

  • 自己相関関数について

     exp(-(x/a)^2)の自己相関関数の相関値が1から1/eまで落ちる長さは2aだと思います。  exp(-(x^2+y^2)/a^2)の自己相関関数の相関値が1から1/eまで落ちる長さはaを使ってどのようにあらわされるのでしょうか?感覚的には2aよりも小さい気がするのですが、、 よろしくお願いします。

  • sin、cosの相互相関係数

    f(t) = sin(t)、g(t) = cos(t) とするとき f(t)、g(t)の相互相関係数を求めたいのですが E{x(t)y(t)} = ∫sin(t)cos(t)dt = (1/2)∫sin(2t)dt (t:0→2π) を計算すると0になりました これが分子にくるので、結果0が答えとなるのですが 周期をずらせば一致する関数なので これはおかしいと思っています 正しい解答を教えてください よろしくお願いします 分母は E{x(t)^2} = π E{y(t)^2} = π となりましたので √π * √πでπとなりました

  • MATLAB 相互相関

    X = zeros(200, 200,3); %%% RGB plane X(95:105, 95:105, 1) = 255; X(95:105, 95:105, 2) = 255; X(95:105, 95:105, 3) = 255; figure imshow(X,[]); Y = zeros(100, 100,3); %%% RGB plane Y(45:55, 45:55, 1) = 255; Y(45:55, 45:55, 2) = 255; Y(45:55, 45:55, 3) = 255; figure imshow(Y,[]); このあと、XとYの相互相関をnormxcorr2でグラフ化したいのですが・・ 色々なサイトを見て、 c=normxcorr2(X,Y)や、c=normxcorr2(X(:,:,1),Y(:,:,1))などを試したのですがエラーが出てしまいます。 何かこれの前に打たなければいけないプログラムがあるのでしょうか・・・

  • 自己相関関数について

    自己相関関数について x(t)= 1(0<=t<D), 0(t<0, t>D) の自己相関関数を次のように求めたのですがあっているかどうかをチェックお願いします。 (1) 0<=τ<=D のとき R(τ)=D-τ (2) D<τ のとき x(t+τ)とx(t)の重なりがないため0 (3) τ<0 のとき R(τ)=D-τ まとめると R(τ)= D-τ |τ|<=D 0 |τ| > D

  • 自己相関関数とパワースペクトル密度関数、フーリエ変換について。

    自己相関関数とパワースペクトル密度関数、フーリエ変換について。 パワースペクトル、パワースペクトル密度と自己相関関数についての質問です。 (tは時間、hは次数、fは周波数として) ある信号x(t)の自己相関関数r(h)をフーリエ変換すると、その信号のパワースペクトル密度関数p(f)になるとネットにあったのですが、パワースペクトル密度関数p(f)と、信号x(t)をそのままフーリエ変換して得たパワースペクトルX(f)はどう違うんでしょうか。 ちなみに数学的な話というよりはコンピュータ上の処理(離散値)で考えています。 もともとパワースペクトルが『自己相関関数の離散フーリエ変換として定義される』と本にはあったのを読みました。 しかし同じ本の中に、『自己相関関数のフーリエ変換は正しくはピリオドグラムと言い、パワースペクトルとはピリオドグラムの平均値で求められる』とも書いてありました。 パワースペクトルとパワースペクトル密度関数はいったいどう違うのか…?とずっと考えているのですが分かりません。 あと(自己、相互)相関関数と(自己、相互)相関係数にはどのような関係があるのですか。回答よろしくお願いします。 前回1つ回答頂いたんですが解決できなかったのですみません、もう一度お願いします。

  • 相互相関関数

    相互相関関数を求めるプログラムを書いたのですが、入力信号のずらし幅が負の時と正の時で処理を分けた場合、結果の波形の整合性が取れずに困っております。計算上の根本的違いやコードのミスも含めて、ご指南頂きたいです。 ちなみに参考にしたのは、 http://wwwdoi.elec.nara-k.ac.jp/html/jisyu/j5ec02/taro12-sjt0_p13_2... の最後の2ページで、入力のサンプルデータは、 http://www.mech.tohoku-gakuin.ac.jp/nken/java/new2/crosscorr7/cross... のデータを用いました。結果画像を添付します。 #include "stdafx.h" #include <stdio.h> #include <math.h> #include <stdlib.h> #define MAX 10000000 #define buf 256 int main (){ FILE *fpi; //_入力ファイル FILE *fpo; //_出力ファイル int num; //_対象波形全点数 double *x,*y; //_波形データ double *cxy1,*cxy2;//結果 char filename[buf]; /////////ファイルオープン/////////// (省略) ////////////データ読み込み/////////////////// if(((x=(double*)malloc(MAX))==NULL)||((y=(double*)malloc(MAX))==NULL)){ fprintf(stderr,"Can't_allocate_memory.\n"); return 0; } num=0; while ( fscanf(fpi,"%lf,%lf",&x[num],&y[num]) != EOF && num < MAX ){ num++; } printf("%d",num); ////////相互相関関数////////////////////////// if(((cxy1=(double*)malloc(MAX))==NULL)||((cxy2=(double*)malloc(MAX))==NULL)){ fprintf(stderr,"Can't_allocate_memory.\n"); return 0; } int i,j,n; double sum; n = num-1; //////////////ずらし幅が負の時の演算//////////// for( i=0; i<n; i++ ) { sum = 0.0; for( j=0; j < n-i; j++ ) { //k = (j + i) % num; sum = sum + x[j] * y[j+i]; } cxy1[i]=double(sum/(n-i)); } ////ずらし幅が正の時の演算/// n = num-1; for( i=0; i<n; i++ ) { sum = 0.0; for( j=0; j < n-i; j++ ) { //k = (j + i) % num; sum = sum + x[i+j] * y[j]; } cxy2[i]=double(sum/(n-i)); } /////////結果の出力/////////////// for(i=0;i<n;i++){ fprintf(fpo,"%f\n",cxy1[i]); } for(i=0;i<n;i++){ fprintf(fpo,"%f\n",cxy2[i]); } ///////////////////////// fclose(fpi); fclose(fpo); free(x); free(y); free(cxy1); free(cxy2); }

  • 統計

    統計 データ(x1,y1),(x2,y2)・・・・,(xn,yn) に大して標本共分散Sxyと標本相関係数Rxy があるとき -1≦Rxy≦1の示し方と、 Yj=aXj+b(a≠0)のとき|Rxy|=1の示し方を教えてください。

  • 自己相関関数の解き方を教えて下さい。

    x(t)=Aexp(-αt) [αは、正]の自己相関関数を0<t<Tの区間で求めよ。τ<<Tとする。 R(τ)=x(t)x(t-τ) =Aexp(-αt)・Aexp(-α(t-τ)) τ<<tなので、 A/T∫T→0[exp(-αt)・exp(-α(t-τ))]dt = 何となく解いてみましたが、誰か、回答を詳しく説明お願いします。