hc528.m

来自「这个是西安电子科技大学的陈怀琛教授编写的《数字信号处理教程——MATLAB释义与」· M 代码 · 共 13 行

%《数字信号处理教程——MATLAB释义与实现》第三章例5.2.8程序hc528
% z反变换求差分方程全解（有初始条件）
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
format long
a=[1,-1.5,0.5];b=1;       % 分子分母系数向量
n=[0:7]; x=(1/4).^n; 
Y=[4,10];X=0;             % 初始条件
xic=filtic(b,a,Y,X)       % 求初始条件和系数
% xic=[1,-2];
y1=filter(b,a,x,xic)        % 用filter函数求差分方程全解
y2=(1/3)*(1/4).^n+(1/2).^n+(2/3)*ones(1,8)  %MATLAB 验算