lin.m

这是一个dsp的matlab包

function test 
b=[1 2 1];
a=[1 -1.58 0.64];
n=50;
x=[1,ones(1,n-1)];
y=filter(b,a,x);
plot(y./70,'*');
hold on;
plot(x,'.');
grid;