一张灰度图是由多个像素点而组成的,同样,这些像素点的是由一个从0(黑)到255(白)的非负数组成的。假设我们现在有一张小的灰度图像。在第一行的灰度值为110,100,120,140,130,100,100.这些灰度图在图3.1画出

图3.1 灰度

我们很很自然的就会问,能不能用一个函数f(t)来表示这一行的数据?我们可以用从例1.2推出的函数Φ(t)来表示。Φ(t)为

那么我们就开始用Φ(t)来构建f(t)吧

Φ(t)以及f(t)如图3.2所示

注意到,所有的f(t)都是对所有的t∈R是连续的,除了在t=1,2,3,4,5。

  • V0空间的定义

我们可以使用上面的示例创建一个向量空间,并且所有的的断点都在整数。我们不妨试着定义这个空间是由Φ(t)以及它的变换函数Φ(t-k)张成的,其中k∈Z,所以在这个空间的元素可以表示为

在式子3.3中我们故意模糊了求和的范围。我们是不是真的要对求和进行一些约束呢?其实,在图像处理中,k只要覆盖到了所有的行列就可以了。在这种情况下,我们称这些像素为紧密支撑元素。

但是,如果把k设定为有限值,我们就无法对一些有用的函数建模。好比如sin和cos线性组合的这种无限长支撑集的函数。当然了,所有的实际运用的函数都是有限的。所以,我们必须要在L2(R)空间里面对函数进行讨论,并且不需要讨论到无穷大,只是需要讨论到函数取值近似于0的时候就可以了,所以我们有如下定义

定义3.1(哈尔空间V0和哈尔函数Φ(t))Φ(t)已经在(3.1)写过了,我们下面就来定义V0空间

(作者注:这么命名是为了纪念匈牙利的数学家n Alfred Haar (1885-1933),这个数学家创立了正交函数论)

这个空间就是在L2(R)里面断点为正数的所有分段常数函数。在练习3.4你需要证明这个空间是L2(R)的子空间。

通过定义,你就可以知道Φ(t)以及Φ(t-k)张成了整个空间。在练习3.1中你就会知道,他们是线性无关的。因此,他们是这个空间的一组基

  • Φ(t)以及它的变换的正交性

我们来选择下标j!=k的函数Φ(t-j)以及Φ(t-k)。Φ(t-j)在区间[j,j+1)不为0,同样Φ(t-k)在[k,k+1)不为0.如果j!=k,他们的乘积就为0,如图3.3所示

因此内积为

所以他们是正交的。同样,我们注意到了Φ(t-k)2=Φ(t-k),所以无论我们怎么平移,它与横轴所围成的面积始终为1

注意我们本章中所讨论都是实函数。所以我们修改了定义1.8的式子中的共轭来计算内积

主题3.1(V0的正交基)在定义3.1中给出的V0,它的正交基为{Φ(t-k)}k∈Z

证明:通过定义以及练习3.1你会知道{Φ(t-k)}k∈Z是线性无关的,之前的计算也很明V0 的内积为

其中δj,k在式子(2.16)给出

下面我们看一些V0的例子

例3.1(V0中元素)判断下列函数是否位于V0

(a)其中

(b)阶梯函数

(c)函数

(d)函数

解:

对于(a),我们作如下化简

不难看出,间断点为2个,而且分段的地方是正数并且

所以f(t)∈V0

对于(b),当t≥1且k=1,2,3,4....的时候,g(t)≥1。然而,我们在练习1.13(b)知道,函数

并不属于L2(R)。又因为在[0,∞)上g(t)≥r(t)。所以。虽然这个函数的间断点也在正数,因为不是绝对可积的,所以并不属于空间V0

对于(c),我们可以看出这个是函数Φ(t-k)的线性结合。但是我们还是要费点功夫来看。所以我们先化简一下

这些严格的证明应该设计到无穷级数求和以及不定积分的顺序的交换。这些应该在分析数学的课上学过,感兴趣的读者可以看Rudin的书。h(t)2的化简如下

我们在R上对它进行积分

从(3.5)我们知道Φ(t-k)和Φ(t-j)的积分为0,只有当j=k的时候为1,所以,式子可以化简为

从积分学,我们知道是p阶收敛的。(作者注:我们可以尝试通过傅里叶级数证明这个东西的求和为∏2/6,详见Kammler的书)。所以我们知道h(t)属于V0空间。

(d)问题症结所在是每一项的系数都涉及了Φ(4t)这个函数。这是一个经过缩放的哈尔函数

我们在图3.4中画出了Φ(4t)。现在Φ(4t)=Φ(4(t-1/4)),所以我们可以认为这个是向右移动了1/4个单位,其余的以此类推,我们在图3.4画出了这个函数

图3.4函数Φ(4t)以及l(t)

  • 从L2(R)映射到V0

由于{Φ(t-k)}k∈Z是一组正交基,因此,我们可以将一个位于L2(R)的函数映射到V0上面,有

我们下面来看几个L2(R)投影到V0的例子

