Giáo trình cơ sở Matlab v5.2-1 - Phần 2 Bài tập ứng dụng II - Ứng dung xử lý tín hiệu số

Tài liệu tham khảo Giáo trình cơ sở Matlab v5.2-1 - Phần 2 Bài tập ứng dụng II - Ứng dung xử lý tín hiệu số

ToolBox - Digital Signal Processing PhÇn øng dơng vỊ xư lý tÝn hiƯu Sè TÝn hiƯu vμ xư lý tÝn hiƯu Kh¸i niƯm vỊ tÝn hiƯu vμ mét khái niệm rộng Nó đợc định nghĩa Websur (1998) nh l đếm đợc lợng vật lý hay xung (nh hiệu điện thế, dòng, từ thông ) tin thông tin m truyền đợc Ví dụ nh thông tin mong mn cã thĨ lμ nhiƯt ®é vμ tÝn hiƯu ®iƯn thÕ tØ lƯ víi nhiƯt ®é nμy NhiỊu s¸ch cho r»ng tÝn hiƯu thay ®ỉi nh− mét hμm cđa thêi gian Mét sè tÝn hiƯu lμ liªn tơc ; nhiệt độ không khí, sóng biển điểm Một số tín hiệu khác l rời rạc, ví dụ nh th truyền (gửi đi) theo mà Morse Các tín hiệu rời rạc chúng nhận đợc lấy mẫu thông tin không liên tục, ví dụ nh nhiệt độ khí v áp suất đợc truyền theo khoảng thời gian định vô tuyến Máy tính xử lý tín hiệu rời rạc Để xử lý máy tính, hầu hết tín hiệu thể theo chuỗi số 1,0 Tín hiệu đợc sinh bëi c¶m biÕn (sensor), vÝ dơ nh− nhiƯt điện trở, hay tốc độ kế sinh Chuỗi số chiều thực đợc lấy mẫu khoảng không đổi Việc số hoá hình ảnh sinh chuỗi số chiều Nh chuỗi số v ma trận đợc thể đờng khác Chúng đợc lọc khỏi nhiễu, đợc modul hoá v đợc xử lý để lm rõ hình ảnh, nén khoảng động Khi m thao tác đợc thực thời gian tần số chủ đạo, chọn tần số chủ đạo v thuật toán thông thờng có kết l tăng hiệu ứng Matlab dùng để giải thuật toán nhanh v chọn tần số chủ đạo PhÇn nμy ta xư dơng Toolbox xư lý tÝn hiƯu (Signal Processing) Nó gồm 70 hm số khác ®Ĩ ph©n tÝch sè tÝn hiƯu sè vμ x©y dùng lọc với đặc Phần - ứng dụng 140 ToolBox - Digital Signal Processing tÝnh cho tr−íc ë ch−¬ng nμy ta chØ xÐt xÐt mét sè vÝ dơ sử dụng hm hay dùng, hm khác bạn đọc tự tham khảo s¸ch h−íng dÉn sư dơng Hμm läc Mét hm số hay dùng để xử lý chuỗi l hμm filter (hμm läc) Läc sè cã ë mäi n¬i trình xử lý tín hiệu Trờng hợp đơn giản tín hiệu x có thêm nhiễu (noise) Chúng ta loại bỏ chúng cách sử dơng läc H·y xem xÐt tr−êng hỵp tÝn hiƯu hiƯn thị thiết bị đo đợc ngời tự đọc Việc đọc khó khăn số l thay ®ỉi sè lμ thay ®ỉi theo thêi gian, có nhiễu cảm biến, nhng điều ny đợc cải thiện chọn hiển thị khoảng định, giá trị thực x m l tổng trọng số trung bình v giá trị cuối thị v đầu vo mới, cã nghÜa lμ yn = k1.yn-1 + k2.xn Trong ®ã yn-1 l giá trị thị cuối, v xn l ®Çu vμo míi VÝ dơ: nÕu ta chän k = 0.9 vμ k2 = 0.1 » t = linspace (0, 10, 100) ; thời gian ằ s = sin (2*pi/5*t) ; tÝn hiÖu » [t,c] = size(t) ; n = 0.1 * rand(r, c) ; nhiÔu »x=s+n ; ®Çu vμo cã nhiƠu » y(1) = x(1) ; ®iỊu kiện đầu ằ for i = : 100 ằ y(i) = 0.9* Y(i-1) + 0.1 * x(i) yn = 0.9.yn-1 + 0.1.xn Mối quan hệ ny đợc xác định với yn m n > Ta không quan tâm đến việc l giá trị y1 đợc xác định nh thÕ nμo ( ta sÏ xem xÐt ®iỊu kiƯn sau) Chän y1 = x7 → sư dơng Matlab B¹n cã thĨ so s¸nh x vμ y b»ng lƯnh chÊm ®iĨm plot » plot(t, x, t, y) KÕt qu¶ nh− hình vẽ (I.1) Phần - ứng dụng 141 ToolBox - Digital Signal Processing Hình 1.1 Lm thẳng tín hiệu nhiƠu- VÝ dơ Ta cã thĨ tÝnh theo giá trị bớc trớc yn = 0.9.yn-1 + 0.05.xn + 0.05.xn-1 Trong trờng hợp ny giá trị tự y1 , v y2 y đợc sinh nh− sau: » t = linspace (0, 10, 100) ; thời gian ằ s = sin (2*pi/5*t) ; tÝn hiÖu » [t,c] = size(t) ; n = 0.1 * rand(r, c) ; nhiễu ằx=s+n ; đầu vo có nhiễu ằ y(1) = x(1) ; điều kiện đầu » for i = : 100 » y(i) = 0.9* Y(i-1) + 0.05 * x(i) + 0.05*x(i-1): » End Cịng nh− tr−êng hỵp tr−íc chóng ta cịng lμm » plot(t, x, t, y) * Läc sè cã thĨ giíi thiệu đơn giản lọc tơng tự đơn giản lọc RC thông thấp theo ví dụ Phần - øng dơng 142 ToolBox - Digital Signal Processing H×nh 1.2 Lm thẳng tín hiệu nhiễu- Ví dụ Mạch xoay chiều nh hình vẽ Với thông số nh sau: V, 10mV, 50 MHz, R, 15Ω C, 100pF pha Điện trở tơng đơng tụ điện đợc tính nh sau Zc = 1 =j jC C Để đa thông số mạch vo ta lm nh sau PhÇn - øng dơng 143 ToolBox - Digital Signal Processing » V= 0.01; %V » R= 15; % Ohm » omega= 2*pi*50*10^6; %rad/s » C= 100*10^(-12); %F » Zc= -j/(omega*C); %Ohm Ta tính đợc tổng trở mạch ằ Z= R+ Zc Z= 16.0000 -31.8310i V dòng điện ằ I = V/Z I= 1.2606e-04+ 2.5097e-04i Điện áp rơi ®iÖn trë » Vr = R*I Vr= 0.0018 + 0.0039i V qua tụ điện đồ thị biªn pha Sư dơng MATLAB subplot » subplot (1, 2, 1) Các thông số 1, 2, có ý nghÜa lμ vμ ma trËn cđa h×nh vÏ sinh v chấm điểm theo điểm đầu Trong MATLAb 3.5 thông số ny đợc viết l (121) Gồm giá trị V, Vr, Vc nh pha cuối Để định nghĩa pha cần đa ®iĨm gèc cđa chóng » VV= [0 V]; VVc = [0 Vc]; VVr = [0 Vr] TiÕp ®Õn ta chÊm ®iĨm ®iƯn thÕ pha » plot (real(VV) , imag(VV) §Ĩ quan sát đợc thang góc pha xác trục thật cần phải nh trục ảo Có nghĩa l khung chấm điểm l hình vuông ằ axis('square') Khung trục đợc sinh nh sau » axis ([0 0.012 -0.006 0.006]) PhÇn - øng dơng 144 ToolBox - Digital Signal Processing PhÇn më réng hai trục l 0.012 Các giá trị ny đợc nhận thử nghiệm v sai số, nh đợc thị đồ thị , đánh giá trục tốt v dùng câu lệnh axis Hoặc ta cã thĨ sư dơng mét hμm kh¸c lμ abs để tìm giá trị biên pha lớn v để ®iỊu chØnh l¹i trơc Chóng ta cịng cã thĨ dïng hm hold để giữ cho đồ thị chồng lên hai hình khác,Vr v Vc ằ hold on ằ plot(real(VVr), imag(VVr)) ằ plot(real(VVc), imag(VVc)) Ton đồ thị đợc viết nh− sau: » title('(a)') » xlabel('Real'), ylabel('Imaginary') B¹n cịng cã thể xác định hình cách tự đánh dấu điểm tức l dùng hm gtext ằ text(real(V), imag(V), 'V') » text(real(Vr), imag(Vr), 'Vr') » text(real(Vc), imag(Vc), 'Vc') Sau ®ã chÊm ®iĨm toμn bé » hold off Ta cịng tính tổng điện áp nh sau(xem thêm h×nhvÏ): » Vr + Vc ans = 0.0100 - 0.0000i » 180*(angle(Vr) - -angle(Vc))/pi ans = 90 PhÇn - øng dơng 145 ToolBox - Digital Signal Processing H×nh 1.3 Diện áp mạch nói tiếp RC Để có đợc đồ thị v tính toán điện áp ta lm theo bớc sau ằ f= 50*10^6; % tần số, Hz » T = 1/f; % chu kú, s » omega= 2*pi*f; %tÇn sè gãc, rad/s » t=0: T/50 : 2*T; % chuỗi giá trị ằ v = V*sin(omega*t); ằ vr = abs(Vr)*sin(omega*t + angle(Vr)); » vc = abs(Vc)*sin(omega*t + angle(Vc)); Để có đợc điểm nh phần tử thứ hai ma trận đồ thị ằ subplot(1, 2, 2) v hcấm điểm theo lệnh sau: ằ plot (t, 1000*v, t, 1000*vr, t, 1000*vc) PhÇn - øng dơng 146 ToolBox - Digital Signal Processing » grid, title('(b)') » xtabel('t,s'), ylabel('mV') »text(t(5), 1000*v(5), 'V') »text(t(20), 1000*v(20), 'Vr') »text(t(50), 1000*v(5), 'Vc') Nh ta đà biết mạch RC l mạch lọc thông thấp Giả sử ta cho vo đầu vo 10mV, ta xem đầu tụ điện Ký hiệu điện áp vo l Vin v điện áp Vout Z = R + Rc = R + jωC Vin = Vin Z R − j / ωC 1 Vout = ZcI = Vin jC R j / C I= Hm đáp ứng tần (dùng hm biến đổi) lọc xác định bëi tØ sè H ( jω ) = Vout = Vin + jRC Giá trị c = 1/RC gọi l tần số cắt dạng không thứ nguyên H ( jω ) = 1 + jω / ω c hc ω/ωc = 2πf/2πfc =f/fc Chóng ta nhËn đợc dạng không thứ nguyên hm biến đổi H( f ) = 1 + jf / f c Chóng ta sÏ vÏ biªn vμ pha cđa hμm biÕn đổi theo hai đồ thị theo tỉ số tần v theo dÃy giá trị biến đổi Phần - ứng dụng 147 ToolBox - Digital Signal Processing Hình 1.4 Đồ thị biên v tần mạch thông thấp RC ằ fratio= : 0.01: 5; » H = ones(size(fratio))/(1+j*fratio); H×nh 1.5 Mạch thông thấp dự tính trớc Phần - øng dơng 148 ToolBox - Digital Signal Processing VÏ ®å thị theo logarithm Bây ta xem xét bi toán theo quan ®iĨm xư lý tÝn hiƯu L di v +i = R dt R Víi τ = L/R - h»ng sè thêi gian cđa m¹ch - y = i v x = v/R Đẳng thức trở thnh dy + y=x dt (1.1) Ta giả thiết y đợc lÊy mÉu kho¶ng TS nhá so víi thêi gian τ Do ®ã ta cã thĨ xÊp xØ dy y n − y n−1 = dt TS Trong ®ã yn - lấy đợc từ lần đo thứ n giá trị y, v yn-1 bớc trớc Đẳng thức (1.1) đợc viết lại nh sau y n y n −1 + y n = xn TS yn = Ta l¹i cã (1.2) b1 a x n − y n −1 a1 a1 (1.3) BiÓu thøc nμy ta có đợc từ lọc Viết cho trờng hợp tổng quát a1 y n + a y n −1 + + a na y n − na +1 = b1 x n + b2 x n −1 + + bnb x n − nb (1.4) hc yn = bn a a b1 b x n + x n−1 + + b x n− nb+1 − y n−1 − − na y n− na +1 a1 a2 a1 a1 a1 (1.5) Sư dơng MATLAB l lọc N tăng, N = max(na ,nb) - Ta có chuỗi Phần - øng dông 149 ToolBox Simulink Tool box Simulink 1.ThÕ no l Simulink? Simulink l phần mềm dùng để mô hình hoá, mô v phân tích hệ thèng ®éng Simulink cung cÊp cho ta hƯ thèng tun tính, hệ phi tuyến, mô hình thời gian liên tục hay gián đoạn hay hệ lai bao gồm liên tục v gián đoạn Hệ thống có nhiều tốc độ khác có nghĩa l phần khác lấy mẫu v cập nhật số liệu tốc độ khác Để mô hình hoá Simulink cung cấp giao diện đồ hoạ để xâu dựng mô hình nh l sơ đồ khối sư dơng thao t¸c “ nhÊn vμ kÐo” cht Víi giao diện ny bạn xây dựng mô hình nh ta xây dựng giấy Đây l khác xa mô trớc m yêu cầu ta đa vo phơng trình vi phân v phơng trình sai phân ngôn ngữ hay chơng trình Simulink bao gồm ton th viện khèi nh− khèi nhËn tÝn hiƯu, c¸c ngn tÝn hiƯu, phần tử tuyến tính v phi tuyến, đầu nối Ta thay đổi hay tạo khối riêng Các mô hình l có thứ bậc, bạn xây dựng mô hình theo cách từ dới lên hay từ xuống Bạn xem hệ thống mức cao hơn, ta nháy kép v khối để xem xét chi tiết mô hình Cách ny cho phép ta hiểu sâu sắc tổ chức mô hình v tác động qua lại phần nh no Sau tạo đợc mô hình, ta mô nã Simulink hay b»ng nhËp lƯnh cưa sỉ lệnh MATLAB Các Menu đặc biệt thích hợp cho công việc có tác động qua lại lẫn nhau, sử dụng dòng lệnh hay đợc dùng để chạy loạt mô Sử dụng Scope v khối hiển thị khác ta xem kết chạy mô Hơn bạn thay đổi thông số v xem có thay đổi cách trực tiếp Kết mô đặt vo MATLAB để xử lý đa máy in hay hiển thị Công cụ phân tích mô hình bao gồm công cụ tuyến tÝnh ho¸ vμ "trimming" mμ ta cã thĨ truy nhËp từ dòng lệnh MATLAB, ta có nhiều công cụ MATLAB v chơng tr×nh øng dơng cđa nã Vμ bëi v× MATLAB vμ Phần - ứng dụng 194 ToolBox Simulink Simulink đà đợc tích hợp nên ta mô phỏng, phân tích v sửa chữa mô hình hai môi trờng điểm no Để xem xét chơng trình cách tốt l ta xem xét vi ví dụ Bi toán thứ 2.1 Đặt bi toán cho mô hình Một ví dụ đáng ý Simulink l mô hình nhiệt động học nh Để chạy mô hình ny ta thực bớc dới đây: Chạy MATLAB Để chạy mô hình ta đánh "Thermo" cửa sổ lệnh cđa MATLAB LƯnh nμy sÏ ch¹y Simulink vμ t¹o cửa sổ chứa mô hình sau (Hình 2.2.1) Hình 2.1 Sơ đồ mô hình mô tả Simulink Phần - øng dơng 195 ToolBox Simulink Khi b¹n xem mô hình,Simulink đa hai khối hiển thị có tên "Indoor vs Outdoor Temp" v " Heat cost" 3.Để bắt đầu mô phỏng, vo menu Simulation v chọn lệnh Start ( Hoặc ấn phím Start công cụ cửa sổ Simulink) Khi chạy mô phỏng, nhiệt độ vμ ngoμi nhμ sÏ hiÓn thi khèi Scope "Indoor vs Outdoor Temp" vμ sè tiỊn nhiƯt ph¶i tr¶ sÏ xuÊt hiÖn khèi Scope " Heat Cost" §Ĩ dõng m« pháng, chän lƯnh Stop menu Simulation ( Hoặc ấn phím Pause công cụ ) 5.Khi bạn đà kết thúc việc chạy mô hình ny, đóng mô hình lệnh Close từ Menu File 2.2 Mô tả mô hình bi toán Mô hình mô nhiệt động nh l mô hình đơn giản Máy điều nhiệt đợc đặt 700 F v bị tác động nhiệt độ bên ngoi biến đổi theo luật hình sin có biên độ l 150 xung quanh nhiệt độ 500 Đây l mô thay đổi nhiệt độ hng ngy Mô hình sử dụng hệ để đơn giản hoá sơ đồ mô hìnhv tạo hệ thống sử dụng đợc Hệ l nhóm khối m đợc đại diện hệ Mô hình náy có hệ con: máy điều nhiệt, nh v hệ biến đổi nhiƯt ®é ( hai hƯ biÕn ®ỉi tõ F sang 0C vμ mét biÕn ®ỉi tõ 0C sang 0F) Nhiệt độ bên v ngoi nh đợc cấp tới hệ "House",v cập nhật nhiệt độ nh.Nhấp kép vo khối 'House" để xem khối hệ phụ ny Hình 2.2 Mô hình nhiệt động nh Phần - ứng dụng 196 ToolBox Simulink Mô hình hệ ổn định nhiệt l hoạt động máy ổn nhiệt, định no hệ thống nhiệt bật hay tắt Nháy kép vo khối để xem khối hệ ny Hình 2.3 Mô hình máy ổn nhiệt Cả nhiệt độ bên v bên ngoi nh đợc biến đổi tõ 0F sang 0C bëi mét hƯ chung H×nh 2.4 Mô hình hệ biến đổi từ độ F sang độ C Khi nhiệt đợc bật,tiền nhiệt phải trả đợc tính toán v hiển thị khối "Heat Cost", nhiệt độ bên nh đợc hiển thị khối "Indoor Temp" 2.3Thử lại số trình Có số trình m ta cần thử lại để xem mô hình đáp ứng nh no thông số khác ã Một khối hiển thị bao gồm vùng hiển thị tín hiệu v điều khiển m nã cho phÐp ta lùa chän kho¶ng tÝn hiƯu hiĨn thị, phóng to phần tín hiệu v thực công việc khác Trục honh biểu diễn thời gian v trục tung biểu diễn giá trị tín hiệu ã Khối số có tên l "Setpoint" đặt nhiệt độ yêu cầu nh Mở khối ny v đặt giá trị tới 800 F chay mô Xem nhiệt độ bên nh v tiền nhiệt thay ®ỉi Cịng nh− vËy ta cịng cã thĨ thay đổi nhiệt độ bên ngoi v xem ảnh hởng mô hình ã Điều chỉnh độ biến ®ỉi nhiƯt ®é h»ng ngμy bëi viƯc më khèi ph¸t sóng sin có tên "Daily Temp Variation" v thay đổi thông số biên độ Phần - ứng dụng 197 ToolBox Simulink 2.4 Hiệu việc mô trình Ví dụ ny lm sáng tỏ vi công việc đà đợc sử dụng để xây dựng mô hình ã Chạy mô bao gồm đặt thông số v bắt đầu mô với lệnh Start ã Bạn gói gọn ton khối có liên quan khối đơn gọi l hệ ã Bạn tạo biểu tợng v thiết kế hộp đối thoại cho khối công việc sử dụng "masking".Trong mô hình nhiệt tất hệ đợc tạo biểu tợng sử dụng "Masking" ã Khối hiển thị hiển thị đồ hoạ nh máy sóng thực Khối hiển thị hiển thị tín hiệu vo 2.5 Các ví dụ cã thĨ sư dơng kh¸c cđa Simulink C¸c vÝ dơ khác lm sáng tỏ khái niệm mô hình đợc sử dụng Bạn xem ví dơ nμy tõ cưa sỉ th− viƯn cđa Simulink §¸nh "simulink" cưa sỉ lƯnh cđa Matlab Cưa sỉ th viện khối xuất Hình 2.5 Cửa sỉ c¸c th− viƯn cđa Simulink NhÊp kÐp vμo biĨu t−ỵng " Demos" Cưa sỉ "Matlab demos" sÏ xt Cửa sổ ny có vi ví dụ đáng quan tâm m lm sáng tỏ đặc điểm sử dơng cđa Simulink PhÇn - øng dơng 198 ToolBox Simulink Phơng pháp Xây dựng mô hình Ví dụ ny trình bầy cho ta cách xây dựng mô hình nh thết no, cách sử dụng lệnh v thao tác bạn sử dụng để xây dựng mô hình Ta xây dựng mô hình tích phân sóng sin v hiển thị kết với sóng sin Sơ đồ khối mô hình nh sau: Hình 2.6 Mô hình tích phân sóng hình sin Đánh lệnh "simulink" từ cửa số lệnh Matlab ®Ĩ hiĨn thÞ cưa sỉ th− viƯn Simulink vμ nÕu cửa sổ mô hình no đợc mở cửa sổ mô hình đợc tạo Cửa sỉ th− viƯn cđa Simulink nh− sau: H×nh 2.7 Cưa sổ th việncủa Simulink Trong mô hình ny bạn lấy c¸c khèi sau tõ c¸c th− viƯn: - Th− viƯn c¸c ngn tÝn hiƯu (Khèi ph¸t sãng sin) - Th− viện khối nhận tín hiệu(Khối hiển thị) - Th viện hm tuyến tính (Khối tích phân) Phần - ứng dụng 199 ToolBox Simulink - Th viên đầu nối ( Khối chuyển mạch) Mở th viện nguồn tín hiệu để vo khối sóng sin Để mở mét th− viƯn ta nh¸y kÐp vμo nã Simulink sÏ hiển thị cửa sổ chứa tất khối cđa th− viƯn ®ã Trong th− viƯn ngn tÝn hiƯu tất khối l nguồn tín hiệu Th viện nguồn tín hiệu đợc thể nh hình 2.8 Hình 2.8 Cửa sổ th viện nguồn tín hiệu Bạn thêm khối vo mô hình bạn cách chép từ th viện hay từ mô hình khác Trong bi tập ny bạn cần chép khối phát sóng hình sin Đặt trỏ khối ấn v giữ phím chuột, kéo khối tới cửa sổ mô hình Phần - øng dơng 200 ToolBox Simulink H×nh 2.9 Copy khèi sin vo mô hình Khi bạn di chuyển khối bạn cã thĨ thÊy khèi vμ tªn cđa nã di chun với trỏ Hình 2.10 Khối v tên khối di chun cïng trá PhÇn - øng dơng 201 ToolBox Simulink Khi trỏ tới nơi bạn cần đặt khối mô hình bạn nhả phím chuột, copy khối phát hình đà mô hình bạn Hình 2.11 Cửa sổ mô hình bạn đà copy khối sóng sin Theo cách ny chép khối lại vo mô hình bạn Bạn di chuyển khối mô hình sử dụng kỹ thuật nh bạn chép khối, di chuyển khối khoảng nhỏ cách chọn khối v ấn phím mũi tên Với tất khối đà chép cửa sổ mô hình nh sau: Hình 2.12 Cửa sổ mô hình với khối đà copy NÕu b¹n xem kü tõng khèi, b¹n thÊy dÊu > bên phải khối sin v dấu bên trái khối MUX Dấu đầu khối l cổng ra, đầu vo khối l cổng vo Tín Phần - ứng dụng 202 ToolBox Simulink hiệu từ đầu khối tới đu vo khối khác theo đờng nối Khi cổng đà đợc nối biểu tợng cổng Hình 2.13 Đầu vo v khối Bạn thấy khèi MUX cã ba cæng vμo nh−ng chØ cã tín hiệu vo Để thay đổi số cổng vo bạn mở khối MUX cách nháy kép khối v thay đổi giá trị thông số " Number of Input" tíi Sau ®ã Ên phÝm Close, Simulink sÏ ®iỊu chỉnh số cổng vo Hình 2.14 Cửa sổ thông số khèi MUX HiƯn ta cã thĨ nèi c¸c khèi Nối đầu khối phát sinh tới đầu vo khối MUX Đặt trỏ tới đầu khèi sin, lóc ®ã trá sÏ thay ®ỉi thμnh chữ thập nhỏ Phần - ứng dụng 203 ToolBox Simulink Hình 2.15 Cửa sổ mô hình trớc nối dây Giữ v kéo chuột tới đầu vo khèi MUX Chó ý ®−êng lμ nÐt ®øt phÝm chuột giữ v trỏ thay đổi thnh chữ thập kép lại gần khối MUX Hình 2.16 Cửa sổ mô hình nối dây Ta nhả phím chuột v khối đà đợc nối Bạn nối cách nhả phím chuột trỏ bên khối Khi đờng nối nối vo cổng gần vị trí trỏ Phần - ứng dụng 204 ToolBox Simulink Hình 2.17 Hình sin nối vo đầu khối MUX Phần lớn đờng nối từ đầu khối tới đầu vo khối khác Có đờng nối từ đờng tới đầu vo khối gọi l đờng rẽ nhánh Vẽ đờng rẽ nhánh có khác biệt nhỏ so với vẽ đờng Để nối đờng đà có ta thực theo bớc sau: Đặt vị trí trỏ đờng cần rẽ nhánh Hình 2.18 Con trỏ đặt vo điểm cần rẽ nhánh ấn vμ gi÷ phÝm Ctrl, Ên vμ gi÷ phÝm chuét kÐo trỏ tới đầu vo khối Phần - ứng dụng 205 ToolBox Simulink Hình 2.19 Nối khối Nhả phím chuột, Simulink vẽ đờng từ điểm bắt đầu tới cổng vo khối Hình 2.20 Đờng nối dây rẽ nhánh Kết thúc việc nối dây, mô hình nh sau: Phần - ứng dụng 206 ToolBox Simulink Hình 2.21 Cửa sổ mô hình ta ®· vÏ xong B©y giê ta më khèi Scope ®Ĩ hiển thị tín hiệu v chạy mô 10s Đầu tiên ta phải đặt thông số mô b»ng lƯnh Parameter menu Simulation, hép héi tho¹i xt Chú ý Stoptime đặt l 10.0s Hình 2.22 Hộp đối thoại Parameter Simulink Phần - ứng dụng 207 ToolBox Simulink Để đóng hộp đối thoại Simulink Parameter ta ấn phím Close Simulink áp dụng thông số ta đặt v đóng cửa sổ hội thoại Chọn Start menu Simulation vμ xem sù thay ®ỉi cđa đầu vo khối hiển thị Hình 2.23 Cửa sổ hiển thị tín hiệu khối Scope Để lu mô hình sư dơng lƯnh Save menu File vμ nhËp tên v vị trí file File ny chứa mô tả mô hình Để kết thúc Simulink v Matlab chän lƯnh exit Matlab menu File hc ta đánh lệnh Quit cửa sổ lệnh Matlab Nếu bạn muốn thoát khỏi Simulink m không thoát khỏi Matlab đóng tất cửa sổ Simulink Phần - øng dông 208 ... 1, 2, 1, 2] ; » B = [0.5]; » A = [1, -0 .25 , -0 .25 ]; Gọi hm lọc nh sau: Các tần số đợc tính toán nh sau Bớc Vo Phần - øng dông Ra Z1 151 Z2 ToolBox - Digital Signal Processing 1 0. 125 0. 125 0 2 1. 125 ... 0. 125 0. 125 0 2 1. 125 0.4063 0 .28 13 0.9063 0.5078 0 .22 66 1.508 0.6035 0.3770 1.104 0.6 528 0 .27 59 1.653 0.6891 0.41 32 1.189 0.7105 0 .29 73 1.71 0. 724 9 0. 427 6 1 .22 5 0.7338 0.30 62 10 1.734 0.7397... nhận đợc mối quan hÖ sau: y1 = Zi1 + b1 + b1x1 (1.6) y2 = Zi2 + [b1x2 + b2x1] - [a2 - y1] (1.7) y3 = Zi3 + [b1x3 + b2x2 + b3x1] - [a2 - y2 + a1y3] (1.8) k yk = Zik k −1 ∑ bi x k +1−i − ∑ +1