高斯函数sigma变量的选择

Question

我是新来的信号处理，并希望应用低通滤波器使用fft。 我发现this post它回答我的问题。在使用它的时候，我有一个问题： 假设我的截止频率是3500Hz，采样率是25600Hz，sigma在产生高斯曲线时将使用什么值，如eigenchris'answer以下代码给出的？

gauss = zeros(size(Y));  
sigma = 8;       % just a guess for a range of 20 

gauss(1:r+1) = exp(-(1:r+1).^ 2/(2 * sigma^2)); % +ve frequencies 
gauss(end-r+1:end) = fliplr(gauss(2:r+1));   % -ve frequencies 
y_gauss = ifft(Y.*gauss,1024); 

下面是我使用功能的代码：

clf; clc; 
Fs = 25600; 
file = '01cKhaitan181015M4_Opp_LeftS1H8_a.dat'; 
signal = dlmread(file); % read file from specified location 
signal = signal - mean(signal); 

N  = size(signal, 1); 
time = 1000*(0 : N-1)/Fs; % in msec 
freq = (-Fs/2 : Fs/N : Fs/2-Fs/N)'; 

Y = fft(signal, 1024); 

r = 141; % range of frequencies we want to preserve 

gauss = zeros(size(Y)); 
sigma = 119.75; 
gauss(1:r+1) = exp(-(1:r+1).^ 2/(2 * sigma^2)); % +ve frequencies 
gauss(end-r+1:end) = fliplr(gauss(2:r+1));   % -ve frequencies 
y_gauss = ifft(Y.*gauss,1024); 

hold on; 
plot(time, signal, 'k'); plot(time, abs(y_gauss), 'c'); 
legend('signal', 'gaussian', 'Location', 'southwest') 
hold off; 

下面是数据文件的链接

https://www.dropbox.com/s/edb8g43j4a54jvq/01cKhaitan181015M4_Opp_LeftS1H8_a.dat?dl=0

Answer 1

截止频率被定义为频率，其中衰减是0.5或〜6dB。您已经知道所需的截止频率为3500Hz。下一步将与获得相应的索引：

cutoff_frequency = 3500; 
sampling_rate = 25600; 
cutoff_index = 1 + cutoff_frequency/sampling_rate * length(Y); 

那个cutoff_index你想收到预期的0.5衰减。求解sigma这将产生在cutoff_index 0.5衰减的高斯公式：

%% Derivation of sigma value 
% 0.5 = exp(-cutoff_index^2/(2 * sigma^2)); 
% log(0.5) = -cutoff_index^2/(2 * sigma^2)); 
% sigma^2 = -cutoff_index^2/(2 * log(0.5)); 
% sigma^2 = cutoff_index^2/(2 * log(2)); 

产生：

sigma = cutoff_index/sqrt(2 * log(2)); 

因此，例如，如果length(Y)是1024，你会得到

sigma = (3500/25600 * 1024)/sqrt(2 * log(2)); % approx. 119.75