Xem mẫu
- TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI KHOA ĐIỆN CÁC BÀI THỰC HÀNH MÔN HỌC TÍN HIỆU VÀ HỆ THỐNG ( 6 bài) H ọ và tên sinh viên: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Mã lớp thí nghiệm: . . . . . . . . . . Mã số sinh viên: . . . . . . . . . . (Sinh viên phải nộp lại vào buổi bảo vệ tuần cuối cùng để chấm điểm) HÀ NỘI - 2010 2-1
- BÀI 2 Khảo sát đáp ứng của hệ thống trên miền thời gian I. Mục đích Làm quen với thư viện Simulink của MATLAB Sử dụng Simulink đ ể xây dựng sơ đồ khối hệ thống và mô phỏng các tín hiệu trong hệ thống Sử dụng MATLAB để soạn thảo các chương trình tính tích chập, tính đáp ứng của hệ thống trên miền thời gian II. Yêu cầu đối với sinh viên Thực hiện trước bài 2 .1 ở nhà. Đọc tài liệu hướng dẫn và trả lời được các câu hỏi của CBHD trước khi làm thực hành Hoàn thành nội dung bài thực hành (kể cả các b ài về nhà) trước khi tham dự buổi tiếp theo III. Nội dung 3.1 Làm quen với Simulink Mục này do CBHD hướng dẫn cho SV (khoảng 10ph). 3.2 Vẽ đồ thị đáp ứng của hệ thống với MATLAB-Simulink Bài 2.1 (tự chuẩn bị và làm ở nhà): Giải phương trình vi phân tuyến tính cấp hai Cho mạch điện trên Hình 2.1, với các sơ kiện là iL (0 ) 0 và vC (0 ) 0.5V, đ iện áp ngu ồn là tín hiệu bước nhảy đ ơn vị u (t ). Hãy tính đ iện áp của tụ điện vC (t ). dv Trong bài này, i iL iC C C (2.1) dt Và theo định luật điện áp của Kirchoff, di RiL (t ) L L vC (t ) u (t ) (2.2) dt Thay (2.1) vào (2.2) ta được Hình 2.1 Mạch điện cho Bài 2.1 d 2v dv RC C LC 2C vC (t ) u (t ) (2.3) dt dt Thay giá trị của các linh kiện vào và sắp xếp lại (2.3), ta có 1 d 2 vC 4 dvC vC (t ) u (t ) . 3 dt 2 3 dt hay d 2vC dv 4 C 3vC 3, t 0. (2.4) 2 dt dt Các em hãy giải phương trình vi phân (2.4) và ghi lời giải vào phần dưới đây. ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... 2-2
- ..................................................................................................................................................... ..................................................................................................................................................... ..................................................................................................................................................... Với kết quả tìm đ ược ở trên, các em hãy sử dụng MATLAB để vẽ đồ thị của vC (t ) vào Hình 2.2 dưới đây 1 0.9 0.8 0.7 0.6 0.5 0 1 2 3 4 5 6 7 8 9 10 Hình 2.2 Đồ thị cho lời giải của phương trình (2.4) Bài 2.2 (thực hiện có hướng dẫn tại PTN): Mô tả hệ thống bằng mô hình trạng thái Với mạch điện ở Bài 2.1, ta có: diL RiL (t ) L vC (t ) u0 (t ) dt Nếu định nghĩa các biến trạng thái là x1 iL , x2 vC và tín hiệu ra là y x2 vC đồng thời thay các giá trị vào, ta sẽ thu đ ược mô hình trạng thái như sau (sinh viên tự kiểm tra): x1 4 4 x1 4 x 3 4 0 x 0 u0 (t ) 2 2 (2.5) x y 0 1 1 x2 Sau đây chúng ta sẽ sử dụng Simulink để thu được và hiển thị tín hiệu ra vC (t ). Cách 1: Mở một file mới trong Simulink và đặt tên file là, ví d ụ mhtt_2_5. Sử dụng các khối trong Simulink Library Browser để tạo được sơ đồ khối như Hình 2.3. Lưu ý: Khối Step nằm trong S ource, khối State-Space trong Continuous, khối Bus Creator trong Common Used Blocks, và khối Scope trong Sink. . Lời chú thích cho các khối có thể đ ược thực hiện bằng tay như sau: kích đúp vào dưới mỗi khối và khi một hình chữ nhật nhỏ xuất hiện cùng với con trỏ, ta có thể gõ các thông tin cho khối đó. Để bắt đầu một dòng mới, ta ấn phím Enter. Hình 2.3 Cá ch thu được vC (t ) cho Bài 2.2 với khối State-Space 2-3
- Tiếp theo, ta kích đúp vào khối State-Space, và trong cửa sổ Function Block Parameters, ta nhập các thông số như Hình 2.4. Hình 2.4 Các thông số cho khối chức năng trong khối State-Space Các sơ kiện [x1 x2]’ được chỉ định trên Cửa sổ lệnh của MATLAB như sau x1=0; x2=0.5; Chúng ta cũng cần chỉ định thời gian mô p hỏng. Điều này được thực hiện b ằng cách, trên menu chọn Simulation Configuration Parameters rồi chọn thời gian mô p hỏng là 10s. Để bắt đầu mô phỏng, ta chọn Simulation Start ho ặc kích vào biểu tượng và để quan sát tín hiệu ra, ta kích đúp vào khối Scope, rồi kích vào biểu tượng tự động co giãn đ ồ thị . Các em hãy vẽ lại đồ thị đó vào Hình 2.5 và so sánh với đồ thị ở Hình 2.2 trước đó. 1 0.8 0.6 0.4 0.2 u (t) 0 v (t) 0 C 0 1 2 3 4 5 6 7 8 9 10 Hình 2.5 Đồ thị của vC (t ) cho Bài 2.2 Cách 2 : Từ phương trình (2.4), ta có d 2 vC dv 4 C 3vC u (t ) (2.6) 2 dt dt Sơ đồ khối biểu diễn quan hệ trong (2.6) đ ược cho trên Hình 2.7. 2-4
- Hình 2.7 Sơ đồ khối của ph ương trình (2.6) Các em hãy mở một file mới trong Simulink, đặt tên file và xây d ựng sơ đồ khối tương tự như Hình 2.7. Sau đó , hãy khai báo tham số cho các khối và thực hiện mô p hỏng tín hiệu ra với thời gian mô p hỏng là 10s (CBHD kiểm tra). Vẽ lại sơ đồ Simulink đó vào khung dưới đây. Lưu ý: Trong Simulink Library Browser, khối Add (hay Sum) và khối Gain nằm trong Math Operations, khối Integrator nằm trong Continuous. Các sơ kiện được khai báo trong các khối Integrator. Bài 2.3 (thực hiện có hướng dẫn tại PTN): Tính toán và vẽ đồ thị tích chập của hai tín hiệu Đọc chương trình MATLAB sau đ ây và cho biết chương trình này dùng để làm gì. Sau đó gõ chương trình vào file soạn thảo và lưu vào với tên file conv.m. Các em hãy tìm hiểu và giải thích từng dòng lệnh, sau đó vẽ lại đồ thị (với chu kỳ lấy mẫu 0.01 giây) trên Hình 2.8. % Program 2.3. Convolution of two signals % Approximate the convolution y(t) of the signals % x(t)=exp(-2(t-2))u(t-2) and h(t)=exp(-3(t+3))u(t+3). % Plot the approximate convolution output. % Plot the exact convolution output also. % Assume that the first nonzero value of the signals start at time t=0 % and find the convolution output. % Then, use the convolution time shift theorem to find the exact starting % time of the convolution output. clear; % clear variables and functions from memory t1=2; % time of the first nonzero value of x(t) t2=-3; % time of the first nonzero value of h(t) t3=t1+t2; % time of the first nonzero value of convolution output for ts=[0.2 0.1 0.01]; % sampling interval clf; % clear current figure 2-5
- n=0:1:4/ts; % sample index x=exp(-2*n*ts);x(1)=0.5*x(1); % sample values of x(t)=exp(-2t)u(t) % sample value is the average of the two limiting values at any % discontinuity h=exp(-3*n*ts);h(1)=0.5*h(1); % sample values of h(t)=exp(-3t)u(t) y=ts*conv(x,h); L=length(y); t=ts*(0:L-1); yf=exp(-2*t)-exp(-3*t); % exact samples of convolution plot(t3+t,y,t3+t,yf,':') % plot with solid and dashed lines ylabel('y(t)') % y-axis label of the graph is y(t) xlabel('t') % x-axis label of the graph is t legend('computed','exact') title('convolution of e^{-2(t-2)}u(t-2) and e^{-3(t+3)}u(t+3)') text(0.6,0.7,['sampling interval = ',num2str(ts)],'units','normalized','fontsize',10) disp('see the graph and press any key to continue') pause % wait until a key is pressed end convolution of e-2(t-2) u(t-2) and e-3(t+3) u(t+3) 0.16 0.14 0.12 0.1 y(t) 0.08 0.06 0.04 0.02 0 -1 0 1 2 3 4 5 6 7 t Hình 2.5. Đồ thị cho Bài 2.3. Bài 2.4 (thực hiện có hướng dẫn tại PTN): Vẽ đáp ứng của hệ thống không liên tục Chương trình sau đây tính toán và vẽ đáp ứng của hệ thống tuyến tính không liên tục đ ược mô tả bởi phương trình sai phân cấp hai: y n y n 1 0.8 y n 2 2 x n 3 x n 1 4 x n 2 với các sơ kiện y 1 2, y 2 1, x 1 x 2 0 và tín hiệu vào x n 0.5n u n . Các em hãy điền các dòng lệnh còn thiếu vào phần chấm chấm đ ể ho àn thiện chương trình (với gợi ý là các dòng chú thích bên cạnh) và vẽ lại đ ồ thị sau khi chạy mô p hỏng lên Hình 2.9. % Program 2.4. System response by iteration % Find the response of the system, for n=0 to n=M, % governed by second-order difference equation % y(n)+a1 y(n-1)+a0 y(n-2)=b2 x(n)+b1 x(n-1)+b0 x(n-2) % by iteration to the input x(n)= (0.5)^n u(n) % with the given initial conditions, y(-1) and y(-2). % Plot the response. clear; % clear variables and functions from memory clf; % clear current figure ym1=2;ym2=-1; % initial output conditions 2-6
- ....................... % initial input conditions b2=2;b1=-3;b0=4; % coefficients of the input terms ....................... % coefficients of the output terms M=10; % iterate up to n=M x=(0:M)*0;y=x; for n=0:M x(n+1)=(0.5)^n; % present input value yc=................................................................................ % present output ym2=ym1;ym1=...........; % update past output values xm2=xm1;xm1=...........; % update past input values y(n+1)=yc; % store the output values in array y end n=0:1:M; % discrete time base sequence plot(n,y,'*',n,x,'o') ; % plot with '*' and 'o' symbols grid on; ylabel('y(n)') % y-axis label of the graph is y(n) xlabel('n') % x-axis label of the graph is n legend('output','input') title('Response of a second-order system by iteration') Response of a second-order system by iteration 5 output input 4 3 y(n) 2 1 0 -1 0 1 2 3 4 5 6 7 8 9 10 n Hình 2.9. Đáp ứng của hệ thống cho Bài 2.4 Bài 2.5 (về nhà tự làm) Cho hệ thống bậc bốn được mô tả b ởi mô hình trạng thái như sau: 0 1 0 0 0 0 0 1 0 x 0 u x 0 0 0 1 0 1 0 2 0 1 y 1 0 0 0 x T trong đó x x1 x2 x3 x4 . Các em hãy xâ y d ựng sơ đ ồ khối Simulink (sử dụng khối State-space) đ ể mô phỏng tín hiệu ra của hệ thống ứ ng với sơ kiện x(0) 0, tín hiệu vào u (t ) sin t , và thời gian mô phỏng là 2 5s. Hãy vẽ lại tín hiệu ra đó vào Hình 2 .10 dưới đây. Lưu ý: Tín hiệu sin có thể được tạo ra từ khối Signal Generator trong Sources của Simulink Library Browser. 2-7
- 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2.5 0 5 10 15 20 25 Hình 2.10. Dạng tín hiệu cho Bài 2.5 Bài 2.6 (về nhà tự làm): Tính toán và vẽ đồ thị đáp ứng của hệ thống không liên tục Hãy sửa lại chương trình trong Bài 2.4 với các giả thiết sau đây: ym1=1;ym2=-2; xm1=0;xm2=0; b2=3;b1=-2;b0=1;a1=0.9;a0=0.8; x(n)=u(n). Vẽ lại đồ thị sau khi thực hiện chương trình vào Hình 2.11. Response of a second-order system by iteration 3 output input 2 1 0 y(n) -1 -2 -3 -4 0 1 2 3 4 5 6 7 8 9 10 n Hình 2.9. Đáp ứng của hệ thống cho Bài 2.6 Bài 2.6 (về nhà tự làm) Tính toán và vẽ đồ thị đáp ứng của hệ thống liên tục Hãy viết chương trình MATLAB để tính toán và vẽ đáp ứng y(t ) x(t ) h(t ) của hệ thống có đáp ứng xung h(t ) e t u (t ) và tín hiệu vào x(t ) u (t 1) u (t 2). 2-8
Chính trị học
Báo chí - Truyền thông
Xã hội học
Giáo dục học
Tâm lý học
Lịch sử - Văn hoá
Triết học
Ngôn ngữ học
Thư viện thông tin
Văn học nước ngoài
Ngư nghiệp
Hành chính - Pháp luật
Địa lý - Địa danh
Văn học Việt nam
Lịch sử Đảng
CNXH - KH
Tư Tưởng HCM
Ngụ ngôn - Cổ tích
Ca dao - Tục ngữ
Hoá học
Sinh học
Y khoa - Dược
Kinh tế học