利用filtic函数 差分方差的初始条件

1. %这里要利用filtic函数 为滤波器的直接II型实现选择初始条件  
2. %求解查分方程 y(n) - 0.4y(n-1) - 0.45y(n-2) = 0.45x(n) +0.4x(n-1) - x(n-2)  
3. %y(-1) = 0 y(-2) = 1 x(-1) = 1 x(-2） = 2  
4. %x(n) = 0.8 ^u(n) 状态方程H(z) = (0.45 + 0.4z^-1 - z^-2)/(1- 0.4z^-1-0.45z^-2)  
5. % x 为滤波前序列 y 为输出 序列  
6. % 总结一下 首先把num 和den 写出来  
7. num = [ 0.45 0.4 -1];  
8. den = [1 -0.4 -0.45]  
9. % 其次把初始条件写来  
10. x0 = [1 2] ;  
11. y0= [0 1];  
12. N = 50 ;  
13. n = [1 :N-1]';  
14. x = 0.8.^n;  
15. %生成初始条件  
16. Zi = filtic(num, den , y0 , x0);  
17. [y , Zf] = filter(num , den ,x, Zi);  
18. plot(n , x ,'R-', n, y, 'b--');  
19. xlabel('n'); ylabel('(n)--y(n)');  
20. legend('Input x' , 'Output y', 1);  
21. grid;