信号奇异性程序

% -----------------------------------------------------------------------－－
%   exa130202.m   例13.2.2
%  求冲激函数、阶跃函数及三角函数的小波变换、模极大线并求它们的李氏指数
%  注：在该程序中，用到了子程序 RWT.m，
%      该程序请读者在如下的网站上下载:
%                      http://www-stat.stanford.edu/~wavelab/
%      因为该网站上的小波分析软件并没有列入MATLAB的工具箱，涉及到知识产权问题，
%      因此不能将其列入国内公开出版物上，但读者可以自由下载。
%------------------------------------------------------------------------
close all;

nvoice = 16;
%   nvoice表示每2倍s之间被分为多少点

wavelet = ‘Sombrero‘; 
% 选择小波：可选择的小波包括： ‘Gauss‘‘DerGauss‘‘Morlet‘and‘Sombrero‘（mexh）
    N = 64*8;
    n = N;
    oct = 5;
    scale = 8;
   	noctave = floor（log2（n））-oct;
	nscale  = nvoice .* noctave;
    ytix   = linspace（2+（oct-floor（log2（scale）））log2（N）+2-floor（log2（scale））nscale）;
%  ytix 从imageRWT函数得到 理论上应等于log2（s）
%  ytix = linspace（log2（scale）-noctavelog2（scale）nscale）;
    xtix   = linspace（0NN）;
for count=1:3
   if count==1
        sig = zeros（1N）;sig（N/2） =sig（N/2） +1;% 冲激函数
        ti=‘冲激函数‘;
    elseif count==2
        sig = zeros（1N）;sig（N/4:3*N/4） = sig（N/4:3*N/4）+1; % 矩形窗
        ti=‘矩形窗‘;
    else
        sig（N/4:N/2） = linspace（01N/4+1）;
        sig（N/2:3*N/4） = linspace（10N/4+1）; % 三角窗
        ti=‘三角窗‘;
    end
    sig = sig+1e-8*randn（1N）;
    r