求问一道阶跃响应(用stepz来求)

Quandom

离散时间系统可由下列差分方程描述：
y[n]-y[n-1]+0.9y=x[n]
计算n=[-20:100]时系统的阶跃响应！
计算n=[-20:100]时系统的阶跃响应！

下面是我的做法
*****************************************************************************************************
h=impz(b,a,n);
subplot(211);stem(n,h,'.');
title('impulse response');xlabel('n');ylabel('h[n]');
a=[1,-1,0.9];b=1;n=[-20:100];
s=stepz(b,a,n);
subplot(212);
stem(n,s,'.');
title('step response');xlabel('n');ylabel('s[n]');
一下是我的疑惑，希望大神帮帮忙
******************************************************************************************************
1为什么这么写的代码，能够把第一个冲击响应求出来，而阶跃响应就求不出来了呢？
2经过我自己的调试，发现在s=stepz(b,a,n);执行后还是空的向量。这是为什么？
3我把s=stepz(b,a,n);写成s=stepz(b,a);貌似能求出s，但是在stem时会出错，它的长度为什么不再和n相等，为什么？
小弟不胜感激！