% 設置小波函數、時間間隔和采樣點數
wname  = 'morl';
A = 0; B = 64; P = 500;
% 計算采樣周期和采樣函數及真實頻率
t = linspace(A,B,P);
delta = (B-A)/(P-1);
tab_OMEGA = [5,2,1];
tab_FREQ  = tab_OMEGA/(2*pi);
tab_COEFS = [5,3,2];
x = zeros(1,P);
for k = 1:3;
    x = x+tab_COEFS(k)*sin(tab_OMEGA(k)*t);
end
% 設置尺度并且使用scal2frq函數來計算準頻率數列
scales = [1:1:60];
tab_PF = scal2frq(scales,wname,delta);
% 計算最近似的準周期和相應的尺度
for k=1:3
  [dummy,ind] = min(abs(tab_PF-tab_FREQ(k)));
  PF_app(k) = tab_PF(ind);
  SC_app(k) = scales(ind);
end
% 進行連續分解并繪圖
str1 = strvcat( ...
   '500 samples of x = 5*sin(5t)+3*sin(2t)+2*sin(t) on [0,64]',...
   ['真實頻率 (Hz): [5 2 1]/(2*pi) =~ [' ...
   num2str(tab_FREQ,3) ']' ]  ...
   );
str2 = ['準周期和尺度:  '];
str3 = [num2str([tab_PF',scales'],3)];
str4 = ['準頻率   = ' num2str(PF_app,3)];
str5 = ['相應尺度 = ' num2str(SC_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',strvcat(str4, str5))
clc; 
disp(strvcat(' ',str1,' ',str2,str3,' ',str4,str5))