q704.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第七章例7.4程序q704
% 求时域离散系统对常用序列的响应
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
%N: 输入数据长度 
%B: 差分方程系数
%A: 差分方程系数
clear;close all;
N=64;n=0:N-1;m=10;
B=[1,6,15,20,15,6,1]*0.0007378;	                    % 设定系统参数A,B
A=[1.0000,-3.1836,4.6223,-3.7795,1.8136,-0.4800,0.0544];
x1=[n==0]; 			                                % 产生输入信号x1(n)
y1=filter(B,A,x1); 	                                % 对x1(n)的响应
x2=[(n-m)==0];		                                % 产生输入信号x2(n)
y2=filter(B,A,x2); 	                                % 对x2(n)的响应
x3=[n>=0];				                            % 产生输入信号x3(n)
y3=filter(B,A,x3); 	                                % 对x3(n)的响应
x4=[(n>=0)&(n<32)];	                                % 产生输入信号x4(n)
y4=filter(B,A,x4); 	                                % 对x4(n)的响应
x5=exp(j*pi*n/8);		                            % 产生输入信号x5(n)
y5=filter(B,A,x5); 	                                % 对x5(n)的响应
subplot(3,2,1);stem(n,y1,'.');line([0,N],[0,0])
set(gcf,'color','w')                                % 设置图形背景色为白色
axis([0,N,min(y1),max(y1)]);ylabel('y1(n)]')
subplot(3,2,2);stem(n,y2,'.');line([0,N],[0,0])
axis([0,N,min(y2),max(y2)]);ylabel('y2(n)]')
subplot(3,2,3);stem(n,y3,'.');line([0,N],[0,0])
axis([0,N,min(y3),max(y3)]);ylabel('y3(n)]')
subplot(3,2,4);stem(n,y4,'.');line([0,N],[0,0])
axis([0,N,min(y4),max(y4)]);ylabel('y4(n)]')
subplot(3,2,5);stem(n,real(y5),'.');line([0,N],[0,0])
axis([0,N,-1,1]);ylabel('Re[y5(n)]')
subplot(3,2,6);stem(n,imag(y5),'.');line([0,N],[0,0])
axis([0,N,-1,1]);ylabel('Im[y5(n)]')