gabor_r = function(x){
exp(-alpha^2*(x-x0)^2)*(cos(2*Pi*beta*(x-x0)+gamma))
}
gabor_i = function(x){
exp(-alpha^2*(x-x0)^2)*(sin(2*Pi*beta*(x-x0)+gamma))
}
alpha=1;
beta=1;
gamma=0.5;
x0=0
x <- seq(-3,3,length=500)
ylim_un=-1
ylim_up=1
plot(x,gabor_r(x),type="l",ylim=c(ylim_un,ylim_up))
par(new=T)
plot(x,gabor_i(x),type="l",lty = 2,ylim=c(ylim_un,ylim_up))
axis(2, pos = 0, at = -3:3, adj = 1, col = 2)
alpha=1;
beta=2;
gamma=2.5;
x0=0
x <- seq(-3,3,length=500)
ylim_un=-1
ylim_up=1
plot(x,gabor_r(x),type="l",ylim=c(ylim_un,ylim_up))
par(new=T)
plot(x,gabor_i(x),type="l",lty = 2,ylim=c(ylim_un,ylim_up))
axis(2, pos = 0, at = -3:3, adj = 1, col = 2)