function ff=frftc2d(f,ax,ay);
ax=ax*pi/2;
ay=ay*pi/2;
N=length(f);
for ax=0.1:-1:1;
    for ay=0.1:-1:1;
        f=input(f,ax,ay);      
x=([0:N-1]-[N-1]/2)/sqrt(N);
Ba=exp(i*pi*x.^2*csc(ax));
Bpa=exp(-i*pi*x.^2*tan(ax/2));
fa=f.*Bpa;
C1=fft2shift(fft2(fft2shift(fa)));
C2=exp(-i*pi*x.^2*sin(ax));
CC=fft2shift(ifft2(fft2shift(C1.*C2)));
ff=Bpa.*CC;
ff=ff*sqrt(1-i*cot(ax))*sqrt(1-i*cot(ay));
这是我编的程序  其中for ax=0.1:-1:1;
    for ay=0.1:-1:1;
        f=input(f,ax,ay);  这里需要是矩阵输入  不清楚这样对不对
ff=ff*sqrt(1-i*cot(ax))*sqrt(1-i*cot(ay));
显示不完整  求高人帮忙
[x,y]=meshgrid(linspace(-3,3,256));
X=zeros(size(x));
X(abs(x)<=1&abs(y)<=1)=1;
f1=frftc2d(X,0.6,0.4);
figure;subplot(121);imshow(abs(X),[]);
xlabel('(a)','Fontsize',14,'Fontname','Times New Roman');
subplot(122);imshow(abs(f1),[]);
xlabel('(b)','Fontsize',14,'Fontname','Times New Roman');
这个是需要仿真的信号的程序