小波总结

小波总结
小波总结

关于小波的总结

作者:madong123 转自:研学论坛

学习小波4个月了,也来写个小结,有两个目的,第一个当然是完成强化课的作业,其次就是整理下思路,理清脉络,进一步对相关知识模型化,定型于脑中,把自己的理解曝于众人,接受反馈,

然后归纳没有搞明白的问题,有的放矢。

(一)傅立叶变换:

(1)感觉学习小波最大的困难就在于充分地去理解那些公式背后的意义,由于知识背景上的差异,同一个公式不同的人得到的信息量是不一样的。就拿最普通的内积公式来说,< f,g>=……,最

初接触是用来刻画信号之间的相关性,两个信号内积归一化后就得到相关系数,也就是两个信号越相似,内积就越大,这是第一个理解;然后我们会在Hilbert空间里看到这个东西,用来刻画两个向量的

夹角,当内积为0时,两个向量正交,若g为Hilbert空间里的正交基的时候,内积为f向基上的正交投影,这是第二个理解;(Hilbert空间是一个很直观的空间,我一直都理解为欧氏空间去理解定义在其

上的东西,L^2和l^2同样为Hilbert空间)。

(2)接下来不得不提的就是基和框架的概念:

基的概念大家都有,也很常见,只是很少涉及到数学的严格性,大家都知道颜色的RGB分解法,其实就是基的一种表示,取定RGB各自的比重就可以条配出任何一种颜色,信号和函数的表达也是如此。人

们一直渴望找到一组基来表示所有信号,从三角基,DCT,小波,脊波,曲波,基越来越复杂,表示能力也越来越强,当然计算复杂度也大幅度增加,似乎这已经逐渐成为一个不可调和的矛盾,好在计算

机的能力也在按摩尔定律增长。

基具有非冗余性,即使基不是正交的,有相关性,但若去掉其中任何一个,则不成为基,这一点也叫完备性;基的表示有唯一性,即给定一族基对一个函数的表达是唯一的;一般情况下基非正交,也称

为为exact frame(Resize basis),这个时候要表示信号可以将基正交化成唯一的正交基(对偶为其自身);也可以求其对偶框架(dual frame),其对应了小波变换中的双正交情形!信号可以依框架

分解,然后用对偶框架重构。

若在基集里添加一些新的向量,并随意调整空间位置,则有可能成为框架。把函数与基或框架作内积,也可以说成是一种函数空间到系数空间的变换。若某种变换后的能量(内积的平方和度量)仍然有

一个大于0的上下界,才可以成为框架,由于框架的冗余性,所以系数的表达也不具有唯一性。若上下界相等,则为紧框架,且界表示冗余度。若上下界相等为且为1,称为pasval identity frame,此时

不一定为正交基(想象把一组正交基中某一个拆成两个同方向的基之和,则pasval identity 仍然成立),此时若加上基的长度均为一的条件,则框架退化为正交基。可能你会问我们用基来表示信号就行

了啊,为什么还要框架呢?其实很多信号表示方法不能构成基,却能构成框架,如短时傅立叶变换中如要求窗函数满足基条件,则可推出该函数有很差的时频局部化性质(事实上退化为了傅立叶变换)

,框架的冗余性为我们提供了很多自由度,如果有时间的话会专门写一篇框架应用的总结。(3)傅立叶变换

傅立叶变换将函数投影到三角波上,将函数分解成了不同频率的三角波,这不能不说是一个伟大的发现,但是在大量的应用中,傅立叶变换的局限性却日趋明显,事实上在光滑平稳信号的表示中

,傅立叶基已经达到了近似最优表示,但是日常生活中的信号却并不是一直光滑的,而且奇异是平凡的,傅立叶在奇异点的表现就着实让人不爽,从对方波的傅立叶逼近就可以看出来,用了大量不同频

率的三角波去逼近其系数衰减程度相当缓慢,而且会产生Gibbs效应。其内在的原因是其基为全局性基,没有局部化能力,以至局部一个小小的摆动也会影响全局的系数。实际应用中很需要时频局部化,

傅立叶显然缺乏此能力了。即使如此,由于其鲜明的物理意义和快速计算,在很多场合仍然应用广泛。

傅立叶变换在从连续到离散的情形是值得借鉴与学习的,大家都知道,时间周期对应频域离散,时间离散对应频域周期,时间离散+周期对应频域离散+周期,DFT其实是将离散信号做周期延拓然后做傅立

叶变换再截取一个周期,反变换同样如此,所以DFT用的是块基的概念,这样如果信号两端的信号连接后不再光滑(即使两边都光滑),同样会在边界上产生大幅值系数(边界效应),延伸到图像中就是

块效应。

当对信号做对称+周期延拓后再做傅立叶变换得到的正弦系数全部为0,也就是任何对称函数可以写成余弦的线性组合,同样按照离散的思路构造得到的是离散块余弦基,即DCT变换,虽然DCT可以通过对

称后周期延拓再变换减少了边界效应(两边信号接上了,但不一定平滑),但任不能消除块效应,尤其是图像变换中人为将图像分成8*8处理后块效应更加明显。但是DCT很好的能量聚集效应让人惊奇,

加之快速计算方法使它替代DFT成为图像的压缩的标准了很长时间(JPEG)。

(二)短时傅立叶变换到小波变换

三角基的局部化性质这么差,但是很多应用场合又要求比较精确的时频定位,傅立叶的缺点就越来越突出了,窗口傅立叶变换将信号乘上一个局部窗,然后再做傅立叶变换,获得比较好的时频定

位特性,再沿时间轴滑动窗口,得到整个时间轴上的频率分布,似乎到这里就应该结束了,因为我们可以把窗设计小点获得较高的时间分辨率,并期望有同样高的频率分辨率,但测不准原理无情的告诉

我们,没有这么好的窗能在时间和频率都任意小的,最优的就是gaussian window了(窗的选取还需满足频率域也为窗函数,并不是每个时窗都满足这个条件的)。

通过短时傅立叶变换我们可以画出时频图,但是存在问题:当我们分析频率较高部分信号时应该用更窄的窗,反之用宽窗,但STFT一旦选定窗过后,分辨率就固定了,若要其他分辨率则需要更换窗。接

下来用于分析窗函数的平移本身不能构成基,没有简化计算的可能性,使得时频分析的计算量一直很大(若为正交基,系数的计算相当方便)。另外一个问题:由于时间和频率都使用连续表达,连续窗

口傅立叶变换具有极大的冗余性,怎样去离散时间和频率参数以减少冗余,而又不导致信息丢失,一个明显的要求就是时频盒子一致时间和频率平移必须完全覆盖整个时频平面。由前面的框架分析可以

得知,离散窗口调制不能成为基,但可构成框架(时频采样密度大于临界值,即盒子的有效铺叠刚好邻接并充满整个时频空间),并当时频采样密度为临界采样率一半的时候(盒子有大量重叠),框架

差不多是紧的,并可证实冗余因子为4。

针对短时傅立叶变换的第一个问题,马上想到我们能不能用一个窗函数的伸缩(S)和平移(以前只有平移b)来分析信号,要求这个变换是可以完全重构的且保持能量守恒,如果存在,这样的窗

函数应满足什么条件呢?答案当然是能,这个窗函数就是连续小波,其应满足容许条件即可达到完全重构。容许条件暗含该函数在0点的傅立叶变换为0,这也解释了小波为什么必须有零值平均。由重构

公式可以看到尺度s的积分范围为(-inf,+inf),而我们分析一般只取(1,So)对于大于So以上的信息(对应实际应用的过高频率),我们认为已经不重要了,予以丢弃也不会对我们的分析造成太大

的影响,对于尺度小于等于1的信息(低频成分),我们需要对其截断保留,以便恢复信号。这可通过引入尺度函数(低通滤波)来做到。这样连续小波变换就得到一个在不同尺度和不

同时间的系数图,

系数较大的地方说明信号在该时间与该尺度的窗有很大的相似性(由变换公式和前面提到的内积的含义可知)。这样我们就可以根据系数图来分析信号的特性了。

但是连续小波变换存在很大的冗余性,事实上由连续小波变换的模极大值就可完全稳定重构信号(mallat的模极大迭代投影重构),很自然就就可以想到对尺度和平移的离散,这一想法导致了小波

框架,离散要求满足时频窗完全覆盖整个时频平面,Daubechies提供了在给定离散参数s和b的情况下估计小波框架的上下界的定理,在满足一定条件的小波函数和一定离散参数选择下,框架界可近似达

到1(当为二进离散正交基时为1),即最大程度上消除了冗余性,此即二进离散小波;但是此举将小波的(任意)平移不变性去掉了(信号平移与采样网格未必一致),即信号移动后再做小波变换后系

