cwt関数の結果から​、振幅の最大値を抽出​する方法を教えて下さ​い。

8 次查看（过去 30 天）
勝也 秋元
勝也 秋元 2022-8-25
评论： 勝也 秋元 2022-8-29
＜実施したい内容＞
cwt関数の結果に対し、探索する領域（時間、周波数）を設定し、その範囲で振幅最大値を抽出する。
振幅最大値と時間、周波数をアウトプットにする。
添付の例）
探索区間：時間：１～２sec、周波数：4Hz以上
アウトプット（時間、周波数、最大振幅）：（1.22、14.32、12.1075）

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

回答（1 个）

takemoto
takemoto 2022-8-29
Ran in:
load mtlb
[wt,f] = cwt(mtlb,Fs);
t = 0 : 1/Fs : numel(mtlb)/Fs - 1/Fs;
tidx = 100:2000; % ROI (time)
fidx = (20:40)'; % ROI (freq)
wtROI = wt(fidx,tidx); % ROI
fROI = f(fidx);
tROI = t(tidx);
[maxwt,idx] = max(abs(wtROI(:)));
[ii, jj] = ind2sub(size(wtROI), idx);
surf(tROI,fROI,abs(wtROI).^2), view(0,90)
shading interp, hold on
plot3(tROI(jj),fROI(ii),maxwt,'ro')
text(tROI(jj)+0.01,fROI(ii),maxwt, ...
{['time: ',num2str(tROI(jj))], ...
['freq: ',num2str(fROI(ii))], ...
['mag: ',num2str(maxwt)]})
hold off
cwt関数は、出力引数を与えることで、変換結果を獲得できますので、max関数等で最大値を抽出してはいかがでしょうか？
  1 个评论
勝也 秋元
勝也 秋元 2022-8-29
ありがとうございます。
内容理解できました。

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 連続ウェーブレット変換 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!