fs=10000;
%采樣頻率
t=（DUT20180413215506040（:1））;
s=（DUT20180413215506040（:2））;
%兩不同頻率正弦信號合成的信號
%%%%小波時頻圖繪制%%%%
wavename=‘cmor3-3‘;
totalscal=256;
%尺度序列的長度，即scal的長度
wcf=centfrq（wavename）;
%小波的中心頻率
cparam=2*wcf*totalscal;
%為得到合適的尺度所求出的參數
a=totalscal:-1:1;
scal=cparam./a;
%得到各個尺度，以使轉換得到頻率序列為等差序列
coefs=cwt（sscalwavename）;
%得到小波系數
f=scal2frq（scalwavename1/fs）;
%將尺度轉換為頻率
figure（1）plot（ts）;
xlabel（‘時間/ms‘）;
ylabel（‘幅值‘）;
title（‘時域波形圖‘）;
%繪制原始信號圖
figure（2）;
%繪制信號的傅里葉變換的頻譜圖
N=length（t）;
y=fft（sN）;
amg=abs（y）;
plot（（1:N/2）*fs/Namg（1:N/2）*2/N）;
xlabel（‘頻率?f/Hz‘）;
ylabel（‘幅值‘）;
title（‘fft頻域圖‘）;
figure（3）;
c=cwt（sawavename‘plot‘）;
%直接使用尺度
figure（4）imagesc（tfabs（coefs））;
%繪制色譜圖
colorbar;
xlabel（‘時間?t/ms‘）;
ylabel（‘頻率?f/Hz‘）;
title（‘小波時頻圖‘）;
figure（5）;
mesh（tfabs（coefs））;?
axis?tight;????
colorbar;
xlabel（‘時間?t/ms‘）;
ylabel（‘頻率?f/Hz‘）;?????
title（[‘小波時頻幅值圖‘‘（‘num2str（wavename）‘）‘]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????431??2018-03-29?11:35??cwt2.m
?????文件????????1109??2018-06-23?10:49??cwt1.m