数发生了很大的变化。如果想保持小波的平移不变性,可以只对尺度进行离散化但平移参数不离散,取尺度序列为二进序列{2^j}即为具有平移不变性的二进小波变换,但得到一个高冗余度的信号表示(

时间未被离散)。怎样继续较少冗余而又保持平移不变性呢?这就要求将采样网格自适应地作与信号相同的平移,也就是信号的平移必须以采样网格为参考。后面我们看到非抽取小波变换在离散了平移

参数情况仍然保持了平移不变性,这是因为信号和小波均离散化后采样网格与信号平移是一致的。

注:在对离散二进小波做重构的时候,需要求其对偶框架,而不是用离散二进小波框架本身来做重构(双正交情形,参看框架部分),但是框架为二进离散正交小波时,对偶为其自身。

到这里,短时傅立叶分析中存在的毛病在这里已经全部解决了,新问题产生了,满足允许条件的离散二进小波怎么构造,其对偶框架又怎么获得,理论上非正交二进小波有无穷多对偶,这带来的是不确

定性还是灵活性?幸而MRA和完全重构滤波器组理论为二进(正交)小波提供了构造方法和快速实现。

在MRA分析出现以前,滤波器组与小波理论是没有融合到一起的,没有人注意到完全重构滤波器中的HG的某种迭代竟是小波的构造和变换的快速实现方法,以至于小波和完全重构滤波器的构造走在

两条不相干的路上。最简单的Harr小波因为其连续性太差得不到好的应用,Meyer试图证明光滑连续的小波不存在,却意外的得到了具有一定衰减性的光滑函数,其二进伸缩和平移构成了L^2空间的标准

正交基!Mallat于同年提出了多分辨分析,为小波的构造提供了统一框架,并提供了函数

分解与重构的快速算法。

(三)MRA(多分辨分析)

我们先来考虑下任何函数的二进网格采样f(2^j*t),随着j的增大,采样网格越来越密,由采样值重构函数的空间也越来越大(频域范围越来越大),可以想象当j趋于+inf 时,L^2中的任何函

数都可以由其采样值重构。j趋于-inf时,这个空间只能包含零空间。直观理解f(2^j*t)就是在分辨率j上的逼近,也可以理解成f(t)在空间Vj上的正交投影。

注:可以这样想象,每个采样点是空间中的一个坐标轴,Vj就是由这些坐标轴张成的一个空间,j每增大一次,坐标轴就越来越多,坐标张成的空间也越来越大,且原来的坐标轴也包含在更大的空间里面

,一般的连续函数f是在一个很大的空间里面的,这样把f向Vj垂直投影,得到的是在Vj 空间上对函数的最佳逼近,可以拿三维的情况想象,再进行扩展。

若Vj空间的有以下几个性质,就称为MRA:

a)属于Vj的函数随网格平移仍属于Vj。(这个条件暗含了平移不变性)

b)空间是随j的增大二进膨胀的且严格包含(一个比一个大,且逐层包含),表明Vj-1的信息完全包含在Vj中且比Vj包含信息要少。

c)若f(t)属于Vj,则f(2t)属于Vj+1,

d)j趋于+inf,Vj能逼近属于L^2的所有函数,j趋于-inf时,Vj丢失所有细节信息,成为零空间。

e)V0存在riesz基{m(t –n)},即属于V0的f都可以由{m(t –n)}线性组合构成。

则称m为尺度函数,并称m生成了一个多分辨分析{Vj}。可以预见,不同的m张成的多分辨分析性质是不一样的,如由Haar函数张成的空间连续性非常差(分段线性函数空间),用这个空间去逼近

连续函数的效果是很差的。理所当然想找到好的m使得张成的空间有更好的性质,更适合逼近一般的函数。

前面的讨论要求我们构造小波函数,使其二进伸缩平移函数族构成L^2的一组基。怎么构造呢?在MRA背景下,我们可以先找到V0的riesz基{m(t –n)}(可以证明Vj空间的riesz基是{2^(j/2)*m(2^j*t

–n)})。MRA框架下,Vj是一系列嵌套且随j不断(二进)膨胀的空间,Vj与Vj-1的差空间记为Wj-1空间,很显然,Wj-1空间的任何函数(包括基)可由Vj的基构成,这样,我们就找到了每个Wj空间的基

,把所有Wj空间的基合起来就构成了整个L^2的Riesz基(而且不同空间的Wj是正交的,空间不相交)。实际上,我们一般从V0开始,把V0的基和的Wj(j>0)空间的基合起来作为L^2空间的Riesz基。

与Vj的基类似,我们当然希望由一个函数g(t)的二进伸缩平移构成整个W空间的基,当然MRA保证了有这样的函数存在,它就是小波函数。如果m(t)二进伸缩平移构成的是V0规范正交基,则可以得到

L^2的规范正交基。

由空间的关系和我们的希望可以写出几个关系:

a)Vj中的基函数可以由Vj+1的基线性组合构成(尺度方程)

b)Wj的基函数可以由Vj+1的基线性组合构成(小波方程)

c)不同尺度的小波基是正交的(不同尺度空间分割关系)

d)同一尺度上的小波基与尺度函数是正交的(同一尺度空间分割关系)

e)同一尺度中的不同的基之间是正交的(正交基关系),若无此条件则为双正交情形。

到这里就有两条路摆在面前,一是找尺度函数,然后得到小波函数,迭代求极限可得到尺度函数和小波函数。对于第一条路,我们根据关系(a)找到h,然后由(d)找到g,再由(b)求得小波函数。但

是一般情况下可供选择的尺度函数都太少,而且求得的小波函数性质都比较差。

第二条路是滤波器组,通过找到满足条件的h,g,迭代求极限可得到尺度函数和小波函数。

可以说滤波器组完美重构理论为我们提供了小波构造的一般方法。

(四)滤波器组完美重构与小波快速算法。

前面的分析可以知道Vj相当于在j分辨率的逼近,Vj-1相当于j-1分辨率的逼近,这样Wj-1相当于两个分辨率逼近的差。在高分辨率下,我们可以用f在(2^j*t)的采样值来代替向Vj空间的投影,

但是这是需要说明的,否则成为“小波的罪恶”,本来在Vj上的投影需要函数对Vj上的基{2^(j/2)*m(2^j*t –n)}投影,用采样值来代替是因为当j足够大的时候,如一般情况下j=7~9时,尺度函数已

经非常窄,以致用delta采样来表示误差可以忽略(其实数学推导中还因为尺度函数的消失矩有关),但是必须理解这个取代是近似的。

得到了第j层的系数,如何求得j-1层的逼近系数和两个逼近层次间的误差系数呢?这就是mallat由MRA得出的快速算法。这个算法由上面列出的几个空间和相应基的关系很容易得出。同样重构算法也可以

推导得出。(具体可以查看任意一本小波书)

到这里是不是就结束了呢,那这样的话MRA也得不到这么大的名声了,它的伟大之处是与完美滤波器重构桥接上了,从而为小波构造提供了一个普遍的方法,个人认为双正交小波也是由滤波器组理论发展

出来。

把快速算法的一级分解和重构结构画出来,这不就是一个完美滤波器重构么?之前对完美滤波器的重构的结论大都可以搬上来了,大家熟知的两个PR方程其实也可以通过MRA下的空间关系推出。当然单由

这两个方程得到的h和g有很多解,并不是每个解都可以收敛到尺度函数和小波函数,必须附加其他条件,其中以Daubechies的p阶消失矩条件构造出的小波应用得最多(Daubechies 系列小波)。

很遗憾,除了Haar小波以外(haar小波可由一阶消失矩条件构造出来),没有正交小波满足对称性条件,也就是不满足线性相位,这样在分解重构后会造成失真,在一些需要对称性的场合(如图像的分

解重构,奇异点的检测等),结果是不能满足要求的。

为了构造具有光滑特性,一定消失矩,对称的小波,就不得不放弃正交条件,也就是前面提到的双正交多分辨分析

在正交情形下,我们只需要知道H0,就可以由共轭镜像滤波器条件推导得出其他滤波器为G0,H1,G1(是H0的逆序及调制),也就是我们只需要知道一个滤波器。在双正交情形下,由完美重构滤波器条

件可从H0,G0推导出H1和G1(通过逆序及调制),这表示我们需要知道两组滤波器。{(H0,G0)(H1,G1)}这两组(双)正交滤波器是可以对调的,就是谁做分解另一组就做重构。由滤波器来造小波

的步骤前面已经提及!四个滤波器之间的关系是相互交叉的,即G1由H0逆序调制,而H1由G0逆序调制。

当然我们希望尺度函数和小波是紧支的(暗含滤波器也是紧支的),否则在计算时需要进行截断,正交紧支小波的对偶为其自身,当然也是紧支的,但是有一个定理:紧支非正交小波,其对偶必然是无限

