通信原理Matlab仿真----周期方波的傅里叶级数

% 周期信号(方波)的展开
clear;clc;close all;
N = 100;            % 取展开式的项数为2N+1项
T = 1;              % 方波周期
fs = 1/T;           % 方波频率
N_sample = 128;     % 为了画出波形，设置每个周期的采样点数
dt = T/N_sample;    % 时间采样间隔
t = 0:dt:4*T-dt;    % 时间变量，这里设置了4个周期的长度
n = -N:N;           % 傅里叶级数项数
Fn = sinc(n/2).*exp(-1i*n*pi/2);    % 傅里叶级数系数
Fn(N+1) = 0;        % 当n=0时，Fn=0
ft = zeros(1, length(t));   % 初始化傅里叶级数
% 得到傅里叶级数
for I = -N:N
    ft = ft +Fn(I+N+1).*exp(1i*2*pi*I*fs*t);
end
% 画图
figure;
plot(t, ft);
xlabel('t'),ylabel('f(t)'),title('周期方波傅里叶级数展开后的波形');