Lecture
load dryer2; ts=0.08; dan=iddata(у2,u2,ts);
dan.outputn='temperature'; dan.inputn='gas flow rate'; dan.inputUnit='m3/h'; dan.outputUnit='deg.C 100'; get(dan);
plot(dan);
zdan=dtrend(dan);
plot(zdan);
cra(zdan);
[ir,R,cl]=cra(zdan,[],[],2);
[g,phiv]=spa(dan,[],[]);
bodeplot([g]);
bodeplot([g],’sd’,3,'fill'); bodeplot([phiv], ‘sd',3,'fill’); zdanv=zdan(1:500); zdane=zdan(501:1000); NN=struc(l:10,1:10,1); v=arxstruc(zdane,zdanv,NN); [nn,vmod]=selstruc(v,'plot'); darx=arx(zdanv,пп);
darmax = armax(zdanv,[2 2 2 1]); zoe=oe(zdanv [2 2 1]); zbj=bj(zdanv,[2 2 2 2 1]); zpem~em(zdanv); zn4s=n4sid(zdanv,[1:10],[1:10]); [A,В]=th2arx(dагх); [g,рhiv]=th2ff(dаrx); [А,В,С,D,К,lan,T]=th2роlу(zpem); [А,В,С,D,К,х0]=th2ss(darmax); [num,dеn]=th2tf(zn4s); zzn4s=tf(num,den,0.08); sn4s=th2thc(zn4s); [num,den]=th2tf(sn4s); sysn4s=tf(num„den); [zepo,k]=th2zp(zn4s);
168
[zero,polus]=getzp(zepo);
zpplot(zpform(zepo)); compare(zdane,zn4s,zpem,zoe,zbj,dаrx,darmax; [е,r]=resid(zdan dаrx);
resid(r); [еl,r1]=resid(zdan,zbj); resid(rl); step(zzn4s,sysn4s); impulse(zzn4s,sysn4s); k=dcgain(sysn4s); bode (zzn4s,sysn4s);
[Gm,Pm,Wcg,Wcp]=margin(sysn4s); [Gm1,Рm1,Wcg1,Wcp]=margin(zzn4s); Gmlog=20*log10(Gm); Gmlog=20*log10(Gm1;) nyquist(zzn4s,sysn4s) ;
plot(e); [А,В,С,D]=ssdata(sn4s); Mу=ctrb(A,В); n=rank(Mу); Mн=obsv(А,С); n=rank(Mн)
Comments