kendi sorumun cevabını kendim veriyim;
ana fonksiyon
sigma_x=7;
sigma_y=7;
theta=0;
lambda=9;
psi=0;
gamma=1;
gabb=gabor_deneme(sigma_x,sigma_y,theta,lambda,psi,gamma);
figure, mesh(gabb);
figure, imagesc(gabb),colormap('gray')
çağırılan fonksiyon
function gb=gabor_deneme(sigma_xx,sigma_yy,thetaa,lambdaa,psii,gammaa )
sz_x=fix(6*sigma_xx);
if mod(sz_x,2)==0,
sz_x=sz_x+1;
end
sz_y=fix(6*sigma_yy);
if mod(sz_y,2)==0, sz_y=sz_y+1;end
[x y]=meshgrid(-fix(sz_x/2):fix(sz_x/2),fix(-sz_y/2):fix(sz_y/2));
% Rotation
x_theta=x*cos(thetaa)+y*sin(thetaa);
y_theta=-x*sin(thetaa)+y*cos(thetaa);
gb=exp(-.5*(x_theta.^2/sigma_xx^2+gammaa^2*y_theta.^2/sigma_yy^2)).*cos(2*pi/lambdaa*x_theta+psii);