WVD时频代码matlab

clear;
clc
N=4000;
T=2.5;
f0=100;
f1=900;
f3=960;
%三个初始频率
F=800/T/T;
t=linspace（0TN）;
t1=zeros（1N）;
t2=zeros（1N）;
t1（1100:1200）=t（1100:1200）;
t2（1250:1350）=t（1250:1350）;
x=sin（2*pi*（f0+F.*t）.*t）+cos（2*pi*f3.*t1）+cos（2*pi*f3.*t2）;
fs=N/T;