Selecting level in Wavelet Multiscale Principal Components Analysis

3 ビュー (過去 30 日間)
Aske Ejdrup
Aske Ejdrup 2020 年 5 月 14 日
コメント済み: Megha S Kumar 2020 年 12 月 19 日
I have a time series of neurotransmitter levels in the brain. I'd like to explore the contribution to the signal at different frequencies, and I've stumbles over the wmspca() function in the guide from https://se.mathworks.com/help/wavelet/ug/wavelet-multiscale-principal-components-analysis.html.
I deconstruct the signal using modwt and input the individual levels as the multivariate signal. However, wmspca prompts me for a level, but I'm unsure what the given level represents in the wmspca function and how it impacts the result. I'm familiar with levels in wavelet analysis, and so far I've been extracting frequency ranges manually, but now I'd like to explore the PCAs of the signal.

回答 (1 件)

Rajani Mishra
Rajani Mishra 2020 年 7 月 31 日
LEVEL input argument accounts for number of detail coefficients in decomposition step. If “LEVEL” argument is not passed, then “DEC” is passed. Where “DEC” is output of mdwtdec which contains detail and approximate coefficients. wmspca works level wise on a signal.
Signal has level wise coefficients (suppose level is 5 then it has 5 detail coefficients (matrix) and 1 approximate coefficient (matrix) ). On passing signal and level wmspca calculates coefficients (this process is called as decomposition).
  1 件のコメント
Megha S Kumar
Megha S Kumar 2020 年 12 月 19 日
Hi, i am trying to implement wmspca inbuilt function step by step and compare it with eigen vectors in Wavelet Analyzer app also. when i input a 1024 x 4 signal matrix and apply 5 level dwt on each column. If i choose 1 principal component how can i reconstruct using inverse dwt as wavelet coefficients length and eigen vector lengths are totally different?.
Step after application of pca to each scale is confusing me. i am working around this for the past few weeks.
Any help is appreciated. Thanks in advance.

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeEEG/MEG/ECoG についてさらに検索

製品


リリース

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by