[求助] 请了解利用匹配傅立叶变换估计LFM信号参数的朋友帮忙看看，谢谢!

clear;
N=1024;   %采样点数
fs=1;   %采样频率（此时经调频傅立叶变换后的成像可观察到峰值，当采样频率扩大时，例如变为100时，最终得到的图像没有明显的峰值）
c=3e8;
h0=10e7;
f0=11.1e6;   %LFM信号中心频率
w0=5.2e5;    %LFM信号调频率

for n=1:N
    s(n)=5*exp(complex(0,4*pi/c*(h0+f0*(n/fs)+w0*(n/fs)^2)));     %LFM信号表达式
end

ff=f0-50*f0/100:f0/100:f0+49*f0/100;
pp=w0-50*w0/100:w0/100:w0+49*w0/100;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %%调频傅立叶变换%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mf_fft=zeros(length(pp),length(ff));
for i=1:length(ff)
    temp=zeros(length(pp),1);
    for k=1:N;
        temp=temp+s(k)*exp(complex(0,-4*pi/c*(ff(i)*(k/fs)+pp'*(k/fs).^2)));
    end
    mf_fft(:,i)=temp;
end

figure,
mesh(abs(mf_fft));