frft_square2.m

分数阶傅里叶变换

clear;
clc;
close all;
f=[zeros(1,23) ones(1,15) zeros(1,23)];
f_2=zeros(61,61);
for i=1:23
    f_2(i,:)=zeros(1,61);
end
for i=24:38
    f_2(i,:)=f;
end
for i=39:61
    f_2(i,:)=zeros(1,61);
end
figure;
imshow(f_2)
p=0.4;
frft_square2=(frft(f_2,p))';
figure(2);
frft_square2_real=real(frft_square2);
frft_square2_imag=imag(frft_square2);
plot(frft_square2_real);
hold on;
plot(frft_square2_imag,'r--');

frft_mat=zeros(61,61);
for k=1:61
   frft_mat(k,:)=frft_square2(1+61*(k-1):61*k);
end
figure(3);
frft_mat_real=real(frft_mat);
frft_mat_imag=imag(frft_mat);
imshow(frft_mat);

p_inv=-p;
frft_inv=(frft(frft_square2,p_inv))';
frft_mat_inv=zeros(61,61);
for k=1:61
   frft_mat_inv(k,:)=frft_inv(1+61*(k-1):61*k);
end
figure(4);
frft_mat_inv_real=real(frft_mat_inv);
frft_mat_inv_imag=imag(frft_mat_inv);
imshow(frft_mat_inv);