chirp信号的frft的峰值搜索仿真

chirp信号的matlab仿真，显示了chirp信号的frft的峰值搜索，参数估计

clear all
clc;
N=511;
fi=0;
f0=-28.9;  %初始频率
fs=97.6;   % 抽样频率
k=-8;      %调频率
f=zeros（1N）;
t=zeros（1N）;
for n=-255:255
    x（n+256）=exp（j*fi+j*2*pi*f0*n/fs+j*pi*k*（（n/fs）*（n/fs）））;%generate signal t=n/fs
    t（n+256）=n/fs;%time                  
    f（n+256）=n*fs/N;%freq
end
figure（1）;
plot（real（x））;
title（‘LFM信号‘）;
X=awgn（x20）;
r=0.01;
a=[0:r:2];%fractional power
G=zeros（length（a）length（X））; %for different power do frft
h=0;
for l=1:length（a）
    T=frft（Xa（l））;
    G（l:）=abs（T（:））;
  if（h<=max（abs（T（:））））
    h=max（abs（T（:）））;
    p1=a（l）                %阶数
  end
end
for m=1:length（a）;
    y（m）=max（abs（G（m:）））;
end
figure（2）;
stem（y） %可以看看投影在阶数轴的曲线
figure（3）;
mesh（G）;

tt=t（N）-t（1）;
kk=tt/（fs）;
kr=-cot（p1*pi/2）/kk % k参数的估计值

FF=frft（xp1）;
[aabb]=max（abs（FF））;
u0=f（bb）;
ff0=u0*csc（p1*pi/2）  % 中心频率f0的估计值

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        911  2014-03-05 18:45  frft峰值搜索\fengzhisousuo.m

     文件       1538  2010-03-11 15:45  frft峰值搜索\frft.m

     文件        978  2014-03-05 19:51  frft峰值搜索\Untitled.m

     文件        250  2014-03-04 15:47  frft峰值搜索\说明.txt

     目录          0  2014-04-05 16:23  frft峰值搜索

----------- ---------  ---------- -----  ----

                 3677                    5