Archive

Posts Tagged ‘Math’

特征值与特征向量(转)

November 13, 2012 Leave a comment

特征值特征向量

特征向量确实有很明确的几何意义,矩阵(既然讨论特征向量的问题,当然是方阵,这里不讨论广义特征向量的概念,就是一般的特征向量)乘以一个向量的结果仍是同维数的一个向量,因此,矩阵乘法对应了一个变换,把一个向量变成同维数的另一个向量,那么变换的效果是什么呢?这当然与方阵的构造有密切关系,比如可以取适当的二维方阵,使得这个变换的效果就是将平面上的二维向量逆时针旋转30度,这时我们可以问一个问题,有没有向量在这个变换下不改变方向呢?可以想一下,除了零向量,没有其他向量可以在平面上旋转30度而不改变方向的,所以这个变换对应的矩阵(或者说这个变换自身)没有特征向量(注意:特征向量不能是零向量),所以一个变换的特征向量是这样一种向量,它经过这种特定的变换后保持方向不变,只是进行长度上的伸缩而已(再想想特征向量的原始定义Ax=cx,你就恍然大悟了,看到了吗?cx是方阵A对向量x进行变换后的结果,但显然cx和x的方向相同),而且x是特征向量的话,ax也是特征向量(a是标量且不为零),所以所谓的特征向量不是一个向量而是一个向量族, 另外,特征值只不过反映了特征向量在变换时的伸缩倍数而已,对一个变换而言,特征向量指明的方向才是很重要的,特征值不是那么重要,虽然我们求这两个量时先求出特征值,但特征向量才是更本质的东西!

Categories: Uncategorized Tags: ,

卷积的通俗理解

November 11, 2012 Leave a comment

光学成像里面经常会用到卷积,对卷积的概念总是有时候清楚有时候模糊。以前看过一篇《大牛很通俗的通俗的介绍信号与系统》的转帖,原帖已经找不到了。转帖是科学网看到的:http://bbs.sciencenet.cn/thread-554282-1-1.html。今天又拿出来读了一遍,感觉又理解得清楚了一些,可是我可以预料到过几天我又会忘了,所以还是用自己的话再加工加工,争取理解正确。

时域信号f(t)与g(t)的卷积在文中有三个非常形象的例子:

例子1:无穷时间连续信号的卷积;

经理让新来的技术员(小王)测试一个器件,他给了技术员1000多个输入函数,让技术员告诉他这些函数的输出是什么。

小王犯憷了,1000多个函数难道一个个测啊,那什么时候能测完啊?

上帝给了小王一个方法:输入一个单位脉冲信号(能力1J)到器件里,得到这个脉冲的输出(响应),然后把输入函数想象为无数个不同幅值脉冲按时间先后顺序排队进器件,输出就是这些脉冲的响应的叠加,但是最先输入的脉冲在时间轴的最左边,因此需要翻转后叠加。

例子2:无穷时间周期信号的傅里叶变换。

经理又给了小王一个任务,这次他告诉小王一个输入波形(不能用数学函数描述的波形),这个波形时间有限长,但是总是不停的周期性的输入,让小王告诉他输出会是什么样的。

这次第一次的方法不能用了,函数表达式不知道,没法直接算出卷积来,小王束手无策。

好心的上帝又告诉了小王一个方法:你可以把这个周期信号转换(投影)到另外一个域(一组不同的基构成的域),然后这个不平整的无穷周期信号在那个域中就不是无穷的了,然后回变得平整,这个变换叫傅里叶变换,时间域的卷积在频率域是简单的相乘,完了再将乘积变化回时间域就得到了输出信号了。

上面这两个例子形象的简介卷积的定义和它与傅里叶变换的关系(其实并没有将为什么卷积在频域是乘积)。

对于光学成像系统(linear imaging system),输入一般是空间域中的pure amplitude or phase or complex amplitude,卷积的kernel就是PSF(point spread function),输出就是输入与PSF的卷积。从频域来看的话,kernel的傅里叶变化叫OTF(optcal tranform function), 输入输出分别是频率谱.

Categories: Uncategorized Tags: