Matlab——周期信号的傅里叶级数

通信系统中利用信号来传递信息，确定信号是时间的确定函数。通信系统中的通信信道及其收发设备中的很多部分可以等效成线性时不变系统（LTI）来建模，因此，确定信号的性质通过线性系统的分析是分析通信系统的数学工具。

周期信号的傅里叶级数

        周期信号定义为随时间变换，取值呈周期变化的信号，即, k为整数，T称为信号的周期。一个正弦型信号源即为一个典型的周期信号。如果周期信号在一个周期内可积，则可以通过傅里叶级数展开该周期信号。傅里叶级数展开如式

                                                               

                            

其中，T为周期信号的最小周期，；为傅里叶展开系数，其物理意义为频率分量的幅度和相位。 

例：设周期信号的一个周期波形为

                                                   

求该周期信号的傅里叶级数展开解析式，并用Matlab画出傅里叶级数展开后的波形。

解：

由上式可以表明，信号可以展开成一系列频率为1/T整数倍的正弦、余弦信号的加权叠加，其中相应的频率分量加权系数即为。下例中采用Matlab程序画出了取 2N+1项近似式的波形（N=100）

 

%周期信号（方波）的展开，fb_jinshi.m
close all;
clear all;
N = 100;      %去展开式的项数为2N+1项

T = 1;
fs = 1/T;
N_sample = 128;  %为了画出波形，设置每个周期的采样点数
dt = T/N_sample;

t = 0:dt:10 * T - dt;  %不管是否减去抽样间隔dt，所画出图是一样的
n = -N:N;
Fn = sinc(n/2) .* exp(-j * n * pi/2);
figure(1)
plot(n,Fn);
Fn(N+1) = 0;
ft = zeros(1, length(t));    %为了使画出的-1、1延续一段时间
figure(2)
plot(t,ft);
for m = -N:N;
    ft = ft +Fn(m+N+1) * exp(j * 2 * pi * m * fs * t);
end
figure(3)
plot(t,ft);

 

傅里叶系数展开

原周期信号波形

 

 

 

 

（资料来源于通信原理基于Matlab的计算机仿真）