支集的,可能你会很奇怪,我们平时用的双正交小波不都是紧支的么,其实这些紧支双正交小波是经过提升的,Daubechies有一个定理,任何双正交滤波器可通过对惰性滤波器不断做提升和对偶提升而

生成!你只需要了解这一事实即可,深入的理解恐怕需要太多的数学知识。至于提升我也希望我能有时间写一个总结,从框架的角度来理解提升,恐怕会容易得多,只是这个愿望可能不太好实现,因为

这必须要大量的推导来表述。

最后说说消失矩这个条件,Haar小波得不到应用是因为它的消失矩为1,也就是对大于一次多项式的函数的“消失”效果不好,所谓消失矩其实就是对多项式的抑制能力,消失矩越高,与信号做内积得到

的系数越少越小,这在度量信号局部正则性和压缩方面是相当重要的。提升就是一种提高小波消失矩和正则性的及其重要的手段。

注:由mallat算法得到的快速算法不具有平移不变性,其中的原因是因为采样因子是不具有平移不变性的。如果需要保持平移不变性,则需要去掉抽取这一步(多孔算法),其实不是去掉,如果去掉了

,就不是小波变换了,而是利用noble identity把抽样算子移到每一分支的分解完全结束之前而已。得到的是每次分解得到原来两倍长信号,而mallat算法是每次抽掉了其中的一部分,这样随着分解层

次的增加,小波系数也越来越稀。多孔也抽,不过移到最后一起抽,然后在重构之前同数目的上采样。

注:实际编程实现的时候由于要做滤波器卷积,每次卷积完后要用wkeep保持原来的长度。

(五)小波的性质及构造

很明显随附加条件不同,我们可以构造出无穷多的正交或双正交小波,Mallat说过,如果没有应用的刺激,小波的构造将变成聊的游戏,这里我们先分析小波具有的性质,然后谈谈如果根据应用

来构造(更多情况是选择)小波。

前面虽然给出了小波构造的统一方法,但感觉太笼统,太灵活,这需要些著名的定理来逐步加强认识,既然是定理,大家就试着接受,记住,再理解好了。这里我们只谈紧支正交小波的构造,因为紧支

双正交小波的构造已经涉及到了提升的概念。

性质1:消失矩(Vanishing Moment),这可以说是小波最具杀伤力的一个性质,压缩,去噪,快速计算等无不希望小波VM越高越好,虽然是通过滤波器卷积来求小波系数,但是思考上仍然用信

号与小波的内积来表示,这样有助于理解小波的性质。由VM的定义可知具有p阶消失矩的小波与小于p次的多项式是正交的,也就是内积为0,这样若函数f是正则的且小波有足够的消失矩,则内积产生很

小的系数。注,讨论函数的局部正则性其实是一个比较复杂的问题,在这里,姑且将f在某点附近想成一个k阶多项式和一个误差函数的逼近,这样当k

f与误差函数产生较小的系数。

小波的消失矩与对应滤波器h的傅立叶变换在pi处的零点重数是等价的,事实上我们也用这一点作为附加条件来构造具有p阶消失矩的小波(Daubechies系列的构造方法),当然我们可以构造出任

意消失矩的小波,但是我们不得不注意到滤波器长度随小波消失矩增加而增加的事实,时的我们同样只能在支集与消失矩之间折衷!

性质2:正则性,小波的正则性(光滑程度),即使理论上的完全重构,在计算机实施的时候会引起由于量化和截断造成的误差,这样在重构的时候,若小波基不够光滑,则引入的误差很容易被人察觉,

如haar小波,若小波基足够正则,引入的光滑误差不容易被察觉。有结论说明,对重要的共轭镜像滤波器族,如样条和Dx系列,小波的正则性随消失矩提高而提高

性质3:紧支性,我们当然希望我们的滤波器越短越好,因为这意味着计算量的大大减少,同时考虑到小波滑动作内积时包含到函数奇异点的时候同样可能造成大幅值系数(支集越长,包含奇异

点的小波次数越多),这样若想大幅值系数数目最小,必须尽可能减小支集长度。前面提到了,紧支与消失矩是矛盾的,为了得到更小的系数,需要高的消失矩,而要得到更少的系数,希望小的支集。

有定理在给定的消失矩情形下,Dx系列小波具有最小支集

在这个意义上,这个系列小波是最优的。

注:紧支性和消失矩可以在多小波构造下得到更好的权衡,在实际情况中,我们更多是根据信号的奇异程度来权衡,若信号奇异点很少,则考虑高的消失矩,若奇异点很平凡,则考虑更短的支集。

性质4:对称性,对正交小波来说,除haar外不存在对称或反对称小波,这点daubechies 已经严格证明了,而且由消失矩条件通过最小相位构造的小波是极不对称的,这在某些应用中是不好的性质,如奇

异点的确定,图像的重构等等,至于symmlet滤波器虽然更对称,但是产生了复小波系数双正交小波的性质与正交小波是类似的,但是可以做到完全对称,上面讨论的性质对双正交小波依然适用,前面提到双正交情形两组滤波器的位置是可以互调的,在这里我们一般选用消失矩高的

小波来做内积,然后用另一个小波做重构(它通常都是最光滑的那个)

在了解了小波的各种性质之后,我们就可以根据不同的应用来构造或选择我们适合的小波了,这在后面具体的应用中会反复提到。

(六)再谈滤波器与小波的关系

大家肯定都熟知了小波构成了L^2的(双)正交基,我们习惯在脑海中把小波系数的幅度看成未被采样的函数和小波之间的相似性度量,这也是我们获得的清楚的物理意义,那么滤波器组的角色仅

仅是提供一个快速计算?那就太浅显了,我们回忆下mallat的快速算法,对函数进行采样后近似表示系数,是不是可以考虑成l^2(Z)中的函数呢?当然可以!那l^2(Z)的基是什么呢?

再回头看看快速分解的公式:第j层的低频系数Aj(n)与滤波器h(n)卷积后做下采样得到j-1层低频系数Aj-1,这个过程可以写成第j层的系数Aj(n)与下采样的滤波器h(2n)做卷积,也可以

写成Aj-1(k)= < Aj(n),h(2k-n)>,这个形式是不是很熟悉呢,仔细看是l^2(Z)中的函数在基h(2k-n)下投影,系数为Aj-1(k);同理可以得到第j-1层细节系数Dj-1(k)= < Aj(n),g

(2k-n)>,同样是l^2(Z)中的函数Aj(n)在基g(2k-n)下投影,那这个h和g是不是就是我们要找的离散小波基呢?

可以由滤波器完美重构的条件推得如下结论:

如果h(-n),g(-n),h1,g1是完全重构滤波器组且傅立叶变换都有界,则{h(2k-n),g(2k-n)|k属于Z}和{h1(2k-n),g1(2k-n)|k属于Z}构成了l^2(Z)的双正交resize基如果上面的h(n)=h1(n),g(n)=g1(n),则{h(2k-n),g(2k-n)|k属于Z}构成了l^2

(Z)的规范正交基。

这个结论的含义是什么呢?我们找到了l^2(Z)的离散小波(双)正交基!而且这些滤波器基按照小波树形分解结构得到的基仍然是l^2(Z)的离散小波(双)正交基!

注,好好理解这两句话,可能你需要对正交基或者双正交基的好处有一些体会!其实同样的方法可以构造离散的小波包的基。

大家都知道我们计算机处理的都是离散信号,以前的快速算法似乎通过将采样数据离散来近似高分辨率数据,然后通过连续小波基之间的尺度关系来推出的快速算法,而现在我们可以完全在离散

情况下来考虑这个问题了,因为我们从大自然采集离散数据本身就很方便,我们只需要满足采样定律来保存原信号的信息;同时我们把滤波器看成基,这对以后的提升理解是有帮助的(对这句话有兴趣

的可以和我交流)。

那滤波器基和小波基有什么关系呢,考虑对尺度方程和小波方程两边做傅立叶变换,可以通过反复迭代取极限求得尺度函数和小波基的傅立叶变换,也就是我们把滤波器通过某种方式无穷迭代最

后会收敛到尺度函数和小波基,这个就叫cascade算法吧,事实上我们经过三五次迭代的结果就与连续尺度函数和小波基非常相似了。

注:并不是所有的h,g都可以最后收敛的,必须加条件,但我认为这样的h,g同样也可以用来做分解的,这似乎决定了构造离散的基比连续情况的基要简单。至于用不能收敛的h,g来分析有什么

后果我还没有深入研究过。

其实以前就接触到很多结论中就有尺度函数和小波与对应滤波器组之间的关系了,如消失矩,滤波器在pi处的零点重数与小波基的消失矩是对应的,这些都是通过尺度方程和小波方程联系起来的

,可以说这两个方程为我们的连续与离散架起了一座桥梁

(完整版)MATLAB常用函数大全

