Bài toán nhận dạng trong lý thuyết điều khiển

TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI KHOA TOÁN HỒNG THỊ THƠ BÀI TỐN NHẬN DẠNG TRONG LÝ THUYẾT ĐIỀU KHIỂN KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC Chuyên ngành: Toán ứng dụng Người hướng dẫn khoa học TS HÀ BÌNH MINH HÀ NỘI - 2013 Lời cảm ơn Em xin bày tỏ lòng biết ơn sâu sắc đến tiến sĩ HÀ BÌNH MINH, người tận tình giúp đỡ bảo cung cấp cho em kiến thức tảng để em hồn thành khóa luận Thầy người giúp em ngày tiếp cận có niềm say mê khoa học suốt thời gian làm việc Thầy Em xin bày tỏ tới anh PHẠM VĂN DUẨN, người nhiệt tình giúp đỡ bảo hướng dẫn em q trình gõ Tex hồn thành khóa luận Anh người cung cấp thêm tư liệu kiến thức giúp em giải đáp điều chưa hiểu băn khoăn Em xin bày tỏ lòng biết ơn tới thầy, cơng tác Khoa Tốn Trường Đại học sư phạm Hà Nội Thầy cô trực tiếp giảng dạy, truyền đạt cho em kiến thức quý báu chuyên môn kinh nghiệm nghiên cứu khoa học thời gian qua Cuối cùng, em xin chân thành gửi lời cảm ơn đến người thân gia đình, bạn bè giúp đỡ, động viên tạo điều kiện cho em suốt trình học tập hoàn thiện luận văn Lời cam đoan Tên tơi là: Hồng Thị Thơ, sinh viên Đại học khóa 2009 – 2013 lớp K35CN Toán Khoa Toán – Trường Đại học Sư phạm Hà Nội Tôi xin cam đoan đề tài: “ Bài toán nhận dạng lý thuyết điều khiển”, kết nghiên cứu thu thập riêng Các luận cứ, kết thu đề tài trung thực, không trùng với tác giả khác Nếu có khơng trung thực luận văn tơi xin hồn tồn chịu trách nhiệm trước hội đồng khoa học Mục lục Mở đầu i Nội dung iii Chương 1: Hệ động lực tuyến tính 1.1 Hệ động lực tuyến tính 1.2 Hàm Truyền 1.2.1 Phép biến đổi Laplace .5 1.2.2 Một số phép toán với ma trận hàm truyền 1.3 Các tham số Markov hàm truyền 1.3.1 Biểu diễn không gian trạng thái hàm truyền 1.3.2 Biểu diễn điều khiển quan sát 10 1.4 Dạng biểu diễn tối thiểu hệ 12 1.4.1 Tính chất .13 Chương 2: Xác định biểu diễn tối thiểu từ tham số Markov 19 2.1 Một số đặc tính ma trận Hankel tham số Markov 19 2.2 Thuật toán SVD cho biểu diễn tối thiểu .24 2.2.1 Thuật toán SVD .24 2.2.2 Ví dụ .25 2.3 Thuật toán SVD sửa đổi cho biểu diễn tối thiểu 29 2.3.1 Ví dụ .30 Tài liệu tham khảo 36 Mở đầu Lý chọn đề tài Điều khiển tốn có ý nghĩa ứng dụng quan trọng đời sống, đặc biệt lĩnh vực điện tử, viễn thông xử lý tín hiệu nói riêng Ta thường xây dựng mơ hình tốn học từ q trình vật lý Có nhiều vấn đề cần nghiên cứu lĩnh vực điều khiển Một số vấn đề có tính chất kinh điển tốn điều khiển Nó có ứng dụng rộng rãi ngành tốn ứng dụng, nên từ trước đến nay, ln đề tài mà nhà khoa học quan tâm nghiên cứu Để hiểu rõ toán em chọn đề tài “Bài toán nhận dạng lý thuyết điều khiển” để làm đề tài nghiên cứu cho khóa luận tốt nghiệp Khái quát nội dung phạm vi nghiên cứu Bài tốn điều khiển tuyến tính phần tảng quan trọng lý thuyết điều khiển nói chung phát triển khái niệm điều khiển nâng cao có gợi ý tư tưởng từ lý thuyết điều khiển tuyến tính Khóa luận tơi chọn đề tài tốn nhận dạng lý thuyết điều khiển Nội dung bao gồm phần sau: • Chương 1: Hệ động lực tuyến tính • Chương 2: Xác định biểu diễn tối thiểu từ tham số Markov Mục đích- Yêu cầu i Hồng Thị Thơ - Tốn K35-CN MỤC LỤC • Đây dịp để tập dượt nghiên cứu (với định hướng giáo viên hướng dẫn) nội dung khoa học • Nắm bắt nội dung lý thuyết (Các khái niệm, tính chất, tốn đặt ra, số ứng dụng, ) • Biết cách thể hiểu biết Đối tượng nghiên cứu Bài tốn nhận dạng hệ tuyến tính thời gian liên tục kiến thức liên quan Phạm vi • Các tài liệu tham khảo cá nhân tự tìm hiểu thu thập thêm • Thời gian thực khóa luận ii Hồng Thị Thơ - Tốn K35-CN Nội dung Tên đề tài Bài tốn nhận dạng lý thuyết điều khiển Kết cấu nội dung Gồm chương: • Chương 1: Hệ động lực tuyến tính - Hệ động lực tuyến tính - Hàm truyền - Các tham số Markov hàm truyền - Dạng biểu diễn tối thiểu hệ • Chương 2: Xác định biểu diễn tối thiểu từ tham số Markov - Một số đặc tính ma trận Hankel tham số Markov - Thuật toán SVD cho biểu diễn tối thiểu - Thuật toán SVD sửa đổi cho biểu diễn tối thiểu Phương pháp nghiên cứu • Thu thập, tra cứu, phân tích,tổng hợp tài liệu • Sử dụng phương pháp nghiên cứu lý thuyết điều khiển • Phương pháp quan sát, đọc sách iii Hồng Thị Thơ - Tốn K35-CN Chương Hệ động lực tuyến tính Chương giới thiệu hệ động lực tuyến tính mơ hình khơng gian trạng thái khái niệm liên quan 1.1 Hệ động lực tuyến tính Định nghĩa 1.1.1 Một hệ động lực tuyến tính liên tục với tham số bất biến mơ tả qua hai phương trình : x(t) = Ax(t) + Bu(t), x(t0) = x0 (1.1) y(t) = Cx(t) + Du(t), (1.2) đó: x(t) vectơ n chiều gọi trạng thái hệ, u(t) vectơ m (m ≤ n) chiều gọi đầu vào hệ, y(t) vectơ r chiều gọi đầu hệ, x(t0) điều kiện ban đầu Ma trận A, B, C D ma trận thực, có số chiều n × n, n × m, r × n, r × m Hồng Thị Thơ - Tốn K35-CN Nếu n=r=1, ta nói hệ thống có đầu vào đầu hay gọi hệ thống SISO Hồng Thị Thơ - Toán K35-CN CHƯƠNG HỆ ĐỘNG LỰC TUYẾN TÍNH Nếu có nhiều đầu vào đầu ta gọi hệ thống MIMO Hình 1.1: Hình minh họa hệ động lực tuyến tính Ví dụ 1.1.2 Ta xét ví dụ [2] Ta xét mạch điện RLC nối tiếp có sơ đồ sau: Hình 1.2: Mạch điện RLC nối tiếp {H1, H2, H3, H4, H5, H6, H7} = {73, 574, 4318, 32816, 249260, 1892808, 14374504} Ta sử dụng Matlab để tính tốn %Thuật tốn SVD cho biểu diễn tối thiểu %Thuật toán 9.3.1 %Đầu vào Tham số Markov : {H1, H2 H2N+1} % (N should be at least equal to the McMillan degree) %Đầu Ma trận A, B, C biểu diễn tối thiểu %Bước O %Bước Tìm SVD khối ma trận Hankel function [A,B,C]=amr(H) N=(size(H,2)-1)/2; %Xay M M=zeros(N+1,N+1); for i=1:N+1 M(i,:)=H(i:i+N); end [U,S,V] = svd(M);%M=USV’ %Tìm p tol=10^(-6); p=0; sizeS=size(S,1); for i=1:sizeS if (S(p+1,p+1) Bước 3: Tính: A = Σ −12 U M n T V Σ−12 n B = Σ 21 V T Er n n rT C = E r R1 n n m U nΣ n Ở Un Vn ma trận n cột U và m Er V , Er xác định 2.3.1 r Ví dụ Ví dụ 2.3.1 Ta xét ví dụ trước cho đầu vào đầu Dùng Matlab để tính tốn %Thuật tốn sửa đổi SVD cho biểu diễn tối thiểu %Đầu vào Tham số Markov { H1, H2 H2R }, R %Đầu A, B, C biểu diễn tối thiểu function [A,B,C]=Mamr(H) R=size(H,2)/2; %Xay MR,MR1 30 Hồng Thị Thơ - Tốn K35-CN > n MR=zeros(R,R); MR1=zeros(R,R); for i=1:R MR(i,:)=H(i:i+R-1); 31 Hoàng Thị Thơ - Tốn K35-CN MR1(i,:)=H(i+1:i+R); end [U,S,V] = svd(MR); % %Tìm n tol=10^(-6); n=0; sizeS=size(S,1); for i=1:sizeS if (S(n+1,n+1)