wavelet波形を3次元にする際，X軸を時間に変更する方法を教えてください。

Question

Ryosuke Takahashi 2017-7-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/348334-wavelet-3-x

编辑： Ryosuke Takahashi 2017-7-13

采纳的回答： Tohru Kikawada

wavelet波形を3次元にする際，X軸を時間に変更する方法を教えていただけないでしょうか？

現在，sin波のモデル波形に対してwavelet変換を行っています。 CWTcoeffsを使用して2次元の図を作成するまではx軸を時間に変更できるのですが， 3次元にする際には規定をする必要があるのでしょうか？

fs=1000; t = 0:1/fs:2-1/fs;

S1 = sin(2*pi*20*t)+sin(2*pi*100*t);

time = (1:fs*2)/fs;

figure(1);

subplot(2,1,1);

plot(time,S1);

subplot(2,1,2);

wname = 'morl';

fc = centfrq(wname); % 中心周波数

fa = 1:200; % 擬似周波数(Hz)

scal2frq(200,wname,1/fs)%周波数分解能

sf = fc./(fa.*1/fs); % スケールファクタ

[CWTcoeffs,frq] = cwt(S1,sf,wname,1/fs);

abs_CWT = abs(CWTcoeffs);

imagesc(time,fa,abs_CWT);

colormap(jet)

axis xy

axis([0,inf, -inf, 200])

title('Scalogram')

ylabel('Hz')

colorbar

figure(2);

mesh(abs(CWTcoeffs));

colormap(jet)

colorbar

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Tohru Kikawada 2017-7-12

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/348334-wavelet-3-x#answer_273812

在 MATLAB Online 中打开

X軸の表示を時間にしたいということでしょうか。

表示だけであれば mesh に軸のラベルデータを指定することで実現できます。

mesh(time,fa,abs(CWTcoeffs));

3 个评论
显示 1更早的评论隐藏 1更早的评论

Tohru Kikawada 2017-7-13

在 MATLAB Online 中打开

掲載されているコードを実行したところ、 time 変数のサイズは 1x2000 のようです。

>> size(time)
ans =
           1        2000

単純に mesh の引数を変えただけで私の環境では表示できました。

サンプリング周波数も1000Hzとありますが何か条件が異なっていますか。

fs=1000; t = 0:1/fs:2-1/fs;
S1 = sin(2*pi*20*t)+sin(2*pi*100*t);
time = (1:fs*2)/fs;
figure(1);
subplot(2,1,1);
plot(time,S1);
subplot(2,1,2);
wname = 'morl';
fc = centfrq(wname); % 中心周波数
fa = 1:200; % 擬似周波数(Hz)
scal2frq(200,wname,1/fs)%周波数分解能
sf = fc./(fa.*1/fs); % スケールファクタ
[CWTcoeffs,frq] = cwt(S1,sf,wname,1/fs);
abs_CWT = abs(CWTcoeffs);
imagesc(time,fa,abs_CWT);
colormap(jet)
axis xy
axis([0,inf, -inf, 200])
title('Scalogram')
ylabel('Hz')
colorbar
figure(2);
mesh(time,fa,abs(CWTcoeffs));
colormap(jet)
colorbar