构造B样条小波和基于小波变换模极大值（WTMM）的奇异点检测的Matlab实现

构造B样条小波和基于小波变换模极大值（WTMM）的奇异点检测的Matlab实现。

clear all;clc;
N=4;        % 设置B样条阶数，对应N-1次B样条小波
dual_N=2;   % 设置对偶阶数  N+dual_N之和必须为偶数，且dual_N>=1
syms z;
i=1;
n_min=-round（N/2）;
n_max=N-round（N/2）;
for n=n_min:n_max
    h（i）=2^（-N）*nchoosek（Nn+round（N/2））;
    i=i+1;
end
h             % 显示低通h（z）的系数，求滤波器系数时需乘以根号2
rat_h=rats（h）        % 分数显示
if rem（N2）==0
    dual_n=-（dual_N+（N-2）/2）:（dual_N+（N-2）/2）;
    L=0;
else
    dual_n=-（dual_N+（N-3）/2）:（dual_N+（N-1）/2）;
    L=1;
end
k=（N+dual_N）/2;
p1=z+z^（-1）-2;
sum=0;
for n=0:k-1
    p2=nchoosek（k-1+nn）*（-1/4）^n*p1^n;
    sum=sum+p2;
end
if rem（dual_N2）==0
   p3=z+z^（-1）+2;
   f=2^（-dual_N）*p3^（dual_N/2）;
else
    p3=z+1;
    p4=z+z^（-1）+2;
    f=2^（-dual_N）*p3*p4^（（dual_N-1）/2）;
end
dual_hz=f*sum;
% dual_hz=expand（dual_hz）    % 显示对偶低通dual_h（z）表达式，从而可以得出系数
if rem（N-32）==0
    a=（N-3）/2;      % 通过观察总结得出
else
    a=（N-2）/2;
end
dual_hz=z^（dual_N+a）*dual_hz;  % 由于dual_hz中含有负幂顶，为了利用sym2poly函数提取多项式系数，所以乘以z^n将最高负幂项变为0
dual_hz=expand（dual_hz）
dual_h=sym2poly（dual_hz）;  % 提取多项式系数
dual_h=fliplr（dual_h）      % 小数显示
rat_dual_h=rats（dual_h）    % 分数显示

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1756  2012-06-03 17:04  任意次样条小波的构造及模极大值检测奇异点例子\myb3.m

     文件        668  2012-06-02 21:53  任意次样条小波的构造及模极大值检测奇异点例子\myb2.m

     文件       1230  2012-06-04 11:13  任意次样条小波的构造及模极大值检测奇异点例子\bs_filter.m

     文件    1019392  2012-06-05 09:44  任意次样条小波的构造及模极大值检测奇异点例子\Cubic B-Spline Wavelets.doc

     文件     142336  2012-06-03 17:12  任意次样条小波的构造及模极大值检测奇异点例子\2-3次B样条小波构造.doc

     文件        347  2012-06-03 10:19  任意次样条小波的构造及模极大值检测奇异点例子\modulus_maxima.m

     文件      36123  2012-06-02 20:19  任意次样条小波的构造及模极大值检测奇异点例子\data.mat

     文件       1663  2012-07-14 20:40  任意次样条小波的构造及模极大值检测奇异点例子\data_test.m

     文件    2756260  2012-06-03 20:55  任意次样条小波的构造及模极大值检测奇异点例子\002 Biorthogonal bases of compactly supported wavelets.pdf

     目录          0  2012-07-14 20:27  任意次样条小波的构造及模极大值检测奇异点例子

----------- ---------  ---------- -----  ----

              3959775                    10