第一个:
x1n=[1,1,1,1];
stem(abs(x1n))
x1k8=fft(x1n,8);
x1k16=fft(x1n,16);
n=8;
f=2/n*(0:n-1);
figure(1);
subplot(1,3,2);stem(f,abs(x1k8),'r','.');
title('8');
n=16;
f=2/n*(0:n-1);
subplot(1,3,3);stem(f,abs(x1k16),'r','.');
title('16');
subplot(1,3,1);stem(f,abs(x1n),'r','.');
title('1');
第二个:
xa=[1,2,3,4,4,3,2,1];
stem(xa);
x1=fft(xa,8);
x2=fft(xa,16);
subplot(2,2,1);stem(xa);
title('1');
subplot(2,2,2);stem(abs(x1));
title('8');
subplot(2,2,3);stem(abs(x2));
title('16');
第三个:
xa=[4,3,2,1,1,2,3,4];
stem(xa);
x1=fft(xa,8);
x2=fft(xa,16);
subplot(2,2,1);stem(xa);
title('1');
subplot(2,2,2);stem(abs(x1));
title('8');
subplot(2,2,3);stem(abs(x2));
title('16');
第四个:
a=8;
b=0:a-1;
x1=cos(pi/4*b);
c=16;
d=0:c-1;
x2=cos(pi/4*d);
x1k8=fft(x1,8);
x1k16=fft(x1,16);
figure(1)
subplot(2,2,1);stem(x1);
title('1');
n=8;
f=8/n*(0:n-1);
subplot(2,2,2);stem(f,abs(x1k8));
title('8');
n=16;
f=16/n*(0:n-1);
subplot(2,2,4);stem(f,abs(x1k16));
title('16');
subplot(2,2,3);stem(x2);
title('16');
第五个:
a=8;
b=0:a-1;
x1=sin(pi/8*b);
c=16;
d=0:c-1;
x2=sin(pi/8*d);
x1k8=fft(x1,8);
x1k16=fft(x1,16);
figure(1)
subplot(2,2,1);stem(x1);
title('8原图');
n=8;
f=8/n*(0:n-1);
subplot(2,2,2);stem(f,abs(x1k8));
title('8');
n=16;
f=16/n*(0:n-1);
subplot(2,2,4);stem(f,abs(x1k16));
title('16');
subplot(2,2,3);stem(x2);
title('16原图');
第六个:
F=64;
T=1/F;
N=16;
n=0:N-1;
x=cos(pi*8*n*T)+cos(pi*16*n*T)+cos(pi*20*n*T);
xk16=fft(x,16);
xk32=fft(x,32);
xk64=fft(x,64);
N=32;
n=0:N-1;
x32=cos(pi*8*n*T)+cos(pi*16*n*T)+cos(pi*20*n*T);
N=64;
n=0:N-1;
x64=cos(pi*8*n*T)+cos(pi*16*n*T)+cos(pi*20*n*T);
figure(1)
subplot(2,3,1);stem(x);title('原图');
N=32;
n=0:N-1;
x=cos(pi*8*n*T)+cos(pi*16*n*T)+cos(pi*20*n*T);
xk32=fft(x,32);
subplot(2,3,2);stem(x32);title('32原图');
N=64;
n=0:N-1;
x=cos(pi*8*n*T)+cos(pi*16*n*T)+cos(pi*20*n*T);
xk64=fft(x,64);
subplot(2,3,3);stem(x64);title('64原图');
N=16;
f=2/N*(0:N-1);
subplot(2,3,4);stem(f,abs(xk16));title('16');
N=32;
f=2/N*(0:N-1);
subplot(2,3,5);stem(f,abs(xk32));title('32');
N=64;
f=2/N*(0:N-1);
subplot(2,3,6);stem(f,abs(xk64));title('64');