求二维离散分数傅里叶变换Matlab仿真程序

Roof0102

求二维离散分数傅里叶变换Matlab仿真程序

Blah

Code：
* ?6 v- y! [( P  cX=zeros(30,30);%本来这里用图片的，图片都太大了，所以创建一个小的。
9 T) S2 o8 ~  f' s+ cfor I=14:16
for J=13:17
1 E) ?2 C% D0 {X(I,J)=255;
end
/ T3 E! }2 K: n; B$ ~# Qend
X=mat2gray(X);
subplot(1,2,1); imshow(X); title 原图;
, \% J1 Z" c2 vb=size(X);
* V2 [5 r' ?  c, V1 ?M=1/b(1)*b(2);
Y=zeros(b);
2 l7 z. S- X+ }" l% }) Gfor(u=1:b(1))
/ a/ `2 \/ }6 X' E: r3 Cfor(v=1:b(2))
$ [. d% Q5 }* Y+ pfor(x=1:b(1))
for(y=1:b(2))
- ^9 }! D. c0 Y9 PY(u,v)=(X(x,y)*exp(-j*2*pi*(u*x/b(1)+v*y/b(2))))+Y(u,v);
9 _( }6 a( ~* F4 T6 t: rend;
end;
* [/ {: h- q" G& `! W( D# P; ?2 |Y(u,v)=M*Y(u,v);
9 h* P, l, ?2 Zend;
end;
2 Y3 h; o" g/ g: m! gfor(u=1:b(1))
+ {/ d) [5 d, ^8 h/ cfor(v=1:b(2))
Y(u,v)=sqrt(real(Y(u,v))*real(Y(u,v))+imag(Y(u,v))*imag(Y(u,v)));
/ s' v7 a2 @, ~. Z) Y4 i3 e( s# [end
end
& d% D# l, E* \" c  r! D9 cY=mat2gray(Y);
: `. s7 W" @) f8 I) Ysubplot(1,2,2); imshow(Y); title 经FDT变换后

Jame33

再看看别人是怎么说的
. h5 N/ s. X; o2 g$ ?

land

试试一楼的方法
2 I( u9 \) U( F7 _; \4 u

理论的

同问，求结果