mat_fast.txt

来自「matlab的代码优化」· 文本 代码 · 共 45 行

close all;
clear all;
clc;
load gindb   %取点的数据库（每图40点）
load imgdb  %图片数据库（20张）

% Creating 40 Gabor Filters
%G = cell(5,8);
Kmax = pi/2   
f = sqrt(2)
sig = 2*pi
m=1:20  %循环运算20张图片
n=1:40  %循环运算每张图片上的40个点
x=B(n,1,m)
y=B(n,2,m)  %x与y为取的点，也用循环实现全部运算
%Z(:,:,n,m)=[x;y]
for s = 0:4
for t = 0:7
K = Kmax/f^s * exp(i*t*pi/8);
Kreal = real(K);
Kimag = imag(K);
NK = Kreal^2+Kimag^2;
G=(sig^(-2)*K^2)*exp((-1)*NK*(x^2+y^2)/(2*sig^2))*(exp(i*[Kreal Kimag]*[x;y])-exp(-(sig^2)/2));
W=zeros(1,1)
W(1)=G
%W(2)=imag(G{t,s,n,m})
C=zeros(1,1)
C(1)=A(round(x),round(y)) 
R=conv(C,W)  
%用J存储幅值,以备计算
%J{:,:,n,m}=sqrt(real(R)^2+imag(R)^2)  %用一个四维细胞数组存储计算结果
end
end
save gabor J




Kmax = pi/2   
f = sqrt(2)
sig = 2*pi
%Z(:,:,n,m)=[x;y]
s = 0:4
t = 0:7
K = Kmax/f^s * exp(i*t*pi/8);