Trang 16 Nhận xét: • Pha dao động có sự thay đổi lệch hoàn toàn so tín hiệu gốc do dịch pha ban đầu.. Tích chập, phép biến đổi Fourier và lọc tín hiệu I.Tích chập và lọc tín hiệu âm th
Trang 1TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI
VIỆN ĐIỆN
- -
BÁO CÁO THÍ NGHIỆM Tín Hiệu Và Hệ Thống
Giảng viên hướng dẫn: Đặng Văn Mỹ
Sinh viên thực hiện:
Họ và tên Cao Hoàng Lân
Mã lớp thí nghiệm 724985
Hà Nội 2022
Trang 2Bài 1 Tín hiệu liên tục
I Hàm bước nhảy đơn vị (unit step) và hàm dốc đơn vị (ramp) Bài 1:
Trang 9Tín hiệu x1 không phải là tín hiệu tuần hoàn b)
Tín hiệu x2 là tín hiệu tuần hoàn với chu kỳ T0=5
IV Năng lượng, công suất của 1 tín hiệu
Trang 12Bài 2 Hàm tuyến tính Bài 1 Tần số và nốt nhạc
Hãy viết một chương trình Matlab để chơi bản nhạc sau: CCGGAAG , FFEEDDC ,
% program to play the ABC song with pure tone
Trang 14Bài 2: Tiến kèn Trumpet
a.Xây dựng note B của kèn Trumpet với số liệu cho trước.
-Th'*ones(1,Fs)); sound(x,Fs);
b Vẽ đồ thị của x(t) trong vòng 3 chu kì đầu Do tần số của B là 494HZ
và được được lưu 44100 samples mỗi giây nên 44100 / 494 là số samples của 1 chu kì
Vậy 3 chu kì có số samples xấp xỉ 268
Trang 17Bài 3 Tích chập, phép biến đổi Fourier và lọc tín hiệu
I Tích chập và lọc tín hiệu âm thanh bằng bộ lọc thông thấp lý
% bo loc thong thap ly tuong voi tan so cat wb
= wb/(2*pi)*sinc(wb*t/(2*pi));
% dau vao: du lieu, dap ung tuyen tinh, bat bien: ht
% dau ra: y = tich chap cua du lieu voi ht
Bài 1 Tự tạo các hàm của bạn
Bạn hãy viết chương trình tạo hàm FourierTransform để tính ảnh Fourier của một tín hiệu như dưới đây Sau đó lưu vào với tên là FourierTransform.m
Trang 18function [f,X]=FourierTransform(t,x)
% computes the Fourier transform of signal x(t)
% ns: length(x)=number signal points
% dt: signal point spacing
Viết chương trình tạo hàm IFourierTransform để tính ảnh Fourier ngược như dưới đây Sau đó lưu vào với tên là IFourierTransform.m
function [t, x] = IFourierTransform(f, X) %
computes the inverse Fourier transform of X(f)
% ns: length(X)=number of transform points
% df: transform point spacing
Bài 2 Lọc tín hiệu điện tim
1 Vẽ tín hiệu điện tim và phổ biên độ của nó
Trang 19% bo loc thong thap ly tuong voi tan so cat wb
= wb/(2*pi)*sinc(wb*t/(2*pi));
% dau vao: du lieu, dap ung tuyen tinh, bat bien: ht
% dau ra: y = tich chap cua du lieu voi ht
Trang 20Hình 1
=> Nhận xét:
• Từ độ thị ta thấy hầu hết năng lượng phân bổ ở vùng có tần số thấp
• Phần có tần số cao được coi là nhiễu và cần loại bỏ
title("Bo loc tin hieu dien tim");
xlabel('t (sec)'); ylabel('h(t)');
Trang 21xlabel('f (Hz)'); ylabel('\angle H(f)');
Hình 2
=> Nhận xét: Bộ lọc được sử dụng ở đây là 1 bộ lọc
thông thấp
• Bộ lọc nhằm mục đích lọc lấy tín hiệu dưới 50Hz tuy nhiên do tính "không
lý tưởng" của bộ lọc nên cho tín hiệu từ dưới 100Hz thông qua
3 Tín hiệu sau khi lọc
Yf=Xf.*Hf; figure(3) subplot(2,1,1)
plot(f,abs(Yf), 'r' , 'linewidth' ,1);
grid; axis([-150 150 0 0.12]);
title( "Tin hieu dien tim sau xư
ly" ); xlabel( 'f (Hz)' ); ylabel( '|
Trang 22Hình 3
=> Nhận xét:
• Sau khi qua bộ lọc, hầu hết nhiễu với tần số cao đều đã được khử đi
• Đồ thị tín hiệu theo thời gian trông rõ ràng hơn so hình 1