% 設置小波函數、時間間隔和采樣點數
wname = 'db10';
A = -64; B = 64; P = 224;
% 計算采樣周期和采樣函數及真實頻率
delta = (B-A)/(P-1);
t = linspace(A,B,P);
omega = 5; x = cos(omega*t);
freq  = omega/(2*pi);
% 設置尺度并且使用scal2frq函數來計算準周期數列
scales = [0.25:0.25:3.75];
TAB_PF = scal2frq(scales,wname,delta);
% 計算最近似的準周期和相應的尺度
[dummy,ind] = min(abs(TAB_PF-freq));
freq_APP  = TAB_PF(ind);
scale_APP = scales(ind);
% 進行連續分解并繪圖
str1 = ['224 samples of x = cos(5t) on [-64,64] - ' ...
        '真實頻率 = 5/(2*pi) =~ ' num2str(freq,3)];
str2 = ['準周期數組和尺度: '];
str3 = [num2str([TAB_PF',scales'],3)];
str4 = ['準頻率 = ' num2str(freq_APP,3)];
str5 = ['對應尺度 = ' num2str(scale_APP,3)];
figure; cwt(x,scales,wname,'plot'); ax = gca; colorbar
axTITL = get(ax,'title');
axXLAB = get(ax,'xlabel');
set(axTITL,'String',str1)
set(axXLAB,'String',[str4,'  -  ' str5])
clc ; 
disp(strvcat(' ',str1,' ',str2,str3,' ',str4,str5))