分数阶傅里叶matlab,怎么用分数阶傅里叶变化(FrFT)求线性调频信号(LFM)的调频率...

各位大神帮帮忙啊，急用

怎么利用分数阶傅里叶变换求线性调频信号的调频率

FrFT程序我已经有了，还需要怎么做才能求LFM信号的调频率

或者看看我的程序有没有什么问题，谢谢

下面是FrFT程序

function Faf = frft(f, a)

% The fast Fractional Fourier Transform

% input: f = samples of the signal

%        a = fractional power

% output: Faf = fast Fractional Fourier transform

error(nargchk(2, 2, nargin));

f = f(:);

N = length(f);

shft = rem((0:N-1)+fix(N/2),N)+1;

sN = sqrt(N);

a = mod(a,4);

% do special cases

if (a==0), Faf = f; return; end;

if (a==2), Faf = flipud(f); return; end;

if (a==1), Faf(shft,1) = fft(f(shft))/sN; return; end

if (a==3), Faf(shft,1) = ifft(f(shft))*sN; return; end

% reduce to interval 0.5 < a < 1.5

if (a>2.0), a