数字信号处理实验

x1=[1 1 1 1 0 0 0 0]; x2=[1 2 3 4 4 3 2 1]; x3=[4 3 2 1 1 2 3 4]; x4=cos*pi*n); N=8;n=0:7;k=0:7; X1k=fft(x1,N);

subplot(2,2,1);stem(n,x1,'.'); xlabel('n');ylabel('|x1(n)|'); subplot(2,2,2);stem(k,abs(X1k),'.'); xlabel('k');ylabel('|X1(k)|'); X2k=fft(x2,N);

subplot(2,2,3);stem(n,x2,'.'); xlabel('n');ylabel('|x2(n)|');

subplot(2,2,4);stem(k,abs(X2k),'.'); xlabel('k');ylabel('|X2(k)|');

1430.5|X1(k)||x1(n)|210024n680024k6843210|X2(k)||x2(n)|20151050024n68024k68

>> x1=[1 1 1 1 0 0 0 0]; x2=[1 2 3 4 4 3 2 1]; x3=[4 3 2 1 1 2 3 4]; x4=cos*pi*n); N=8;n=0:7;k=0:7; figure(2) X3k=fft(x3,N);

subplot(2,2,1);stem(n,x3,'.');

xlabel('n');ylabel('|x3(n)|'); subplot(2,2,2);stem(k,abs(X3k),'.'); xlabel('k');ylabel('|X3(k)|'); X2k=fft(x4,N);

subplot(2,2,3);stem(n,x4,'.'); xlabel('n');ylabel('|x4(n)|'); subplot(2,2,4);stem(k,abs(X2k),'.'); xlabel('k');ylabel('|X4(k)|');

43210|X3(k)||x3(n)|20151050024n68024k6810.50-0.5-1|X4(k)||x4(n)|43210024n68024k68

>> %ex3(2)main

fs=64;N=16; n=0:N-1;k=n;

x5=cos(n*pi/4)+cos(n*pi/8);

x6=cos(8*pi*n/fs)+cos(16*pi*n/fs)+cos(20*pi*n/fs); X5k=fft(x5,N); X6k=fft(x6,N); figure(3)

subplot(2,2,1);stem(n,x5,'.'); xlabel('n');ylabel('|x5(n)|'); subplot(2,2,2);stem(abs(X5k),'.'); xlabel('k');ylabel('|X5(k)|'); subplot(2,2,3);stem(n,x6,'.'); xlabel('n');ylabel('|x6(n)|'); subplot(2,2,4);stem(abs(X6k),'.');