信号的Fouirer分解与合成MATLAB实现

用振幅为0.8的方波进行Fourier分析，并用分析得到的系数求解当k取不同值时的合成信号

程序如下;

clear all;
close all;%关闭所有打开的所有图形窗口
N=200;dt=4/N;%数据点数和采样间隔
for n=1:N%得到方波信号
    if(n*dt>=2)
        x(n)=0.8;
    else
        x(n)=-0.8;
    end
end
figure(1) %打开第一个图形窗口，绘制原始信号及其Fourier分析函数
subplot(211),plot((1:N)*dt,x),hold on;
plot((1:N)*dt,zeros(1,N),'k'),xlabel('时间/s')
title('原始信号');
a=zeros(1,N);b=zeros(1,N);
nn=floor(N/2)+1;
for k=0:nn-1
    a(k+1)=0;
    b(k+1)=0;
    for ii=0:N-1
        a(k+1)=a(k+1)+2/N*x(ii+1)*cos(2*pi*k*ii/N);%求解Fourier系数
        b(k+1)=b(k+1)+2/N*x(ii+1)*sin(2*pi*k*ii/N);
    end
    c(k+1)=sqrt(a(k+1).^2+b(k+1）.^2);
end
subplot(212),plot((0:nn-1)/(N*dt),c);title('Fourier变换')
xlab