(TIỂU LUẬN) báo cáo THÍ NGHIỆM bài 1 GIỚI THIỆU MATLAB cơ bản

Thực hành 5 Thực hành 5: Viết một file-M vẽ đáp ứng tần số của hệ thống có hàm truyền trên theo 2 cách: i tính toán đáp ứng tần số và ii sử dụng hàm Matlab, và lưu với tên Bai1_NHOMx_fre

ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH

TRƯỜNG ĐẠI HỌC BÁCH KHOA

-⁂ -BÁO CÁO THÍ NGHIỆM

BÀI 1: GIỚI THIỆU MATLAB CƠ BẢN

Họ và tên sinh viên: Trần Trung Kiên

MSSV: 1913863

GVHD: Huỳnh Văn Phận

x = 0:pi/20:2*pi; % Define vector x from 0 to 2pi with step size pi/100

y = sin(x); % Define vector y

subplot(1,2,1); % Create a graph with two sub-graphs in 1 row 2 columnand with subplot position 1

plot(x,y,'b-'); % Plot x versus y

axis([0 2*pi -1 1]);

xlabel('x (pi)');

ylabel('y=sin(x)'); %Label of x-axis and y-axis

title('Graph of Continuous Sine from 0 to 2pi'); %Title of graph

subplot(1,2,2);stem(x,y,'r-') %Plot p versus q with subplot position 2axis([0 2*pi -1 1]);

Kết quả :

stem(t,Xk_magnitude,'bo-'); hold on;

stem(t,Xk_2_magnitude,'r* '); hold on;

ylabel ('Amplitude'); xlabel ('K');

title('Magnitude Response');

Kết quả:

4 Thực hành 4

Thực hành 4: Viết một file-M tìm ngõ ra y(n) của hệ thống nhân quả với ngõ vào x(n) = [1, 3, 5, 3, 6, 3] và đáp ứng xung h(n) = [1, 4, 7, 2, 8] theo 2 cách, và lưu với tên Bai1_NHOMx_conv.m

Kết quả:

5 Thực hành 5

Thực hành 5: Viết một file-M vẽ đáp ứng tần số của hệ thống có hàm truyền trên theo 2 cách: i) tính toán đáp ứng tần số và ii) sử dụng hàm Matlab, và lưu với tên Bai1_NHOMx_freqz.m

Kết quả:

6 Thực hành 6

Thực hành 6: Viết các chương trình Matlab thực hiện các yêu cầu sau:

1 Tạo 5 chu kỳ mẫu tín hiệu với Hz, tần số lẫy mẫu 8000 Hz.

title('Đồ thị tín hiệu theo thời gian);

%Đồ thị tín hiệu theo mẫu

Kết quả:

2) Tạo mẫu tín hiệu được lấy mẫu với 20 trong 1 chu kỳ T Các lệnh matlab:

title('Đồ thị tín hiệu theo thời gian');

%Đồ thị tín hiệu theo mẫu

Kết quả:

3) Tạo tín hiệu sinc với Hz, (s) và được lấy mẫu ở tần số 500 Hz.

title('Đồ thị tín hiệu theo thời gian');

%Đồ thị tín hiệu theo mẫu

Kết quả:

7 Thực hành 7

Thực hành 7: Thực hiện việc thiết kế bộ lọc FIR chắn dải với các thông số như trong ví dụ 1 Lấy các hệ số của bộ và viết chương trình vẽ đáp ứng tần số của bộ lọc trên Lưu lại với tên BAI_1_NHOMx_bs2700_freqz.

Sử dụng hàm semilogy thay cho plot đưa ra nhận xét về sự khác biệt.

Kết quả:

8 Thực hành 8

Thực hành 8: Thực hiện việc thiết kế bộ lọc FIR chắn dải với các thông số như trong ví dụ 2 Lấy các hệ số của bộ và viết chương trình vẽ đáp ứng tần số của bộ lọc trên Lưu lại với tên BAI_1_NHOMx_bs2700_freqz.

Sử dụng hàm semilogy thay cho plot.

Thiết kế bộ lọc:

Các lệnh matlab:

[z,p,k] = tf2zp(filt2.tf.num, filt2.tf.den);sos = zp2sos(z,p,k);

Kết quả:

10.Thực hành 10

Thực hành 10: Thiết kế bộ lọc FIR chắn đa dải 1000-1500 và 2500-3000,

có bậc 62, tần số lấy mẫu là 10 kHz Sau đó vẽ đáp ứng tần số của bộ lọc Thiết kế bộ lọc:

F N=1

2 10 kHz=5000 Hz

chuẩn hóa f/FN

11.Thực hành 11

Thực hành 11: Thiết kế bộ lọc IIR chắn đa dải 1000-1500 và 2500-3000, có bậc

62, có tần số lấy mẫu là 10 kHz Sau đó vẽ đáp ứng tần số của bộ lọc

Kết quả: