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

Roof0102

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

Blah

Code：
X=zeros(30,30);%本来这里用图片的，图片都太大了，所以创建一个小的。
+ j$ E* t! T4 A2 T* @7 g+ rfor I=14:16
. x0 Q7 f0 x. M8 x( i- D' E9 Z/ ^- Vfor J=13:17
4 f" k3 G4 Y4 LX(I,J)=255;
% B+ Y. f8 B) h* }' H& k9 Yend
9 c8 j$ N( m( o2 [  Kend
- ^" n" K3 c+ T( F4 V7 `X=mat2gray(X);
$ K( K: R4 P6 `! R* Psubplot(1,2,1); imshow(X); title 原图;
b=size(X);
M=1/b(1)*b(2);
Y=zeros(b);
for(u=1:b(1))
9 j0 w! V$ L% `2 R% hfor(v=1:b(2))
for(x=1:b(1))
8 C( j' H1 D# d; _for(y=1:b(2))
9 l) o1 u" p/ D7 w" }Y(u,v)=(X(x,y)*exp(-j*2*pi*(u*x/b(1)+v*y/b(2))))+Y(u,v);
end;
, W+ Y" d1 x" Z1 e) |" {1 H2 }end;
" `' }1 E# e) ]. j* _, l2 pY(u,v)=M*Y(u,v);
end;
end;
for(u=1:b(1))
for(v=1:b(2))
Y(u,v)=sqrt(real(Y(u,v))*real(Y(u,v))+imag(Y(u,v))*imag(Y(u,v)));
# }& _! a8 z& ^& x  ?end
! _  g: W- E8 J7 t% ^; w0 n+ dend
1 ]  ^/ Q9 ]' W0 p6 CY=mat2gray(Y);
3 z* G/ b2 x' o$ |6 R# lsubplot(1,2,2); imshow(Y); title 经FDT变换后

Jame33

再看看别人是怎么说的

land

试试一楼的方法

理论的

同问，求结果