一、MATLAB常用的基本数学函数 abs(x):纯量的绝对值或向量的长度 angle(z):复数z的相角(Phase angle) sqrt(x):开平方 real(z):复数z的实部 imag(z):复数z的虚部 conj(z):复数z的共轭复数 round(x):四舍五入至最近整数 fix(x):无论正负,舍去小数至最近整数 floor(x):地板函数,即舍去正小数至最近整数ceil(x):天花板函数,即加入正小数至最近整数rat(x):将实数x化为分数表示 rats(x):将实数x化为多项分数展开 sign(x):符号函数(Signum function)。 当x<0时,sign(x)=-1; 当x=0时,sign(x)=0; 当x>0时,sign(x)=1。 rem(x,y):求x除以y的馀数 gcd(x,y):整数x和y的最大公因数 lcm(x,y):整数x和y的最小公倍数 exp(x):自然指数 pow2(x):2的指数 log(x):以e为底的对数,即自然对数或 log2(x):以2为底的对数 log10(x):以10为底的对数 二、MATLAB常用的三角函数 sin(x):正弦函数 cos(x):余弦函数

tan(x):正切函数 asin(x):反正弦函数 acos(x):反馀弦函数 atan(x):反正切函数 atan2(x,y):四象限的反正切函数 sinh(x):超越正弦函数 cosh(x):超越馀弦函数 tanh(x):超越正切函数 asinh(x):反超越正弦函数 acosh(x):反超越馀弦函数 atanh(x):反超越正切函数 三、适用於向量的常用函数有: min(x): 向量x的元素的最小值 max(x): 向量x的元素的最大值 mean(x): 向量x的元素的平均值 median(x): 向量x的元素的中位数 std(x): 向量x的元素的标准差 diff(x): 向量x的相邻元素的差 sort(x): 对向量x的元素进行排序(Sorting)length(x): 向量x的元素个数 norm(x): 向量x的欧氏(Euclidean)长度sum(x): 向量x的元素总和 prod(x): 向量x的元素总乘积 cumsum(x): 向量x的累计元素总和cumprod(x): 向量x的累计元素总乘积 dot(x, y): 向量x和y的内积 cross(x, y): 向量x和y的外积 四、MATLAB的永久常数

小波变换图像去噪综述

科技论文写作大作业小波变换图像去噪综述 院系: 班级: 学号: 姓名:

摘要小波图象去噪已经成为目前图象去噪的主要方法之一.在对目前小波去噪文献进行理解和综合的基础上,首先通过对小波去噪问题的描述,揭示了小波去噪的数学背景和滤波特性;接着分别阐述了目前常用的3类小波去噪方法,并从小波去噪中常用的小波系数模型、各种小波变换的使用、小波去噪和图象压缩之间的联系、不同噪声场合下的小波去噪等几个方面,对小波图象去噪进行了综述;最后,基于对小波去噪问题的理解,提出了对小波去噪方法的一些展望 关键词:小波去噪小波萎缩小波变换图象压缩 1.前言 在信号数据采集及传输时,不仅能采集或接收到与所研究的问题相关的有效信号,同时也会观测到各种类型的噪声。在实际应用中,为降低噪声的影响,不仅应研究信号采集的方式方法及仪器的选择,更重要的是对已采集或接收的信号寻找最佳的降噪处理方法。对于信号去噪方法的研究可谓是信号处理中一个永恒的话题。传统的去噪方法是将被噪声污染的信号通过一个滤波器,滤除掉噪声频率成分。但对于瞬间信号、宽带噪声信号、非平稳信号等,采用传统方法具有一定的局限性。其次还有傅里叶(Fourier)变换也是信号处理中的重要手段。这是因为信号处理中牵涉到的绝大部分都是语音或其它一维信号,这些信号可以近似的认为是一个高斯过程,同时由于信号的平稳性假设,傅立叶交换是一个很好的信号分析工具。但也有其不足之处,给实际应用带来了困难。 小波变换是继Fourier变换后的一重大突破,它是一种窗口面积恒定、窗口形状可变(时间域窗口和频率域窗口均可改变)的时频局域化分析方法,它具有这样的特性;在低频段具有较高的频率分辨率及较低的时间分辨率,在高频段具有较高的时间分辨率及较低的频率分辨率,实现了时频窗口的自适应变化,具有时频分析局域性。小波变换的一个重要应用就是图像信号去噪。将小波变换用于信号去噪,它能在去噪的同时而不损坏信号的突变部分。在过去的十多年,小波方法在信号和图像去噪方面的应用引起学者广泛的关注。本文阐述小波图像去噪方法的原理,概括目前的小波图像去噪的主要方法,最后对小波图像去噪方法的发展和应用进行展望。 2小波图像去噪的原理 所谓小波变化,即:

营销战略开题报告

营销战略开题报告 引导语:开题报告是指开题者对科研课题的一种文字说明材料。这是一种新的应用写作文体,这种文字体裁是随着现代科学研究活动计划性的增强和科研选题程序化管理的需要而产生的。接下来是为你带来收集整理的营销战略开题报告,欢迎阅读! 摘要:在国际化的发展过程中,为了满足市场需求,本土化发 展战略已成为跨国公司开拓海外市场纷纷采取的主要措施。本文以沃尔玛中国营销策略为背景,从适应环境、降低成本、满足消费者需求、树立良好形象等层面对沃尔玛在中国实施营销策略的必要性进行解析,并针对其目前在华的管理团队、分销渠道管理和市场战略本土化的现状,分析沃尔玛公司在中国经营中遇到的困难及存在的主要问题,最后提出沃尔玛在华经营的几点建议,即采用“外包”和“自营”相结合的物流配送、协助供应商加强信息化管理、跨文化沟通的社会责任以及政府公关。 关键词:沃尔玛、跨国企业、中国战略、问题策略 1研究背景 1.1研究背景 沃尔玛百货有限公司是由美国零售业的传奇人物山姆·沃尔顿 先生于1962年在阿肯色州成立。经过四十多年的发展,公司已经成 为美国本土最大的私人雇主和世界上最大的连锁零售企业。到目前为止,沃尔玛在全球15个国家开设了超过8,000家商场,下设53个品牌,员工总数210多万人,每周光临沃尔玛的顾客2亿人次。1991

年,沃尔玛年销售额突破400亿美元,成为全球大型零售企业之一。据1994年5月美国《财富》杂志公布的全美服务行业分类排行榜,沃尔玛1993年销售额高达673.4亿美元,比上一年增长118亿多,超过了创始人:山姆·沃尔顿1992年排名第一位的西尔斯(Sears),雄居全美零售业榜首。1995年沃尔玛销售额持续增长,并创造了零售业的一项世界纪录,实现年销售额936亿美元,在《财富》杂志95美国最大企业排行榜上名列第四。至今,沃尔玛己拥有2133家沃尔玛商店,469家山姆会员商店和248家沃尔玛购物广场,分布在美国、中国、墨西哥、加拿大、英国、波多黎各、巴西、阿根廷、南非、哥斯达黎加、危地马拉、洪都拉斯、萨尔瓦多、尼加拉瓜等14个国家。它在短短几十年中有如此迅猛的发展,不得不说是零售业的一个奇迹。 从沃尔玛角度来看,其拥有先进技术、专利和充足的资金、先进的管理经验等优势,实施“本土化”经营可以把产品的生产、采购转移到中国,充分利用中国本土市场的资源,利用中国市场制造成本低、生产成本低、人力资源成本低的有利条件,生产和采购产品,使企业迅速融入中国市场[],享受本土化带来的好处,为沃尔玛“天天平价”的经营理念创造有力条件,从而在激烈的市场争夺中获得竞争优势。 沃尔玛进入中国后,势必面临与其在美国、在世界其他国家不同的市场环境,这主要表现在社会文化、企业文化、语言、消费观念、政府的法律制度和倾向、地理条件、原料供应、市场替代者等很多方

1关于小波基函数选择的相关研究

