Báo cáo thực hành tín hiệu hệ thống 2019 2020

Báo cáo thực hành tín hiệu hệ thống Sinh viên Nguyễn Tiến Hoàng 20173900 DK TDH 05Bài 1 Tín hiệu liên tụcI.Hàm bước nhảy đơn vị (unit step) và hàm dốc đơn vịBài 1 %Hàm ustepfunction y=ustep(t,t0)N=length(t);y=zeros(1,N);for i=1:N if t(i)>=t0 y(i)=1; end;end; Bài 2 %Hàm urampfunction y=uramp(t,t0)N=length(t);y=zeros(1,N);for i=1:N if t(i)>=t0 y(i)=t(i); end;end;Bài 4 %Hình 1>> t=10:0.01:10;>> y1=uramp(t+4,0)2uramp(t,0)+uramp(t4,0)2;>> plot(t,y1,linewidth,2) %Hình 2>> y2=ustep(t+8,0)+uramp(t+4,0)2uramp(t,0)+uramp(t4,0)2ustep(t8,0);>> plot(t,y2,linewidth,2) II. Tín hiệu chẵn, lẻ Bài 1: Xây dựng hàm số trả về kết quả là phần chẵn và phần lẻ của một tín hiệu như sau:function ye,yo=evenodd(y)yr=fliplr(y);ye=0.5(y+yr);yo=0.5(yyr);endBài 2: Sử dụng hàm số ở bài một tìm phần chẵn lẻ của các tín hiệu liên tục và vẽ đồ thị của tín hiệu chính cũng như phần chẵn và lẻ của nó trong cùng một đồ thị sử dụng các dạng đường thẳng và màu sắc khác nhau( miền t là 10,10): y(t)=2r(t+2,5)5r(t)+3r(t2)+u(t4).>> t=linspace(10,10,500);>> y= 2uramp(t,2.5)5uramp(t,0)+3uramp(t,2)+ustep(t,4);>> evenodd(y);>> yr=fliplr(y);>> ye=0.5(y+yr);>> yo=0.5(yyr);>> plot(t,yo,b); hold on; >> plot(t,ye,r); >> xlabel(t); >> ylabel(y); >> title(do thi ham so y(t)); >> plot(t, y, k, linewidth, 1)>> legend(yo,ye,y); >> grid on; III)Tồng các tín hiệu tuần hoàna.>> t=10:0.01:10;>> O=pi10; %omega>> x1=1+1.5cos(2piOt)0.6cos(4Ot);>> plot(t,x1)>> grid on =>Hàm không tuần hoànb. x2=1+1.5cos(6pit)0.6cos(4Ot);>> plot(t,x2)>> grid on =>Hàm tuần hoàn với chu kì T=5IV. Năng lượng, công suất của một tín hiệu:syms tT=20;u(t)=heaviside(t);x(t)=exp(t).cos(2pit).u(t);f=(abs(x(t))).2;E=int(f, t, T2, T2)P=int(f, t, T2, T2)T>> BaiIV E = (exp(20)(2pi2 + 1)(exp(20) 1))(2(4pi2 + 1)) P = (exp(20)(2pi2 + 1)(exp(20) 1))(40(4pi2 + 1))V. Phép dịch, phép co dãn và phép đảo tín hiệuBài 1: Tín hiệu bị dịcht=10:0.01:10;x=exp((abs(t)));x1=exp((abs(t2)));x2=exp((abs(t+2)));plot(t,x,b,linewidth,2);hold on;plot(t,x1,k,linewidth,2);plot(t,x2,r,linewidth,2);legend(x(t),x(t2),x(t+2));set(gca,fontsize,24);title(x(t)=exp(|t|)); Bài 2: Co dãn tín hiệut=10:0.01:10;x=exp((abs(t)));x1=exp((abs(2t)));x2=exp((abs(0.5t)));plot(t,x,b,linewidth,2);hold on;plot(t,x1,k,linewidth,2);plot(t,x2,r,linewidth,2);legend(x(t),x(2t),x(0,5t));set(gca,fontsize,24);title(x(t)=exp(|t|)); Bài 3: Phép đảo tín hiệut=10:0.01:10;x1=exp((abs(t)));x2=exp((abs(t)));hold on;plot(t,x1,k,linewidth,2);plot(t,x2,y,linewidth,2);legend(x(t),x(t));set(gca,fontsize,24);title(x(t)=exp(|t|));

