数字信号处理实验例c2.txt

基于matlab的数字信号处理经典实验代码及所有完整实验报告

% 例C二
%x1(n)=R4(n)
%x2(n)=cos(pi/4*n)+cos(pi/8*n)
%X1:存放x1(n)的向量,Y1:存放X1(K)
%X2:存放x2(n)的向量,Y2:存放X2(K) 
%
 x1=[1,1,1,1,0,0,0,0]; 
n=0:15;
x2=cos(pi*n/4)+cos(pi*n/8);
i=0:7;
subplot(3,2,1);
stem(i,x1,'.');
axis([0 7 0 1]);
xlabel('n');
ylabel('x1(n)');

y1=fft(1,8);
subplot(3,2,3);
stem(i,abs(y1),'.');
xlabel('(N=8,wk=2pik/N)k');
ylabel('/X1(k)/');
y1=fft(x1,64);
i=0:63;
subplot(3,2,5);
stem(i,abs(y1),'.');
axis([0 63 0 4]);
xlabel('(N=64 wk=2pik/N)k');
ylabel('/X1(k)/');

y2=fft(x2);
figure;
subplot(2,2,1);
stem(n,x2);
title('x2(n)的时域序列');
xlabel('n');
ylabel('x2(n)');
subplot(2,2,3);
stem(n,abs(y2));
title('x2(n)的幅频特性');
xlabel('(N=16,wk=2pik/N)k');
ylabel('/X2(k)/');