Gaussian filter - with Matlab example

July 11, 2021

After fourier transformation:

fv = fftshift(fft(x));
semilogy(((1:length(fv))-length(fv)/2)/length(fv).*60,real(fv).^2)
xlabel('freqency(herz)');

Gassian filter

sigma = 5;
k = -30:.1:30;
plot(k,exp(-0.5*abs(k).^2/sigma^2));xlabel('freqency(herz)');

clf
gaussianF(y,1e2); xlabel('time(s)');ylabel('position(cm)')

Invere gaussian filter

c = 5;
sigma = 1;
plot(k,(1-exp(-0.5*abs(k+c).^2/sigma^2)).*(1-exp(-0.5*abs(k-c).^2/sigma^2)))
ylim([0 1.2]);xlabel('freqency(herz)');

clf
invergaussF(x,[0.357 5.5],1,0.1e3);

spectrogram

open("signalProcessing.mldatx")