Converting code from old CWT to new CWT

6 次查看（过去 30 天）
Michael
Michael 2020-3-20
评论： Michael 2020-5-1
I am trying to understand and re-implement an algorithm using someone else's code, which uses the old cwt function. Here is the portion of code in question:
fb = 13.5;
fc = 0.5;
scale = 3.7;
output = cwt(signal,scale,['cmor' num2str(fb) '-' num2str(fc)]);
The variable signal is time series data at 100Hz. The goal of this is to return the scale parameters of a bandpassed version of the signal at approximately 10-16Hz.
Here are my questions:
  1. Is this actually getting the proper bandpass range? Presumably this is using scale parameter. However, Fb and Fc seem mislabeled in this code and seem redudant to the scale parameter.
  2. How can I properly implement this using the new CWT functions? That is, given a fixed center frequency, bandwidth, and sampling frequency what is the best way to implement the appropriate single-wavelet transformation using a complex morlet wavelet using the new framework?
Thanks!

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

回答（1 个）

Prabhan Purwar
Prabhan Purwar 2020-3-28
编辑：Prabhan Purwar 2020-3-28
Hi,
Variables fb and fc represents properties of the Morse Wavelet used to determine the Continuous Wavelet Transform (CWT) of the signal.
Following link illustrates the use of filterbank in cwt
To verify your results take help of Magnitude Scalogram (as illustrated in the above example)
For further insights refer to the following links:
  1 个评论
Michael
Michael 2020-5-1
Thank you.
If I just wanted to have the wavelet coefficient values for a single wavelet with a center frequency of Fc and a bandwidth of Fb at sampling frequency Fs, what would be the most direct way to go about it using the updated functions?
I cannot see how to directly do this with these functions, which seem more designed ot give you the scaleogram or filtered estimate from multiple wavelets, rather than the coefficients for a single wavelet of specific properties.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Continuous Wavelet Transforms 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by