hc531.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第四章例5.3.1程序hc531
% 用逆z变换求离散系统的输出
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=[2,4,3,5];		% 输入序列及初始时间
nfx=length(x)-1;		% 计算序列终止时间
b=-3;a=[2,-2.2,0.5]; 		% 系统函数的分子分母系数,
B=conv(-3,x);A=a;			% 输入与分子z变换的多项式乘积,分母不变
[r,p,k]=residuez(B,A)		% 求留数r,极点p及直接项k		
nf=input('终点时间nf= ');	% 要求键入终点时间
n = 0:nf;		            % 生成总时间数组
% 求无限序列yi和直接序列yd
yi=(r(1)*p(1).^n+r(2)*p(2).^n).*stepseq(0,0,nf);
yd =k(1)*impseq(0,0,nf)+k(2)*impseq(1,0,nf);
y=yi+yd;										 	% 合成输出
subplot(2,1,1),
stem(n,[x,zeros(1,nf-nfx)],'.')             % 画输入序列
xlabel('n'),ylabel('x(n)')
subplot(2,1,2),stem(n,y,'.')                % 画输出序列
xlabel('n'),ylabel('y(n)')
set(gcf,'color','w')                        % 置图形背景色为白