某信号展开成傅里叶级数，再恢复

傅里叶级数

  

对比，1) 频域法没有直流分量，f=0的傅里叶变换幅值为0 ；2）an 是 hn，bn是gn。 

clc,clear all,close all
%% 正确
dFs = 50e3;
M = 80;
Fs = dFs*M*2 + 1;% 采样频率
T = 1/Fs;	% 采样时间
L = 80+1;	% 信号长度
t = (0:L-1)*T; % 时间
x = sin(2*pi*1e5*t ).*sin(2*pi*1e6*t );	%cos为底原始信号 + 
y = x ;%+ randn(size(t));	%添加噪声
figure;
plot(t,y)
Fy = fft(y,163)/163*2;
Fy2 = fftshift(abs(Fy));
figure
% plot(abs(Fy2));
An = Fy2;
% An(An<0.1) = 0;
ag = fftshift(angle(Fy));
yy = 0;
% plot(t,sin(Fs*t))
for i = -81:81
%     abc = sin(2*i*dFs*t) ;
%     plot(abc)
    temp = An(i+82)/2.*(cos(ag(i+82)+pi/2)*sin(2*i*pi*dFs*t) + sin(ag(i+82)+pi/2)*cos( 2*i*pi*dFs*t)); 
%     temp = An(i