ĐẠIHỌC BÁCHKHOAHÀNỘITRƯỜNGĐIỆN –ĐIỆNTỬ BÁOCÁOTHỰCHÀNH TÍN HIỆU HỆ THỐNG HỌC KỲ 20221 Họvàtên NguyễnHữuThái MSSV 20212962 Mãhọcphần EE2000 Mãlớp 724990 Hànội2/2023 Bài1 Tínhiệuliêntục I hàm bước nhảyđơ[.]
ĐẠIHỌC BÁCHKHOAHÀNỘITRƯỜNGĐIỆN –ĐIỆNTỬ BÁOCÁOTHỰCHÀNH TÍN HIỆU HỆ THỐNG HỌC KỲ 20221 Họvàtên:NguyễnHữuThái MSSV:20212962 Mãhọcphần:EE2000 Mãlớp:724990 Hànội2/2023 Bài1:Tínhiệuliêntục I hàmbướcnhảyđơnvị(unitstep)vàhàmdốcđơnvị(ramp) Bài 1:Viếthàmustep(t) functiony=ustep(t,m); n = length(t); y=zeros(1,n);fori=1:n ift(i)>=-m y(i)=1; end end Bài2:viếthàmuramp(t) functiony=uramp(t,m);n=length(t) ; y=zeros(1,n); fori=1:n ift(i)>=-m y(i)=t(i)+m; end end Bài3:Sửdụngcác hàmvừaviếtđểvẽđồthịcủa cáctínhiệuliêntụcsau trênđoạn -10≤t≤10 +5u(t-2) >>t=linspace(-10,10,1000); >>y=5*ustep(t,-2); >>plot(t,y) >>gridon >>title('5u(t-2)') >>axis([-1010-16]) +3r(t+5) >>t=linspace(-10,10,1000); >>y=3*uramp(t,5); >>plot(t,y); >>title('3r(t+5)'); >>axis([-1010-550]) >>gridon +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); >>gridon >>axis([-1010-26]) >>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); >>gridon >>title('y=sin(t)*[u(t+3)+u(t-3)]') Bài4:Sửdụnghaihàmtrênđểtạo rađồthị: Đồ thị1 >>t=linspace(-10,10,1000); >>y=0.5*uramp(t,4)+0.5*uramp(t,-4)-uramp(t,0); >>plot(t,y) >>gridon >>axis([-1010-0.52]) Đồ thị2 >>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) >>gridon >>axis([-101003]) II Tínhiệuchẵn,lẻ Bài1.Xâydựnghàmsốtrảvềkếtquảlàphầnchẵnvàphầnlẻcủamộttínhiệunhưsau: function[ye,yo]=evenodd(y)yr=fliplr( y); ye=0.5*(y+yr); yo=0.5*(y-yr); Bài 2.Sửdụnghàmsốtrênđểtìmphầnchẵnvàphầnlẻcủacáctínhiệuliêntụcsauvà 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'); >>holdon; >>plot(t,yo,'red '); >>plot(t,ye,'green-.'); >>legend('y','yo','ye'); >>gridon III tổngcủacáctínhiệutuầnhồ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); >>gridon -Tínhiệukhơngphảitín hiệutuầnhồn vì: x1(t)=1+1.5cos(2πΩ0t)-0.6cos(4Ω0t) Hayx1=1+f1-f2 Trongđó: f1 =1.5cos(2πΩ0t)v i T1 = f =0.6cos(4Ωt) =π2Ω0 =2 T2 π Khônglàsốhữutỷ 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); >>gridon Ω0 vớiT T1 Bài2:Hàmtuyếntí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=[CCGGAAGSiSiFFEEDDCSiSi]; sound(song,Fs) a, t=linspace(0,1,44100);F=494; C=[0.11550.34170.17890.12320.06780.04730.0260.00450.002]; Th=[-2.131.67-2.5450.661-2.0392.16-1.04671.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))axistigh t,girdon; c, t=linspace(0,1,44100);F=494; C=[0.11550.34170.17890.12320.06780.04730.0260.00450.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)) gridon Nhậnxét:Saukhiđổicácphathành0,tínhiệmthanhthuđượccóđồthịâm thay đổi, nhiên âm nghe tai người khơng đổi Bài3:Tíchchập,phépbiếnđổiFourier vàlọctínhiệu I Tíchchậpvàlọc tínhiệmthanhbằngbộlọcthơngthấplý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)); ... >>legend(''y'',''yo'',''ye''); >>gridon III tổngcủacáctínhiệutuầnhồ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); >>gridon -Tínhiệukhơngphảitín hiệutuầnhồn vì: x1(t)=1+1.5cos(2πΩ0t)-0.6cos(4Ω0t)... -Tínhiệul? ?tín hiệutuầnhồn vì: x1(t)=1+1.5cos(6πt)-0.6cos(4Ω0t) Hayx1=1+f1-f2 Trongđó: f1 =1.5cos(6πt) f2=0.6cos(4Ω0t) vớiT1 =13s vớiT2=5s T1 = làsốhữutỷ=>T T2 15 IV Nănglượngcơngsuấtcủamộttínhiệu... gridon Nhậnxét:Saukhiđổicácphathành0,tínhiệmthanhthuđượccóđồthịâm thay đổi, nhiên âm nghe tai người khơng đổi Bài3:Tíchchập,phépbiếnđổiFourier vàlọctínhiệu I Tíchchậpvàlọc tínhiệmthanhbằngbộlọcthơngthấplýtưởng