a/ Vi t hai chu i fế ỗ 0 và m0 đ xác đ nh đ c tính c u b l c trên d a vào ể ị ặ ả ộ ọ ự remez.
b/ Thay đ i cho đúng v i l i gi i (a), b ng cách ch m đi m mổ ớ ờ ả ằ ấ ể 0 theo f0
c/ S d ng hàm ử ụ remez tìm h s c a b l c cho 10 đi m, 20, 30 (g i chúng t ng ngệ ố ủ ộ ọ ể ọ ươ ứ
v i bFIR10, bFIR20, bFIR30) và x p x b l c đ c đ a ra.ớ ấ ỉ ộ ọ ượ ư
d/ So sánh đ ho đ c tính c a b l c nh n đ c v i Fồ ạ ặ ủ ộ ọ ậ ượ ớ
6) Hàm Bilinear cùng v i tính toán kho ng l y m uớ ả ấ ẫ
Hàm chuy n đ i c a b l c thông th p đ c xác đ nh theo m t s nh sau:ể ổ ủ ộ ọ ấ ượ ị ặ ư
H(s) =
1
12
2
2 2
++ +
s p
w n
w w s s
n n
/.
ξ
v i P = 6 rad/s; wớn = 15rad/s và ξ = 0.6
a/ Hãy vi t H(s) nh là t s c a 2 đa th c ế ư ỷ ố ủ ứ num và den tính nhanh chúng nh s d ngờ ử ụ
chuy n đ i MATLAB (chu i các h s c a vi c gi m năng l ng c a s)ể ổ ỗ ệ ố ủ ệ ả ượ ủ
Yêu c u : S d ng hàm ầ ử ụ cour .
b/ Ch m các đi m biên đáp ng c a b l c đ c xác đ nh b i hàm t s H(s), trongấ ể ứ ủ ộ ọ ượ ị ở ỉ ố
kho ng 0.1 đ n 100 rad/s, b ng bi n th các t n s theo thang logarithm và biên theo dB.ả ế ằ ệ ị ầ ố
c/ Ch m đi m đáp ng pha trên cùng m t kho ng, b ng bi n th t n s theo rad/s và phaấ ể ứ ộ ả ằ ệ ị ầ ố
theo đ (s d ng ộ ử ụ unwarp n u có kh năng cài đ t vào máy c a b n).ế ả ặ ủ ạ
d/ Tìm các h s c a s t ng đ ng c a nó, l y t ng các t n s l y m u c a 50Hzệ ố ủ ố ươ ươ ủ ấ ổ ầ ố ấ ẫ ủ
Chú ý: S d ng hàm ử ụ bilinear.
7) X p x kho ng l y m u cho hàm bilinearấ ỉ ả ấ ẫ
Hàm bi n đ i c a b l c thông th p, xác đ nh theo s nh sau:ế ổ ủ ộ ọ ấ ị ư
H(s) =
1
12
2
2 2
++ +
s p
w n
w w s s
n n
/.
ξ
v i P = 8rad/s, wớn = 20rad/s và ξ = 0.65. B n mu n có b l c s t ng đ ng và mu n xác đ nhạ ố ộ ọ ố ươ ươ ố ị
t n s l y m u v i s d ng help c a MATLAB. Vi t m u Ch ng trình b ng MATLAB sau:ầ ố ấ ẫ ớ ử ụ ủ ế ẫ ươ ằ
a/ Nh c ng i s d ng hàm s l y m u fs;ắ ườ ử ụ ố ấ ẫ
b/ S d ng hàm ử ụ bilinear, c u t o b l c s t ng đ ng ng v i fs1ấ ạ ộ ọ ố ươ ươ ứ ớ
c/ Hi n th ra đ c tính và pha c a b l c t ng t và b l c s t ng đ ng.ể ị ặ ủ ộ ọ ươ ự ộ ọ ố ươ ươ
d/ L p l i b c (a), (b), (c) cho đ n khi ng i dùng tho mãn, lúc đó thoát ra kh i vòngặ ạ ướ ế ườ ả ỏ
l p.ặ
e/ D a trên vòng l p trên hi n th t n s l y m u và các thông s c a b l c s .ự ặ ệ ị ầ ố ấ ẫ ố ủ ộ ọ ố
f/ So sánh đ th c a đáp ng c a b l c s đ c xây d ng (d) v i b l c t ng tồ ị ủ ứ ủ ộ ọ ố ượ ự ở ớ ộ ọ ươ ự
đ u tiên.ầ
8)FFT và bins.
X là 256 đi m FFT, có đ c b ng DFT đ n chu i x có 256 đi m, có đ c b i t n s c aể ượ ằ ế ỗ ể ượ ở ầ ố ủ
tín hi u l p l i t i 64Hz. Nhìn vào ph công su t, chúng ta th y có nh y bin th 33. Đi u đóệ ặ ạ ạ ổ ấ ấ ả ở ứ ề
ch ra kh năng thành ph n c a t n s lên tín hi u g c (có nghĩa là d u hi u có th đ c th cỉ ả ầ ủ ầ ố ệ ố ấ ệ ể ượ ự
hi n)ệ
dùng công th c hình 1.12ứ ở
9)Aliasing (ký t )ự
Gi s có tín hi u đ c xác đ nh bài 8. N u ký t là hi n t i, v i các t n s khác cóả ử ệ ượ ị ở ế ự ệ ạ ớ ầ ố
cùng b c nh y không?ướ ả
10) Ph công su t c a tín hi u tam giác.ổ ấ ủ ệ
Tín hi u l p l i hình tam giác đ c l y m u t i 256Hz theo th i gian kho n 0.5 giây,ệ ặ ạ ượ ấ ẫ ạ ờ ả
sinh ra chu i x có 4 đi m c a giá tr 0 v i 16 đi m giá tr 1.ỗ ể ủ ị ớ ể ị
a/ Ch m đi m x nh m t hàm theo th i gian, m i kho ng 0.5 giây.ấ ể ư ộ ờ ỗ ả
b/ Tính toán ch m đi m ph t n c a xấ ể ổ ầ ủ
c/ Ch ra 5 t n s đ u tiên nó g m nh ng công su t l n nh t .ỉ ầ ố ầ ồ ữ ấ ớ ấ
Chú ý: S d ng hàm ử ụ sort
d/ X p x x b i 5 giao đ ng (g i k t qu x ppr5) và ch m đi m x và ấ ỉ ở ộ ọ ế ả ấ ể xapprs.
11) L c và tín hi uọ ệ
B l c thông th p F đ c xác đ nh v i chu i c a các h s a và b có đ c b ng cách sộ ọ ấ ượ ị ở ỗ ủ ệ ố ượ ằ ử
d ng l nh MATLAB.ụ ệ
» [b, a] = butter (5, 0.5)
Tín hi u x, đ c xác đ nh nh sau:ệ ượ ị ư
» Ts = 1/100 ; t = Ts * (1 : 500);
» f = 25 ; x = sin 92 * pi * f * t);
Có sin 2πft và đ c l y m u l i 100Hz trong đo n [0, 5] s d ng hàm ượ ấ ẫ ạ ạ ử ụ filter mà không
xác đ nh gi i h n tr ng thái zi (có nghĩa là MATLAB t xác đ nh 0), tìm tín hi u y, có đ c b iị ớ ạ ạ ự ị ệ ượ ở
cho x qua b l c F, và nh n d ng hi u qu u tiên b ng cách ch m đi m y nh hàm c a trungộ ọ ậ ạ ệ ả ư ằ ấ ể ư ủ
đo n [0,0.2).ạ
12) B l c v i mô t tr ng thái t i h n.ộ ọ ớ ả ạ ớ ạ
B l c thông th p F đ c xác đ nh b ng 1 chu i c a các h s a và b có đ c nh sộ ọ ấ ượ ị ằ ỗ ủ ệ ố ượ ờ ử
d ng dòng l nh c a MATLABụ ệ ủ
» [b, a] = butter (5, 0.5)
Tín hi u x đ c xác đ nh nh sau.ệ ượ ị ư
» Ts = 1/100 ; f = 25 ; x = sin (2 + pi * z * Ts * (1 : 500) ;
Tìm vector Z; nh l nh sauư ệ
»y = filter (b,a,x,zi)
Sinh ra vector y b t đ u v i 5 zeros. Ki m tra câu tr l i b ng cách ch m đi m b t đ u 30 nhómắ ầ ớ ể ả ờ ằ ấ ể ắ ầ
c a y.ủ
Chú ý: Dùnghàm filteric.