Giao trinh matlab v5.2 P16

Tương tự với phương pháp đặt tên kho nhớ, là phương pháp đặt tên những nhóm của các chỉ thị. Trong hầu hết các ngôn ngữ lập trình, đều có cho phép gọi đến các macro hay các chương trình con như là các câu lệnh để thi hành nội dung mô tả trong các macro hay chương trình con này thông qua tên.
Nội dung trích xuất từ tài liệu:
Giao trinh matlab v5.2 P16ToolBox - Digital Signal Processing » [b, a] = butter (5 , 30/50) b= 0.1084 0.5419 0.0837 1.0837 0.5149 0.1084 a= 1.000 0.9853 0.9738 0.3864 0.1112 0.0113 Vμ tr¹ng th¸i ®Çu cña kh©u trÔ » Zi = filteric (b, a, x, yi) ; Chóng ta liªn hÖ víi hμm läc » y = filter(b, a, x, Zi) ;Chóng ta cã thÓ hiÖn thÞ 5 nhãm ®Çu cña chuçi vμo vμ ra vμ thay ®æi cho chóngb»ng nhau » [ X(1:5), Y(1:5) ] ans= 0.2190 0.2190 0.0970 0.0470 0.6789 0.6789 0.6793 0.6793 0.9347 0.93474. thiÕt kÕ c¸c bé läc sè 4.1. C¸c ®Þnh nghÜa : Mét trong nh÷ng vÊn ®Ò chung nhÊt xuÊt hiÖn trong xö lý sè tÝn hiÖu lμ cÊu tróccña bé läc víi c¸c ®Æc tÝnh biªn t¹i c¸c tÇn sè kh¸c nhau. Mét trong tools (c«ng cô)trong toolbox xö lý tÝn hiÖu lμ 2 hμm yulewalk vμ remez. Chóng ta gäi chóng víi bé läc sè H cña N ®iÓm, ®Æt tÇn sè lÊy mÉu sè liÖu x,sinh ra tÇn sè míi y, quan hÖ víi x theo ®¼ng thøc. a1yn + a2yn-1 + ... aN+1 yn-N = b1xn + b2xn-1 + ... bN+1 xn-N (1.12) C¸c hÖ sè B = [b1, b2, ... , bN+1+] vμ A = [a1, a2,..., aN+1] ®Òu x¸c ®Þnh c¸c hÖ sè ≠ 0.Song chóng ta cã thÓ gi¶ thiÕt chóng chuÈn theo a1. H¬n thÕ n÷a, t¹i c¸c hÖ sè cuèi 155PhÇn 2 - øng dôngToolBox - Digital Signal Processingcïng aN +1 hoÆc bN+1 cã thÓ kh¸c 0, trong c¸c tr−êng hîp kh¸c bÞ läc cÇn x¸c ®Þnh vectorthu gän A vμ B, vμ chóng cÇn nhá h¬n N. Hμm trong toolbox MATLAB sinh ra c¸c hÖsè cña bé läc . (yulewalk, cheb1 vμ c¸c hμm kh¸c ...) lu«n sinh ra c¸c hÖ sè qui chuÈn,thμnh phÇn cña hμm läc (filter). Khi sö dông thao t¸c dÞch thêi gian nh− x¸c ®Þnh phÇn tr−íc, bé läc H ®−îcbiÓu diÔn b»ng hμm ph©n thøc sau. b1 + b2 Z −1 + bN Z − ( N −1) B (x) H (z ) = = (1.13) H ( z ) a1 + a 2 Z −1 +...+ a N Z − ( N −1) Víi a1 = 1 vμ hÖ sè lín nhÊt aN vμ bN ≠ 0 RÊt tiÕc lμ trong tÊt c¶ c¸c version cña MATLAB ta khi dïng help yulewalkkh«ng ®−îc ®óng l¾m, chóng sÏ hiÖn ra ®o¹n v¨n b¶n nh− sau yulewalk Recusive filter design using a least-squaresmethod. [B,A] = yulewalk(N,F,M) finds the N-th order recursive filter coefficients B and A such that the filter: -1 -(n-1) B(z) b(1) + b(2)z + .....+ b(n) z ------ = ----------------------------------------- -1 -(n-1) A(z) 1+ a(1)z + ........+ a(n)z Trong ®ã chØ sè cña A bÞ sai dÞch nh− sau n = N + 1 , song cho vÝ dô, bé läccña 4 sÏ ®−îc x¸c ®Þnh bëi vectors B = [b1, b2 , b3, b4, b5] A = [a1, a2, a3, a4, a5] Víi a1 = 1 vμ hÖ sè cuèi b5 hoÆc a5 ≠ 0 D¹ng nh− ®· nãi ë trªn bÞ giíi h¹n v× lçi ë trong s¸ch sö dông ®óng cña t−¬ngøng víi hμm sè nh− bé läc. NÕu c¸c nhãm a2, a3,..., aN ®Òu b»ng 0 th× bé läc sÏ gäi FIR (bé läc ®¸p øngxung h÷u h¹n). yulewalk ®−îc dïng ®Ó tæng hîp bé läc IIR , khi hμm remez ®−îc södông cho FIR 4.2 X¸c ®Þnh ®Æc tÝnh tÇn cña bé läc. 156PhÇn 2 - øng dôngToolBox - Digital Signal Processing MATLAB cho phÐp ta ®Þnh nghÜa sè cña tÇn sè fr1 , fr2 , ..., fri ,..., frk vμ biÕnt−¬ng øng mag1, mag2 , ..., magi,..., magk vμ m« t¶ bé läc sè gÇn ®óng (xÊp xØ) víi ®¸pøng cña bé läc t−¬ng tù. TÇn sè ®¸p øng cña bé läc sè phô thuéc vμo tÇn suÊt lÊy mÉu . Mét nöa tÇn s许p øng cña lÊy mÉu ®−îc gäi lμ Nyquist .Mét sè hμm cña MATLAB ®−îc nhanhchãng ®©û vμo vïng tÇn sè kh«ng thø nguyªn , cã nghÜa lμ b»ng viÖc ®Þnh nghÜa tÇn sèkh«ng thø nguyªn 50/(1000/2) = 0.1 vμ 150 Hz sÏ t−¬ng øng víi 150/(1000/2) = 0.3. §Ó x¸c ®Þnh ®Æc tÝnh cña bé läc, chóng ta cÇn cã 2 chuçi : 1 lμ tÇn sè kh«ng thønguyªn, f=[f1, f2 ,..., fk ] vμ 1 t−¬ng øng víi biªn m = [m1, m2, ..., mk] MATLAB ®ßi háif1 = 0 vμ fk = 1. Nh− trong vÝ dô , chóng ta gi¶ thiÕt r»ng tÝn hiÖu x ®−îc lÊy mÉu ë 500Hz vμchóng ta muèn x©y dùng bé läc víi tÇn sè biªn nh− sau: Tõ §Õn Biªn 0 100 1.0 100 150 gi¶m ®Õn tõ 1 ®Õn 0.5 150 200 lu«n lμ h»ng 0.5 200 225 t¨ng ®Òu tõ 0.5 ®Õn 1 200 250 - 157PhÇn 2 - øng dôngToolBox - Digi ...