function ft(N)
%x(n)為輸入信號;N為所取點的個數
y=zeros(1,N);
for n=1:N
    x(n)=cos(0.48*pi*n)+cos(0.52*pi*n);
end
for k=1:N
    for m=1:N
        y(k)=y(k)+x(m)*exp(-j*2*pi*(k-1)*(m-1)/N);
        z(k)=norm(y(k));
    end
end
plot(z)