fft适应goertzel写到算法
2视图(30天)
显示旧的评论
我想做的就是利用FFT和Goertzel写到算法寻找罪两个信号之间的相位差。我相信我与FFT代码运行得很好。它可以输出正确的相位差值,但Goertzel写到算法不能输出预期的结果。据我所知,Goertzel写到算法之间的差异和FFTis Goertzel写到算法有一个感兴趣的区域。我不知道为什么当我做Goertzel写到算法,它只有一个值为每个犯罪的功能。结果截图和代码如下所示。
清晰;
clc;
f = 10000;
p = 1 *π;
t = linspace(0, 10 *π,100000);
s =[罪(2 *π* f * t);罪(2 *π* f * t + p)];%设置信号
Ts =意味着(diff (t));
Fs = 1 / Ts;
Fn = f / 2;
% fts = fft (s) /长度(t);%做转换
indxs =找到(> = 9500 & f < = 10500);
fts = goertzel写到(年代,indxs);
阵线= linspace(0, 1,修复(长度(t) / 2) + 1) * Fn;
4 = 1:长度(艘);
(pk, ix) = max (abs (fts (Iv,:)));
phsmin = [min(角(fts(九(1):)));马克斯(角(fts(九(2):))));
phsdif = diff (phsmin);%计算相位差
PhaseDiff = sprintf('相位差= %。3 f rad, phsdif)
