信号与系统频域分析(1)——傅里叶级数

        三角函数形式的傅里叶级数:在直流分量(常数项)上叠加不同频率的正弦波来近似周期函数。

      一、基本形式:

        函数f(t)的周期为T_1,角频率为\omega_1=\frac{2\pi}{T_1},频率f_1=\frac{1}{T_1},傅里叶级数展开表达式为:

                f(t)=a_0+a_1cos(\omega{t})+b_1sin(\omega{t})+a_2cos(2\omega{t})+b_2sin(2\omega{t})+...+a_ncos(n\omega{t})+b_nsin(2\omega{t})+... =a_0+\sum_{n=1}^{\infty}[a_ncos(n\omega_1t)+b_nsin(n\omega_1t)]   ....(1)

        其中,n为正整数,各次谐波成分的幅度值按以下各式计算(三角函数的正交性):

        直流分量(等式两边在一个周期上积分,数值相等,确定系数):

                         a_0=\frac{1}{T_1}\int_{t_0}^{t_0+T_1}f(t)dt

        余弦分量的幅度(等式两边同时乘以cos(n\omega{t}),然后在一个周期上积分,数值相等,确定系数):

                         a_n=\frac{2}{T_1}\int_{t_0}^{t_0+T_1}f(t)cos(n\omega_1t)dt                                     ...(3)

        正弦分量的幅度(等式两边同时乘以sin(n\omega{t}),然后在一个周期上积分,数值相等,确定系数):

                         b_n=\frac{2}{T_1}\int_{t_0}^{t_0+T_1}f(t)sin(n\omega_1t)dt                                      ...(4)

 

       二、合并形式

       将(1)中同频率项加以合并,可以写成另一种形式

                         f(t)=c_0+\sum_{n=1}^{\infty}c_ncos(n\omega_1t+\phi_n)    或   f(t)=d_0+\sum_{n=1}^{\infty}d_nsin(n\omega_1t+\theta_n)

        其中

                         a_0=c_0=d_0      c_n=d_n=\sqrt{a_n^2+b_n^2}   a_n=c_ncos\phi_n=d_nsin\theta_n   b_n=-c_nsin\phi_n=d_ncos\theta_n

                         tan\theta_n=\frac{a_n}{b_n}     tan\phi_n=-\frac{b_n}{a_n}

    

       三、指数形式(从三角波叠加到螺旋波叠加!!!)

       利用欧拉公式,将三角函数表达成指数形式。

       由欧拉公式 e^{j\theta}=cos\theta+jsin\theta 

       得到  cos(n\omega_1{t})=\frac{1}{2}(e^{jn\omega_1{t}}+e^{-jn\omega_1{t}})      sin(n\omega_1{t})=\frac{1}{2j}(e^{jn\omega_1{t}}-e^{-jn\omega_1{t}})

        将上式带入,可得:

        f(t)=a_0+\sum^{\infty}_{n=1}(\frac{a_n-jb_n}{2}e^{jn\omega_1t}+\frac{a_n+jb_n}{2}e^{-jn\omega_1t})

         令    F(n\omega_1)=\frac{1}{2}(a_n-jb_n)                                                          ...(9)

         考虑到a_n是n的偶函数,b_n是n的奇函数,可得

                 F(-n\omega_1)=\frac{1}{2}(a_n+jb_n)     

         将上式带入,得到:

                 f(t)=a_0+\sum_{n=1}^{\infty}[F(n\omega_1)e^{jn\omega_1{t}}+F(-n\omega_1)e^{-jn\omega_1{t}}]

         令F(0)=a_0,考虑到

                 \sum_{n=1}^{\infty}F(-n\omega_1)e^{-jn\omega_1{t}}=\sum_{n=-1}^{-\infty}F(n\omega_1)e^{jn\omega_1{t}}

         得到f(t)的指数形式傅里叶级数,它是:

                  f(t)=a_0+\sum_{n=1}^{\infty}F(n\omega_1)e^{jn\omega_1{t}}

         将(3)(4)带入式(9)得到指数形式傅里叶级数的系数F(n\omega_1)(或简写为F_n),如下:

         (太美妙了,和谐的一匹!!!)

                 F_n=\frac{1}{T_1}\int_{t_0}^{t_)+T_1}f(t)e^{-jn\omega_1{t}}dt

          因为F_n一般是复函数,所以称这种频谱为复数频谱!!!根据F_n=|F_n|e^{j\phi_n},可以画出复数幅度谱|F_n|-\omega与复数相位谱\phi_n-\omega

          关于负频率的讨论:负频率的出现完全是数学运算的结果,并没有任何物理意义(很多问题都是从物理意义的角度入手,最后通过抽象的数学得到纯数学量的结论,数学的魅力啊!!),只有把负频率项与相应的正频率项成对地合并起来,才是实际的频谱函数。

 

       PS: 盗图http://blog.sina.com.cn/s/blog_57ad1bd20100txgs.html

                              

                 

 

 

首先,需要知道三角波的数学表达式为: $$f(t)=\frac{4A}{T} \cdot \left( \frac{t}{T} - \left\lfloor \frac{t}{T} +\frac{1}{2}\right\rfloor \right) $$ 其中 $\lfloor \cdot \rfloor$ 表示向下取整函数。 根据傅里叶级数的公式,可以推导出三角波的频谱表达式: $$F(k)=\frac{2A \cdot (-1)^k}{\pi k}$$ 其中 $k$ 为频率,$F(k)$ 为频率 $k$ 对应的频谱值。 根据 Nyquist-Shannon 采样定理,需要对该三角波进行采样,使得采样频率至少是三角波的最高频率的两倍。因为三角波最高频率为 $f_{max}=1/2T$,所以采样频率应该为 $f_s \geq 2f_{max}=1/T$。 假设采样频率为 $f_s=10/T$,则可以使用离散傅里叶变换来求解该三角波的频谱。离散傅里叶变换的公式为: $$X(k)=\sum_{n=0}^{N-1} x(n)e^{-i 2 \pi kn/N}$$ 其中 $x(n)$ 为时域信号,$X(k)$ 为频域信号,$N$ 为采样点数,$k$ 为频率。 因为三角波为周期信号,所以这里只需要记录一个周期的波形即可。取 $N=10$,则每个周期内有10个采样点。 根据三角波的数学表达式,可以计算出每个采样点的幅值,如下表所示: | n | t | f(t) | |---|---|-----| | 0 | 0 | 0 | | 1 | T/10 | 2A/5 | | 2 | 2T/10 | 4A/5 | | 3 | 3T/10 | 2A/5 | | 4 | 4T/10 | 0 | | 5 | 5T/10 | -2A/5 | | 6 | 6T/10 | -4A/5 | | 7 | 7T/10 | -2A/5 | | 8 | 8T/10 | 0 | | 9 | 9T/10 | 2A/5 | 将其进行离散化,得到时域信号 $x(n)$: $$x(n)=[0, 2A/5, 4A/5, 2A/5, 0, -2A/5, -4A/5, -2A/5, 0, 2A/5]$$ 利用离散傅里叶变换可以求得该信号的频域信号 $X(k)$。具体计算过程可以使用 FFT 算法进行计算,得到的结果如下图所示: ![image.png](attachment:image.png) 从图中可以看出,三角波的频谱是一系列的谐波,其中基波的频率为 $f_1=1/T$,幅值为 $8A/\pi$,其他谐波频率的幅值为基波幅值的 $\frac{2}{n\pi}$ 倍。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值