Time-frequency spectrum for a low frequency signal

Question

Océane ESPINASSE 2022-5-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1712060-time-frequency-spectrum-for-a-low-frequency-signal

回答： Vinay 2024-9-25，6:02

在 MATLAB Online 中打开

Hi,

I'm a begginer with matlab and I search to do the time-frequency spectrum for a low frequency signal.

I tried the cwt function on a "normal" frequency signal :

t = 0:0.001:20-0.001;
x = cos(2*pi*50*t).*(t > 0 & t<=5)+ ...
    cos(2*pi*50*t).*(t >= 10 & t < 15)+ ...
    cos(2*pi*100*t).*(t >= 10 & t < 15)+ ...
    cos(2*pi*100*t).*(t >= 15 & t < 20)+ ...
    0.5*randn(size(t));
[y, f] = cwt(x,30);
figure;
subplot(2,1,1)
plot(x)
subplot(2,1,2)
imagesc(abs(y).^2)

And I get this correct result :

But when I try on a low frequency signal, it doesn't work at all...

t = 0:.001:20-0.001;
x = cos(2*pi*0.005*t).*(t > 0 & t<=5)+ ...
    cos(2*pi*0.005*t).*(t >= 10 & t < 15)+ ...
    cos(2*pi*0.01*t).*(t >= 10 & t < 15)+ ...
    cos(2*pi*0.01*t).*(t >= 15 & t < 20)+ ...
    0.5*randn(size(t));
[y, f] = cwt(x,30);
figure;
subplot(2,1,1)
plot(x)
subplot(2,1,2)
imagesc(abs(y).^2)

How could I make it work? Do I need to change functions or settings?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Vinay 2024-9-25，6:02

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1712060-time-frequency-spectrum-for-a-low-frequency-signal#answer_1522180

Hi @Océane ESPINASSE,

The continuous wavelet transform requires few cycles of a signal to accurately determine its frequency content within a given window. When a signal with a frequency of 0.01 Hz is used, its time period is 100 seconds. However, the signal duration is only 20 seconds, therefore cycles are not fully captured.

The issue can be resolved by extending the signal duration to exceed the signal's time period.

Additionally, for low-frequency signals, noise can introduce frequency components ranging from 0 to (fs/2 ), resulting in a wide range of low-amplitude frequency detections due to noise.

Kindly refer to the below documentation of “cwt” for more details:

https://www.mathworks.com/help/wavelet/ref/cwt.html

I hope this helps!