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