clear allsyms z n p T syms z n p T h=0.1; A=[exp(-h) int(exp(-T),0,h);-0.5 1]; b=[0;0.5]; c=[exp(-p) int(exp(-T),0,p)]; G=c/(z*eye(2)-A)*b; R=z/(z-1); Y=R*G; y=iztrans(Y); t=1; for m=0:9, v(1,m+1)=subs(y,n,m); for l=0:9, w(1,l+t)=subs(v(1,m+1),p,l*1e-2); end t=t+10; end figure plot(w) grid on -- (たむむ) 2007-10-23 14:39:18