一张灰度图是由多个像素点而组成的,同样,这些像素点的是由一个从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. DELL vostro V5460 装mSATA 半高SSD硬盘后装win8.1系统

    1.介绍 DELL vostro V5460的超极本性能虽然不是那么好,但的确是我觉得样子最好看的一款.满足了,轻,薄.如果是红色的,则很适合女生用. 2.缺点 这款笔记本最大的缺点,就是有的无线网卡 ...

  2. 推荐一些不错的计算机书籍(php c mysql linux等等)

    推荐一些不错的计算机书籍. # PHP<PHP程序设计>(第2版)  --PHP语法和入门最好的书<PHP5权威编程>  --PHP入门后升级书<深入PHP:面向对象.模 ...

  3. JS中直接从java后台获得对象的值(数组的值)

    这里举得例子是:JS直接从后台Contorller中(SpringMVC中的model中)获得数值的值 Contorller 此处将 talentIntegralRecordsDay talentIn ...

  4. java 读写properties

    网速不好:先贴上资料: Java配置文件Properties的读取.写入与更新操作 [Spring] - Property注入 http://www.360doc.com/content/14/073 ...

  5. RNN 入门学习资料整理

    建议按序阅读 1. RNN的一些简单概念介绍 A guide to recurrent neural networks and backpropagation Deep learning:四十九(RN ...

  6. 2016年中国大学生程序设计竞赛(合肥)-重现赛1008 HDU 5968

    异或密码 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submis ...

  7. 20160405互联网新闻<来自涛涛大产品>

    1.滴滴或将收购腾讯地图,打造“滴滴地图”(滴滴与神州.uber之间的争斗,归根到底还是BAT的代理人之战)2.优信二手车否认合并传言 并谴责58同城仿冒优信网站(商战无所不用其极)3.京东旗下的拍拍 ...

  8. mySql中IFNULL的使用说明

    IFNULL(expr1,expr2) 如果expr1不是NULL,IFNULL()返回expr1,否则它返回expr2.IFNULL()返回一个数字或字符串值 具体用法如:现有学生表(tbl_stu ...

  9. Angularjs 服务注册

    $injector: (When you request a service, the $injector is responsible for finding the correct service ...

  10. Android驱动开发前的准备(三)

    Git使用入门 3.1安装Git 3.2查看Git文档 3.3源代码的提交与获取 3.1安装Git # apt-get install git # apt-get install git-doc gi ...