1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
| function [CWTcoeffs,scales] = analyse_adapted_wavelet(S,wavelet_sig,scalogramme)
%ANALYSE_ADAPTED_WAVELET_FUNC Summary of this function goes here
% Detailed explanation goes here
delete('adp_FRM1.mat')
X = wavelet_sig(1,:);
Y = wavelet_sig(2,:);
save adp_FRM1 X Y
wavemngr('add','AdapF1','adpf1',4,'','adp_FRM1.mat',[0 1]);
% Ondelette = load('adp_FRM1.mat');
long = length(wavelet_sig);
%Trace coefficients
stepSIG = 1/1000; %meilleure solution
%stepSIG = 1/10;
stepWAV = 1/long;
% wname = 'adpf1';
%scales = (1:1.5*long)*stepSIG; %meilleure solution
scales = (1:1.5*long)*stepSIG;
WAV = {'adpf1',stepWAV};
SIG = {S,stepSIG};
if scalogramme ==1
figure;
disp(['//adapted wavelet analysis //start cwt' datestr(now,'yyyy-mm-dd/HH:MM:SS')])
CWTcoeffs = cwt(SIG,scales,WAV,'scalCNT'); colormap jet
else
disp(['//adapted wavelet analysis //start cwt' datestr(now,'yyyy-mm-dd/HH:MM:SS')])
CWTcoeffs = cwt(SIG,scales,WAV);
end |
Partager