%生成信號(hào)
N=128;
sig=gdpower(N,-1);
figure(1);
plot(1:N,real(sig));
xlabel('時(shí)間 t');
ylabel('幅值 A');
%計(jì)算D-Flandrin 分布
[tfr,t,f]=tfrunter(sig,1:N,'A',0.01,0.22,172,1);

%顯示
figure(2);
contour(t,f,tfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');
figure(3)
mesh(t,f,tfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');
zlabel('幅值 A');