%生成信號
N=128;
sig=gdpower(N);
figure(1);
plot(1:N,real(sig));
xlabel('時間 t');
ylabel('幅值 A');
%計算bertrand分布
[tfr,t,f]=tfrbert(sig,1:128,0.01,0.22);
%顯示
figure(2);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
figure(3)
mesh(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');