wcoherence
子波一致性和互谱
语法
描述
例子
两个正弦波的小波相干
使用默认wcoherence设置获得正弦波之间的小波相干性和随机噪声调频信号的频率降低。
t = linspace (0, 1, 1024);x = sin(8 *π* t) + 0.4 * randn (1024);x = x / max (abs (x));y = wnoise (“多普勒”10);wcoh = wcoherence (x, y);
默认的相干计算使用分析Morlet小波,12每八度的声音,抚平12尺度。默认的八度=地板(log2(元素个数(x))) 1在这种情况下,这是9。
采样间隔对小波相干性的影响
获得两个信号的小波相干数据,指定一个采样间隔为0.001秒。信号由两个正弦波(10赫兹和50 Hz)在白噪声。正弦波有不同的时间支持。金宝app
设置随机数发生器为再现性其默认设置。然后创建了两个信号。
rng默认的;t = 0:0.001:2;x = cos(2 *π* 10 * t) * (t > = 0.5 & t < 1.1) +…因为(2 *π* 50 * t)。*(t>= 0.2 & t< 1.4)+0.25*randn(size(t)); y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+…罪(2 *π* 50 * t)。*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t)); subplot(2,1,1) plot(t,x) title(“X”次要情节(2,1,2)情节(t, y)标题(“Y”)包含(的时间(秒))
获得两个信号的一致性。
[wcoh, ~,期间,coi] = wcoherence (x, y,秒(0.001));
使用pcolor命令情节连贯性和锥的影响力。
图周期=秒(时期);coi =秒(coi);h = pcolor (t, log2(时期),wcoh);h。EdgeColor =“没有”;甘氨胆酸ax =;ytick =圆(pow2 (ax.YTick), 3);ax.YTickLabel = ytick;ax.XLabel.String =“时间”;ax.YLabel.String =“时间”;ax.Title。字符串=“小波相干”;hcol = colorbar;hcol.Label。字符串=“平方一致性”;持有在;情节(ax, t, log2 (coi),“w——”,“线宽”,2)
使用wcoherence (x, y,秒(0.001))没有任何输出参数。这个情节包括阶段箭头和锥的影响力。
wcoherence (x, y,秒(0.001));
采样频率对小波相干性的影响
获得两个信号的小波相干,指定采样频率为1000赫兹。信号由两个正弦波(10赫兹和50 Hz)在白噪声。正弦波有不同的时间支持。金宝app
设置随机数发生器为再现性其默认设置,创建两个信号。
rng默认的t = 0:0.001:2;x = cos(2 *π* 10 * t) * (t > = 0.5 & t < 1.1) +…因为(2 *π* 50 * t)。*(t>= 0.2 & t< 1.4)+0.25*randn(size(t)); y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+…罪(2 *π* 50 * t)。*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t));
得到小波相干。情节连贯性对情节翻面在前面的例子中,它指定一个采样间隔采样频率。
wcoherence (x, y, 1000)
获得scale-to-frequency转换输出f。
[wcoh wcs f] = wcoherence (x, y, 1000);
的平滑尺度小波相干性的影响
获得两个信号小波相干。信号由两个正弦波(10赫兹和50 Hz)在白噪声。使用默认的鳞片光滑。这个值的数量相当于每八度的声音。两个值默认为12。
设置随机数发生器为再现性其默认设置。然后,创建两个信号,获得一致性。
rng默认的;t = 0:0.001:2;x = cos(2 *π* 10 * t) * (t > = 0.5 & t < 1.1) +…因为(2 *π* 50 * t)。*(t>= 0.2 & t< 1.4)+0.25*randn(size(t)); y = sin(2*pi*10*t).*(t>=0.6 & t<1.2)+…罪(2 *π* 50 * t)。*(t>= 0.4 & t<1.6)+ 0.35*randn(size(t)); wcoherence(x,y)
天平光滑的数量设置为18。平滑的增加会导致低频率分辨率降低。
wcoherence (x, y,“NumScalesToSmooth”,18)
阶段显示阈值对小波的影响天气数据的一致性
比较使用不同阶段的影响显示阈值小波相干。
情节厄尔尼诺现象之间的小波相干时间序列和所有印度平均降雨量指数。每月的数据采样。指定采样间隔1/12年多年来显示时间。使用默认阶段显示阈值为0.5,显示阶段的箭只在相干大于或等于0.5。
负载ninoairdata;wcoherence(尼诺、空气、年(1/12));
将相位显示阈值设置为0.7。阶段的箭的数量减少。
wcoherence(尼诺、空气、年(1/12),“PhaseDisplayThreshold”,0.7);
输入参数
输出参数
更多关于
兼容性的考虑
引用
[1]Grinsted, J。,C。Moore, and S. Jevrejeva. “Application of the cross wavelet transform and wavelet coherence to geophysical time series.”在地球物理非线性过程。11卷,期5/6,2004年,页561 - 566。
[2]Maraun D。,J. Kurths, and M. Holschneider. "Nonstationary Gaussian processes in wavelet domain: Synthesis, estimation and significance testing.”75年《物理评论。2007年,页016707-1-016707-14。
[3]一如,C。,P. Webster. "Interdecadal changes in the ESNO-Monsoon System."杂志的气候。12卷,1999年,页2679 - 2690。
