Transformada wavelet
Publicado por Maria (14 intervenciones) el 07/09/2017 23:34:16
Buen día
Necesito aplicar a una señal la transformada Wavelet con el comando 'cwt', pero debo crear el tipo de wavelet que voy a usar dado que ninguna de las ya definidas en MATLAB cumple mi propósito. Mi pregunta es: ¿cómo puedo crear una wavelet nueva (wavelet quadratic spline - WQS) si tengo los filtros de descomposición y reconstrucción?
He probado con el comando 'wavemngr' pero cuando coloco la wavelet creada en los argumentos del comando cwt me regresa un error.
Anexo el código.
Necesito aplicar a una señal la transformada Wavelet con el comando 'cwt', pero debo crear el tipo de wavelet que voy a usar dado que ninguna de las ya definidas en MATLAB cumple mi propósito. Mi pregunta es: ¿cómo puedo crear una wavelet nueva (wavelet quadratic spline - WQS) si tengo los filtros de descomposición y reconstrucción?
He probado con el comando 'wavemngr' pero cuando coloco la wavelet creada en los argumentos del comando cwt me regresa un error.
Anexo el código.
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
%% Esta es la wavelet quadratic spline
function [psi, = quaspliwavf_tipoIV()
Dlow=[0.2500 0.7500 0.7500 0.2500];
Rlow =Dlow;Dhigh = [-0.5 0.5];
Rhigh =[-0.5 0.5];
[PHIS,PSIS,PHIA,PSIA,XVAL] = bswfun(Dlow,Dhigh,Rlow,Rhigh,'plot');
psi=PHIS;close all;
%% Aqui creo la wavelet
function crear_wavelet_nueva
wavemngr('del','quaspli');
FN='quadratic spline';
FSN='quaspli';
WT=4;NUMS='';
FILE='quaspliwavf_tipoIV';
wavemngr('add',FN,FSN,WT,NUMS,FILE,[-1,1]);
wavemngr('read')
%% Uso de mi wavelet
function [pxR, pyR]=detectar_R(senal_filtrada)
[S_ECG]=Prueba1_v3();
% senal_filtrada=filtrar_senal()
close all
senal_filtrada = S_ECG;cwsib = cwt(senal_filtrada,1:1:9,'quaspli');
Valora esta pregunta
0