小波基函数的选择 1011208041 材料学院冯梦楠 多分辨率分析方法使小波分析成为一种实用的信号分析工具。同传统的Fourier 变换相比,理论上来说小波变换可以刻画信号的任意细节,但在实际应用中,信号分析的好坏很大程度上依赖小波基波的选择。因为与Fourier变换不同,小波基不具有惟一性,它是不规则的,不同的小波基波波形差别很大,其支撑长度和规则性也有很大的差别。因此,对同一个信号选用不同的小波基进行处理所得的结果往往不尽相同。同时,小波变换又是一种在基波可变的情形下其尺度仍可变的信号分析方法,它可在不同尺度下对信号进行分析处理。因此这也意味着即使小波基选定,如尺度选择不当,对信号分析的效果仍然会有一定的影响。因此,最优小波基函数的选择就成为了小波分析在工程应用中的一个重要问题。 小波基函数的选择是一个重要而复杂的问题,它受到测不准原理、小波基函数的性质和具体应用的特点等多方面的综合制约。因此如何选择小波基函数,到目前为止还没有一个统一的理论标准。在实际工程应用中,通常是根据具体问题的特点,结合小波基函数的性质和时频测不准原理进行经验性的选择。如Morlet小波一般用于信号的表示和分类图像识别、特征提取;Mallat小波多用于系统辨识;样条小波则常用于材料探伤;Shannon正交基用于差分方程求解:对于数字信号则往往选择Haar或Daubechies作为小波基。同时我们也通过用小波基函数处理信号的结果与理论结果的误差,来判断小波基函数的好坏,并由此选定小波基函数。 1.关于小波基函数选择的相关研究 通过查阅国内外相关文献,我了解到目前还是有一些学者在这方面做了一定的探索性工作。宋国乡等人提出了根据小波消失矩来选择小波基的思想。消失矩的定义为:小波Ψ(x)称为具有n阶消失矩,如果对于所有非负整数k,0≤h≤n, k??(x)dxx?0,选择方法是如果被检测信号的奇异度为α,n-1均有<α<n,R则需要具有n阶以上的消失矩的紧支撑的小波。其理论依据是假设小波具有紧支撑和n阶消失矩,且n次连续可微(其中n是正整数),设x为突变点,如果f(x)0在x的Lipschitz奇异度为α(α<n),而在x附近n次连续可微,则可以证明 00Wf(x)在x达到极大值。由此可以检测缓变信号的奇异点。0S周小勇等提出了采用小波规则性系数相似性来选择小波基的方法,其思想来源于Fourier变换。Fourier变换是以正弦信号为基波,用其各次谐波来近似某一函数或信号,其Fourier系数代表了各次谐波分量和原信号的相似性。小波系数的大小也反映了小波和函数某段的相似程度。同时函数和小波的规则性均表示着各自的可微和平滑程度。因此由相似性,可以用平滑的小波,即规则性系数大的小波表示平滑的函数;用不平滑的小波,即规则性系数小的小波表示非平滑的函数。当然这里所说的相似并不是绝对的相似或相近,而只是一种趋势。 2. 小波基函数的选择标准

二维各向同性不可分小波变换特性分析

文章编号:167320291(2006)0520024204 二维各向同性不可分小波变换特性分析 章春娥,裘正定 (北京交通大学计算机与信息技术学院,北京100044) 摘 要:从滤波器设计、采样方面与标准二维可分小波变换的比较、分析二维各向同性不可分小波变换在图像处理中表现出来的特性,并针对不可分小波提出了新的父子关系定义和改进型零树结构.二维不可分小波变换相对标准可分小波变换而言,尺度函数和小波函数不可分且各向同性,具有更细的渐进尺度,更好的紧支撑特性,各个子带有清晰的频率特征及重建特性.本文通过实验统计分析了二维不可分各向同性小波特性,利用改进型零树结构大大提升了其在图像压缩中的性能;对不可分小波变换在图像处理中的应用具有指导性的意义.关键词:各向同性;二维不可分小波;零树结构;图像处理中图分类号:TP391.4 文献标识码:A Analysis on 2-D Nonseparable and Isotropic W avelet T ransform ZHA N G Chun-e ,Q IU Zheng-di ng (School of Computer and Information Technology ,Beijing Jiaotong University ,Beijing 100044,China ) Abstract :C ompared with 2-D separable wavelet the properties of 2-D nonseparable wavelet are analyzed from the aspects of filter design ,sampling and others.For s pecial nonseparable wavelet ,a new father-s on relation 2ship and modified zero-tree structure are proposed.For 2-D nonseparable wavelet trans form ,its scale and wavelet functions are nonseparable and is otropic ,and it has more progressive scales and more tight energy sup 2port.Each subband has more s pecific frequency feature and could gain better reconstruction.This paper ana 2lyzes the properties by experimental data ,increasing greatly the image compression performance by the modi 2fied zerotree ,and als o providing guide for 2-D nonseparable wavelet in image processing.K ey w ords :isotropic ;2-D nonseparable wavelet ;zerotree structure ;image processing 小波变换凭借良好的时频局部化特性,其理论和方法在语音分析,模式识别,数据压缩,图像配准、数据融合、数字水印等信号处理方面得到了广泛的应用.一维小波的理论研究比较深入,多维小波通常自然而然地由一维小波的张量积得到,属于多维可分小波,即通常意义上的标准小波.多维标准小波基于成熟的一维小波理论,构造和实现基本上都基于一维小波进行.一维小波张量积形成的多维小波分析同时带来了一些问题[1].当用张量积表示多维尺度空间时,不同基函数的数目将增加为2d (其中d 表示空间维数).由于基函数是不同尺度下一维小波和尺度函数的积,一些多维小波将会高度偏离坐标方向,失去与一维小波一致的紧支撑特性.张量积形成的基函数还混淆了尺度的运算表达,降低了多维信号小波表示的稀疏性.换而言之,这类由一维小波张量积形成的多维标准小波引入了各向异质性,且难以从物理上解释相应尺度子带代表的频率信息.一维小波能很好地描述一维序列点的奇异性,却很难将对点序列奇异性的描述通过张量积的形式拓展到二维以上的信号处理. 收稿日期:2005212217 基金项目:国家中小型企业创新基金资助项目(04C26213301189);北京交通大学创新基金资助项目(2005SM009) 作者简介:章春娥(1976— ),女,四川达州人,博士生.em ail :z-ce @https://www.360docs.net/doc/b71882205.html, 裘正定(1944— ),男,浙江嵊县人,教授,博士生导师.第30卷第5期2006年10月 北 京 交 通 大 学 学 报JOURNAL OF BEI J IN G J IAO TON G UN IV ERSIT Y Vol.30No.5 Oct.2006

(完整版)matlab函数大全(非常实用)

信源函数 randerr 产生比特误差样本 randint 产生均匀分布的随机整数矩阵 randsrc 根据给定的数字表产生随机矩阵 wgn 产生高斯白噪声 信号分析函数 biterr 计算比特误差数和比特误差率 eyediagram 绘制眼图 scatterplot 绘制分布图 symerr 计算符号误差数和符号误差率 信源编码 compand mu律/A律压缩/扩张 dpcmdeco DPCM(差分脉冲编码调制)解码dpcmenco DPCM编码 dpcmopt 优化DPCM参数 lloyds Lloyd法则优化量化器参数 quantiz 给出量化后的级和输出值 误差控制编码 bchpoly 给出二进制BCH码的性能参数和产生多项式convenc 产生卷积码 cyclgen 产生循环码的奇偶校验阵和生成矩阵cyclpoly 产生循环码的生成多项式 decode 分组码解码器 encode 分组码编码器 gen2par 将奇偶校验阵和生成矩阵互相转换gfweight 计算线性分组码的最小距离 hammgen 产生汉明码的奇偶校验阵和生成矩阵rsdecof 对Reed-Solomon编码的ASCII文件解码rsencof 用Reed-Solomon码对ASCII文件编码rspoly 给出Reed-Solomon码的生成多项式syndtable 产生伴随解码表 vitdec 用Viterbi法则解卷积码 (误差控制编码的低级函数) bchdeco BCH解码器 bchenco BCH编码器 rsdeco Reed-Solomon解码器 rsdecode 用指数形式进行Reed-Solomon解码 rsenco Reed-Solomon编码器 rsencode 用指数形式进行Reed-Solomon编码 调制与解调

外文翻译小波变换在图像处理中的仿真及应用

