Giáo trình matlab v5.1 P15

Vẽ trường vec tơ: Cũng như đồ thị lớp màu, việc hiển thị trường vec tơ rất cần thiết trong các ngành khoa học - vật lí. Để vẽ trường véc-tơ hai chiều của các ma trận u và v, dùng lệnh:
Nội dung trích xuất từ tài liệu:
Giáo trình matlab v5.1 P15Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoänghaøng cuûa ma traän u töông öùng vôùi moät thôøi ñieåm môùi. Neáu duøng theâm ñoái soá x0 ôû veáphaûi thì leänh lsim(a,b,c,d,u,x0) seõ chæ ra ñieàu kieän ban ñaàu cuûa caùc traïng thaùi. lsim(num,den,u) veõ ra ñaùp öùng thôøi gian cuûa haøm truyeàn ña thöùc: G(z) = num(z)/den(z) trong ñoù num vaø den chöùa caùc heä soá ña thöùc theo chieàu giaûm daàn soá muõ cuûa s. Neáu giöõ laïi caùc ñoái soá ôû veá traùi thì: [y,c] = dlsim(a,b,c,d,u) [y,c] = dlsim(a,b,c,d,u,x0) [y,c] = dlsim(num,den,u) seõ khoâng veõ ra caùc ñoà thò ñaùp öùng maø taïo ra caùc ma traän y vaø x, trong ñoù matraän y laø ñaùp öùng ngoõ ra vaø ma traän x laø ñaùp öùng traïng thaùi cuûa heä thoáng. Ma traän y coùsoá coät baèng soá ngoõ ra y vaø moãi haøng öùng vôùi moät haøng cuûa ma traän u. Ma traän x coù soácoät baèng soá traïng thaùi x vaø moãi haøng öùng vôùi moät haøng cuûa ma traän u.d) Ví duï: Moâ phoûng ñaùp öùng cuûa heä thoáng giaùn ñoaïn coù haøm truyeàn: 2 z 2 − 3.4 z + 1.5 H ( z) = 2 z − 1.6 z + 0.8 vôùi 100 maãu cuûa nhieãu ngaãu nhieân. num = [2 -3.4 1.5]; den = [1 -1.6 0.8]; rand(‘nomal’) u = rand(100,1); dlsim(num,den,u) title(‘Dap ung nhieu’) vaø ta ñöôïc ñoà thò ñaùp öùng cuûa heä nhö sau:Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 45 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng7. Leänh STEPa) Coâng duïng: Tìm ñaùp öùng naác ñôn vò.b) Cuù phaùp: [y,x,t] = step(a,b,c,d) [y,x,t] = step(a,b,c,d,iu) [y,x,t] = step(a,b,c,d,iu,t) [y,x,t] = step(num,den) [y,x,t] = step(num,den,t)c) Giaûi thích: Leänh step tìm ñaùp öùng naác ñôn vò cuûa heä tuyeán tính lieân tuïc. Neáu boû qua caùc ñoái soá ôû veá traùi cuûa doøng leänh thì leänh step veõ ra ñaùp öùng naáctreân maøn hình. step(a,b,c,d) veõ ra chuoãi ñoà thò ñaùp öùng naác, moãi ñoà thò töông öùng vôùi moái quanheä giöõa moät ngoõ vaøo vaø moät ngoõ ra cuûa heä lieân tuïc LTI: . x = Ax + Bu y = Cx + Du vôùi vector thôøi gian ñöôïc xaùc ñònh töï ñoäng. step(a,b,c,d,iu) veõ ra ñoà thò ñaùp öùng naác töø moät ngoõ vaào duy nhaát tôùi taát caû caùcngoõ ra cuûa heä thoáng vôùi vector thôøi gian ñöôïc xaùc ñònh töï ñoäng. Ñaïi löôïng voâ höôùng iuThöïc hieän: PHAÏM QUOÁC TRÖÔØNG - 46 - GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoänglaø chæ soá ngoõ vaøo cuûa heä thoáng vaø noù chæ ra ngoõ vaøo naøo ñöôïc söû duïng cho ñaùp öùngxung. step(num,den) veõ ra ñoà thò ñaùp öùng naác cuûa haøm truyeàn ña thöùc: G(s) =num(s)/den(s) trong ñoù num vaø den chöùa caùc heä soá ña thöùc theo chieàu giaûm daàn soá muõ cuûa s. step(a,b,c,d,iu,t) hay step(num,den,t) cuõng veõ ra ñaùp öùng naác cuûa heä khoânggian traïng thaùi hay haøm truyeàn vôùi vector thôøi gian t do ngöôøi söû duïng xaùc ñònh.Vector t chæ ra nhöõng thôøi ñieåm maø taïi ñoù ñaùp öùng naác ñöôïc tính vaø vector t phaûi ñöôïcchia thaønh nhöõng ñoaïn ñeàu nhau. Neáu giöõ laïi caùc ñoái soá ôû veá traùi cuûa doøng leänh thì: [y,x,t] = step(a,b,c,d) [y,x,t] = step(a,b,c,d,iu) [y,x,t] = step(a,b,c,d,iu,t) [y,x,t] = step(num,den) [y,x,t] = step(num,den,t) khoâng veõ ra caùc ñoà thò ñaùp öùng maø taïo ra caùc ma traän ñaùp öùng ngoõ ra y vaø matraän ñaùp öùng traïng thaùi x cuûaheä thoáng ñöôïc xaùc ñònh taïi nhöõng thôøi ñieåm t. Ma traän ycoù soá coät baèng soá ngoõ ra vaø moãi haøng öùng vôùi moät thaønh phaàn trong vector t. Ma traänx coù soá coät baèng soá traïng thaùi vaø moãi haøng öùng vôùi moät thaønh phaàn trong vector t.d) Ví duï: Veõ ñoà thò ñaùp öùng naác cuûa heä khoâng gian traïng thaùi baäc 2 sau: ⎡ . ⎤ ⎡− 0.5572 − 0.7814 ⎡ x ⎤ ⎡1⎤ ⎤1 ⎢ x1 ⎥ = ⎢ +u 0 ⎥ ⎢x2 ⎥ ⎢0⎥ . ⎢x2 ⎥ ⎣ 0.7814 ⎦⎣ ⎦ ⎣ ⎦ ⎣⎦ ⎡x ⎤ y = [1.9691 6.4493] ⎢ 1 ⎥ + [0] u ⎣ x2 ⎦ a = [-0.5572 -0.7814 ; 0.7814 0]; b = [1 ; 0]; c = [1. ...