num=1.064;
den=[2 -3.685 1.791];
u1=[ones(1,50),-1*ones(1,50)];
u=[u1,u1,u1];
figure(1)
dlsim(num,den,u);
grid on;
title('離散系統仿真')
xlabel('時間')
ylabel('振幅')