Giáo trình matlab v5.1 P8

Năm 2000 MATLAB 6 cho đổi mới môi trường làm việc MATLAB, thêm LAPACK và FFTW (Fastest Fourier Transform in the West - "biến đổi Fourier nhanh nhất của phương Tây").Năm 2002 MATLAB 6.5 phát hành đã cải thiện tốc độ tính toán, sử dụng phương pháp dịch JIT (Just in Time) và tái hỗ trợ MAC.Năm 2004 MATLAB 7 phát hành, có khả năng chính xác đơn và kiểu nguyên, hỗ trợ hàm lồng nhau, công cụ vẽ điểm, và có môi trường phân tích số liệu tương tác....
Nội dung trích xuất từ tài liệu:
Giáo trình matlab v5.1 P8Khaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng d = [0 0 ; 0 0]; Gc = gram(a,b) Ta nhaän ñöôïc ma traän: Gc = 1.2016 -0.0318 -0.0318 1.0708 Tìm haïng ma traän baèng leänh: r = rank(Gc) ta ñöôïc r = 2 vaø baèng kích thöôùc cuûa ma traän ñaùnh giaù. Vaäy heä thoáng naøy coù theå ñieàukhieån ñöôïc.6. Leänh DSORT, ESORTa) Coâng duïng: Saép xeáp caùc giaù trò rieâng theo thöù töï phaàn thöïc hoaëc bieân ñoä soá phöùc.b) Cuù phaùp: s = dsort(p) [s,ndx] = dsort(p) s = esort(p) [s,ndx] = esort(p)c) Giaûi thích: s = esort(p) xeáp caùc giaù trò rieâng phöùc trong vector p theo thöù töï giaûm daàn cuûa phaàn thöïc.Ñoái vôùi caùc giaù trò rieâng lieân tuïc, caùc giaù trò rieâng khoâng oån ñònh xuaát hieän tröôùc. s = dsort(p) xeáp caùc gí trò rieâng phöùc trong vector p theo thöù töï giaûm daàn cuûa bieân ñoä. Ñoáivôùi caùc giaù trò rieâng giaùn ñoaïn, caù giaù trò rieâng khoâng oån ñònh xuaát hieän tröôùc. [s,ndx] = dsort(p) hay [s,ndx] = esort(p) cuõng taïo ra vector ndx chöùa caùc chæ soá duøng theothöù töï.d) Ví duï: Xeáp caùc phaàn töû cuûa vector p = [2+3j -3+j 1-9j 3-7j 5+2j 6-j] theo thöù töï giaûm daàncuûa phaân thöïc vaø ñoä lôùn soá phöùc. p = [2+3j -3+j 1-9j 3-7j 5+2j 6-j] % Xeáp theo thöù töï giaûm daàn cuûa ñoä lôùn soá phöùc: s = dsort(h) s= 1.0000 + 9.0000j 3.0000 + 7.0000j 6.0000 + 1.0000j 5.0000 – 2.0000j 2.0000 + 3.0000j -3.0000 + 1.0000j % Xeáp theo thöù töï giaûm daàn cuûa phaàn thöïc: s’ = esort(h) 6.0000 + 1.0000j 5.0000 – 2.0000j 3.0000 + 7.0000jThöïc hieän: PHAÏM QUOÁC TRÖÔØNG -7- GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng 2.0000 – 3.0000j 1.0000 + 9.0000j -3.0000 – 1.0000j7. Leänh EIGa) Coâng duïng: Tìm caùc giaù trò rieâng vaø caùc vector rieâng cuûa heä thoáng.b) Cuù phaùp: E = eig(X) [V,D] = eig(X) [V,D] = eig(X) [V,D] = eig(X,’nobalance’) E = eig(A,B) [V,D] = eig(A,B)c) Giaûi thích: E = eig(X) laø moät vector chöùa caùc giaù trò rieâng cuûa ma traän vuoâng X. [V,D] = eig(X) taïo ra moät ma traän ñöôøng cheùo D cuûa caùc giaù trò rieâng vaø ma traän ñuû vôùicaùc coät laø caùc vector rieâng töông öùng ñeå cho X*V = V*D. [V,D] = eig(X,’nobalance’) gioáng nhö [V,D] = eig(X) nhöng boû qua söï caân baèng. Caùchnaøy ñoâi khi cho keát quaû chính xaùc hôn. E = eig(A,B) laø vector chöùa caùc giaù trò rieâng phoå bieán cuûa caùc ma traän vuoâng A vaø B. [V,D] = eig(A,B) taïo ra ma traän ñöôøng cheùo D cuûa caùc giaù trò rieâng phoå bieán vaø caùc matraän ñuû V vôùi caùc coät laø caùc vector rieâng töông öùng ñeå cho A*V = B*V*D.d) Ví duï: Cho X = [2+3j -3+j 1-9j ; 3-7j 5+2j 6-j ; 0+7j 6-8j 2+5j]. tìm caùc giaù trò rieângcuûa X. X = [2+3j -3+j 1-9j ; 3-7j 5+2j 6-j ; 0+7j 6-8j 2+5j]; [V,D] = eig(X) V= 0.4158 + 0.3442j 0.5455 + 0.4929j 0.4344 – 0.2255j -0.3275 + 0.3580j 0.1837 – 0.2659j 0.5974 + 0.1368j 0.1209 – 0.6772j -0.5243 + 0.2831j 0.4954 + 0.3734j D= -9.3743 + 4.7955j 0 0 0 9.2099 + 0.2831j 0 0 0 9.1644 – 2.2542j8. Leänh PRINTSYSa) Coâng duïng: In ra caùc tham soá cuûa heä thoáng tuyeán tínhb) Cuù phaùp: printsys(a,b,c,d) printsys(a,b,c,d,ulabels,ylabels,xlabels) printsys(num,den,‘s’)Thöïc hieän: PHAÏM QUOÁC TRÖÔØNG -8- GVHD: PHAÏM QUANG HUYKhaûo saùt öùng duïng MATLAB trong ñieàu khieån töï ñoäng printsys(num,den,‘z’)c) Giaûi thích: printsys in caùc tham soá cuûa heä khoâng gian traïng thaùi vaø haøm truyeàn theo daïng ñaëc bieät.Ñoái vôùi heä khoâng gian traïng thaùi, caùc ngoõ vaøo, ngoõ ra vaø traïng thaùicuûa heä ñöôïc ñaët teân vaøhaøm truyeàn ñöôïc hieån thò döôùi daïng tyû soá cuûa hai ña thöùc. printsys(a,b,c,d) in ra heä khoâng gian traïng thaùi (a,b,c,d) vôùi teân tham soá ôû phía treân vaøphía beân traùi cuûa ma traän heä thoáng. printsys(a,b,c,d,ulabels,ylabels,xlabels) in ra heä khoâng gian traïng thaùi vôùi teân tham soáñöôïc chæ ñònh bôûi caùc vector ulabels, ylabels vaø xlabels. ulabels, ylabels vaø xlabels chöùa teânngoõ vaøo, ngoõ ra vaø traïng thaùi cuûa heä thoáng. printsys(num,den,‘s’) hoaëc printsys(num,d ...