[MATLAB信号处理] 程序求助

[langyun]

K=[1;7;9];
numc=1;
denc=conv([1 0],[0.67 1]);
T=0.1;
t=0:T:100;
for i=1:length(K)
[numd,dend]=c2dm(numc,denc,T,'zoh');
[num1,den1]=cloop(numd*K(i),dend,-1);
[yd,xd,td]=dstep(num1,den1,t);
subplot(1,2,1),plot(td,yd);title('离散系统');axis([0 100 -10 10]);hold on
[num,den]=cloop(numc*K(i),denc,-1);
[y,x,t1]= step(num,den,t);
subplot(1,2,2),plot(t1,y);title('连续系统');axis([0 20 0 1.7]);hold on
end
这是我写的一段程序，可是那个离散系统的就是出不来，请大家帮个忙，麻烦诸位了