w1=0.23;
w2=1;
f0=3.4*10^11;
fs=100;
c=3*10^8;
r=0.00025;
r1=0.005;
r2=0.0008;
t=0:0.01:30;
a=0.25-r*w2;
t3=w2.*t;
t4=floor(t3);
y=abs(t3-0.5-t4);
aaa=1/(1-2*a).*(y-a+abs(y-a));
s=r1.*sin(2.*pi.*w1.*t)+r2.*aaa;
sr=exp(1j*2*pi*f0.*(t-2*s/c));
L=length(sr);
%[tfr,t,f] = tfrstft(sr.',1:L,L,hamming(63),0);
[tfr,t,f] = tfrstft(sr.',1:L,L-2000,hamming(63),0); % 此处修改了频率的数量,通过tfr的矩阵大小可以分清楚哪个是时间,哪个是频率
figure();
%contour(t/fs,f(1:length(f))*fs,abs(tfr(1:length(f),:)));
mesh(t/fs,f(1:length(f))*fs,abs(tfr(1:length(f),:)));
xlabel('time/s')
ylabel('frequency/Hz')
view([0,90]); % 时间和瞬时频率
xlim([0,30])
tfr=abs(tf