Báo cáo thực hành tín hiệu hệ thống

Sinh viên Nguyễn Tiến Hoàng 20173900 DK

TDH 05

Bà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ị

Bài 1 %Hàm ustep

function y=ustep(t,t0)

N=length(t);

y=zeros(1,N);

for i=1:N

if t(i)>=t0

y(i)=1;

end;

end;

Bài 2

%Hàm uramp

function y=uramp(t,t0)

N=length(t);

y=zeros(1,N);

for i=1:N

if t(i)>=t0

y(i)=t(i);

end;

end;

Bài 4

-10 -8 -6 -4 -2 0 2 4 6 8 10 -0.5

0

0.5

1

1.5

%Hình 2

>> y2=ustep(t+8,0)+uramp(t+4,0)/2-uramp(t,0)+uramp(t-4,0)/2-ustep(t-8,0);

>> plot(t,y2,'linewidth',2)

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

0.5

1

1.5

2

2.5

II Tín hiệu chẵn, lẻ

Bài 1: Xây dựng hàm số trả về kết quả là phần chẵn và phần lẻ của một tín hiệu như sau:

function [ye,yo]=evenodd(y)

yr=fliplr(y);

ye=0.5*(y+yr);

yo=0.5*(y-yr);

end

Bài 2: Sử dụng hàm số ở bài một tìm phần chẵn lẻ của các tín hiệu liên tục và vẽ đồ thị của tín hiệu chính cũng như phần chẵn và lẻ của nó trong cùng một đồ thị sử dụng các dạng đường thẳng

và màu sắc khác nhau( miền t là [-10,10]):

y(t)=2r(t+2,5)-5r(t)+3r(t-2)+u(t-4).

>> t=linspace(-10,10,500);

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

>> evenodd(y);

>> yr=fliplr(y);

>> ye=0.5*(y+yr);

>> yo=0.5*(y-yr);

>> plot(t,yo,'b '); hold on;

III)Tồng các tín hiệu tuần hoàn

a >> t=-10:0.01:10;

>> O=pi/10; %omega

>> x1=1+1.5*cos(2*pi*O*t)-0.6*cos(4*O*t);

>> plot(t,x1)

>> grid on

-10 -8 -6 -4 -2 0 2 4 6 8 10 -1.5

-1

-0.5

0

0.5

1

1.5

2

2.5

=>Hàm không tuần hoàn

b x2=1+1.5*cos(6*pi*t)-0.6*cos(4*O*t);

>> plot(t,x2)

>> grid on

-1.5

-1

-0.5

0

0.5

1

1.5

2

2.5

3

3.5

=>Hàm tuần hoàn với chu kì T=5

>> BaiIV

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, phép co dãn và phép đảo tín hiệu Bài 1: Tín hiệu bị dịch

t=-10:0.01:10;

x=exp(-(abs(t)));

x1=exp(-(abs(t-2)));

x2=exp(-(abs(t+2)));

plot(t,x,'b','linewidth',2);

hold on;

plot(t,x1,'k','linewidth',2);

plot(t,x2,'r','linewidth',2);

legend('x(t)','x(t-2)','x(t+2)');

set(gca,'fontsize',24);

title('x(t)=exp(-|t|)');

Bài 2: Co dãn tín hiệu

t=-10:0.01:10;

x=exp(-(abs(t)));

x1=exp(-(abs(2*t)));

x2=exp(-(abs(0.5*t))); plot(t,x,'b','linewidth',2); hold on;

plot(t,x1,'k','linewidth',2); plot(t,x2,'r','linewidth',2); legend('x(t)','x(2t)','x(0,5t)'); set(gca,'fontsize',24);

title('x(t)=exp(-|t|)');

Bài 3: Phép đảo tín hiệu

t=-10:0.01:10;

x1=exp(-(abs(t)));

x2=exp(-(abs(-t)));

hold on;

plot(t,x1,'k','linewidth',2); plot(t,x2,'y ','linewidth',2); legend('x(t)','x(-t)');

set(gca,'fontsize',24);

title('x(t)=exp(-|t|)');