Báo cáo môn học môn học thí nghiệm xử lý tín hiệu số

TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP TP... P.2.2 Tạo các chuỗi tuần hoàn và vẽ đồ thị của chúng sử dụng hàm stem dựa theo chu kì tuần hoàn... Sir dung ham filter tinh toán và vẽ đáp ứng xung của h

TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP TP HÒ CHÍ MINH

Họ và tên SV

KHOA CÔNG NGHỆ ĐIỆN TỬ

BÁO CÁO MÔN HỌC MON HOC: THÍ NGHIỆM XỬ LÝ TÍN HIỆU SÓ

NHÓM SỐ: 3

Nguyễn Minh Khang _ Lé Van Sy

'_ Phạm Đỗ Đức Tài Ngô Trịnh Quốc Thái

Tên lớp học phân : ĐHDTVTI7C

Mã lớp học phân : 420300381414

CIáo viên giảng

dạy :_ Nguyễn Tiến Tùng

Thành phố Hô Chí Minh, tháng 1 năm 2024

Bai 2: Tin hiéu va hé thong roi rac, Hé thong LTI

Ham xung don vi(impseq)

Tao function su dung ham xung don vi

function [x,n] = impseq(n0,nl,n2)

nl = -S;

n2 = 15;

n= [nl:n2]; x = [(n-n0) == 0];

en

Ham budc nhay don vi(stepseq)

Tao functiom stepseq

1 function [x,n] = stepseq(n0,nl,n2)

li E= n=[nl:n2]; x = [(n-n0) >= 0];

SS end

P.2.1 Tạo các chuỗi sau sử dụng các hàm tín hiệu cơ bản của MATLAB được đề

cập trong chương này Sử dụng hàm stem đê biêu diễn

x, (n) = 30(n + 2) + 26(n) — ð(n — 3) + Sð(n — 7) » amas _Š* ots,

Š x,(n)= >) Zo o(n 2k) -10<n<10 x,(n) = 10u(n) —Su(n—5) —10u(n -10) + Su(n —15)

d_ xứ0=€" '*[u(n + 20)—u(n —10)]

x;(n) = 5[cos(0.49Zm) + cos(0.Š lzm)] ~200<n<200

f x,(n)=2sin(0.01m)cos(0.5m) , —200<n<200

za

n = [-5:15];

x = 3*impseq(-2,-5,15) + 2*impseq(0,-5,15) - impseq(3,-5,15) + S*impseq(7,-5,15)

subplot (4,2,1); stem(n,x); title('Sequence in Problem 2.1a')

stem(n,x); title('Sequence in Problem 2.1a')

xlabel('n'); viabel('x(n) ');

4L 4

œÈ—©—> e ; oO a oO S OOO

5 0 5 10 15

n ) = cau b

n = -10:1:10;

temp = 0;

[-)for k = -5:5

temp = temp + exp(-l*abs(k)) *impseq(2*k,-190,10);

end

xn=temp;

subplot (4,2,2); stem(n,xn); title('2.1b')

t——e—>———P*®

05Ƒ 1

-10 8 4 ot 2 0 2 4 6 8 10

c)

3 cau c

n = [0:15];

xl = 10.*(stepseq(0,0,15));

x2 = -5* (stepseq(5,0,15));

x3 = -10* (stepseq(10,0,15));

x4 = 5*(stepseq(15,0,15));

x = x1+x2+x3+x4;

subplot (4,2,3); stem(n,x); title('2.1c')

xlabel('n'); ylabel('x(n)');

2" Li tyi

+

d)

3® cau d

n = [-20:10];

xl = exp(0.1*n).* (stepseq(-20,-20,10) -stepseq(190,-20,10));

x = xl;

subplot (4,2,4); stem(n,x); title('2.1d')

xlabel('n'); ylabel('x(n)');

©)

3 cau e

n = [-200:2001; x = 5 * (cos(0.49*pi*n) + cos(0.51*pi*n));

subplot (4,2,5); stem(n,x); title('2.1e')

xlabel('n'); viabel('x(n) ');

100

lệ S) b

Ăn

i |

| Cà,

0

0

% cau f

n = [-200:200]; x = 2 * (sin(0.01*pi*n) + cos(0.5*pi*n));

subplot (4,2,6); stem(n,x); title('2.1f')

xlabel('n'); ylabel('x(n)');

g)

$ cau g

n = [0:100];

xl = exp(-0.05*n) * sin(0.1*pi*n + pi/3);

subplot (4,2,7); stem(n,xl); title('2.1¢')

xlabel('n'); ylabel('x(n)')?

te 21g _ —_ <, 8 t=]tJ Cứ]

h)

3® cau h

= [0:100];

xl = exp(0.01*n) * sin(0.1*pi“n);

x = xl;

subplot (4,2,8); stem(n,x); title('2.ih')

xlabel('n'); ylabel('x(n)');

0 to 20 30 40

x(n} Z

P.2.2 Tạo các chuỗi tuần hoàn và vẽ đồ thị của chúng (sử dụng hàm stem) dựa

theo chu kì tuần hoàn

bộ 5@)=€”“k@)—u@—20) | Chu ki = 3

c 5ứ)=sin(0.1zml(w)~w(x~10)] _ Chụ kì = 4

Fm) = fob ot bob 2d AB hilt Bla „

nhiéu?

a) bs - Len

