%生成信號
N=128;
sig=gdpower(N,1/2);
figure(1);
plot(1:N,real(sig));
xlabel('時間 t');
ylabel('幅值 A');
%計算D-Flandrin 分布
[tfr,t,f]=tfrdfla(sig,1:N,0.01,0.22,N,1);

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