https://www.zhihu.com/question/21686447?nr=1
对于初学者,我推荐用复利的例子来理解卷积可能更直观一些:
小明存入100元钱,年利率是5%,按复利计算(即将每一年所获利息加入本金,以计算下一年的利息),那么在五年之后他能拿到的钱数是,如下表所示:用求和符号来简化这个公式,可以得到:
在上式中,为小明的存钱函数,而为存入银行的每一笔钱的复利计算函数。在这里,小明最终得到的钱就是他的存钱函数和复利计算函数的卷积。
为了更清晰地看到这一点,我们将这个公式推广到连续的情况,也就是说,小明在从到的这一段时间内,每时每刻都往银行里存钱,他的存钱函数为,而银行也对他存入的每一笔钱按复利公式计算收益:,则小明到时间将得到的总钱数为:
这也就是卷积的表达式了,上式可以记为。
相信通过上面这个例子,大家应该能够很清晰地记住卷积公式了。下面我们再展开说两句:
如果我们将小明的存款函数视为一个信号发生(也就是激励)的过程,而将复利函数视为一个系统对信号的响应函数(也就是响应),那么二者的卷积就可以看做是在时刻对系统进行观察,得到的观察结果(也就是输出)将是过去产生的所有信号经过系统的「处理/响应」后得到的结果的叠加,这也就是卷积的物理意义了。一个关于卷积的血腥的例子:
比如说你的老板命令你干活,你却到楼下打台球去了,后来被老板发现,他非常气愤,扇了你一巴掌(注意,这就是输入信号,脉冲),于是你的脸上会渐渐地(贱贱地)鼓起来一个包,你的脸就是一个系统,而鼓起来的包就是你的脸对巴掌的响应,好,这样就和信号系统建立起来意义对应的联系。下面还需要一些假设来保证论证的严谨:假定你的脸是线性时不变系统,也就是说,无论什么时候老板打你一巴掌,打在你脸的同一位置(这似乎要求你的脸足够光滑,如果你说你长了很多青春痘,甚至整个脸皮处处连续处处不可导,那难度太大了,我就无话可说了哈哈),你的脸上总是会在相同的时间间隔内鼓起来一个相同高度的包来,并且假定以鼓起来的包的大小作为系统输出。好了,那么,下面可以进入核心内容——卷积了!
如果你每天都到地下去打台球,那么老板每天都要扇你一巴掌,不过当老板打你一巴掌后,你5分钟就消肿了,所以时间长了,你甚至就适应这种生活了……如果有一天,老板忍无可忍,以0.5秒的间隔开始不间断的扇你的过程,这样问题就来了,第一次扇你鼓起来的包还没消肿,第二个巴掌就来了,你脸上的包就可能鼓起来两倍高,老板不断扇你,脉冲不断作用在你脸上,效果不断叠加了,这样这些效果就可以求和了,结果就是你脸上的包的高度随时间变化的一个函数了(注意理解);如果老板再狠一点,频率越来越高,以至于你都辨别不清时间间隔了,那么,求和就变成积分了。可以这样理解,在这个过程中的某一固定的时刻,你的脸上的包的鼓起程度和什么有关呢?和之前每次打你都有关!但是各次的贡献是不一样的,越早打的巴掌,贡献越小,所以这就是说,某一时刻的输出是之前很多次输入乘以各自的衰减系数之后的叠加而形成某一点的输出,然后再把不同时刻的输出点放在一起,形成一个函数,这就是卷积,卷积之后的函数就是你脸上的包的大小随时间变化的函数。本来你的包几分钟就可以消肿,可是如果连续打,几个小时也消不了肿了,这难道不是一种平滑过程么?反映到剑桥大学的公式上,f(a)就是第a个巴掌,g(x-a)就是第a个巴掌在x时刻的作用程度,乘起来再叠加就ok了,大家说是不是这个道理呢?我想这个例子已经非常形象了,你对卷积有了更加具体深刻的了解了吗?
数学上:
好了,这时候水波就是g(t),扔进去的石头就是一个冲击函数。
物理上:
那如果连着扔好多块呢?产生的水波就是一个个石头产生的水波的叠加吧。(图不画了自己脑补吧╮(╯_╰)╭)只是这时候因为扔进去的石头有先后顺序,最近扔进去的那块石头对现在的水波影响比较大咯。(这就是g(t)要反一下的原因,g(t)在t=0处就是石头刚扔进去时候水的反应)
数学上:
这时候这一个个石头就是f(t),然后现在的水波,就是这一块块石头激起的水波的叠加,也就是f(t)与g(t)的卷积,只是这个时候f(t)是一个离散的函数罢了。看,因为我们扔石头是有先后的咯,所以产生的水波叠加实际上就是一个个有时间差的g(t)的简单叠加咯。
物理上:
那那...如果我们不扔石头了,我们...尿尿呢 (/ω?\) ...尿尿也会激起水~~波~~荡~~漾~~吧 _(:3 」∠)_ ...你看我们尿尿的时候,水波会不停地产生呢...
数学上:
这时候呢,我们尿的尿就是一个f(t)了...如果我们尿尿的力道忽大忽小...那么f(t)就会一会变大一会变小。(所以我们可以尿出我们想要的函数f(t)...嗯 (??`ω´?))然后这时候g(t)是什么呢。。。你可以理解为一滴尿液滴到水里产生的水波呢~~然后这忽大忽小的尿液,在湖里产生了连续不断的水波就是f(t)*g(t),卷积啦!
物理意义…不知楼主问的是不是卷积在信号系统中输入输出计算的意义。我觉得要想理解好卷积首先要理解好一个和线性代数相关的概念:线性映射。
什么是线性映射呢?
设想有个向量v
通过某种映射可以从v得到w,记为
w=T(v)
如果对于任意v,T都满足线性代数里两条巨基本的性质:
1) T(u+v)=T(u)+T(v)
2) aT(v)=T(av)
那么就称T是一个线性映射,其实就是简单的线性代数嘛
所以如果我们把输入信号和输出信号分别和v和w对应起来的话,似乎就是个线性映射系统啊。当然实际上没有这么简单,还多了那么一点东西,所以在信号与系统里,叫做线性时不变(Linear Time-Invariant, LTI)系统。
什么叫时不变呢?以冲击函数为例子说明一下,假设在t=0的位置上有个单位冲击,在某个LTI system中,造成的响应是在t=1到t=2区间上的某种曲线。
时不变的意思就是如果把这个t=0上的单位冲击移动到任意时间点,比如t0去,那么响应的曲线也会响应地移动到t=t0+1到t=t0+2的位置上去。
结合之前线性系统的定义,整个系统是线性的,所以是线性+时不变。
其实从不太严谨的线性代数角度也可以理解线性时不变系统,本质还是线性映射。比如考虑对如下对一个向量的变换:
可以看到,x移动了两个单位,这似乎就是时间轴上的移动操作啊。那如果变换矩阵再换换样子呢?
是不是觉得眼熟,其实这种斜边元素相等形式的矩阵就是一个不严谨的简单的线性时不变系统。更不严谨地说,如果我们把每个时间点都理解成一个向量中的一维的话,矩阵乘法操作的本质就是将这个向量投射到对应的维上,并且矩阵中的每个向量的Norm就是投影之后的系数。
总之,这就是线性时不变,说白了还是线性两个字最重要,于是,卷积求响应就很清晰了,就是简简单单的加法而已:已知冲激信号在输出中的响应的情况下,给定一个输入信号,只要把输入信号看成一个个时间点上的冲激信号的叠加,直接在输出的t时间点上把这所有冲激信号的响应简单加起来,就可以了。
写出来的话就是:
(注意这里t和τ∈N)
对应前面解释来看,其中f(τ)是输入信号在τ位置的值,也就是和单位冲激信号的比值,作为系数,而h(t-τ)则是距离t时刻τ那么多时间点上的一个冲激信号在t时刻所对应的响应。我这里说的都是离散情况举例,连续情况可以脑补推广。总之,不严谨地说,都是线性代数。
有一个挺有意思的例子,我以前在万门大学童哲校长的公开课上看到的:
我们在研究一个人一天24小时肚子里还剩多少东西,记为函数(小时)图:“今儿的吃法”
好了,我们现在看考察一下,比如说下午15点,这人肚子里还剩多少东西:加起来就是:
考察任意时间,就有卷积
你可以仔细看一下,上面t=15的情况可以从这个式子里面导出[注1]。
图:(上)吃法儿; (下)肚里食物量
其实在前面的讨论里面,我们把“吃饭的量”看作输入,“肚子里还剩的食物”看作输出。任何一个线性时不变系统都可以由冲激响应完全描述。
它的输入输出关系正是卷积:
/*-----------------------------------------------------------*/
上面有几个参数可以变化,我们来玩一玩:
比如说:这兄弟最近消化不大好,也就是冲激响应里面的衰减系数比较小,那么会怎么样呢?(最近消化不大好)上:吃法儿; 下:肚里食物量
你可以看到前一顿儿还没消化完,就到吃下一顿了。如果按原计划吃,实在要把肚子撑坏。。(峰值大概12,比上面的结果大了很多)
(花半小时吃饭,而且在正常消化状态)上:吃法儿; 下:肚里食物量
你可以看到,峰值(大概8)进一步减小了,细嚼慢咽就是好。
(花一个小时吃饭,而且在正常消化状态)上:吃法儿; 下:肚里食物量
嗯,对肚子压力更小一点了。
(过年的时候,每顿花三个小时,吃得多而且肠胃状态不那么好)上:吃法儿; 下:肚里食物量
可以看到,这个峰值和最前面那种瞬间吃好饭的状况差不多了,而且过了半夜还没消化完。
所以,不要暴饮暴食噢!
最后这一部分有点娱乐性质,主要是画图的时候玩儿上了。。哈哈
并且用线性时不变的模型来模拟人吃饭消化,可能有不妥,比如:人的胃不是线性的,因为有总容量的限制,而且估计也不是时不变,因为吃的太多时肚子消化能力大概会变弱。拿这个仅仅用于举例子咯。
希望有帮到你。也请各位朋友赐教!
/*-----------------------------------------------------------*/
[注1]t=15的情况:
激励(输入)
冲击响应(系统)
求响应(输出)
信号反褶是这样的:信号波形是我们按着信号脉冲出现的先后顺序记录下来的,所以信号波形的左边先出现,而信号反褶恰恰说明哪个信号脉冲先来,则这个脉冲就先进入系统,但是,在计算机上用软件进行科学计算时可以将系统单位冲激响应进行反褶,相应的,此系统可以看成非因果系统。
上述内容是摘自我写的《信号与系统分析和应用》书中的几节内容,希望能对大家有所帮助,并请提出意见和建议。简单来说就是卷积
就是x点所受的其它所有ξ点的影响的叠加。
比如信号与系统里,LTI系统的输出就是每个时刻的输入所产生的输出的叠加。
把x(t)送入冲激响应为h(t)的LTI系统得到y(t)。
τ时刻的输入是x(τ);τ时刻附近dτ时间的输入在t时刻产生的输出,就是把h(t)乘一个x(τ)dτ,再平移τ,就是x(τ)h(t-τ)dτ;把所有τ产生的输出叠加就是y(t)
写作
再如概率论中的全概率公式。
比如已知两个独立的随机变量X和Y的概率密度分布fx(x)和fy(y),求Z=X+Y的密度分布。
考虑X=x的条件下Y=z-x的处概率密度,正是fy(z-x),由全概率公式,此时Z的概率密度就是所有x处fx(x)fy(z-x)叠加。
卷积可以解决这个问题:
已知 如果输入是 冲击函数,输出是 函数 f()。
问 如果输入是 给定的函g(),输出是什么。
输出就是f g卷积一下。
1.
e是一个重要的常数,但是我一直不知道,它的真正含义是什么。
它不像π。大家都知道,π代表了圆的周长与直径之比3.14159,可是如果我问你,e代表了什么。你能回答吗?
维基百科说:
"e是自然对数的底数。"
但是,你去看"自然对数",得到的解释却是:
"自然对数是以e为底的对数函数,e是一个无理数,约等于2.718281828。"
这就构成了循环定义,完全没有说e是什么。数学家选择这样一个无理数作为底数,还号称这种对数很"自然",这难道不是很奇怪的事情吗?
2.
昨天我读到一篇好文章,它把这个问题解释得非常清楚,而且一看就懂。
它说,什么是e?简单说,e就是增长的极限。
下面就是它的解释。
3.
假定有一种单细胞生物,它每过24小时分裂一次。
那么很显然,这种生物的数量,每天都会翻一倍。今天是1个,明天就是2个,后天就是4个。我们可以写出一个增长数量的公式:
上式中的x就表示天数。这种生物在x天的总数,就是2的x次方。这个式子可以被改成下面这样:
其中,1表示原有数量,100%表示单位时间内的增长率。
4.
我们继续假定:每过12个小时,也就是分裂进行到一半的时候,新产生的那半个细胞已经可以再次分裂了。
因此,一天24个小时可以分成两个阶段,每一个阶段都在前一个阶段的基础上增长50%。
当这一天结束的时候,我们一共得到了2.25个细胞。其中,1个是原有的,1个是新生的,另外的0.25个是新生细胞分裂到一半的。
如果我们继续修改假设,这种细胞每过8小时就具备独立分裂的能力,也就是将1天分成3个阶段。
那么,最后我们就可以得到大约2.37个细胞。
很自然地,如果我们进一步设想,这种分裂是连续不断进行的,新生细胞每分每秒都具备继续分裂的能力,那么一天最多可以得到多少个细胞呢?
当n趋向无限时,这个式子的极值等于2.718281828...。
因此,当增长率为100%保持不变时,我们在单位时间内最多只能得到2.71828个细胞。数学家把这个数就称为e,它的含义是单位时间内,持续的翻倍增长所能达到的极限值。
这个值是自然增长的极限,因此以e为底的对数,就叫做自然对数。
5.
有了这个值以后,计算银行的复利就非常容易。
假定有一家银行,每年的复利是100%,请问存入100元,一年后可以拿多少钱?
回答就是271.828元,等于100个e。
但是,实际生活中,银行的利息没有这么高,如果利息率只有5%,那么100元存一年可以拿到多少钱呢?
为了便于思考,我们取n等于50:
我们知道,在100%利息率的情况下,n=1000所得到的值非常接近e:
因此,5%利息率就相当于e的20分之一次方:
20分之一正好等于5%的利率率,所以我们可以把公式改写成:
上式的rate就代表增长率。这说明e可以用于任何增长率的计算,前提是它必须是持续不断的复合式增长。
6.
再考虑时间因素,如果把钱在银行里存2年,可以得到多少钱?
在时间t的情况下,通用公式就是:
上式就是计算增长量的万能公式,可以适用于任何时间、任何增长率。
7.
回到上面的例子,如果银行的利息率是5%的复利,请问100元存款翻倍需要多少时间?
计算结果是13.86年:
上式最后一个等号,表明用72除以增长率,可以得到翻倍的大致时间,这就是72法则的来源。
source: 数学常数e的含义
**************************************************************************************
https://www.zhihu.com/question/22298352?rf=21686447
有那么麻烦吗?
不推荐用“反转/翻转/反褶/对称”等解释卷积。好好的信号为什么要翻转?导致学生难以理解卷积的物理意义。
这个其实非常简单的概念,国内的大多数教材却没有讲透。
直接看图,不信看不懂。以离散信号为例,连续信号同理。
已知下面通过演示求的过程,揭示卷积的物理意义。
第一步,乘以并平移到位置0:简单吧?无非是平移(没有反褶!)、叠加。
====================================================
从这里,可以看到卷积的重要的物理意义是:一个函数(如:单位响应)在另一个函数(如:输入信号)上的加权叠加。
重复一遍,这就是卷积的意义:加权叠加。
对于线性时不变系统,如果知道该系统的单位响应,那么将单位响应和输入信号求卷积,就相当于把输入信号的各个时间点的单位响应 加权叠加,就直接得到了输出信号。
通俗的说:
在输入信号的每个位置,叠加一个单位响应,就得到了输出信号。
这正是单位响应是如此重要的原因。
在输入信号的每个位置,叠加一个单位响应,就得到了输出信号。
这正是单位响应是如此重要的原因。
看了好多关于卷积的答案,看到这个例子才彻底地理解了这个过程~
关于卷积的一个血腥的讲解
比如说你的老板命令你干活,你却到楼下打台球去了,后来被老板发现,他非常气愤,扇了你一巴掌(注意,这就是输入信号,脉冲),于是你的脸上会渐渐地(贱贱地)鼓起来一个包,你的脸就是一个系统,而鼓起来的包就是你的脸对巴掌的响应,好,这样就和信号系统建立起来意义对应的联系。下面还需要一些假设来保证论证的严谨:假定你的脸是线性时不变系统,也就是说,无论什么时候老板打你一巴掌,打在你脸的同一位置(这似乎要求你的脸足够光滑,如果你说你长了很多青春痘,甚至整个脸皮处处连续处处不可导,那难度太大了,我就无话可说了哈哈),你的脸上总是会在相同的时间间隔内鼓起来一个相同高度的包来,并且假定以鼓起来的包的大小作为系统输出。好了,那么,下面可以进入核心内容——卷积了!
应题主邀,把我在卷积的物理意义是什么? - 数学的答案搬过来。
——————————————————————————————————————
对于初学者,我推荐用复利的例子来理解卷积可能更好理解一些:
用求和符号来简化这个公式,可以得到:
在上式中,为小明的存钱函数,而为存入银行的每一笔钱的复利计算函数。在这里,小明最终得到的钱就是他的存钱函数和复利计算函数的卷积。
为了更清晰地看到这一点,我们将这个公式推广到连续的情况,也就是说,小明在从到的这一段时间内,每时每刻都往银行里存钱,他的存钱函数为,而银行也对他存入的每一笔钱按复利公式计算收益:,则小明到时间将得到的总钱数为:
这也就是卷积的表达式了,上式可以记为。
最近在看Feedback Control of Dynamic Systems,趁此机会复习了一下卷积。
先看下图,左边是输入信号,右边是输出信号:(a)中,输入信号p(t)经过系统后得到输出信号h(t);
(b)中,输入信号较之于(a)延迟了τ,表示为p(t-τ),由于是LTI(线性时不变系统),输出信号也延迟τ,变为h(t-τ);
(c)、(d)两图阐释了LTI的叠加原理:若以p(t)+p(t-τ)为输入,则输出为h(t)+h(t-τ);
但我们仍有一些不太满意的地方,用p(t)表示的u(t)并不是精确的u(t)啊,那些小长条的面积比u(t)的面积可少了不少呢。除非Δ尽可能的小,长条尽可能的窄。
诶,这不就是积分么?
所以:这就是卷积,与其理解成翻转,不如理解成延迟后叠加。
20170708更新并纠正部分细节。
补充一下
和 的答案。
想要形象地理解卷积这个数学概念,那么思路当然是找它的“实际应用”,最容易想到的例子就是它对动态系统的输入输出关系的描述了。许多答主提到了“打脸”的例子,非常形象,就不复述了。
以下略纠结,嫌麻烦可以跳过直接看图。
之前我对
答案的有部分理解是不准确的;现在看来他对 的定义确实有问题(感谢 和 指出)。顺便提一下, 的答案引用了Franklin的那本经典教材,但是我越来越感觉它讲卷积的那部分并没有写得很清楚,比如它和配图和文字就对不上号,徒增了理解的难度。我尽量梳理一遍吧。
首先是单位脉冲和单位脉冲响应 的定义,大家应该都懂我就不啰嗦了。(注意 的“面积”为 。)再是定义 和 ,当 时趋近于 和 。(同样注意 的高度为 ,面积为 。)
一个输入信号可以拆成很多 的和,即 (所以
的答案确实少了个 )
输出即
取极限写成积分形式即
注意不要忘了 ,因为我们是要跟单位脉冲比较面积,即 是 的多少倍。
纠结部分完。
有了这个基础,再去理解卷积的公式就没那么intimidating了。放张图,心血来潮随手画的,但愿能顶个卵用。
另外,对于实际系统(因果系统),由于未来的打脸( )不会造成现在的脸肿,所以积分上限设为 就行了;此外我们还一般假设 之前没人打脸并且脸也不肿,所以积分下限设为 就行了。因此
从信号处理的角度讲,首先将上面定义中所有的x改为t并理解为时间。这样就是用一个激励信号f(t)去不完美地获取系统g(t)的响应信息。这种不完美体现在f(t)不是delta函数,具体说是两点,其一时宽度不为零,这样就不能简单确定信号发射的时间点(这不是一个完美的瞬时时刻),其二是强度也存在分布。卷积的问题是:当理论上的完美响应g(t)和这个不完美的激励信号f(t)都已知的情况下如何获取这个不完美的结果h(t)?这个不完美的结果h(t)就称之为f(t)和g(t)的卷积。
为了回答卷积为什么写成定义的样子,具体分析如下:
如果f(t)是完美宽度的,那么h(t)=f(0)*g(t)(没错,这里*是普通的乘积)就应该是结果了。也就是说仅仅是把g(t)乘以一个常数f(0)而已,而g(t)的形状严格不变。而且注意,这里激励的时间是没有误差的,从激励信号开始给出的时候开始计时(定义了完美的时间零点),得到的h(t)相对于g(t)就是完美的,没有时间延迟,只有一个f(0)的常数造成绝对强度的f(0)倍,而f(0)也是已知的,所以g(t)= h(t)/f(0) 或者h(t)=f(0)*g(t)。
现在的问题是激励信号f(t)不完美,那么就把激励信号f(t)分解为一系列的delta函数,如上图。这些delta函数的宽度都是完美的(宽为微分dt,是瞬时的),但是他们的高度不同(或者说有不同的权重),称之为“加权delta函数”。
下面,定义f(0)为时间零点,那么当f(0)这个加权delta函数激励的时候,扫描是完美的,时间上精确地定义了,强度上只是乘以了f(0)这个常数而已,这和上面说的完美扫描完全一模一样。但是f(1)这个加权delta函数激励的时候,扫描是不同而且不完美的,强度变了,变成f(1)倍而不是f(0)倍了,而且时间上因为f(1)自己出发的时候就晚了1秒,所以得到的这个h(t)分量的结果肯定也不完美,应该正好晚了1秒,就是说应该是g(t)这个函数平移1秒,也就是g(t-1),再乘以f(1)就对了,就是h(1)这个分量了。这样以此类推,每个加权delta函数激励的时候都有一个时间差,这个时间差就定义为tau。对于每一个tau,都有一个强度加倍因子叫做f(tau),还有一个时间差,叫做tau,所以最后得到的这个分量(或者说贡献)h(tau)就是等于f(tao)*g(t-tau)。这个式子就模拟了f(tau)这个加权delta函数扫描g(t)的过程以及得到的结果h(tau)。
最后,把这些h(tau)加起来,就是积分,也就是定义式了。信号反褶是这样的:信号波形是我们按着信号脉冲出现的先后顺序记录下来的,所以信号波形的左边先出现,而信号反褶恰恰说明哪个信号脉冲先来,则这个脉冲就先进入系统,但是,在计算机上用软件进行科学计算时可以将系统单位冲激响应进行反褶,相应的,此系统可以看成非因果系统。
上述内容是摘自我写的《信号与系统分析和应用》书中的几节内容,希望能对大家有所帮助,并请提出意见和建议。很喜欢唐常杰老师的解释:科学网—辐射、服碘、补盐、空袭和卷积-----教学难点讨论之一
公式有错位,推荐在源链接中看
另外推线唐老师的系列文章:科学网—从这里开始 -----系列博文的入口及其DIY
幼童背古诗文的感觉,来自数学系的同学觉得卷积是小菜一碟,随手就写出卷积定义
F(t)= ∫ f(τ)g(t-τ)dτ (积分限从-∞ 到+∞)
并指出这是含参积分,t是参数,觉得浅而又显,无须解释。而部分(例如来自工科和医学专业的)选修数据挖掘的学生,还是觉得稍有点难,说:相关公式能默写、能推导、能通过考试,自己还是觉得不踏实,觉得没有真正理解;发明者是怎样想出来的?有何直观背景?用在哪些场合?
一言以蔽之,在逻辑上认可,而直观上迷茫。好像很小的时候背诵古诗文那种感觉。
鉴于数学老师已经讲解过理论推导,作为一种补充,这里用生活实例做一些直观解释,给出一个大框架和物理直观,为叙述简单,忽略一些细节。需要说明,直观的解释仅用于辅助理解,不能取代严格的描述和证明。
几个时髦(但可能不很贴切)的例子.
辐射:设某核电站事故中,某工作人员每天到抢险现场工作T分钟,接受一定剂量的辐射,辐射会自然地衰减,如此工作N天,总的辐射量用什么计算工具来(粗略地)估计?回答:可以用卷积。
服碘:某人为了防辐射,自己找来碘片,每天口服若干,体内碘残量会随人体代谢衰减,N天后体内积累的碘残量如何(粗略地)估计?还是卷积;(后面科普部分将给出简单的推导过程);
补盐:某人为了反辐射,抢购来碘盐,每餐口服若干,体内盐残量会随人体代谢衰减。N天后体内积累的盐量和碘残量如何(粗略地)估计?可以用卷积;
空袭:某多国部队每隔N小时对桀骜不驯的某地区或国家实行间歇性空中打击,每次打击后,其物理破坏和心理震慑作用会随时间衰减(例如,被打方会组织抢修,心里承受度增加等等),如此进行M天后,累积的打击总效果如何(粗略地)估计?还是可以用卷积。
还有其他例子,如长期服药的血药浓度,长期吸入污染物在人体内的积累,吸烟或喝咖啡的积累效应,多次喷洒农药的残留量,等等,也可以用卷积来估计。
上面的有些例子可能不很贴切,有几个原因:,
(a)卷积是积分运算,处理对象要求是可以积分的函数,在工程中,一般对应于连续现象而不是离散对象;把离散对象当做连续的现象处理,只能粗略估计。
(b)社会问题,政治问题比较复杂,即使加上很多假定,也只是框架性的估算。
但是,有计算、有依据的估计总比算命先生的神仙数字可信。
难懂之因:为了数学美,拆卸了脚手架。 教科书书常用“定义—定理”的体系,先给出数学定义,然后给出若干性质, 从公式 到 公式,逐步推导。有的教科书采用用信号“反褶、平移、相乘、积分”给出几何解释,属于用数学解释数学,提问者不满足这种解释。
这不是当年发明卷积的大师们的“需求–猜想—发现—证明—应用”的路径,大师们建设好“卷积”大厦后,为了数学美,拆卸了脚手架,现在人们看到的是炼成的钢铁,看不出钢铁是怎样炼成的。造成了部分非数学专业学生的一个难点。
一次输液引出的班门弄斧 一次偶感风寒,服药未愈,转作静脉滴注,无聊地望着那药液慢腾腾地滴,忽然灵感一闪:
(1)这是一个可离散观察的连续过程。透明玻璃管构成了可视化的界面,能离散地对药滴计数,而下面是相对稳定的液柱高度,保证了药液连续(有点脉动)地注入静脉,比较适合积分处理;(口服和注射,就相对离散,结果就更粗略一些)。
(2)药动学有个术语血药浓度,怎样来保证血药浓度在安全阈值之下,又在有效阈值之上呢?
立刻在草稿本上写划,哇噻,原来可以用卷积!而且只需要简单的积分知识。于是,对此常问难点,有了一个易懂的直观解释。正是:小恙滴注,焉知非福?
下面将叙述这次双重的(数学与医学)的班门弄斧,疏漏之处,请专家指正。
静脉滴注与体内药物浓度 为简单又不失一般性,给出下列符号和假定:
从t=0开始,每隔τ秒,输入药物一次(离散化是为了简单);药量随时间变化, 在时刻t时的那次给药量为f(t),关注的时刻点为 t=0, τ,2τ,3τ,…
一滴药液的在体内衰减规律 药物以多种方式代谢(衰减),按假设,在τ1时的那滴药液含药量f(τ1),当时间流逝到t时刻,假设那一滴药物在体内的残量是f(τ1) g(t,τ1),其中g(t,τ1)称为衰减因子函数,怎么找出衰减因子的具体结构呢?药动学中有两种衰减方式 :
(a)零级动力学消除,即恒速消除,如乙醇血浓>0.05 mg/ml时,较简单;
(b)一级动力学消除,即恒比消除,消除速度与血药浓度成正比,如乙醇血浓<0.05 mg/ml时的衰减规律,这也类似于简单热传导中散热速度与温差成正比。
设在τ时刻 ,输入一滴药,药量为f(τ), 根据一级动力学消除,建立最简单的微分方程 ;
dg/dt =-kg
考虑t=τ时不衰减的初始条件,容易求得 g=e-k(t-τ)
,
为下面方便,把衰减因子改写为
g(t-τ)= e-k(t-τ)
于是,在τ>0时,给药一次,药量为f(τ),
当t为2τ时,血药浓度降到 f(τ)*g(t-τ)= f(τ)*g(τ)= f(τ)(1/ek)
当t为3τ时,血药浓度降到 f(τ)*g(t-τ) = f(τ)*g(2τ)= f(τ)(1/e2k)
当t为4τ时,血药浓度降到 f(τ)*g(t-τ) = f(τ)*g(3τ)= f(τ)(1/e3k )
可见,只给药一滴,血药浓度衰减很快,难以治疗那种要与病毒或细菌打持久战的疾病。
多次密集给药 或连续给药 上面是只在τ>0时,给药一次,现在,让τ动起来,n次给药,给药时间依次取τ1, τ2,….. τn ;则n次密集给药后,当时间流逝到t时的血药浓度是
∑j f(τj)*g(t-τj) ( 对 j=1,2,….. n 求和)
前面说过,静脉滴注是一个可离散观察的连续过程。,所以,上面的和式可写为积分形式,即卷积
F(t)= ∫ f(τ)g(t-τ)dτ
曲线光滑工具 当f(τ)是脉冲函数时(例如考察一滴药引发的血药浓度),曲线显得不够光滑,而卷积F(t)是多次脉冲的(平均)累积效应,或可视为是一种加权平均,所以,F(t)的曲线就光滑一些,所以,医生要考察N小时的滴注效果,而不察几分钟或一滴药的效果。选择适当的g(t-τ)函数,(例如,3/2次方衰减型、平方衰减性、指数衰减型、周期兼指数衰减型,...),可用卷积作为突出不同加权方式的曲线光滑工具。
比较光滑、不是陡升陡降的血药浓度曲线表明,静脉滴注能较好地控制血药浓度;这大概也是有些医生和病人喜欢它的原因;当然,如果过分依赖静脉滴注,则减少了免疫系统的锻炼机会,所以很多医生主张,如果服药能解决问题,就不要滴注。
更多的应用实例 卷积的结果可辅助人们定量地协调脉动式输入 f(τ) 和 衰减g(t-τ) 这一对矛盾,使得累积效应F(t)= ∫ f(τ)g(t-τ)dτ 在控制范围内。
例如,研究干预规则,(例如,叶酸干预新生儿脑畸形缺陷),干预为f(τ),复杂的衰减g(t-τ),总的干预效果可否用卷积来粗略描述?
再例如,制定正确的给药剂量和周期,例如照医嘱摄入碘或盐;
又例如,制定空中打击方案的强度和频度,常识告诉人们,足够的强度和密度才能有效打击。卷积作为工具,或许可定量计算出最经济打击强度和密度。而被打击的一方,可计算出足够的衰减因子,使得能在被轰炸后有效恢复;战争是铁血与智慧的较量,当双方的铁与血差不多时,如《孙子.计篇》所说,“多算胜,少算不胜”,而卷积只不过在众多的计算方法基础上,增加了一个算法,仅此而已,
最近,在这个不平静的世界上,有一场空袭和反空袭的较量,不知持续多久?10天,100天,还是200天?研究军事的专家或许可用卷积做个模型。
武侠小说中,有时候看见一方逐步投入兵力,使用添油战术,好像是多次服药,每次都没有服够量,血药浓度低于有效门限。被逐次歼灭。
卷积并不神秘,它有其退化版,例如水池一面进水,一面放水,求瞬时水量。当进水匀速且放水速度服从零级或一级动力学消除规律时,偶尔也作为中小学生的数学奥赛题,基础好的聪明学生能用初等方法计算。但当进水是1+sin(t)这样的脉动函数,或更复杂的函数时,就只能用卷积了。
原文:http://www.cnblogs.com/xfzhang/p/7261180.html