A. Đặt vấn đề Nhận dạng đối là một trong những bước đầu tiên và quan trọng để thực hiện quá trình thiết kế điều khiển cho đối tượng. Sau khi thu thập được dữ liệu vào ra theo thời gian hoặc là phổ tín hiệu của đối tượng nhiệm vụ của việc nhận dạng đối tượng là tìm mô hình toán học, hàm truyền đạt thích hợp mô tả gần đúng nhất đối tượng thực. Để dễ dàng cho quá trình nhận dạng thì trong matlab có toolbox : system identification toolbox để giúp chúng ta thực hiện dễ dàng trực quan , nhanh chóng. Tuy nhiên để sử dụng tốt công cụ này chúng ta cần phải hiểu rõ các phương pháp nhận dạng, phạm vi sử dụng và ưu nhược điểm của từng phương pháp cộng với khả năng về phân tích hệ thống thông qua các đặc tính thu được. Để giúp cho quá trình tìm hiểu về toolbox này ở đây chúng tôi sẽ giới thiệu cho các bạn các bước làm cơ bản để thực hiện quá trình nhận dạng. Với mong muốn rằng các bạn có thể sử dụng nó để hiểu được trình tự nhận dạng trong matlab. Và chúng tôi hy vọng rằng sẽ thực hiện tìm hiểu thật cụ thể về toolbox này trong thời gian tới. 1 B. Các bước thực hiện nhận dạng bằng system identification toolbox I. Chuẩn bị dữ liệu để nhận dạng : 1. Dữ liệu được nhập trực tiếp từ cửa sổ matlab ta nhập giá trị vào u1,y1 rồi lưu nó vào trong file .mat >>u1= [1.2;1.3;…]; >>y1=[3;6;; ]; >>save identdata u1 y1 2. Dữ liệu lưu trong exel Để đọc dữ liệu từ trong exel ta dùng lệnh xlsread A=xlsread(‘path\file.xls’,1,’A1:A10’) A= xlsread(‘path\file.xls’,2) A=xlsread(‘path\file.xls’) Trong đó các đặc tính path\file.xls là đường dẫn đến file exel chứa số liệu mình cần lấy ví dụ :C:\Documents and Settings\phong\My Documents\identdata.xls 1 có nghĩa là sheet 1 trong exel A1:A10 mình giới hạn số liệu cần lấy để khảo sát từ A1:A10 Để sử dụng có hiệu quả mình phải biết được đặc điểm của dữ liệu lưu trong exel. Ví dụ như trong exel với tên file và đường dẫn của nó là : C:\Documents and Settings\phong\My Documents\identdata.xls nếu A2:B6 sẽ lấy ra ma trận từ cột A2 đến B6 Để lấy dữ liệu tín hiệu vào ra với cột A là tín hiệu vào u1 và cột B là tín hiệu ra y1 ta dùng đoạn lệnh >> X=xlsread('C:\Documents and Settings\phong\My Documents\identdata.xls',1); >> u1=X(:,1); >> y1=X(:,2); Và để dùng để mô phỏng ta phải lưu giá trị u1,y1 vào file identdata.mat để sử dụng ta dùng lệnh >> save identdata u1 y1 3. Nếu dữ liệu lưu trong text 2 C = textscan(fid, 'format') C = textscan(fid, 'format', N) C = textscan(fid, 'format', param, value, ) C = textscan(fid, 'format', N, param, value, ) C = textscan(str, ) [C, position] = textscan( ) Fid là tên file cũng như xlsread ta cần có cả đường dẫn của nó. Phần format thì tham khảo bảng định dạng dữ liệu sau. 3 Có thể tham khảo các hàm sau : dlmread , dlmwrite, csvread, csvwrite II. Các loại dữ liệu cho việc mô phỏng Gõ ident vào workpace của matlab 4 1. Time-DomainData dữ liệu trong miền thời gian Sau đó vào Import data > Time domain data Input : biến vào Output :biến ra Data name : là tên là tên file chứa dữ liệu vào/ra (session.mat) Starting time : nhập 0 nếu muốn khảo sát tại thời gian gốc Sampling interval : khoảng thời gian lấy mẫu. Sau đó kích vào more trên cửa sổ Import data Sau đó để đặt thêm tên biến đơn vị biến ta click vào More (đặc tính chung của các quá trình thực hiện vào dữ liệu) 5 InterSample : Lựa chọn phương pháp chuyển từ liên tục sang rời rạc (zoh, foh,bl ) Period : là thời gian khảo sát.ở đây mình chọn là inf Channel Names • Input : nhập tên của tín hiệu vào (ở đây là power) • Output: nhập tên của tín hiệu ra (ở đây là temperat ure) Physical Units of Variables (Đơn vị vật lý của các tín hiệu vào/ra) • Input : tương ứng đây sẽ là W • Output: tương ứng là độ C 2. Frequency-Domain Data dữ liệu trong miền tần số 6 Các dạng của dữ liệu : a. Freq.Function(Complex) • Freq.Func.: nhập tên biến chứa dữ liệu hoặc biểu thức của đáp ứng tần số G ( exp(jw)) • Frequency : nhập tên biến chứa mảng tần số khảo sát; định dạng của nó là vector cột • Frequency unit : đơn vị tần số (Hz hay rad/s) b. Amplitudeand Phase Frequency-Response Data Thì các đặc tính cần cài đặt là biên độ, pha, và tần số ở đây là các biến có giá trị là các vector. 7 III. Các mô hình nhận dạng I. Nhân dạng tham số mô hình ARMA ( Autoregressive moving average): 1.1. Lý thuyết: Mô hình ARMA là mô hình rời rạc, nhận dạng mô hình ARMA là phương pháp nhận dạng tham số K, a , b cho mô hình rời rạc 1 1 1 1 1 ( ) ( ) ( ) 1 nb nb na na b z b z Y z G z K U z a z a z − − − − + + + = = + + + (1) Trên cở sở quan sát, đo tín hiệu vào u(t) và ra y(t) sao cho sai lệch giữa mô hình và đối tượng là nhỏ nhất. Với những điều kiện mô tả sai lệch khác nhau sẽ có các phương pháp nhận dạng khác nhau. Các phương pháp này được chia làm 2 loại chính: - Loại nhận dạng chủ động ( active). Tín hiệu đầu vào u(t) được chọn là tín hiệu ồn trắng có giá trị mật độ phổ bằng 1, tức là: m u =0 và S u (w)=1 8 - Loại nhận dạng bị động (passive). Đặc biệt, khi n b =0 thì mô hình trở thành: 1 1 ( ) 1 ( ) ( ) 1 na na Y z G z K U z a z a z − − = = + + + (2) và được gọi là mô hình AR (Autoregressive) của đối tượng. Khi n a =0 thì mô hình trở thành: 1 1 1 ( ) ( ) ( ) 1 nb nb b z b z Y z G z K U z − − + + + = = (3) có tên gọi là mô hình MA (Moving average). Chuyển (1) sang miền thời gian, ta có mô hình tương đương sang sai phân: 1 1 ( ) na nb n k n k n k n k k k y a y K u b u − − = = + = + ∑ ∑ Hoàn toàn tương tự, dạng thời gian của mô hình AR là: 1 na n k n k n k y a y Ku − = + = ∑ của mô hình MA là: 1 ( ) nb n n k n k k y K u b u − = = + ∑ Để đơn giản trong quá trình nhận dạng ta có thuật toán: - Chuyển bài toán nhận dạng bị động sang bài toán nhận dạng chủ động. - Nhận dạng từng khâu AR và MA trong mô hình ARMA. 1.2. Linear parametric models . 1. Trong System Identification Tool GUI, kiểm tra dữ liệu đã được nhập vào chưa. Nếu sử dụng một dạng dữ liệu khác ta có thể thay thế một mô hình mẫu tính toán khác. Kiểm tra dữ liệu trong Working Data. 2. Chọn Estimate > Linear parametric models để mở hộp Linear Parametric Models, nếu hộp này chưa được mở. 3. Trong hộp Linear Parametric Models chọn By Initial Model từ danh sách Structure. 4. Nhập tên mẫu trong trường Initial model, nhấn Enter. Tên này phải là tên của mẫu trong Model Board trong System Identification Tool GUI hoặc là biến trong MATLAB Workspace. 5. Thay đổi thuật toán cài đặt, hiển thị trong Linear Parametric Models nếu cần thiết. 6. Chọn Estimate để lọc mẫu. 7. Chuẩn hoá mẫu mới này. Để tiếp tục lọc mẫu chọn Continue iter. 9 2. Spectral models: 2.1. Lý thuyết: Giống như h(t), thông qua g(t) ta luôn có được u(t) và do đó g(t) cũng có thể được xem như một mô hình không tham số của đối tượng. Việc nhận dạng mô hình không tham số sẽ đồng nghĩa với việc nhận dạng hàm quá độ h(t) hay hàm trọng lượng g(t). Một trong những phương pháp nhận dạng mô hình không tham số đơn giản là xác định hàm quá độ h(t) bằng cách kích thích đối tượng với tín hiệu Heaviside tại đầu vào rồi đo tín hiệu đầu ra. Từ g(t) ta cũng có thể suy ra h(t) do đó bài toán sẽ được coi là giải quyết xong nếu xác định được g(t) hay ảnh Fourier của nó. Nhận dạng đường đặc tính tần G(jw) cũng như các giá trị {G n } với G n =G(jnΩ) cho một đối tượng tuyến tính cần nhận dạng trên cơ sở quan sát cả 2 tín hiệu vào và ra bao gồm: - Nhận dạng mật độ phổ tín hiệu bằng kỹ thuật DFT. Đánh giá sai số và các kỹ thuật làm giảm sai số. - Xác định hàm trọng lượng, hàm quá độ từ phổ tín hiệu theo thuật toán cực tiểu sai lệch. 2.2 . Nhận dạng trong toolbox Matlab: Trước khi sử dụng Toolbox nay, cần xuất dữ liệu vào trong GUI như đã nêu trong bước chuẩn bị dữ liệu. Các bước thực hiện thuật toán nhận dạng mô hình tham số nhờ phân tích tín hiệu phổ trong System Identification Tool GUI. 1. Trong System Identification Tool GUI, chọn Estimate > Spectral models mở hộp Spectral Model. 2. Trong danh sách Method, chọn thuật toán phân tích tín hiệu phổ cần thiết. Có thể lựa chọn các dạng: • etfe (Empirical Transfer Function Estimate). • spa (SPectral Analysis). • spafdr (SPectral Analysis • with Frequency Dependent Resolution). 3. Nhập một số thông số cần thiết cho tính toán thuật toán: • Trong trường Frequencies, nhập giá trị của vector hoặc tên vector đã được nhập trong cửa sổ Matlab.Ví dụ: logspace(-1,2,500). • Sử dụng kết hợp Frequency Spacing và Frequencies để xây dựng tần số của vector. • Trong danh sách Frequency Spacing chọn Linear hoặc Logarithmic cho không gian tần số. • Trong trường Frequencies nhập tần số của các điểm ( tần số nằm trong khoảng từ tần số nhỏ nhất dến tần số lớn nhất của dữ liệu). 4.Trong trường Frequency Resolution, nhập tần số giải. Nếu sử dụng giá trị mặc định, chọn default hoặc để trống trường đó. 5. Trong trường Model Name, nhập tên tương thích với mẫu phân tích. Tên này chỉ xuất hiện duynhất trong Model Board. 10 [...]... lý của hệ System Identìication Toolbox cung cấp rất nhiều mô hình mẫu, theo thời gian, theo tần số với số điểm cực khác nhau, có thể thêm khâu tích phân, thêm hoặc bỏ thờì gian trễ, điểm không Trong toolbox này ta có các mô hình bậc 1, 2 hoặc 3 với nhiều điểm cực thực hoặc ảo Trước khi sử dụng Toolbox cần phải chuẩn bị dữ liệu ( theo các cách đã trình bày) Các bước sử dụng Toolbox: 1 Trong System Identìfication...6 Kích Estimate, nhập mẫu từ Model Board trong System Identification Tool GUI 7 Trên hộp Spectral Model kích close 8 Muốn hiển thị hình ảnh của nhiễu quá trình, chọn hộp thoại Noise spectrum trong System Identification Tool GUI Để xuất mẫu ra của sổ Matlab ta sử dụng lệnh idfrd hoặc bode 3 Process models: Hàm truyền đạt của... lại trong thời gian sớm nhất 19 Chọn < Preprocess > Remove means để trừ tín hiệu vào/ra từ dữ liệu vào/ra mục đích là đưa về gốc Lúc này time plot sẽ là 20 Chọn < Preprocess > Select range Trong windown bạn nhập dữ liệu thành 2 phần phần thứ nhất là phần dữ liệu của tín hiệu đo và phần thứ hai là phần dữ liệu của mô hình mẫu • ở đây samples là số lần đo (ở đây chọn 500) • Data name : nhập vào mô hình... hộp Trace Trong quá trình tính toán, thông tin được hiển thị trên màn hình: • Loss command • Parameter values • Search direction • Fit improvements 12 Chọn Estimate để nhập mẫu vào Model Board trong System Identification Tool GUI 13 Để ngừng quá trình tính toán và lưu kết quả lại, ta chọn Stop Iterations Để tiếp tục quá trình ta chọn Continue iter 14 Để hiển thị đồ thị của mẫu ta chọn Model Views 15... Estimate Có thể lặp lại các công việc trên đến khi xác định được một cách tương đối chính xác mô hình đối tượng IV Ví dụ tham khảo Ta chọn khâu quán tính bậc 2 G(s) = 56/((1+45s)*(1+4s)); Để thu thập dữ liệu dạng tần số ta làm như sau : syms s w1=[0:0.1:10]; G=56/((1+45*s)*(1+50*s)); w=j*[0:0.1:10]; G=subs(G,w); G=double(G); 14 A=abs(G); B=angle(G); Save data1 G w1; Save data2 A B w1; Sau đó click vào... Response 60 G 50 G1` Amplitude 40 30 20 10 0 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 Time (sec) Kết luận Do chưa tìm hiểu hết các phương pháp nhận dạng mà ở đây chỉ mang tính hướng dẫn sử dụng toolbox này của matlab nên chúng tôi chưa có được cái nhìn tổng quan để rút ra kinh nghiệm nhận dạng cho từng đối tượng và đánh giá được hết các ưu nhược điểm của các phương pháp nhận dạng.Với sự tiếp... có các mô hình bậc 1, 2 hoặc 3 với nhiều điểm cực thực hoặc ảo Trước khi sử dụng Toolbox cần phải chuẩn bị dữ liệu ( theo các cách đã trình bày) Các bước sử dụng Toolbox: 1 Trong System Identìfication Toolbox, chọn Estimate > Process models để mở Process models 11 2 Nếu mô hình có nhiều đầu ra, lựa chọn kênh trong danh sách Input Danh sách này chỉ xuất hiện khi mô hình có nhiều đầu vào 3 Trong Model... số lần đo (ở đây chọn 500) • Data name : nhập vào mô hình ước lượng sau đó cick vào insert 21 Để lưu file nay chọn File > Save session as Why Identify Models Using Default Settings? Sau khi chuẩn bị dữ liệu cho quá trình ước lượng (Plotting and Preprocessing Data),bạn có thể dung Quick Start để có kết quả nhận dạng,và so sánh với một số cấu trúc của mô hình mẫu Bạn có thể dùng các mô hình mẫu này để . tìm hiểu thật cụ thể về toolbox này trong thời gian tới. 1 B. Các bước thực hiện nhận dạng bằng system identification toolbox I. Chuẩn bị dữ liệu để nhận dạng : 1. Dữ liệu được nhập trực tiếp. . 1. Trong System Identification Tool GUI, kiểm tra dữ liệu đã được nhập vào chưa. Nếu sử dụng một dạng dữ liệu khác ta có thể thay thế một mô hình mẫu tính toán khác. Kiểm tra dữ liệu trong. đúng nhất đối tượng thực. Để dễ dàng cho quá trình nhận dạng thì trong matlab có toolbox : system identification toolbox để giúp chúng ta thực hiện dễ dàng trực quan , nhanh chóng. Tuy nhiên để