gabor wavelet 初步代码 已修改

function [gb_R,gb_I]=gabor_wt(lambda,theta,m,gamma,a,psi,b,delta)
%[gb_R,gb_I]==gabor_wt(lambda,theta,m,gamma,a,psi,b,delta)
%lambda 为函数波长
%theta  为gabor 核函数的方向
%psi    为相位偏移
%sigma  为高斯函数的标准差
%gamma  This parameter, called more precisely the spatial aspect ratio, specifies the ellipticity of the support of the Gabor function. 
% For γ = 1, the support is circular. For γ < 1 the support is elongated in orientation of the parallel stripes of the function. 
%Default value is γ = 0.5
%m      尺度参数
%a      尺度参数 默认为2
%以h(x,y)=exp(-.5*(x_theta.^2/sigma_x^2+y_theta.^2/sigma_y^2)).*exp(1i*(2*pi/lambda*x_theta+psi));
% 为母小波，对它进行尺度收缩和旋转 hmn(x,y)=a^(-m).*h(x_theta,y_theta);a >1  
%a^(-m) 为尺度因子 x_theta=a^(-m)*(x*cos(theta)+y*sin(theta));
% y_theta=-a^(-m)*(x*sin(theta)+y*cos(theta)); theta=n*pi/k; S,K 为尺度和方向的数目，m=0,1,…,S?1,n=0,1,…,K