ĐẠI HỌC BÁCH KHOA HÀ NỘI TRƯỜNG ĐIỆN – ĐIỆN TỬ Hà nội 2/2023 BÁO CÁO THỰC HÀNH TÍN HIỆU HỆ THỐNG HỌC KỲ 20221 Họ và tên Nguyễn Hữu Thái MSSV 20212962 Mã học phần EE2000 Mã lớp 724990 Bài 1 Tín hiệu li[.]
ĐẠI HỌC BÁCH KHOA HÀ NỘI TRƯỜNG ĐIỆN – ĐIỆN TỬ BÁO CÁO THỰC HÀNH TÍN HIỆU HỆ THỐNG HỌC KỲ 20221 Họ tên : Nguyễn Hữu Thái MSSV :20212962 Mã học phần : EE2000 Mã lớp : 724990 Hà nội 2/2023 Bài 1: Tín hiệu liên tục I.hàm bước nhảy đơn vị (unit step) hàm dốc đơn vị (ramp) Bài 1: Viết hàm ustep(t) function y= ustep(t,m); n = length(t); y = zeros(1,n); for i=1:n if t(i) >=-m y(i) = 1; end end Bài 2: viết hàm uramp(t) function y=uramp(t,m); n=length(t); y=zeros(1,n); for i=1:n if t(i)>=-m y(i)=t(i)+m; end end Bài 3: Sử dụng hàm vừa viết để vẽ đồ thị tín hiệu liên tục sau đoạn -10 ≤t≤10 + 5u(t-2) >>t=linspace(-10,10,1000); >>y=5*ustep(t,-2); >>plot(t,y) >>grid on >>title('5u(t-2)') >>axis([-10 10 -1 6]) +3r(t+5) >>t=linspace(-10,10,1000); >>y=3*uramp(t,5); >>plot(t,y); >>title('3r(t+5)'); >>axis([-10 10 -5 50]) >>grid on + y(t)= 2r(t+2,5)-5r(t)+3r(t-2)+u(t-4) >>t=linspace(-10,10,1000); >>y=2*uramp(t,2.5)-5*uramp(t,0)+3*uramp(t,-2)+ustep(t,-4); >>plot(t,y); >>grid on >>axis([-10 10 -2 6]) >>title('y=2r(t+2.5)-5r(t)+3r(t-2)+u(t-4)') + y(t)=sin(t)*[u(t+3)-u(t-3)] >>t=linspace(-10,10,1000); >>y=sin(t).*(ustep(t,3)+ustep(t,-3)); >>plot(t,y); >>grid on >>title('y=sin(t)*[u(t+3)+u(t-3)]') Bài 4: Sử dụng hai hàm để tạo đồ thị : Đồ thị >>t=linspace(-10,10,1000); >>y=0.5*uramp(t,4)+0.5*uramp(t,-4)-uramp(t,0); >>plot(t,y) >>grid on >>axis([-10 10 -0.5 2]) Đồ thị >>t=linspace(-10,10,1000); >>y=ustep(t,8)-ustep(t,-8)+0.5*uramp(t,4)+0.5*uramp(t,-4)-uramp(t,0); >>plot(t,y) >>grid on >>axis([-10 10 3]) II Tín hiệu chẵn , lẻ Bài Xây dựng hàm số trả kết phần chẵn phần lẻ tín hiệu sau: function [ye,yo]=evenodd(y) yr=fliplr(y); ye=0.5*(y+yr); yo=0.5*(y-yr); Bài Sử dụng hàm số để tìm phần chẵn phần lẻ tín hiệu liên tục sau vẽ đồ thị tín hiệu phần chẵn phần lẻ đồ thị sử dụng dạng đường thẳng màu sắc khác nhau: (giả sử -10≤t≤10) y(t)=2r(t+2,5)-5r(t)+3r(t-2)+u(t-4) >>linspace(-10,10,1000); >>y=2*uramp(t,2.5)-5*uramp(t,0)+3*uramp(t,-2)+ustep(t,-4); >>[ye,yo]=evenodd(y); >>plot(t,y,'black'); >>hold on; >>plot(t,yo,'red '); >>plot(t,ye,'green-.'); >>legend('y','yo','ye'); >>grid on III tổng tín hiệu tuần hồn a, >>w=pi/10; >>t=-10:pi/100:10; >>x1=1+1.5*cos(2*pi*w*t)-0.6*cos(4*pi*t); >>plot(t,x1); >>grid on -Tín hiệu khơng phải tín hiệu tuần hồn : x ( t )=1+1.5cos(2πΩ0t)-0.6cos(4Ω0t) Hay x1 =1 +f1 -f Trong đó: f =1.5cos(2πΩ0t) với T1 = Ω f =0.6cos(4Ω0t) T1 = Không số hữu tỷ T2 π b, >>w=pi/10; >>t=-10:pi/100:10; >>x2=1+1.5*cos(6*pi*t)-0.6*cos(4*w*t); >>plot(t,x2); >>grid on π với T2 = 2Ω Bài 2: Hàm tuyến tính 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) a, 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 0.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, subplot(211),plot(t(1:200),x(1:200)) axis tight,gird on; c, 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 Nhận xét: Sau đổi pha thành 0, tín hiệu âm thu có đồ thị âm thay đổi, nhiên âm nghe tai người không đổi Bài : Tích chập, phép biến đổi Fourier lọc tín hiệu I Tích chập lọc tín hiệu âm lọc thông thấp lý tưởng [data,Fs]=audioread('female_v oice.wav'); data=data(:,1)'; Ts=1/Fs; sound(data,Fs); t=[-10:Ts:10]; wb=1500*2*pi; ht=wb/(2*pi)*sinc(wb*t/ (2*pi)); y=conv(data,ht,'same'); y=y/max(abs(y)); sound(y,Fs); [data,Fs]=audioread('female_v oice.wav'); data=data(:,1)'; Ts=1/Fs; sound(data,Fs); t=[-10:Ts:10]; wb=1500*2*pi; ht=wb/(2*pi)*sinc(wb*t/ (2*pi)); y=conv(data,ht,'same'); y=y/max(abs(y)); sound(y,Fs); [data,Fs]=audioread('female_v oice.wav'); data=data(:,1)'; Ts=1/Fs; sound(data,Fs); t=[-10:Ts:10]; wb=1500*2*pi; ht=wb/(2*pi)*sinc(wb*t/ (2*pi)); y=conv(data,ht,'same'); y=y/max(abs(y)); sound(y,Fs); [data,Fs]=audioread('female_v oice.wav'); data=data(:,1)'; Ts=1/Fs; sound(data,Fs); t=[-10:Ts:10]; wb=1500*2*pi; ht=wb/(2*pi)*sinc(wb*t/ (2*pi)); ... Nhận xét: Sau đổi pha thành 0, tín hiệu âm thu có đồ thị âm thay đổi, nhiên âm nghe tai người khơng đổi Bài : Tích chập, phép biến đổi Fourier lọc tín hiệu I Tích chập lọc tín hiệu âm lọc thơng thấp... >>legend(''y'',''yo'',''ye''); >>grid on III tổng tín hiệu tuần hồn a, >>w=pi/10; >>t=-10:pi/100:10; >>x1=1+1.5*cos(2*pi*w*t)-0.6*cos(4*pi*t); >>plot(t,x1); >>grid on -Tín hiệu khơng phải tín hiệu tuần hồn : x ( t )=1+1.5cos(2πΩ0t)-0.6cos(4Ω0t)... >>y=ustep(t,8)-ustep(t,-8)+0.5*uramp(t,4)+0.5*uramp(t,-4)-uramp(t,0); >>plot(t,y) >>grid on >>axis([-10 10 3]) II Tín hiệu chẵn , lẻ Bài Xây dựng hàm số trả kết phần chẵn phần lẻ tín hiệu sau: function [ye,yo]=evenodd(y) yr=fliplr(y); ye=0.5*(y+yr);