báo cáo thí nghiệm ee2000 tín hiệu và hệ thống

III.Tổng các tín hiệu tuần hoàn...6IV.Năng lượng và công suất của một tín hiệu...7V.. Phép dịch,co giãn và phép đảo tín hiệu...8.

ĐẠI HỌC BÁCH KHOA HÀ NỘITRƯỜNG ĐIỆN – ĐIỆN TỬ

BÁO CÁO THÍ NGHIỆM EE2000-TÍN HIỆU VÀ HỆ THỐNG

HỌC KỲ 20231Họ và tên: Trần Lê Huy Bình

MSSV: 20222106

Lớp: Kỹ thuật ĐK-TĐH 13 K67Mã lớp: 736431

III.Tổng các tín hiệu tuần hoàn 6

IV.Năng lượng và công suất của một tín hiệu 7

V Phép dịch,co giãn và phép đảo tín hiệu 8

BÀI 1: TÍN HIỆU LIÊN TỤC

I, Hàm bước nhảy, hàm dốc

Bài 1:

function y = ustep(t,m)y = (sign(t)+1)/2;end

Bài 2:

function y = uramp(t,m)y = t.*ustep(t)end

Bài 3:a)

>> y=5*ustep(t,-2);>> plot(t,y);>> axis([-10 10 -1 6])

>> y=3*uramp(t,5); >> plot(t,y)

>>

y=2*uramp(t,2.5)-5*uramp(t,0)+3*uramp(t,-2)+ustep(t,-4); >> plot(t,y)

>> y=sin(t).*(ustep(t,3)-ustep(t,-3)); >> plot(t,y)

Bài 4:

function [y,z] = bai4a(t)t=linspace(-10,10,5000);

plot(t,y,'linewidth',2,'linewidth',0.5);grid on

axis([-10 10 -0.5 2])subplot(2,1,2)

plot(t,z,'linewidth',2,'linewidth',0.5);grid on

axis([-10 10 0 3])end

Bài 2:

[ye, yo] = evenodd(y);hold on;

grid on;

plot(t,y,'r -');plot(t,ye,'m ');plot(t,yo,'b -.');end

III.Tổng các tín hiệu tuần hoàn

a) >> t=linspace(-10,10,5000);>> x=1+1.5*cos(2*pi*pi/10*t)-0.6*cos(4*pi/10*t);

>> plot(t,x)

Tín hiệu này không tuần hoàn vì T1=π2

5 và T2=2∗π5T 1T 2=

2 không là số hữu tỷ

b) >> t=linspace(-10,10,4000);

x=1+1.5*cos(6*pi*t)-0.6*cos(4*pi/10*t);>> plot(t,x)

Tín hiệu này tuần hoàn với chu kì T=5(s)

IV.Năng lượng và công suất của một tín hiệu

>> syms t

>> u(t)=heaviside(t);>> syms x

>> x(t)=exp(-t).*cos(2*pi*t).*u(t);>> T=20;

>> E=int((abs(x(t)))^2,-T/2,T/2)>> P=E/T

Kết quả:E =

(exp(-20)*(2*pi^2 + 1)*(exp(20) - 1))/(2*(4*pi^2 + 1))

P =

(exp(-20)*(2*pi^2 + 1)*(exp(20) - 1))/(40*(4*pi^2 + 1))

V Phép dịch,co giãn và phép đảo tín hiệu

Bài 1:

>> t=linspace(-10,10,5000);>> x=exp(-abs(t));>> x1=exp(-abs(t-2));>> x2=exp(-abs(t+2));>> hold on;>> plot(t,x,'g');>> plot(t,x1,'r');>> plot(t,x2,'b');

>> legend('x(t)','x(t-2)','x(t+2)');>> grid on;

Bài 2:

>> t=linspace(-10,10,5000);>> x=exp(-abs(t));>> x1=exp(-abs(2*t));>> x2=exp(-abs(0.5*t));>> plot(t,x,'g');>> hold on>> plot(t,x1,'r');>> plot(t,x2,'b');

>> legend('x(t)','x(2t)','x(0.5t)')>> grid on

Bài 3.

>> t=linspace(-10,10,5000);>> x=exp(-abs(t));>> x1=exp(-abs(-t));>> plot(t,x,'g ');>> hold on>> plot(t,x1,'b');>> legend('x(t)','x(-t)')>> grid on

BÀI 2 HÀM TUYẾN TÍNHBài 1: Tần số và nốt nhạc

Viết chương trình Matlab để chơi bản nhạc sau: CCGGAAG , FFEEDDC—>> Fs=8000;

>> t=[0:1/Fs:0.5];

>> Fc=262;Fd=294;Fe=330;Ff=349;Fg=392;Fa=440;Fb=494;Fj=0;>> C=cos(2*pi*Fc*t);

>> D=cos(2*pi*Fd*t);>> E=cos(2*pi*Fe*t);>> F=cos(2*pi*Ff*t);>> G=cos(2*pi*Fg*t);>> A=cos(2*pi*Fa*t);>> B=cos(2*pi*Fb*t);>> Si=cos(2*pi*Fj*t);

>> song=[C C G G A A G Si Si F F E E D D C Si Si];>> sound(song,Fs)

Bài 2: Chuỗi Fourier củaTrumplet

a) Tổng hợp tín hiệu kèn từ phép tính xấp xỉ chuỗi Fourier hữu hạn 9 số hạng:

>> t = linspace(0,1,44100);>> F = 494;

>> C =[0.1155 0.3417 0.1789 0.1232 0.0678 0.0473 0.026 0.0045 0.002];

>> Th = [-2.13 1.67 -2.545 661 -2.039 2.16 -1.0467 1.858 -2.39];>> x = C*cos(2*pi*F*[1:9]'*t-Th'*ones(1,44100));

>> sound(x,44100);

b) Vẽ tín hiệu trong 3 chu kỳ của nó:

>> subplot(211)>> plot(t(1:200),x(1:200)) >> axis tight; grid on;

c) Lặp lại phần a và phần b với việc thay đổi pha có ảnh hưởng đến âm thanh của tín hiệu không?

>> t = linspace(0,1,44100);>> F = 494;

>> C = [0.1155 0.3417 0.1789 0.1232 0.0678 0.0473 0.026 0.0045 0.002];

>> Th = zeros(1,9);

>> x = C*cos(2*pi*F*[1:9]'*t-Th'*ones(1,44100)); >> sound(x,44100);

>> subplot(211), plot(t(1:200),x(1:200)) >> grid on;

Việc thay đổi pha không làm ảnh hưởng đến âm thanh của tín hiệu.