例3.2(投影到V0我们从问题。练习1.13(a)知道函数,我们来寻找他在V0上的投影

解:

对于k∈Z,我们必须计算。对于k≥0,在[k,k+1)上有Φ(t-k)=1,其余区间上Φ(t-k)=0所以我们在[k,k+1)上积分就可以了。又由于k>0,因此,e-∣t=e-t,从而内积为

同理k<0也有类似的做法,不再赘述,直接给出式子

我们在图3.5画出P(t)逼近的效果图

图3.5 函数g(t)以及它在V0的投影

3.1 哈尔空间 V0的更多相关文章

  1. 3.2 一般的哈尔空间Vj

    例3.2给予我们继续往下面做的动力.很明显的我们对于g(t)的逼近还是太粗糙了.很自然的,我们会想到,如果继续细分我们的短点,比如每1/2取一个值,甚至每1/4取一个值,那么就会有更好的逼近效果. 不 ...

  2. 3.3 哈尔小波空间W0

    在3.2节我们学习了关于(3.8)定义的Vj的性质.特别的,我们可以乘以系数从一个Vj空间变换到另一个.我们这节学习V0和V1的关系. 将f1(t)∈V1投影至V0 我们考虑一个属于V1的函数f1(t ...

  3. [算法]检测空间三角形相交算法(Devillers & Guigue算法)

    #pragma once //GYDevillersTriangle.h /* 快速检测空间三角形相交算法的代码实现(Devillers & Guigue算法) 博客原地址:http://bl ...

  4. [PaPaPa][需求说明书][V0.3]

    PaPaPa软件需求说明书V0.3 前   言 不好意思,本文是没有前言的. 别说是前言了,其实关于界面的内容我也是不打算写!!   因为我知道你们想要的界面是这样的: 再不济也应该是这样的: 但是我 ...

  5. VMware 虚拟机Red Hat 5.9 交换区及硬盘空间调整

    首先要通过VMware设置简单实现内存扩大.但是系统中的/swap应该如何设置呢? 1. 创建swap 文件 使用如下命令: #dd if=/dev/zero of=/swap/swapfile bs ...

  6. DBA_Oracle Sort排序处理空间耗用(概念)

    2014-12-18 Created By BaoXinjian

  7. 转:内核空间与用户空间数据交换的方式之一 --ioctl(通过字符设备演示)

    对于linux而言,内核程序和用户程序分别运行在内核空间和用户空间,要实现两者的数据交换,主要有以下几种方式:系统调用,读写系统文件(procfs,sysfs, seq_file,debugfs等), ...

  8. 3D空间中射线与三角形的交叉检測算法

    引言 射线Ray,在3D图形学中有非常多重要的应用.比方,pick操作就是使用射线Ray来实现的,还有诸如子弹射线的碰撞检測等等都能够使用射线Ray来完毕. 所以,在本次博客中,将会简单的像大家介绍下 ...

  9. Python巡检Oracle表空间并邮件告警

    最近,自学了Python基础,突发奇想,把以前通过shell自定义通过nagios实现Oracle表空间以及ASM以及备份的脚本改进下,首先感叹的是Python脚本看上去确实挺好的,效率还不错. 这是 ...

随机推荐

  1. CSS教程:div垂直居中的N种方法[转]

    在说到这个问题的时候,也许有人会问CSS中不是有vertical-align属性来设置垂直居中的吗?即使是某些浏览器不支持我只需做少许的CSS Hack技术就可以啊!所以在这里我还要啰嗦两句,CSS中 ...

  2. linux 下 C语言显示中文

    例如:tset.c int main() { printf("你好,世界\n"); retuen 0; } 编译时应该这样: iconv -f gb2312 -t utf8 tes ...

  3. 让dedecms autoindex,itemindex 从0到1开始的办法! 多重样式输出(借鉴)

    原网址:http://my.oschina.net/lyx2012/blog/55888 [field:global name=autoindex runphp="yes"]@me ...

  4. HDU 5446 中国剩余定理+lucas

    Unknown Treasure Time Limit: 1500/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Other ...

  5. js switch 扩展

    //demo var num=99 switch(n){ case 80<n: document.write("优秀");break; case 70<n: docum ...

  6. android WeakReference(弱引用 防止内存泄漏)与SoftReference(软引用 实现缓存机制(cache))

    在Android开发中,基本上很少有用到软引用或弱引用,这两个东东若用的很好,对自己开发的代码质量的提高有很大的帮助.若用的不好,会坑了自己.所以,在还没有真正的去了解它们之前,还是慎用比较好. 下面 ...

  7. pyqt4:连接的一个带有参数的方法

    一般在pyqt4中的信号连接很少连接带参数的方法,很多时候连接带参数的方法节约不少代码量. self.s5_thread=scene.Worker5() self.log_get=QtCore.QTi ...

  8. kali linux 安装nvidia

    开始安装之前需要说明一下几点: 1.安装闭源显卡驱动有一定风险,比如黑屏或者无法进入图形界面什么的.如果您很害怕折腾,而且不在Linux系统中玩开源游戏(比如Nexuiz)或看高清电影,默认的nouv ...

  9. Next Permutation

    Implement next permutation, which rearranges numbers into the lexicographically next greater permuta ...

  10. Docker搭建Java Web运行环境

    1. 前提条件 安装了Docker的64位Linux 操作系统 Linux操作系统镜像 Linux版本的JDK压缩包 Linux版本的Tomcat压缩包 2. 启动容器 容器是在镜像的基础上来运行的, ...