wcoherenceで作成されたグラフについて

2 ビュー (過去 30 日間)
NAOHIRO KIMURA
NAOHIRO KIMURA 2020 年 2 月 12 日
回答済み: Naoya 2020 年 2 月 12 日
wcoherenceで作成されたグラフについて教えて下さい。
wcoherenceについて、理解を深める目的で以下のような例題を自作してwcoherenceでグラフを
作成してみました。
clc;clear;
t=0:0.1:300;
z1=sin(2*pi*t/50);
z2=sin(2*pi*t/50);
x1=sin(pi*t/50);
x2=sin(pi*t/50);
y1=z1+x1;
y2=z2+x2;
subplot(2,1,1)
plot(t,y1);
grid on
subplot(2,1,2)
plot(t,y2,'r');
grid on
y1とy2.jpg
figure;
wcoherence(y1,y2,10);
ウェーブレットコヒーレンス.jpg
すると、最初(時刻0)から最後(時刻5分)まで一様に黄色で、右向きの→となっている
グラフが出現しました。
私は、周波数が0.01Hzと0.02Hzの部分に黄色い帯が出来るのでは無いかと考えていました。
なぜならば、y1もy2も0.01Hzと0.02Hzの正弦波を足し合わせて出来たものだったからです。
どうして全体が黄色い表示となるのかを教えて下さい。
よろしくお願い致します。

採用された回答

Naoya
Naoya 2020 年 2 月 12 日
コヒーレンス対象の2種類の信号が全く同一ということで、全ての周波数において同じ特性を持っているとみなされ、コヒーレンス強度が一定になっていると思われます。
例えば、下記コードのようにホワイトノイズを加えることにより、信号に含まれる正弦波成分以外は相関がないとみなすことができ、相関がある周波数と相関がない周波数の色分け表示が確認できると思います。
t=0:1/1000:2;
z1=sin(2*pi*t*64);
z2=sin(2*pi*t*64);
x1=sin(2*pi*t*8);
x2=sin(2*pi*t*8);
y1=z1+x1+1.0*randn(1,2001);
y2=z2+x2+1.0*randn(1,2001);
subplot(2,1,1)
plot(t,y1);
grid on
subplot(2,1,2)
plot(t,y2,'r');
grid on
figure;
wcoherence(y1,y2,1000);
wcohre.png

その他の回答 (0 件)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!