Gaborfilter=zeros(num,num); x=[-(num-1)/2:1:(num-1)/2]; y=[-(num-1)/2:1:(num-1)/2]; xx=zeros(num,num); yy=zeros(num,num); Gaborfilter=zeros(num,num); for i=1:num for j=1:num xx(i)=x(i).*cos(theta)+y(j).*sin(theta); yy(j)=-x(i).*sin(theta)+y(j).*cos(theta); Gaborfilter(i,j)=exp((xx(i)^2+0.25*yy(j)^2)/(-50))*exp(1i*((2*pi*xx(i)/10))); end end end

这段代码是在生成 Gabor 滤波器，其中 num 表示滤波器的大小，theta 表示滤波器的方向。具体来说，代码首先初始化了一个 num × num 的矩阵 Gaborfilter，并定义了 x 和 y 分别表示从 -(num-1)/2 到 (num-1)/2 的序列。然后，循环遍历 Gaborfilter 矩阵中的每个位置，计算 xx 和 yy 分别表示该位置在经过旋转变换后的坐标，最后根据 Gabor 滤波器的公式计算该位置的值。其中，exp 表示自然指数函数，1i 表示虚数单位。