Supongamos que tengo la siguiente señal:
fs = 20; % Sampling rate [Hz]
Ts = 1/fs; % Sampling period [s]
duration = 25; % Duration [sec]
t = 0: Ts : duration-Ts; % Time vector
ss1 = (160.74*exp(-0.15*t)).*cos((2*pi*0.5*t)); %signal 1%
ss4 = zeros(1,length(t)); %signal 2%
tz = 0: Ts : 12-Ts;
for i = 1:length(tz)
ss4(i) = 30.54*exp(-0.25*tz(i)).*cos((2*pi*4*tz(i)));
end
xp = ss1+ss4;
La longitud de mi ventana es de 10 segundos. En los primeros 10 segundos, tiene 200 muestras y calcula dos frecuencias, es decir, 4 Hz y 2 Hz. En la estimación recursiva, esto continúa hasta 240 muestras y, a partir de las 241 muestras, el cálculo de la frecuencia se desvía debido a las muestras mixtas. Por ejemplo, si uno elige su ventana actual como 41: 241 (200 muestras o 10 segundos), la estimación de frecuencia es de .5001 y 4.001 Hz y el segundo valor sigue aumentando mientras que en realidad debería disminuir. Todo está sucediendo debido a puntos de muestreo mixtos. Algunos de ellos son componentes de doble frecuencia, mientras que otros son únicos a medida que uno avanza. Quiero saber cuál fue el motivo real detrás de Prony recursivo? Si alguien necesita el código de Prony por favor mencione. Espero que estés al tanto de esto.