用MATLAB语言辅助《自动控制原理》的教学.PDF

　1999 5 () May, 1999　 　22 　3 Jo rnal of Sich an Normal University(Nat ral Science) Vol.22, No.3　 用MATLAB 语言辅助《自动控制原理》的教学 赵仕良 (, 610066) 　　:MATLAB 《》, 、、、. :MATLAB ;;;; :TP312MA 　　:A 　　:1001-8395(1999)03-0341-03 , Ba- [1] sic、Fortran C (CACSD) .Matlab 1984 ,、、、、、 “”.Matlab ,. Matlab ,《》 .: 1 　 ., [2] , Ro th J ry .Matlab , 2 .K (s )=1/s , G(s )=1/s (s +1)(s +2), K (s )G (s )=1/s2 ＊1/s (s +1)(s +2). ,. %MATLAB PROGRAM-1 nc=[zeros(1, dimdgk-dimngk), ngk]; f nction[ngk, dgk, n m, den]=cltf(ng, dg, nk, dk) den=nc+dgk ngk=conv(ng, nk);dgk=conv(dg, dk); n m=ngk dimngk=length(ngk);dimdgk=length(dgk); 　　cltf ,., . %MATLAB PROGRAM -2 i=find(real(p)0);n1=length(i); ng=1;dg=poly([0,-1,-2]);nk=1;dk=[1,0,0]; j=find(real(z)0);n2=length(j); [ngk, dgk, n m, den]=cltf (ng, dg, nk, dk); if n10, disp(`[System is nstable, and nstable poles printsys(n m, den) have]);disp(n1); [z, p]=tf2zp(n m, den) elseif n20, disp(`System is a Nominimal Phase One); - - 　　:1999 03 22 　　:, , 27 , 　　 () 22 342 else disp(`System is a Minimal Phase One); disp(`please enter any key ) end pa se 　　,PROGRAM-2 (ng 、dg 、nk 、dk ), , . 2 　 [2] k ()0 +∞, . K (s )=k , , 1 . %MATLAB PROGRAM-3 pa se ng=1;dg=poly([0,-1,-2]); clpole=rloc s(ng, dg,[0.3:0.1:7]) rloc s(ng, dg), grid range=[0.3:0.1:7];[range, clpole] title(`rloc s