生体信号のwavelet解析について

12 次查看(过去 30 天)
cho hunseki
cho hunseki 2023-9-22
現在とある生体信号のwavelet解析を行おうとしています。
データはcsvデータで、1行n列のものです。
下記codeを試しました。
data = readmatrix('ファイル名.csv'); %
fs = 1000; %
wavelet_name = 'morl'; %
scales = 1:64; %
cwt_coeffs = cwt(data, 'wavelet', wavelet_name, 'scales', scales, 'samplingfrequency', fs);
figure;
imagesc((1:size(data, 2))/fs, scales, abs(cwt_coeffs));
set(gca, 'YDir', 'normal');
title('ウェーブレット解析');
xlabel('時間(秒)');
ylabel('ウェーブレットスケール');
colorbar;
上記に対して、下記のエラーがでました。
wavelet toolboxは購入しております。
どのように修正すればよいかご教授頂けますでしょうか。
[WARN] StAXDialectDetector - Unable to determine dialect of the StAX implementation at jar:file:/C:/Program%20Files/MATLAB/R2023b/java/jarext/axis2/woodstox-core-asl.jar!/
次を使用中のエラー: cwt>parseinputs
スカラー入力引数の数が無効です。
エラー: cwt (297)
fbcell = parseinputs(Norig,TTable,varargin{:});
エラー: untitled (12)
cwt_coeffs = cwt(data, 'wavelet', wavelet_name, 'scales', scales, 'samplingfrequency', fs);

采纳的回答

Atsushi Ueno
Atsushi Ueno 2023-9-23
現在の cwt 関数では引数の形式が多数変更されましたが、新旧どちらの構文にも対応します。
> 上記に対して、下記のエラーがでました
  • 「名前と値の引数」は旧構文では用いられず、文法が合致しません
  • samplingfrequency は旧構文に存在せず、代わりに samplingperiod があります
> どのように修正すればよいかご教授頂けますでしょうか
load vonkoch % これはサンプルデータの読み込み
data = vonkoch(1:510); % data = readmatrix('ファイル名.csv'); %
fs = 1000; %
wavelet_name = 'morl'; %
scales = 1:64; %
%cwt_coeffs = cwt(data, 'wavelet', wavelet_name, 'scales', scales, 'samplingfrequency', fs);
cwt_coeffs = cwt(data, scales, wavelet_name, 1/fs);
figure;
imagesc((1:size(data, 2))/fs, scales, abs(cwt_coeffs));
set(gca, 'YDir', 'normal');
title('ウェーブレット解析');
xlabel('時間(秒)');
ylabel('ウェーブレットスケール');
colorbar;
  1 个评论
cho hunseki
cho hunseki 2023-9-24
Atsushi Uenoさん
いつもご回答頂きありがとうございます。
ご教授頂いたcodeで再現できました。
ありがとうございました。
引き続き宜しくお願い致します。

请先登录,再进行评论。

更多回答(0 个)

类别

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

Community Treasure Hunt

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

Start Hunting!