基于frft估计线性调频信号参数

包含线性调频信号的参数估计，两个：初始频率和调制斜率。

clc;
close all;
clear;
% T=10e-6;
% B=30e6;
% K=B/T;
% fs=2*B;
% ts=1/fs;
% N=T/ts;
% t=linspace（-T/2T/2N）;
% st=exp（j*pi*K*t.^2）;

% t=-1*ts/2:1/fs:ts/2;
% S=sqrt（ts/fs）;
% dx=sqrt（ts*fs）;
% tt=-1*dx/2:1/dx:dx/2;%归一化之后的数据

fs=256;T=2;
f0=10;k=10;
S=sqrt（T/fs）;
dx=sqrt（T*fs）;

ts=1/fs;
point=ts*fs+1;

t=-1:0.01:1;
tt=-1*dx/2:1/dx:dx/2;
f1=20;k1=20;

x=exp（1i*2*pi*f0*tt*S+1i*pi*k*tt.*tt*S*S）+exp（1i*2*pi*f1*tt*S+1i*pi*k1*tt.*tt*S*S）;%exp（j*2*pi*（（f0/fs）*t+0.5*k*t.*t））;


p=0.9:0.01:1.8;
u=-dx/2:1/dx:dx/2;
Pp3=zeros（length（u）length（p））;
for k=1:1:length（p）
    p3=frft（xp（k））;
    %vx3=frft（vxa）;
    %vy3=frft（vya）;
    Pp3（:k）=p3.*conj（p3）/point;
end
plot（pPp3）;
figure;
surf（puPp3）;
[CXCI]=max（max（Pp3））;
p（CI）;

Pp31=max（max（Pp3））;
[x1 y1] = find（Pp3==Pp31）;
u = u（x1）;
p = p（y1）;
[ab]=jisuan（pufsT）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         888  2019-05-30 11:41  工程\ce.m
     文件        1638  2019-05-29 09:24  工程\frft.m
     文件         268  2019-05-30 11:04  工程\jisuan.m
     目录           0  2019-10-11 08:57  工程\