Giáo trình matlab v5.1 P10

Trước đây MatLab không phân biệt chữ in, chữ thường (giống như Fortran). Các phiên hơn gần đây lại có sự phân biệt này (theo ngôn ngữ C). Các từ khóa đều viết chữ thường.
Nội dung trích xuất từ tài liệu:
Giáo trình matlab v5.1 P10Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng CAÙC BAØI TAÄPBaøi 1: Trích töø Ví duï 3.13 saùch … taùc giaû Nuyeãn Vaên Giaùp%Ket NOI 2 HE thong SONG SONGa=[1 2 3;4 5 6;7 8 9];b=[3 4;4 5;7 9];c=[0 0 1];d=[0 0];e=[1 9 3;4 5 6;7 8 7];f=[2 4;4 6;7 9];g=[0 1 1];h=[0 0];[A,B,C,D]= parallel(a,b,c,d,e,f,g,h)Keát quaû:A= 1 2 3 0 0 0 4 5 6 0 0 0 7 8 9 0 0 0 0 0 0 1 9 3 0 0 0 4 5 6 0 0 0 7 8 7B= 3 4 0 0 4 5 0 0 7 9 0 0 0 0 2 4 0 0 4 6 0 0 7 9Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 27 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoängC= 0 0 1 0 0 0 0 0 0 0 1 1D= 0 0 0 0 0 0 0 0A= 1 2 3 0 0 0 4 5 6 0 0 0 7 8 9 0 0 0 0 0 0 1 9 3 0 0 0 4 5 6 0 0 0 7 8 7B= 3 4 4 5 7 9 2 4 4 6 7 9Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 28 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoängC= 0 0 1 0 1 1D= 0 0Baøi 2: Keát noái hai haøm truyeàn noái vôùi soá lieäu nhaäp töø baøn phím (vieát chöông trình trong m_file)%Bai tap tong quat ket noi 2 he thong noi tiep%Cu phap SYS=series(SYS1,SYS2,OUTPUTS1,INPUTS2)%Vi du ta se ket noi 2 ham truyennum1=input( Nhap num1= );den1=input(Nhap den1= );num2=input(Nhap num2= );den2=input(Nhap den2= );[num,den]=series(num1,den1,num2,den2)Baøi 3: (Trích trang 11-14 saùch Control System Toolbox)sys1=tf(1,[1 0])Transfer function:1-ssys2=ss(1,2,3,4)a= x1 x1 1b= u1 x1 2Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 29 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoängc= x1 y1 3d= u1 y1 4Continuous-time model.sys=append(sys1,10,sys2)a= x1 x2 x1 0 0 x2 0 1b= u1 u2 u3 x1 1 0 0 x2 0 0 2c= x1 x2 y1 1 0 y2 0 0 y3 0 3Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 30 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoängd= u1 u2 u3 y1 0 0 0 y2 0 10 0 y3 0 0 4Continuous-time model.Baøi 4: moät heä thoáng bieåu dieån nhö hình sau vôùiG0(s)=1;G1(s)=1/(s+1);G2(s)=1/(s+2);G3(s)=1/(s+3);H1(s)=4;H2(s)=8;H3(s)=12. Y(s)R(s) + G0(s) + G1(s) G2(s) + G3(s) 1 2 3 4 H2(s) H3(s) 6 7 H1(s) 5 n1=1; d1=1; n2=1; d2=[1 1]; n3=1; d3=[1 2]; n4=1; d4=[1 3]; n5=4; d5=1; n6=8; d6=1; n7=12; d7=1; nblocks=7; blkbuild q=[1 0 0 0 0Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 31 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng 2 1 -5 0 0 3 2 -6 0 0 4 2 -6 3 -7 53000 63000 7 4 0 0 0]; input=1; output=4; [aa,bb,cc,dd]=connect(a,b,c,d,q,input,output); [num,den]=ss2tf(aa,bb,cc,dd); printsys(num,den,s)Giaûi thích: Ta phaûi ñaùnh soá trong moãi heä thoáng phuï nhö hình treân. Baûy caâu leänh ñaàu tieân bieåu dieån haømtruyeàn cuûa baûy khoái, qui ñònh teân töông öùng vôùi töû vaø maãu laø n1,d1,n2,d2,...trong tröôøng hôïp neáucho daïng laø kieåu bieán traïng thaùi trong töøng heä thoáng phuï thì teân cuûa chuùng töông öùng laøa1,b1,c1,d1,a2,b2,c2,d2,... Ñaët bieán nblock=7 (baèng vôùi soá cuûa heä thoáng phuï). Sau ñoù laø leänh blkbuild duøng nhöõng bieán cuûa nblock ñeå baét ñaàu xaây döïng heä thoáng. Bieánblkbuild chuyeån ñoåi taát caû caùch thöùc dieån taû haøm truyeàn cuûa töøng heä thoáng phuï thaønh kieåu bieántraïng thaùi nhö duøng leänh tf2ss vaø ñöa chuùng vaøo moät khoái lôùn cuûa ma traän trang thaùi goïi laø a, b,c, d. Taïo ra ma traän q ñeå nhaän bieát moái lieân heä giöõa caùc heä thoáng phuï. (Moãi haøng cuûa ma traän qtöông öùng vôùi moät heä thoáng phuï khaùc nha ...