基于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)/');