program_09_08.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

% 设置采样周期和小波函数
delta = 0.1; wname = 'coif3';
% 定义尺度
amax = 7;
a = 2.^[1:amax];
% 计算相关的准频率
f = scal2frq(a,wname,delta); 
% 计算相关的准周期
per = 1./f; 
% 绘图尺度－准周期
subplot(211), plot(a,per)
title(['小波: ',wname, ', 采样周期: ',num2str(delta)])
xlabel('尺度')
ylabel('准周期')
% 对于每个尺度 2^i：
% - 产生周期为 per(i)的正弦函数
% - 进行小波分解
% - 确定最大能量
% - 计算检测出的准周期
for i = 1:amax
   % 以采样周期delta，产生周期为的per(i)周期信号Generate sine function of period
   t = 0:delta:100;
   x = sin((t.*2*pi)/per(i));
   % 对 x 做第9级分解
   [c,l] = wavedec(x,9,wname);
   % 估计系数的标准差
   stdc = wnoisest(c,l,[1:amax]);
   % 计算检测周期
   [y,jmax] = max(stdc);
   idper(i) = per(jmax);
end
% 比较检出准周期和计算准周期
subplot(212), plot(per,idper,'o',per,per)
title('检出准周期和计算准周期')
xlabel('计算准周期')
ylabel('检出准周期')