论文翻译 通信102 吴志昊 译文: 小波变换在图像处理中的仿真及应用 一、课题意义 在传统的傅立叶分析中, 信号完全是在频域展开的, 不包含任何时频的信息, 这对于某些应用来说是很恰当的, 因为信号的频率的信息对其是非常重要的。但其丢弃的时域信息可能对某些应用同样非常重要, 所以人们对傅立叶分析进行了推广, 提出了很多能表征时域和频域信息的信号分析方法, 如短时傅立叶变换, Gabor 变换, 时频分析, 小波变换等。而小波分析则克服了短时傅立叶变换在单分辨率上的缺陷, 具有多分辨率分析的特点, 使其在图像处理中得到了广泛应用。 传统的信号理论,是建立在Fourier分析基础上的,而Fourier变换作为一种全局性的变化,其有一定的局限性。在实际应用中人们开始对Fourier变换进行各种改进,小波分析由此产生了。小波分析是一种新兴的数学分支,它是泛函数、Fourier分析、调和分析、数值分析的最完美的结晶;在应用领域,特别是在信号处理、图像处理、语音处理以及众多非线性科学领域,它被认为是继Fourier分析之后的又一有效的时频分析方法。小波变换与Fourier变换相比,是一个时间和频域的局域变换因而能有效地从信号中提取信息,通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析(Multiscale Analysis),解决了Fourier变换不能解决的许多困难问题。 小波变换是一种快速发展和比较流行的信号分析方法, 其在图像处理中有非常重要的应用, 包括图像压缩, 图像去噪, 图像融合, 图像分解, 图像增强等。小波分析是傅立叶分析思想方法的发展与延拓。除了连续小波(CWT)、离散小波(DWT), 还有小波包(Wavelet Packet)和多维小波。 小波分析在图像处理中有非常重要的应用, 包括图像压缩, 图像去噪, 图像融合, 图像分解, 图像增强等。小波变换是一种新的变换分析方法,它继承和发展了短时傅立叶变换局部化的思想,同时又克服了窗口大小不随频率变化等缺点,能够提供一个随频率改变的时间一频率窗口,是进行信号时频分析和处理的理想工具。它的主要特点是通过变换能够充分突出问题某些方面的特征,因此,小波变换在许多领域都得到了成功的应用,特别是小波变换的离散数字算法已被广泛用于许多问题的变换研究中。从此,小波变换越来越引进人们的重视,其应用领域来越来越广泛。 二、课题综述 (一)小波分析的应用与发展 小波分析的应用是与小波分析的理论研究紧密地结合在一起的。现在,它已经在科技信息产业领域取得了令人瞩目的成就。电子信息技术是六大高新技术中重要的一个领域,它的重要方面是图象和信号处理。现今,信号处理已经成为当代科学技术工作的重要部分,信号处理的目的就是:准确的分析、诊断、编码压缩和量化、快速传递或存储、精确地重构(或恢复)。从数学地角度来看,信号与图象处理可以统一看作是信号处理(图象可以看作是二维信号),在小波分析的许

文献综述组织发展

文献综述之组织发展 众所周知,企业中唯一不变的就是变革。组织变革已经广泛应用在企业的经营管理当中,用以谋求企业在竞争中的优势地位。1962 年美国学者钱德勒在其代表作《战略与结构:美国工商企业成长的若干篇章》中就提出了著名的“结构跟随战略”假设。文章首次分析了“环境——战略——组织结构”间的相互关系,提出企业的经营战略应该适应环境,满足市场需求,而组织结构又必须适应企业的战略,随着战略变化而变化,即“结构跟随战略”假设。 面对经济全球化,技术快速进步,信息的迅速传播,企业所承受的竞争压力和经营风险比过去的任何时候都更为严峻。 圣吉在1990年出版的《第五项修炼》中提出学习型组织的概念,他认为未来最成功的企业将会是学习型组织,因为未来唯一持久的优势,是有能力比你的竞争对手学习得更快。 1993 年,美国的钱皮和哈默提出“公司再造”的概念,组织或公司的再造被认为是继财务管理、战略计划管理和全面质量管理之后的又一场管理革命。 统计表明,90 年代初的几年中,85%的美国公司都投入到了各种各样的组织再造活动当中。这些公司当中,有60%报道说他们并没有获得所希望的生产力,并且有44%的公司声明他们事实上比以前更糟!在80%的公司当中,员工工作满意程度有所降低。面对这些变革方面的失败,以及新指导方针的缺乏,68%的组织机构在一年之内又不得不重新构建组织。从中我们可以看出,组织变革本身是一项复杂的系统工程,其间充满着风险和不确定性,失败的比例非常高,但由于它是组织迎接时代挑战,获取竞争优势的必由之路,所以通过研究现有的理论,分析实践中的成败,并把两者有机地结合起来,探讨组织变革的内在逻辑关系,用来指导组织变革的实践就显得是一件非常有意义和必需的工作。 本文研究的主要问题可以归结为:如何重建组织战略、组织结构和组织文化,达到提高组织运作效率,实现对环境变化的适应,打造组织核心竞争能力,赢得市场竞争的目的。 通过这些研究主要解决现实组织变革中的两大问题: 1、研究如何通过组织变革使组织增强灵活性与敏捷性,从而有效增强对复杂环境的适应能力与应变能力; 2、探讨如何通过组织变革实现组织内部运作机制与组织成员的心理、行为方式的协调,使组织增强新的活力,体现出团队优势。 由此可见,组织发展可以从以下方面进行讨论:

小波的几个术语及常见的小波基介绍

小波的几个术语及常见的小波基介绍 本篇是这段时间学习小波变换的一个收尾,了解一下常见的小波函数,混个脸熟,知道一下常见的几个术语,有个印象即可,这里就当是先作一个备忘录,以后若有需要再深入研究。 一、小波基选择标准 小波变换不同于傅里叶变换,根据小波母函数的不同,小波变换的结果也不尽相同。现实中到底选择使用哪一种小波的标准一般有以下几点: 1、支撑长度 小波函数Ψ(t)、Ψ(ω)、尺度函数φ(t)和φ(ω)的支撑区间,是当时间或频率趋向于无穷大时,Ψ(t)、Ψ(ω)、φ(t)和φ(ω)从一个有限值收敛到0的长度。支撑长度越长,一般需要耗费更多的计算时间,且产生更多高幅值的小波系数。大部分应用选择支撑长度为5~9之间的小波,因为支撑长度太长会产生边界问题,支撑长度太短消失矩太低,不利于信号能量的集中。 这里常常见到“紧支撑”的概念,通俗来讲,对于函数f(x),如果自变量x在0附近的取值范围内,f(x)能取到值;而在此之外,f(x)取值为0,那么这个函数f(x)就是紧支撑函数,而这个0附近的取值范围就叫做紧支撑集。总结为一句话就是“除在一个很小的区域外,函数为零,即函数有速降性”。 2、对称性 具有对称性的小波,在图像处理中可以很有效地避免相位畸变,因为该小波对应的滤波器具有线性相位的特点。 3、消失矩 在实际中,对基本小波往往不仅要求满足容许条件,对还要施加所谓的消失矩(Vanishing Moments)条件,使尽量多的小波系数为零或者产生尽量少的非零小波系数,这样有利于数据压缩和消除噪声。消失矩越大,就使更多的小波系数为零。但在一般情况下,消失矩越高,支撑长度也越长。所以在支撑长度和消失矩上,我们必须要折衷处理。

浅析中国现当代文学的研究现状与发展

【论文关键词】现当代文学思想情感语言 【论文摘要】中国现当代文学是用现代的文学语言和形式表达现代中国思想情感的文学。本文研究了中国现当代文学的现状,并对现当代文学的发展提出了几点展望。 中国现当代文学通过现代社会和人的意识情感的加入,来改变中国古典文学造成的封闭和隔绝,使文学在内容和表达上与当代中国人的实际有更多的联系和契合。近年来,现当代文学研究过程中,在一种无孔不入的话语的渗透之下,一些像“全球化”之类的词语,成了神圣的词语,只要说出来就具有天然合法性,而不论其所指何物。与此同时,一些诸如“革命”、“救亡”之类的词语,却获得了截然不同的命运,一说出来就成了天然的“非法”,甚至成了过街老鼠,人人喊打。在当下的语境中,“纯文学”、“新时期文学”、“先锋文学”等成了人见人爱的时髦话语,而“十七年文学”、“革命历史小说”、“底层写作”等则成了人见人怕的危险思想。在中国现当代文学的研究过程中,现当代文学究竟如何发展变化,又处在何种现状呢? 一、现当代文学研究目的发生了本质的变化 在中国现当代文学史上,上世纪八十年代成为一个重要的时期。在以和平经济建设为中心时期,文学队伍主要是由两代作家构成,一代是在上世纪五十年代成长起来的作家,他们大部分是在1957年的“反右”运动中遭到不公的批判和打击,并在社会底层度过了一段苦难日子,他们的创作里充满对现实政治生活的干预精神和对人性的赞美。另一代是在“文革”中成长起来的作家,其中大多数人曾在“上山下乡”中感受了民间生活和民间文化的熏陶,因此他们在写作的时候会自然而然地从农村经验中汲取创作素材,由最初的知青题材到稍后的寻根文学,呈现出新的民间化创作趋向。两代作家们胸中涌动着知识分子新生的对现实生活的热情与自信,他们在揭露社会弊病的同时,把希望寄托于批判的社会效果,在这种希望之中正滋生着已经消失了近三十年的知识分子的主体意识。这使上世纪八十年代的文学充满了生机勃勃的创新精神和活跃气氛。 而今天的当代中国作家,一部分忙着追逐“诺贝尔”文学奖,期望以此来奠定自己在文学史上的地位;另一部分则忙于追逐市场。由于读者的欣赏情趣和阅读趋向主导了市场,市场又左右着作家的文学创作,导致现在的文学市场化、快餐化、娱乐化、庸俗化、视觉化倾向。在现当代文学中,以两种文学思想表现最为突出。 (1)精神溃败。在这种精神溃败之中,有的作者干脆放弃信仰,放弃寻找精神出路的企图。如余华就在《活者》中传达了“活着就是一切,活着就是胜利”这样一种人生哲学。正如余杰所说,只有肉体的“活着”而没有灵魂的求索。 (2)性文学进入文坛。90年代末期,性文学卷土重来进入了中国文坛,在这些描写性的文学中声势最响的是贾平凹的《废都》和王小波的《黄金时代》。它们都描写了男人在生存困境中表现出的依赖性和脆弱性。 二、现当代文学存在“现状批评”现象 现状批评是现当代文学研究中最具诱惑力的部分,变化无穷的文学现象和不断推出的新人和作品,既联系着社会生活、文化思潮的搏动,又提供了不重复的新鲜刺激,批评的创造性在现当代文学研究中也最为活跃。因此,许多人把当代文学研究与“现状批评”画上等号。这一观点是值得怀疑的,它可能迟滞了对“十七年文学”、“‘文革’文学”的研究,以至于这三十年的文学已从一些人的研究视野中消失。这些部分文学研究,往往构成现当代文学研究最脆弱的环节。很显然,这三十年的文学创作总体水准的普遍性缺陷,是难以吸引更多注意力的重要原因。但是,从另一角度看,这一曲折的文学进程,又蕴含着某些有关中国现代文学的重要问题和矛盾,这涉及现代文学传统、作家精神结构、现代文学艺术形态等。

