q708.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第七章例7.8程序q708
% 求z的多项式分式的逆z变换
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear, close all
x=[2,4,3,5];nsx=-2;		    % 输入序列及初始时间
nfx=nsx+length(x)-1;		% 计算序列终止时间
Bw=-3; nsbw=-1;				% 系统函数的分子系数,及z的最高次数
Aw=[2,-2.2,0.5]; nsaw=0;	% 系统函数的分母系数,及z的最高次数
B=conv(-3,x);				% 输入与分子z变换的多项式乘积
A=Aw;						% 分母不变
nsy=nsaw-(nsbw-nsx)			% 分子比分母z的次数高nsy
[r,p,k]=residuez(B,A)		% 求留数r,极点p及直接项k		
nf=input('终点时间nf= ');	% 要求键入终点时间
n = min(nsx,nsy):nf;		% 生成总时间数组
% 求无限序列yi和直接序列yd
yi=(r(1)*p(1).^(n-nsy)+r(2)*p(2).^(n-nsy)).*stepseq(nsy,n(1),nf);
yd =k(1)*impseq(nsy,n(1),nf)+k(2)*impseq(-1-nsy,n(1),nf);
y=yi+yd;								% 合成输出
xe = zeros(1,length(n));				% 初始化,将x延拓为xe,
xe(find((n>=nsx)&(n<=nfx)==1))=x;		% 在对应的n处把xe赋值x
subplot(2,1,1),stem(n,xe,'.'),
line([min(n(1),0),nf],[0,0])            % 绘图
set(gcf,'color','w')                    % 设置图形背景色为白色
subplot(2,1,2),stem(n,y,'.'),
line([min(n(1),0),nf],[0,0])
subplot(2,1,1),stem(n,xe,'.'),
line([min(n(1),0),nf],[0,0])            % 绘图
subplot(2,1,2),stem(n,y,'.'),
line([min(n(1),0),nf],[0,0])