Chương 2 trình bày những nội dung chính sau: Khái niệm về tín hiệu rời rạc, các tín hiệu cơ sở, dãy xung đơn vị (tiếp), dãy nhảy bậc đơn vị, dãy tín hiệu hình sin, dãy e-mũ phức, các phép toán trên tín hiệu,. . | CHƯƠNG II Xử lý tín hiệu nâng cao Tín hiệu rời rạc Khái niệm về tín hiệu rời rạc Trong DSP, tín hiệu thời gian rời rạc, được biểu thị bằng một dãy rời rạc: x[n]={-3 , 2, 4, -4, 0, 1 } Quá trình rời rạc hóa còn gọi là quá trình lấy mẫu tín hiệu Các tín hiệu cơ sở Dãy xung đơn vị: hay còn gọi là hàm Delta, có giá trị bằng đơn vị khi đối số = 0 và có giá trị bằng 0 trong các trường hợp còn lại: Một tín hiệu thời gian rời rạc bất kỳ có thể được khai triển từ các dãy xung đơn vị Dãy xung đơn vị (tiếp) Trong Matlab ta có thể biểu diễn như sau: function[x,n]=impseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0]; Ví dụ: Tạo dãy xung đơn vị trong khoảng [-5:5] n=[-5:5] x=impseq(0,-5,5) stem(x) Kết quả: Các tín hiệu cơ sở (tiếp) Dãy nhảy bậc đơn vị: Dãy nhảy bậc đơn vị có giá trị bằng đơn vị khi đối số lớn hơn hoặc bằng 0, và bằng 0 khi đối số nhỏ hơn 0. Một dãy tín hiệu rời rạc theo thời gian bất kỳ x[n] để có thể khai triển thành một tổng các dãy xung nhảy bậc đơn vị Dãy nhảy bậc đơn vị (tiếp) Trong Matlab để tạo ra dãy xung nhảy bậc đơn vị ta xây dựng hàm stepseq: function [x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)>=0]; Ví dụ: tạo dãy nhảy bậc đơn vị trong khoảng [-5:5] x=stepseq(0,-5,5) stem(x) Kết quả: Thực hành Vẽ đồ thị tín hiệu: x[n]=2*δ[n-5]-4*δ[n+7] trên đoạn [-10:10] Vẽ đồ thị tín hiệu: x[n]=(n-10){u(n-5)-u(n-15)} trên đoạn [0:20] Các tín hiệu cơ sở (tiếp) Dãy tín hiệu hình sin: Dãy tín hiệu hình sin được biểu thị bằng hàm số sin (hoặc cos). Trong Matlab, hàm sin (hoặc cos) được sử dụng để tạo ra dãy tín hiệu này. Ví dụ tạo dãy tín hiệu: trên đoạn [0:50] Kết quả: Dãy tín hiệu hình sin (tiếp) Ví dụ 2: Tín hiệu (trong ví dụ trên) bị ảnh hưởng bởi nhiễu Gauss: y[n]=x[n]+*w[n] Kết quả: Các tín hiệu cơ sở (tiếp) Dãy e-mũ phức: được định nghĩa bởi hệ thức: Trong Matlab ta sử dụng hàm exp để tạo ra các dãy e-mũ phức. Ví dụ với dãy trên đoạn [-10:30] clc n=[-10:30]; x=exp(*n*pi); subplot(221);stem(real(x));title('Phan thuc'); subplot(222);stem(imag(x));title('Phan ao'); subplot(223);stem(abs(x));title('Bien do'); subplot(224);stem(angle(x));title('Pha cua x'); Dãy e-mũ phức (tiếp) Hàm e-mũ phức là một hàm tuần hoàn Các phép toán trên tín hiệu Phép dịch chuyển Phép nhân, cộng tín hiệu Phép nhân chập 2 tín hiệu Thực hành Tạo hàm dịch chuyển sigshift, thực hiện dịch tín hiệu xung đơn vị trễ pha 5 mẫu Thực hiện nhân chập 2 tín hiệu sau: x[n]=[1, 2, 3, 4, 5] và h[n]=[2, 2] Phương trình sai phân Với kích thích đầu vào x(n) và đáp ứng ra y(n) của hầu hết các hệ thống tuyến tính thỏa mãn: Phương trình sai phân Biểu diễn trong Matlab, người ta sử dụng một hàm filter, hàm này có định dạng y=filter(b,a,x), Ví dụ: y(n)-y(n-1)+(n-2)=x(n) cần tính h(n) a=[1,-1,]; b=[1]; x=impseq(0,-10, 100); h=filter(b,a,x); stem(h); Phương trình sai phân Tín hiệu hai chiều (ảnh số) Biểu diễn ảnh Một ảnh được biểu diễn dưới dạng một hàm f(x,y) Đối với ảnh số đơn sắc, giá trị f tại (x,y) được gọi là mức xám Kết quả của quá trình lấy mẫu và lượng tử hóa là một ma trận số liệu. Một ảnh có kích thước M x N là một ma trận có M hàng và N cột, mỗi một giá trị trên ma trận gọi là một điểm ảnh (pixel). Biểu diễn ảnh Một ma trận tọa độ ảnh trong Matlab được biểu diễn Đọc ảnh, hiển thị ảnh Trong Matlab sử dụng hàm imread(‘filename’) để đọc một ảnh: Hàm imshow(image) trong được sử dụng để hiện ảnh. Lưu một ảnh từ ma trận I ra một file, sử dụng hàm imwrite(image,’filename’): Cải thiện ảnh Biểu đồ histogram của ảnh I figure, imhist(I) Cân bằng Histogram I2 = histeq(I); Hiệu chỉnh Histogram I = imread(''); J = imadjust(I,[ ],[]); imshow(I), figure, imshow(J); Tạo ảnh nhị phân level = graythresh(I); (graythresh: xác định ngưỡng mức xám của ảnh sử dụng phương pháp Otsu) bw = im2bw(I,level); figure, imshow(bw)
