www.gusucode.com > MATLAB仿真技术与应用__的配套例题和习题的matlab源代码 > 第四章/example3.m

    %MATLAB program4-3
%Solve difference equation
clear all;clc;
b=[1 1 1];    
a=[1,-0.95,0.9025];
Y=[-2,-4];                   %输入初始条件
X=[1,2];
xic=filtic(b,a,Y,X)            %filter()函数选择初始条件
n=[0:50];
x=cos(n*pi/6);
y=filter(b,a,x,xic);
plot(n,y);