a kẻ ones ` q| PỊÊ EÁt View Insert Tools Desktop Window Help ”

YSHES|B\0H\kG r

tinhieuroiracm 32 | impseqm ™ | stepseqm % | lap.m 2 r 7 r 7

BC x2 = xl'*ones(1,5): x2 = (x2(:))': 15

3- stem (n1, x2)

4T sten(nl,x2); titLe|['bai 2.22") 4

08

eC

So

eI

a 0.5 +1 1.5

15 410 5 0 5 10 16

>> bt22a 7

b)

d Indent |E] p= [zp x Y Advance File Edit View Insert Tools Desktop Window Help »

= = Se su '3c7dq|8 ñng|aàzm

tinhieuroiracm % impseqm 3Í | siepseqrm ? | lap.m

1- n= (0:19);

3- & = exp(O.i*n) * (stepseq(o, 0, 18) - ste

6- stem(m, a):

aL

zL

1

0

———— 0 10 2n 30 40 50 60

QC s-::a n

>> " WH double

File Edit View Insert Tools Desktop Window Help

3- mtilde = x’ * ones{1,4); xtilde § (xtiLde(:])+* 4

4— stem(m,Xrilde): tCitle('bai 2.2C')

08

07

_—06

= o

2 05

x

04 0.3

02

=

‘Command Window VUkUHBi3 Z0 ti LUuUV4 ou 0.1

a

m

d)

7 | Comment 9% % 21 Breakpoints Run Runand 9

indent [=] 3 [fe x w Advance File Edit View Insert Tools Desktop Window Help x

1— R~ [0:11]; xL — [L 2 3];

BI 8 “ xi' * ones(1,3); a = (a(:))"; 8

3— n= [0:11]; x2 = [1 2 3 4]:

4- b= x2" * ones(1,3)7 b— (b(:))"z

€- c= x' * ones(1,4): c= (c(:))*:

$- xliabell'n'); *

11— tátle|'2.24!];

1

9.8090 9.5878 9.3080 ° 0.3

a 1x12 doubie

b Ixl2doubie

Columns 37 throuch 40

wee rang tpteeets wee see yey steeweepe

K2

y(n) — 0.5y(n — 1)+0.25y(n — 2) = x(n) + 2x(n — 1) + x(n - 3)

Cho hệ thống tuyến tính bất biến được mô tả bởi phương trình vi phân

a Sir dung ham filter tinh toán và vẽ đáp ứng xung của hệ thống trên trong khoảng 0<n<100

Dựa vào đáp ứng xung xác định hệ thống có ôn định không?

c Nếu tín hiệu đầu vào của hê thống là x(n) = [5 + 3Cos(0.2zn) + 4Sin(0.6mn)]u(n) , xác định đáp ứng xung y(n) trong khoảng 0< n < 200 bằng cách sử dụng hàm filter

Câu a

Câu b hệ thông ôn định tại ans = 5,7 la so thực

p22d.m p225.m p21Hm | p22am p22B.m p22cm

2 XOA MAN HINH, BIEN

4- clear all

S- close all

7

s~

10 -

ils

16

Command Window

ans =

8571

fe ~

2.5

cle

clear all

close all

%

% cau a

b = [1,2,0,1]; a=[1,-0.5,0.25];

n= [0:100];

h =impz(b,a,n);

subplot(1,1,1);

stem(n,h);

title(‘Impulse Response’); xlabel(‘n'); ylabel(‘h(n)')

% cau b

sum(abs(h))

Cau c

¡| sĩ fe v ¥ Advance Time

Editor - C\MATLAB\budi 2\p225c.m

p22d.m p225.m p21H.m p22am_ * | p22B.m p22c.m p226cm | File Edit View Insert Tools Desktop Window Help

— close all Aa il l

a= [1, 0, 0]; -0.5 ¡

n = 0:100; “a8

x = sin((pi*n)/(25) * (stepseq(0,0,100)-stepseq(100,0,100))}) 4

ylabel('y(n)')z £ LIIIE, UU là đi Li

mmand Window

clc

clear all

close all

0

% cau Cc

b=[I, -1];

a=[l1, 0, 0];

n=0:100;

x = sin((pi*n)/(25) * (stepseq(0,0,100)-stepseq(100,0,100)));

y = filter(b, a, x);

subplot(2,1,1);stem(n,x);

subplot(2,1,2);stem(n,y);

title(‘Output Response y(n)');

xlabel('n’);

ylabelCy(n));

P.2.6 Một bộ vi sai số đơn giản được cho bởi

y(n) = x(n) — x(n — 1)

10

Thực hiện bộ vi sai nay cho tin hiéu x(n) sau, vé đồ thị, nhận xét từng kết

quả

a X(n)= Š [ u(n) - u(n - 20)]

b x(n) =n [ u(n) — u(n — 10)] + (20 — n)[ u(m — 10) — u(n — 20)]

c x(n)= sin(=) [ u(n) — u(n — 100)]

Cau a

clc

clear all

close all

0

b=l1.-I]:

a=[1.0, 0];

n=-5:40;

x =5 * (stepseq(0,-5,40)-stepseq(20,-5,40));

y = filter(b, a, x);

subplot(2,1,1);stem(n,x);

subplot(2,1,2);stem(n,y);

tile(Output Response y(n)’);

xlabel(‘n');

ylabel(y(n));

11