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;%歸一化之后的數(shù)據(jù)

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??工程\