Code covered by the BSD License

# Continous Fractional Fourier Transform FrFT

### Ashish Soni (view profile)

16 Apr 2013 (Updated )

The FrFT, which is a generalization of the ordinary Fourier transform.

Untitled4.m
```%Program to Calculate CONTIONOUS FRACTIONAL FOURIER TRANSFORM
% By Ashish Soni,ATC Indore
tic
clc;
clear all;
close all;

T=input('enter the value of T - ');
t=-T:0.01:T;
X=input('enter input function - ');
disp('Program to Calculate Continous Fractional Fourier Transform')
U=input('enter the range of frequency - ');
a=input('Enter the rotation angle a (except 0) - ');

p=0;q=0;
for t=-T:0.01:T
p=p+1;
for u=-U:0.1:U;
q=q+1;
%         x1(p,q)=(sqrt((1-j*cot(a))/(2*pi)))*exp(j*(((0.5*(u^2)+0.5*(t^2))*cot(a))-u*t*(1/sin(a))));
x1(p,q)=(sqrt(1-j*cot(a*pi/2)))*exp(j*pi*((u^2)*cot(a*pi/2)+(t^2)*cot(a*pi/2)-2*u*t*(1/sin(a*pi/2))));
end
q=0;
end

u=-U:0.1:U;
x2=0.01*(X*x1);

% x3=20*log10(abs(x2)./max(x2));
t=-T:0.01:T;

figure,
plot(t,X);xlabel('Time');ylabel('Amplitude');title('Input Function');

figure,
plot(abs(x2));
xlabel('frequency');ylabel('Amplitude');title('Magnitude Response of FRFT');

figure,
plot(imag(x2));
xlabel('frequency');ylabel('Amplitude');title('Imaginary part of FRFT');

toc```