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