%產生暫態信號
transig=amexpo1s(64).*fmconst(64);
sig=[zeros(100,1);transig;zeros(92,1)];
%加噪聲
sign=sigmerge(sig,noisecg(256),-5);

%繪制時域波形
figure(1);
plot(real(sign),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');

%計算功率譜
nfft=1024;
dsp=abs(fft(sign,nfft)).^2;
%繪制頻譜圖
figure(2);
plot((0:nfft/2-1)/nfft,dsp(1:nfft/2));
xlabel('頻率 f');
ylabel('功率譜 PSD');

%計算信號的譜圖
[tfr,t,f]=tfrsp(sign);
%繪圖
%%繪制等高線圖
figure(3);
contour(t,f(1:128),abs(tfr(1:128,:)).^2);
xlabel('時間 t');
ylabel('頻率 f');
%%三維圖

figure(4);
mesh(t,f(1:128),abs(tfr(1:128,:)).^2);
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');