世界的本质是旋转(2) 旋转叠加与傅里叶变换

1.旋转的叠加

上一次,我们观察了单一余弦震动可以表示为两个等模共轭矢量的旋转运动叠加。其实,早在高中时,课本里就有波的叠加的概念。假设现在在复平面上有两个旋转分量,其在某个时刻的状态(姑且叫做0时刻吧)如下图所示:
复平面的两个旋转
从数学上说,在t时刻,两对(共4个)矢量的加成为:
s ( t ) = a 1 ⋅ e 2 π j f 1 t + j θ 1 + e − 2 π j f 1 t − j θ 1 2 + a 2 ⋅ e 2 π j f 2 t + j θ 2 + e − 2 π j f 2 t − j θ 2 2 s(t) = a_1\cdot \frac{e^{2\pi jf_1t+j\theta_1} +e^{-2\pi jf_1t - j\theta_1}}{2} + a_2\cdot \frac{e^{2\pi jf_2t+j\theta_2} +e^{-2\pi jf_2t - j\theta_2}}{2} s(t)=a12e2πjf1t+jθ1+e2πjf1tjθ1+a22e2πjf2t+jθ2+e2πjf2tjθ2
根据上一篇的欧拉公式,
s ( t ) = a 1 ⋅ c o s ( 2 π j f 1 t + θ 1 ) + a 2 ⋅ c o s ( 2 π j f 2 t + θ 2 ) s(t) = a_1\cdot cos(2\pi jf_1t+\theta_1) + a_2\cdot cos(2\pi jf_2t+\theta_2) s(t)=a1cos(2πjf1t+θ1)+a2cos(2πjf2t+θ2)
其中,
$f_1=300Hz, f_2 = 500 Hz, \theta_1 = \frac{\pi}{6}, \theta_2 = \frac{3}{4}\pi $
$a_1 = 1, a_2 = 0.45 $
画出来是酱紫的:
时域

t = 0:0.00025:1-0.00025;
f1=300;
f2=200;
theta1=pi/6;
theta2=pi*3/4;
s = cos (2*pi*f1* t + theta1) + 0.45*cos(2*pi*f2*t+theta2);
>> plot (t,s);

当然了,我们可以任意安排成对的共轭旋转,获得千奇百怪的图案,比如这个由6组幅度、起始角度、转速(频率)组成的运动:

 t = 0:0.00025:1-0.00025;
 a = [7,1300,7,7,535,7];
 f = [ 5 ,501,997, 1005,1501, 1997]';
 theta =[-1.4182,-1.1831,1.4182,-1.4182,-0.3977,1.4182]'; 
 s = a * cos(2 * pi * f * t + theta * ones(1,4000)); 
 plot(t,s)

6个运动

已经和余弦不太像了吧!

如果我们找到几千组呢?
几千个旋转合成

可以想象,使用N多组旋转,就可以组合出任意的形状来。
太好了!搜索互联网,发现在这篇博文中,已经把它说的非常清楚啦!不过要说明的是,下面动画中合成的并不是实信号,而是取的虚半轴投影。
“知乎专栏”
看看方波的合成:
维基百科
看看三角波:
维基百科
上面的图片是维基百科的哦!

2.傅里叶变换

上面的概念如果反过来,给定一个波形,我如何构造一组旋转,可以恰好使我们的复平面叠加点上的小钢笔画出同样的形状呢?显然我们需要知道每个旋转的转速(频率)、棍子长度(幅度)和起始0时刻的角度(相位)。

2.1 解旋转与均值

傅里叶变换是解决这个问题的工具,他的核心理念很简单。想知道这个波形是否含有某个频率f的旋转?? 那就按照f 反着转,看看有没有棍子保持不动吧!先看两个图,
a: 频率为f的旋转以-f解旋转,先向前转了几度,再向后倒退了相同度数,因此无论何时观察,都保持不动.
图中,红色为t时刻的棍子,绿色为 t+时刻其应该在的位置;但是,由于我们施加了一个和f相反的旋转-f,结果t+时刻棍子被转回去了(蓝色)。
被以相同频率反着转的棍子
如果我们不断的对这个图片拍照,该棍子永远停留在相同的位置。如果对这个位置的坐标求重心(均值),则其与上图位置重合。

b: 频率不等于f的旋转以-f解旋转,尽管又叠加了f,但仍旧为旋转状态,各个时刻位置不断变化
频率不动,旋转依旧
如果我们对它拍照,底片是这样的
均值为0
如果对它求均值或者求和,恒为0.

2.2 施加旋转的数学描述

在复平面上,把一个矢量 α \alpha α旋转一个角度 θ \theta θ ,旋转后的矢量 α ′ \alpha ' α为:
α ′ = α ⋅ e j θ \alpha ' = \alpha \cdot e^{j\theta} α=αejθ
对应到2.1节的需求,对任一个时刻t, 波形s(t)按照给定频率-f旋转后的新波形
s ′ ( t ) = s ( t ) ⋅ e 2 π ( − f ) t ⋅ j = s ( t ) ⋅ e − 2 π j f t s'(t) = s(t) \cdot e^{2\pi (-f)t \cdot j} = s(t) \cdot e^{-2\pi jft} s(t)=s(t)e2π(f)tj=s(t)e2πjft

2.3求和/求均值的数学描述

求和在数学上是进行积分

S ( f ) = ∫ − ∞ ∞ s ′ ( t ) d t = ∫ − ∞ ∞ s ( t ) ⋅ e − 2 π j f t d t S(f) =\int_{-\infty}^\infty s'(t) dt = \int_{-\infty}^\infty s(t) \cdot e^{-2\pi jft} dt S(f)=s(t)dt=s(t)e2πjftdt

这正是傅里叶变换的基本公式。

sjf

3What’s the next?

太多问题啦!详见下一章(世界的本质是旋转3-拍照与采样
为什么实信号的傅里叶变换是对称的?
课本里模拟波形的采样,也就是 A/D,怎么好好理解呢?

下一篇在这里

  • 6
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

丁劲犇

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值