Optimization of Wavelet Analysis - Code Execution Time

Question

Velu Prabhakar Kumaravel 2020-7-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/557686-optimization-of-wavelet-analysis-code-execution-time

回答： Chaitanya Mallela 2020-7-13

Hello,

Could someone help me to understand why does it take so much time to execute the following piece of code:

tic
frange = [2 20]; % Frequency range I am interested
wname = 'mexh'; % Wavelet Name
C = 124; % number of sensors 
S = 67534; % number of samples
data = randn(C,S); % Input Data 
energy_mean = zeros(1,C);
f = scal2frq(1:500,wname,1/250); % Computing the Pseudo Frequency
indices = f >= frange(1) & f <=frange(2); % Finding the indices of frange
    
    for c = 1:C % For each sensor data
        
        x = data(c,:);
        CWTcoeffs = cwt(x,1:500,wname); % compute coefficients
        cfs = CWTcoeffs(indices,:).^2;
        energy_mean(ch) = 10*log10(mean(mean(cfs,2))); % average power in db
        
    end
toc
Elapsed time is 918.298398 seconds.