如何解决如何从 Lanczos 滤波器获取频率窗口
dT = 1 % sampling interval
Cf = 1/40 % cutoff frequency
fl = 100 % ?
M = 100 % number of coefficients ? not sure about number
LoH = 1 % low pass
Nf=1/(2*dT); %Nyquist frequency
% normalize the cut off frequency with the Nyquist frequency:
Cf = Cf/Nf;
% lanczos cosine coeficients:
coef = lanczos_filter_coef(Cf,M); coef = coef(:,LoH);
% Filter in frequency space:
[window,Ff] = spectral_window(coef,length(vel)); Ff = Ff*Nf;
% Filtering:
[y,Cx] = spectral_filtering(vel,window);
function coef = lanczos_filter_coef(Cf,M)
% Positive coeficients of lanczos [low high]-pass.
hkcs = lowpass_cosine_filter_coef(Cf,M);
sigma = [1 sin(pi*(1:M)/M)./(pi*(1:M)/M)];
hkB = hkcs.*sigma;
hkA = -hkB; hkA(1) = hkA(1)+1;
coef = [hkB(:) hkA(:)];
end
function coef = lowpass_cosine_filter_coef(Cf,M)
% Positive coeficients of cosine filter low-pass.
coef = Cf*[1 sin(pi*(1:M)*Cf)./(pi*(1:M)*Cf)];
end
function [window,N)
% Window of cosine filter in frequency space.
Ff = 0:2/N:1; window = zeros(length(Ff),1);
for i = 1:length(Ff)
window(i) = coef(1) + 2.*sum(coef(2:end).*cos((1:length(coef)-1)'*pi*Ff(i)));
end
end
function [y,window)
% Filtering in frequency space is multiplication,(convolution in time
% space).
Nx = length(vel);
Cx = fft(vel(:)); Cx = Cx(1:floor(Nx/2)+1);
CxH = Cx.*window(:);
CxH(length(CxH)+1:Nx) = conj(CxH(Nx-length(CxH)+1:-1:2));
y = real(ifft(CxH));
end
我需要从中绘制时间窗口和频率窗口。我已经得到了来自 coef 的时间窗口,但我无法弄清楚哪个输出会给我频率窗口。我已经绘制了输出变量的所有可能组合,并尝试对其中一些变量进行傅立叶变换,但我尝试的任何方法都没有给出预期的数字。
解决方法
频率窗口在输出“窗口”中,因此需要绘制(Ff,window)。您会得到一个在 Cf 附近大幅下降的图表(因此为 1/50),这是您选择的截止频率,它将用于低通滤波器的频率与用于高通滤波器的频率分开。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。