% Fourier.m
% 输入变量：t--时间区间
%          f--输入函数在相应时间区间的取值
%          T--输入函数周期
% 输出变量：F_coefficient--输入函数前四项傅里叶展开系数

function F_coefficient=Fourier(f,T,t)
    omg0 = 2*pi./T; % 基波频率
    N = length(t);
    k1=-8;
    k2=8;
    k=(k1:k2)';
    F=1/N*exp(-1i*kron(k*omg0,t.'))*f;%求指数形式的傅里叶级数系数
    F_coefficient=F(10:13)';
end