0

我想在信号窗口上找到瞬时频率。我正在提取一小部分信号并试图在该窗口中找到瞬时频率。但是该频率与信号的实际频率不匹配。下面是我的代码。

close all; clear all; clc

fs = 25000;
T = 0.5;
t = 0:1/fs:T;
n = length(t);

m = floor(n/4);
f1 = 100;
z1 = cos(2*pi*f1*t(1:m));

f1 = 200;
z2 = cos(2*pi*f1*t(1:m));

f1 = 300;
z3 = cos(2*pi*f1*t(1:m));

f1 = 400;
z4 = cos(2*pi*f1*t(1:(n-3*m)));

z = [z1 z2 z3 z4];

window = 100;
wStart = 1;
wEnd = wStart + window;
freqs = [];

while wEnd < length(z)
    x = z(wStart:wEnd);
    y = t(wStart:wEnd);

    h=hilbert(x);   
    unrolled_phase = unwrap(angle(h));

    dx = diff(unrolled_phase);
    dy = diff(y);
    p = dx./dy;
    inst_freq = p/(2*pi) + 2*pi;
    freqs = [freqs inst_freq];

    wStart = wEnd;
    wEnd = wStart + window;
end

'freqs' 的情节是这样的： 情节（频率）

我认为我计算频率的方式有问题，但我不确定。有人可以帮忙吗？从代码中可以看出，我只需要将频率设为 100,200,300 和 400。我要做的就是找到每个频率的起点（时间）

4

1 回答 1

0

第一个，理论上的答案是，一个信号实际上是几个频率的组合。此外，如果您将信号分解为 100 Hz 范围内的 100Hz、200Hz、300Hz 等，那么您将无法将信号重新创建为与转换后的信号完全相同的信号。也就是说，逆变换只能使用 4 个频率，并且会以输入信号和恢复信号听起来不一样的方式组合！

于 2014-03-24T17:50:38.493 回答