实验三用FFT做谱分析.doc

（1） x(n)={1 0<=n<=5 构造dft 函数计算的10点和20点的
0 其他 dft并画出图形。
function[xk]=dft(xn,N);
n=0:N-1;
k=0:N-1;
w=exp(-j.*2.*pi./N);
nk=n.*k;
wnk=w.^nk;
xk=xn.*wnk;

N1=10;
N2=20;
n1=0:N1-1;n2=0:N2-1;
xn=ones(1,6);
xn1=[xn zeros(1,4)];
xn2=[xn zeros(1,14)];
[xk1]=dft(xn1,N1);
[xk2]=dft(xn2,N2);
subplot(2,1,1);
stem(n1,xk1);
subplot(2,1,2);
stem(n2,xk2);
（2）a：n=0:7;x=[ones(1,4),zeros(1,4)];
subplot(2,2,1);stem(n,x);xlabel('N=8');ylabel('x(n)');
y=fft(x);q=n/8;
subplot(2,2,2);stem(q,abs(y));ylabel('fft N=8');
n=0:15;x=[ones(1,4),zeros(1,12)];
subplot(2,2,3);stem(n,x);xlabel('N=16');ylabel('x(n)');
q=n/16;y=fft(x)
subplot(2,2,4);stem(q,abs(y));ylabel('fft N=16');
b，c在一起：N=8;
n=0:N-1;
q=n/N;
x=cos(pi.*n./4);
y=fft(x,N);
subplot(2,2,1);stem(q,abs(y));grid on;ylabel('fft N=8');
N=16;
n=0:N-1;
q=n/N;
x=cos(pi.*n./4);
y=fft(x,N);
subplot(2,2,2);stem(q,abs(y));grid on;ylabel('fft N=16');
N=8;
n=0:N-1;
q=n/N;
x=sin(pi.*n./8);
y=fft(x,N);
subplot(2,2,3);stem(q,abs(y));grid on;ylabel('fft N=8');
N=16;
n=0:N-1;
q=n/N;
x=sin(pi.*n./8);
y=fft(x,N);
subplot(2,2,4);stem(q,abs(y));grid on;ylabel('fft N=16');

（3）：f1=2,f2=2.05;fs=10;
N=128;
n=0:N-1;
q=n/N;
x=sin(2.*pi.*n.*f1/fs)+sin(2.*pi.*f2.*n/fs);
y=fft(x,N);
subplot(3,1,1);stem(q,abs(y));grid on;ylabel('fft N=128');

f1=2,f2=2.05;fs=10;
N=512;
n=0:N-1;
q=n/N;
x=sin(2.*pi.*n.*f1/fs)+sin(2.*pi.*f2.*n/fs);
y=fft(x,N);
subplot(3,1,2);stem(q,abs(y));grid on;ylabel('fft N=512');

f1=2,f2=2.05;fs=10;
N=512;
n=0:N-1;
q=n/N;
n1=0:128;
x=sin(2.*pi.*n1.*f1/fs)+sin(2.*pi.*f2.*n1/fs);
y=fft(x,N);
subplot(3,1,3);stem(q,abs(y));grid on;ylabel('fft N=128-128');