小波变换相关实验

小波变换相关实验 实验目的 1、通过观察小波变换系数建立对小波变换及其有关性质的感性认识。 2、掌握小波变换及重构方法;了解小波变换基本应用。 实验内容 1、图像二维离散小波变换及其重构; 2、小波变换在去噪、压缩、图像增强上的应用。 实验原理 1、“小波”就是小区域、长度有限、均值为0的波形。所谓“小”是指它具有衰减性;而称之为“波”则是指它的波动性,其振幅正负相间的震荡形式。与 Fourier变换相比,小波变换是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了Fourier变换的困难问题,成为继Fourier变换以来在科学方法上的重大突破。 小波转换分成两个大类:离散小波变换 (DWT) 和连续小波转换 (CWT)。两者的主要区别在于,连续转换在所有可能的缩放和平移上操作,而离散转换采用所有缩放和平移值的特定子集。 小波变换的公式有内积形式和卷积形式,两种形式的实质都是一样的。它要求的就是一个个小波分量的系数也就是“权”。其直观意义就是首先用一个时窗最窄,频窗最宽的小波作为尺子去一步步地“量”信号,也就是去比较信号与小波的相似程度。信号局部与小波越相似,则小波变换的值越大,否则越小。当一步比较完成后,再将尺子拉长一倍,又去一步步地比较,从而得出一组组数据。如此这般循环,最后得出的就是信号的小波分解(小波级数)。 当尺度及位移均作连续变化时,可以理解必将产生大量数据,作实际应用时并不需要这么多的数据,因此就产生了离散的思想。将尺度作二进离散就得到二进小波变换,同时也将信号的频带作了二进离散。当觉得二进离散数据量仍显大时,同时将位移也作离散就得到了离散小波变换。 2、二维离散小波变换常用函数

Matlab作图函数的总结与分析

高等理科教育2005年第6期(总第64期)Matlab作图函数的总结与分析+ 黄琼湘那斯尔江?吐尔逊 (tfi疆大学数学与系统科学学院,新疆乌鲁木齐830046) 摘要Matlab(MATrixLABoratory的简称)是CleverMoler博士用Fortran语言开发的科学计算工具。它已成为科学研究、工程计算、应用开发的重要工具。国外已将它作为理工科大学的必修课程,国内各大学也开始开设这门课程。Matlab有强大的作图功能,有兴趣的读者可参考文献【卜4’。本文对Matlab的作图函数进行分析和总结,以供教学参考和学生学习之用。 关键词Matlab数据可视化作图函数 中图分类号G642.0文献标识码A 一、Matl如作图函数的总结 Matlab提供了丰富的作图函数,有100个之多。在教学和学习中显得有点杂乱。我们先对它们进行总结和分类,并提炼出它们的共性和特性。 Matlab的作图函数从视角的维数上分有三类:一维作图函数、二维作图函数和三维作图函数。它们的代表分别是line、plot和plot3等函数。从类型上分大致有四类:通用作图函数(如plot函数等);专业作图函数(如contour函数、quiver函数等);动画制作函数(如movie、comet3等函数);图形修饰函数(如view等函数)。 Matlab所有的作图函数都可以通过查帮助获得它的功能和用法。这里我们把作图函数按类型分类,列出一些主要和常用的作图函数(见表1),以抓住重点。 作图函数虽然功能不同,但它们的调用格式是一致的。我们用GraphF来表示一般的作图函数,它们的调用格式如下: 1.GraphF(X,Y,S) 这是一、二维函数的作图格式。x和Y表示图形的数据点,s表示图形修饰参数组(可以缺省)。当x,Y都是顶点坐标时,GraphF(X,Y,S)画出以x,Y为端点,s为参数的线;当x是一组顶点坐标,而Y对应于X的函数值时,GraphF(X,Y,S)画出函数Y的二维图形。 2.GraphF(X,Y,Z,S) 这是三维函数的作图格式。z是x和Y的函数。x,Y以二维坐标形式表示函数值z的作图区域D,s表示图形修饰参数组(可以缺省)。GraphF(x,Y,Z,S)画出定义域为D的函数z的三维图形。 值得注意的是,Matlab的作图函数总是描绘数据点(X,Y)(在平面上)或(x,Y,z)(在空间中)的图形。前者视为Y的函数,而后者视为z的函数。函数GraphF在作图前数据点必须事先给定,在作图时函数GraphF将各数据点用光滑的曲线连接成图形。另外,X,Y,Z还 ÷收稿日期2004—02—19 资助项目新疆大学校基金“应用软件程序设计”重点课程建设项目资助 作者简介黄琼湘(1958)男,湖南衡阳人,教授,主要从事组合数学与图论、计算机算法研究

小波相关性与相干性

前言 时域指标参数 1. 均值 当观测时间T 趋于无穷时,信号在观测时间T 内取值的时间平均值就是信号()x t 的均值。均值定义为 ()dt t x T T T x ? ∞→=01lim μ (1) 式中:T 是信号的观测区间。实际T 不可能为无穷,算出的x μ必然包含统计误差,只能作为真值的一种估计。 2. 均方值和方差 当观测时间T 趋于无穷时,信号在观测时间T 内取值平方的时间平均值就是信号()x t 的均方值,定义为: ()dt t x T T T x ?∞→=0221lim φ (2) 如果仅对有限长的信号进行计算,则结果仅是对其均方值的估计。均方值的正平方根,为均方根值(或有效值)m ax x 。 方差定义为 ()[]dt t x T T x T x ?-=∞→0221lim μσ (3) 方差反应了信号()x t 中的动态部分。方差的正平方根x σ称为标准差。若信号()x t 的均值为零,则均方值等于方差。若信号()x t 的均值不为零时,则有下列成立 222x x x μφσ-= (4) 3. 概率密度函数 随机信号()x t 的取值落在区间内的概率可用下式表示 ()[]T T x x t x x P T prb ?=?+≤<=∞→lim (5) 式中:T ?为信号()x t 取值落在区间(]x x x ?+,内的总时间;T 为总观察时间。 当0→?x 时,概率密度函数定义为 ()?? ??????=∞→∞→?T T x x p T x lim 1lim (6) 随机信号()x t 的取值小于或等于某一定值δ的概率,称为信号的概率分布函数。常用()x P 来

小波奇异值检测

