Gabor function

###########################
#Gabor function real part
gabor_r = function(x){
    exp(-alpha^2*(x-x0)^2)*(cos(2*Pi*beta*(x-x0)+gamma))
}
#Gabor function imaginary part
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) # 赤で X 軸を描く

###########################
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) # 赤で X 軸を描く