信号与系统(二)傅立叶变换

函数正交的概念

\left(t_1, t_2\right) 区间的两个函数 \varphi_1(t) \varphi_2(t) 可以看成是由n个点组成的,映射到n维向量空间可得两个向量:

\vec{\varphi_1}= \left [ \begin{array}{} a_1\\a_2\\\dots\\a_n \end{array} \right ] \quad \vec{\varphi_2} = \left [ \begin{array}{}b_1\\b_2\\\dots\\b_n \end{array} \right ]

若这两个向量垂直,有:

\vec{\varphi_1}\cdot \vec{\varphi_2} = \sum_{i=1}^n a_i b_i =0

上式中的 \sum_{i=1}^n a_i b_i 不就是 \int_{t_1}^{t_2} \varphi_1(t) \varphi_2(t) dt 嘛!!事实上,绝大多数满足正交的函数都是复函数,所以我们将正交的概念扩展到复平面:

【定义】在 \left(t_1, t_2\right) 区间的两个函数 \varphi_1(t) \varphi_2(t) ,若满足:

\int_{t_1}^{t_2} \varphi_1(t) \varphi_2(t) ^*dt = 0

\varphi_1(t) \varphi_2(t) \left(t_1, t_2\right) 区间正交。其中 \varphi_2(t)^* \varphi_2(t) 的共轭函数。

广义傅立叶级数

函数映射到向量空间,一定可以分解成两两相互垂直的向量的线性组合:

\vec{V}=c_1 \vec{V}_1+c_2 \vec{V}_2+...+c_n \vec{V}_n

参数也可以用几何方法计算出来:

c_n=\frac{|V| \cos \theta_n}{\left|V_1\right|}=\frac{\vec{V} \cdot \vec{V}_n}{\vec{V}_n \cdot \vec{V}_n}

既然函数映射到向量空间可以分解,那么在函数空间一定也可以分解为两两相互正交的函数的线性组合:

f(t)=C_1 \varphi_1(t)+C_2 \varphi_2(t)+\cdots+C_n\varphi_n(t)+\cdots=\sum_{n=1}^{\infty} C_n \varphi_n(t)

这个公式被称为广义傅立叶级数。其中, \left\{\varphi_1(t),\varphi_2(t),...,\varphi_n(t)\right\} 是完备正交函数集。参数 C_n 也可以类别向量的计算方法得到:

C_n=\frac{\int_{t_1}^{t_2} f(t) \varphi_n(t) \mathrm{d} t}{\int_{t_1}^{t_2} \varphi_n^2(t) \mathrm{d} t}

傅立叶变换

我们找一个具体的完备正交函数集的例子:

\left\{e^{j \omega t}, \omega \in (-\infty,+\infty)\right\}

根据广义傅立叶级数,可得:

f(t)=\frac{1}{2 \pi} \int_{-\infty}^{\infty} F(j \omega) e^{j \omega t} d \omega

它的参数 C_n 是一个关于 j \omega 的函数:

F(j \omega)=\int_{-\infty}^{\infty} f(t) \mathrm{e}^{-j \omega t} \mathrm{~d} t

有趣的地方来了, f(t) 分解后的式子中, F(j\omega) 是参数,只要确定了 F(j\omega) f(t) 就确定了。也就是说函数 F(j\omega) 与函数 f(t) 存在一一对应的关系。

f(t) 是关于t的函数,我们说这是时域里的函数。 F(j\omega) 是关于 j\omega 的函数,我们说这是频域里的函数。现在一个时域函数 f(t) ,不仅在向量域有与之对应的向量 \vec{f} ,还在频域有与之对应的函数 F(j\omega)

任意信号的响应

基本信号 e^{j \omega t} ,直接代入指数形式的傅立叶级数,可求出响应为 H(j \omega) e^{j \omega t} 。与时域类似,由线性系统的性质,可推导出任意信号的响应:

e^{j \omega t} \longrightarrow H(j \omega) e^{j \omega t}

\frac{1}{2 \pi} F(j \omega) e^{j \omega t} \longrightarrow \frac{1}{2 \pi} F(j \omega) H(j \omega) e^{j \omega t}

\frac{1}{2 \pi} \int_{-\infty}^{\infty} F(j \omega) e^{j \omega t} ~d \omega \longrightarrow \frac{1}{2 \pi} \int_{-\infty}^{\infty} H(j \omega) F(j \omega) e^{j \omega t} ~d \omega

f(t) \longrightarrow y(t) = F^{-1} [H(j \omega) F(j \omega) ]

注:因为t是从 -\infty 开始的,而 t=-\infty ,我们认为初始状态为0,所以这里的 y(t) 是全响应。

求解任意信号f(t)的响应y(t)的步骤:

  1. f(t) 的傅立叶变换 F(j \omega)
  2. 求线性时不变系统的傅立叶变换 H(j \omega)
  3. 计算 Y(j \omega) = F(j \omega) H(j \omega)
  4. Y(j \omega) 作傅立叶逆变换得到 y(t)

周期函数的傅立叶级数

上面我们把f(t)分解为 e^{j \omega t} ,这是个复函数,实际运用比较难,能不能把f(t)分解为实数函数呢?当f(t)复合特定条件的时候,是可以将f(t)分解为实函数的线性组合的:

设周期信号f(t),其周期为T,角频率为 \Omega=2 \pi / T ,当满足狄里赫利条件是,可向三角函数集分解:

\left\{\varphi_1(t),\varphi_2(t),...,\varphi_n(t)\right\} = \{1, \cos (n \Omega t), \sin (n \Omega t), n=1,2, \ldots\}

f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty} a_n \cos (n \Omega t)+\sum_{n=1}^{\infty} b_n \sin (n \Omega t)

对应的参数为:

\frac{a_0}{2}=\frac{1}{T} \int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) \mathrm{d} t

a_n=\frac{2}{T} \int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) \cos (n \Omega t) \mathrm{d} t

b_n=\frac{2}{T} \int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) \sin (n \Omega t) \mathrm{d} t

posted @ 2024-08-11 16:41:34
评论加载中...
发表评论