基于连续小波变换的奇异性检测与故障诊断 林京 振动工程学报2000 基于小波变换的奇异性检测方法可以实现对信号局部奇异性的刻画。因此,自该方法提出以来,便获得了广泛的应用。仅在机械故障诊断领域中,它就被用来进行超声无损探伤、柴油机的压力波形识别、纲丝绳断丝检测、切削颤振分析等。从这些应用中可以看到,它们都无一例外地采用二进离散小波变换来做奇异性检测。尽管此时的二进离散只对尺度区间进行,各尺度上信号的时间间隔等于采样间隔,这种小波分解也足以大大缩减了计算量,因此许多研究人员乐于采用。但是,尺度上的二进分割会使奇异性定量过于粗糙,尤其是在低尺度区间的信号,信号中的奇异点往往无法考察,从而出现漏检或定量不准。如果采用连续小波变换,则这些缺陷可迎刃而解。 Lip指数表明了函数f(x)与n次多项式作比较时,其光滑程度是多少。传统的计算Lip指数的计算方法是采用Fourier变换,它只能得到信号整体的Lip指数,所以只能反映信号的全局奇异程度。如果要求得到信号在某点的奇异性,需要借助小波方法来实现。 综合起来,选择小波函数的原则就是在满足能够检测到最大Lip指数值的前提下,选择具有最少消失矩的小波函数。目前比较多的采用墨西哥草帽小波函数和Morlet小波函数。 实际应用中,要根据具体情况选择不同阶的导函数作小波函数。在机械测试信号中,奇异点通常为一些峰值点和突变尖点,这些点的Lip指数总小于1,所以小波函数只具有一阶消失矩即可。 模极大值线的存在要求尺度必需连续变化,当尺度作二进离散后,不存在模极大值线。 尺度区间的二进离散在许多情况下显得过于粗糙,它对Lip指数的定量描述不够精确,采用连续小波变换可以准确对Lip指数作定量计算。 机械状态监测中,信号的突变点往往携带着故障信息,机器运行过程中所产生的撞击、振荡、摩擦、转速突变、结构变形和断裂等都可反映在信号的突变点中,信号突变点的奇异性检测可以有效地揭示机器的故障信息,为机器故障诊断提供有力工具。 用小波做奇异性检测可以对信号的局部奇异程度做定量描述,但目前通用的用二进小波做奇异性检测的定量描述在许多应用场合下仍显得粗糙,不能很好地区分正常点还是异常点。基于连续小波变换的奇异性检测方法可以对局部奇异性做精确的定量描述,具有广泛的应用前景。

matlab函数名称总结

一、常用对象操作:除了一般windows窗口的常用功能键外。 1、!dir 可以查看当前工作目录的文件。!dir& 可以在dos状态下查看。 2、who 可以查看当前工作空间变量名, whos 可以查看变量名细节。 3、功能键: 功能键快捷键说明 方向上键Ctrl+P 返回前一行输入 方向下键Ctrl+N 返回下一行输入 方向左键Ctrl+B 光标向后移一个字符 方向右键Ctrl+F 光标向前移一个字符 Ctrl+方向右键 Ctrl+R 光标向右移一个字符 Ctrl+方向左键 Ctrl+L 光标向左移一个字符 home Ctrl+A 光标移到行首 End Ctrl+E 光标移到行尾 Esc Ctrl+U 清除一行 Del Ctrl+D 清除光标所在的字符 Backspace Ctrl+H 删除光标前一个字符 Ctrl+K 删除到行尾 Ctrl+C 中断正在执行的命令 4、clc可以命令窗口显示的内容,但并不清除工作空间。 二、函数及运算 1、运算符: +:加,-:减, *:乘, /:除,\:左除 ^:幂,':复数的共轭转置,():制定运算顺序。 2、常用函数表: sin( ) 正弦(变量为弧度) Cot( ) 余切(变量为弧度) sind( ) 正弦(变量为度数) Cotd( ) 余切(变量为度数) asin( ) 反正弦(返回弧度) acot( ) 反余切(返回弧度) Asind( ) 反正弦(返回度数) acotd( ) 反余切(返回度数) cos( ) 余弦(变量为弧度) exp( ) 指数 cosd( ) 余弦(变量为度数) log( ) 对数 acos( ) 余正弦(返回弧度) log10( ) 以10为底对数 acosd( ) 余正弦(返回度数) sqrt( ) 开方 tan( ) 正切(变量为弧度) realsqrt( ) 返回非负根 tand( ) 正切(变量为度数) abs( ) 取绝对值 atan( ) 反正切(返回弧度) angle( ) 返回复数的相位角

女大学生就业中性别歧视分析文献综述

女大学生就业中性别歧视分析文献综述 摘要:随着我国高等教育改革的发展,各个高校的招生规模也在逐年的增加。这也就导致了毕业生的数量也在逐年增加,竞争压力也慢慢增大。与此同时它将我国大学毕业生就业歧视问题推向了一个新的局面。本文就搜集众多学者观点及部分文献信息对我国大学生就业歧视问题进行简要分析综述,并提出一些建议及对策。 关键词:大学毕业生;就业歧视;区别对待; 大学生就业环境日趋激烈,就业中的性别歧视问题就成为了一个社会问题。有调查显示,目前我国高等学校中女生占的比例已达44 %,在相同条件下女生就业机会只有男生87% ,女毕业生初次就业率仅为63.4% ,比男生低8.7个百分点。这一现象已引起了社会的广泛关注。本文搜集相关材料,从大学毕业生就业歧视定义、国内外相关研究、就业歧视成因和就业歧视对策三个方面进行分析讨论,力求能得出一些有实质意义的合理化建议和对策。 一、女大学生就业中性别歧视的界定 分析女大学生就业现象,首先应该明确一下概念。 (一)性别歧视 国际劳工组织1958 年通过的《就业与职业歧视公约》(第ll1 号公约)对“就业歧视”的界定是:基于种族、肤色、性别、宗教、政治见解、民族、血统或社会出身的任何区别、排斥或特惠,其效果为取消或损害就业或职业方面的机会平等或待遇平等。 湖南师范大学教育科学学院覃伟丽认为研究我国大学生就业歧视,首先必须将大学生就业歧视与一般的就业歧视进行区分。形式逻辑学“种差加属”的定义方法能够对大学生就业歧视进行较为准确的定义区分。根据其分析得出“大学毕业生就业歧视”是指,基于种族、肤色、性别、宗教、政治见解、民族、血统或社会出身等原因,雇主方对于获得大学文凭就业者的初次或初期就业,采取任何区别、排斥或特惠,其效果为取消或损害就业或职业方面的机会平等或待遇平等。[3] (二)女大学生就业中的性别歧视

小波分析-经典解读

时间序列-小波分析 时间序列(Time Series )是地学研究中经常遇到的问题。在时间序列研究中,时域和频域是常用的两种基本形式。其中,时域分析具有时间定位能力,但无法得到关于时间序列变化的更多信息;频域分析(如Fourier 变换)虽具有准确的频率定位功能,但仅适合平稳时间序列分析。然而,地学中许多现象(如河川径流、地震波、暴雨、洪水等)随时间的变化往往受到多种因素的综合影响,大都属于非平稳序列,它们不但具有趋势性、周期性等特征,还存在随机性、突变性以及“多时间尺度”结构,具有多层次演变规律。对于这类非平稳时间序列的研究,通常需要某一频段对应的时间信息,或某一时段的频域信息。显然,时域分析和频域分析对此均无能为力。 20世纪80年代初,由Morlet 提出的一种具有时-频多分辨功能的小波分析(Wavelet Analysis )为更好的研究时间序列问题提供了可能,它能清晰的揭示出隐藏在时间序列中的多种变化周期,充分反映系统在不同时间尺度中的变化趋势,并能对系统未来发展趋势进行定性估计。 目前,小波分析理论已在信号处理、图像压缩、模式识别、数值分析和大气科学等众多的非线性科学领域内得到了广泛的应。在时间序列研究中,小波分析主要用于时间序列的消噪和滤波,信息量系数和分形维数的计算,突变点的监测和周期成分的识别以及多时间尺度的分析等。 一、小波分析基本原理 1. 小波函数 小波分析的基本思想是用一簇小波函数系来表示或逼近某一信号或函数。因此,小波函数是小波分析的关键,它是指具有震荡性、能够迅速衰减到零的一类函数,即小波函数)R (L )t (2∈ψ且满足: ? +∞ ∞ -=0dt )t (ψ (1) 式中,)t (ψ为基小波函数,它可通过尺度的伸缩和时间轴上的平移构成一簇函数系: )a b t ( a )t (2 /1b ,a -=-ψψ 其中,0a R,b a,≠∈ (2) 式中,)t (b ,a ψ为子小波;a 为尺度因子,反映小波的周期长度;b 为平移因子,反应时间上的平移。 需要说明的是,选择合适的基小波函数是进行小波分析的前提。在实际应用研究中,应针对具体情况选择所需的基小波函数;同一信号或时间序列,若选择不同的基小波函数,所得的结果往往会有所差异,有时甚至差异很大。目前,主要是通过对比不同小波分析处理信号时所得的结果与理论结果的误差来判定基小波函数的好坏,并由此选定该类研究所需的基小波函数。 2. 小波变换 若)t (b ,a ψ是由(2)式给出的子小波,对于给定的能量有限信号)R (L )t (f 2 ∈,其连续小波变换(Continue Wavelet Transform ,简写为CWT )为: dt )a b t ( f (t)a )b ,a (W R 2 /1-f ? -=ψ (3) 式中,)b ,a (W f 为小波变换系数;f(t)为一个信号或平方可积函数;a 为伸缩尺度;b 平移参数; )a b x ( -ψ为)a b x (-ψ的复共轭函数。地学中观测到的时间序列数据大多是离散的,设函数)t k (f ?,

相关文档
最新文档