BÀI PHÚC TRÌNH BÀI 4. MÔ PHỎNG CÁC TÍN HIỆU CƠ BẢN BẰNG MATLAB

Phần II: Lý thuyết tín hiệu BÀI PHÚC TRÌNH BÀI MƠ PHỎNG CÁC TÍN HIỆU CƠ BẢN BẰNG MATLAB Nhóm: Tiểu nhóm: Thứ: 5, ngày 12 tháng 09 năm 2019 STT Họ tên sinh viên MSSV Trần Anh Tuấn B1708987 Vũ Bá Nhân Tô Thiệu Huy B1708963 B1708943 Chữ ký Cam đoan: Chúng ký tên cam đoan nội dung báo cáo (bao gồm đoạn code, đồ thị) kết chúng tơi thực hiện, kết chép từ báo cáo tác giả khác Chúng chấp nhận hình thức kỷ luật Bộ mơn/Giảng viên phụ trách học phần đưa vi phạm điều nêu -Câu Chương trình sau mô hàm bước nhảy đơn vị u(t) clear; % xóa tất biến nhớ Matlab clf; t=-10:0.001:10; u=heaviside(t); plot(t,u); title('Ham buoc nhay don vi u(t)'); %tên tiêu đề đồ thị xlabel('Thoi gian lien tuc t'); %chú thích trục x ylabel('Bien do'); %chú thích trục y axis([-10 10 1.2]); grid on; %hiễn thị lưới cho đồ thị **Yêu cầu a Hãy gõ lại chương trình file soạn thảo Matlab lưu lại với tên step.m b Hãy ghi lời thích vào sau dấu % số dòng lệnh (một số dòng làm mẫu) c Vẽ lại đồ thị tập lên hình 4.1 -1- Phần II: Lý thuyết tín hiệu Hình 4.1 d Hãy sửa lại chương trình Matlab mơ hàm bước nhảy đơn vị u(t) để thực yêu cầu sau đây: - Mô hàm bước nhảy đơn vị dịch chuyển phía trái đơn vị clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t=-10:0.001:10; %tạo vector thời gian u=heaviside(t+3); % tạo hàm bước nhảy đơn vị? plot(t,u); title('Ham buoc nhay don vi u(t)'); tiêu đề đồ thị - Mô hàm bước nhảy đơn vị u(t – t0), u(t + t0) với t0 nhập từ bàn phím Hàm bước đơn vị u(t – t0) Hàm bước đơn vị u(t + t0) clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t0=input('nhap vao t'); t=-10:0.001:10; %tạo vector thời gian u=heaviside(t-t0); %tạo hàm bước nhảy đơn vị plot(t,u); title('Ham buoc nhay don vi u(t)');% tên tiêu đề đồ thị clear; % xóa tất biến nhớ Matlab clf; %xóa đồ thị t0=input('nhap vao t') t=-10:0.001:10; % tạo vector thời gian u=heaviside(t+t0); % t?o hàm bước nhảy đơn vị? plot(t,u); title('Ham buoc nhay don vi u(t)'); -2- Phần II: Lý thuyết tín hiệu xlabel('Thoi gian lien tuc t'); %giải thích trục x ylabel('Bien do'); %giải thích trục y axis([-10 10 1.2]); grid on; xlabel('Thoi gian lien tuc t'); %giải thích trục x ylabel('Bien do'); %giải thích trục y axis([-10 10 1.2]); grid on; Câu Cho đoạn lệnh sau đây: >> syms t >> u=heaviside(t); >> int(u,t) **Yêu cầu a Hãy gõ đoạn lệnh vào cửa sổ Command Window ghi lại kết Kết quả: (t*(sign(t) + 1))/2 b Viết chương trình Matlab mơ hàm dốc đơn vị r(t) với tỉ lệ thời gian hình 4.2 vẽ lại đồ thị hàm lên hình 4.2 clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t=-10:0.001:10; %tạo vector thời gian r=t.*heaviside(t); plot(t,r); title('Ham doc don vi r(t)'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Hình 4.2 c Hãy sửa lại chương trình Matlab mơ hàm dốc đơn vị r(t) để thực yêu cầu sau đây: - Thay đổi màu đồ thị thành màu đỏ đường vẽ đường chấm chấm clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t=-10:0.001:10; %tạo vector thời gian r=t.*heaviside(t); %tạo hàm đơn vị? plot(t,r,'r:'); title('Ham doc don vi r(t)'); %tên tiêu đề đồ thị? xlabel('Thoi gian lien tuc t'); %giải thích cho trục x ylabel('Bien do'); %giải thích cho trục y grid on; - Mô hàm dốc đơn vị r(t – t0), r(t + t0) với t0 nhập từ bàn phím -3- Phần II: Lý thuyết tín hiệu Hàm dốc đơn vị r(t – t0) Hàm dốc đơn vị r(t + t0) clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t0=input('nhap vao t0'); t=-10:0.001:10; %tạo vector thời gian r=(t-t0).*heaviside(t-t0); %tạo hàm doc don vi? plot(t,r); title('Ham doc don vi r(t)'); %tên tiêu de cua thi? xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; clear; %xóa tất biến nhớ Matlab clf; %xóa đồ thị t0=input('nhap vao t0') t=-10:0.001:10; %tạo vector thoi gian r=(t+t0).*heaviside(t+t0); %tao hàm doc don vi? plot(t,r); title('Ham doc don vi r(t)'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Câu Cho đoạn lệnh sau đây: >> syms t >> u=heaviside(t); >> diff(u,t) **Yêu cầu a Hãy gõ đoạn lệnh vào cửa sổ Command Window ghi lại kết Kết quả: dirac(t) b Viết chương trình Matlab mơ hàm phân bố Delta δ(t) với tỉ lệ thời gian hình 4.3 vẽ lại đồ thị hàm lên hình 4.3 clear; %xóa tất biến nhớ Matlab clf; %xóa thi hien thi t=-10:1:10; %tao vector thoi gian delta =dirac(t) %tao hàm phân bo dirac plot(t,delta); title('Ham doc don vi r(t)'); %tên tiêu de cua thi? xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Hình 4.3 c Hãy sửa lại chương trình Matlab mơ hàm phân bố Delta δ(t) để thực yêu cầu sau đây: - Độ rộng cho đường thẳng vẽ đồ thị point clear; clf; t=-10:1:10; %tao vector thoi gian delta =dirac(t) %tao ham phan bo dirac plot(t,delta,'linewidth',2); title('Ham doc don vi r(t)'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; - Mô hàm phân bố Delta δ(t - t0) , δ(t  t0) với t0 nhập từ bàn phím -4- Phần II: Lý thuyết tín hiệu Hàm xung Delta δ(t - t0) clear; clf; t0=input('nhap vao t0') t=-10:1:10; %tạo vector thời gian delta =dirac(t-t0); %tạo hàm phân bo dirac plot(t,delta,'linewidth',2); title('Ham doc don vi r(t)'); %tên tiêu đề đồ thị? xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); Hàm xung Delta δ(t  t0) clear; clf; t0=input('nhap vao t0') t=-10:1:10; delta =dirac(t+t0); plot(t,delta,'linewidth',2); title('Ham doc don vi r(t)'); %tên tiêu đề đồ thị xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); Câu a Viết chương trình Matlab mơ tín hiệu xung vuông xung tam giác clear; đồ thị clf; t=-2:0.001:2; % tạo vector thoi gian x=1*rectpuls(t,1) y=1*tripuls(t,2) subplot(2,1,1) plot(t,x) title('Xung vuong co ban'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); subplot(2,1,2) plot(t,y) title('Xung tam giac co ban'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; b Vẽ lại đồ thị xung với tỉ lệ thời gian hình 4.4 -5- Phần II: Lý thuyết tín hiệu Hình 4.4 c Hãy sửa lại chương trình Matlab để mơ tín hiệu xung vuông xung tam giác dịch chuyển phía trái đơn vị clear; clf; t=-2:0.001:2; x=1*rectpuls(t+1,1) y=1*tripuls(t+1,2) subplot(2,1,1) plot(t,x) title('Xung vuong co ban'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); subplot(2,1,2) plot(t,y) title('Xung tam giac co ban'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Câu a Viết chương trình Matlab mơ hàm Gausse với tỉ lệ thời gian hình 4.4 vẽ lại đồ thị hàm lên hình 4.5 clear; clf; t=-5:0.001:5; x=gaussmf(t,[0.5,0]) plot(t,x) title('Ham Gausse'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Hình 4.5 b Hãy sửa lại chương trình Matlab mơ hàm Gausse để thực yêu cầu sau đây: - Mơ hàm Gausse dịch chuyển phía phải đơn vị độ rộng xung giảm clear; nửa clf; t=-5:0.001:5; x=gaussmf(t,[0.25,3]); plot(t,x); title('Ham Gausse'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; -6- Phần II: Lý thuyết tín hiệu - Mơ hàm Gausse dịch chuyển phía phải đơn vị hàm Gausse dịch chuyển phía trái đơn vị đồ thị clear; clf; t=(-5:0.001:5); x=gaussmf(t,[0.5,2]) y=gaussmf(t,[0.5,(-3)]) plot(t,[x y]) title('Ham Gausse'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; Câu a Viết chương trình Matlab để tạo tín hiệu điều hịa x(t) = Asin(0t + φ) vẽ tín hiệu hàm plot Các thơng số sau nhập từ bàn phím: biên độ (A), tần số góc (ω0), pha ban đầu (φ) Với điều kiện ≤ ω0 ≤ π ≤ φ ≤ 2π Lưu ý: Nếu tần số góc pha ban đầu nhập vào không thỏa mãn điều kiện chương trình báo lỗi clear; clf; t=-5:0.001:5; A=input('Nhap vao bien A=') omega=input('Nhap vao omega=') pha=input('Nhap vao pha ban dau =') if((omegapi)), error('Tan so goc khong hop le'); end; if ((pha2*pi)), error('Pha ban dau khong hop le'); end; x=A*sin(omega*t+pha) plot(t,x) title('Tin hieu dieu hoa'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; b Hãy ghi lời thích vào sau dấu % cho tất dòng lệnh đoạn lệnh -7- Phần II: Lý thuyết tín hiệu clear; clf; t=-5:0.001:5; A=input('Nhap vao bien A=') omega=input('Nhap vao omega=') pha=input('Nhap vao pha ban dau =') if((omegapi)), error('Tan so goc khong hop le'); end; if ((pha2*pi)), error('Pha ban dau khong hop le'); end x=A*sin(omega*t+pha) plot(t,x) title('Tin hieu dieu hoa'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); grid on; c Vẽ lại đồ thị tín hiệu lên hình 4.6 với A = 1, ω0 = π φ = Hình 4.6 Câu a Viết chương trình Matlab dùng hàm square sawtooth để tạo dãy xung vuông, dãy xung cưa Các thông số sau thay đổi người lập trình: tần số (f), biên độ (A) -8- Phần II: Lý thuyết tín hiệu clear; clf; t=-0.05:0.0001:0.05; f=input('Nhap tan so f=') A=input('Nhap bien A=') x=A*square(2*pi*f*t) y=A*sawtooth(2*pi*f*t) subplot(2,1,1) plot(t,x) title('Xung Vuong'); xlabel('Thoi gian lien tuc t') ylabel('Bien do') subplot(2,1,2) plot(t,y) title('Xung rang cua'); xlabel('Thoi gian lien tuc t') ylabel('Bien do') grid on; b Hãy ghi lời thích vào sau dấu % cho tất dòng lệnh đoạn lệnh clear; clf; t=-0.05:0.0001:0.05; f=input('Nhap tan so f='); A=input('Nhap bien A='); x=A*square(2*pi*f*t); %Tao xung Vuong y=A*sawtooth(2*pi*f*t); % Tao Xung Rang Cua subplot(2,1,1); plot(t,x); title('Xung Vuong'); xlabel('Thoi gian lien tuc t'); ylabel('Bien do'); subplot(2,1,2); plot(t,y); title('Xung rang cua'); % ten tieu de cua tin hieu xlabel('Thoi gian lien tuc t'); % giai thich cho truc x ylabel('Bien do'); %giai thich cho truc y grid on; c Vẽ lại đồ thị dãy xung lên hình 4.7 với f = 50Hz A = -9- Phần II: Lý thuyết tín hiệu Hình 4.7 Câu a Viết chương trình Matlab mơ tín hiệu có phương trình tổng hợp sau đây: f(t) t(u(t)  u(t  2))  (6  2t)(u(t  2)  u(t  3)) Với  t 5 clear; % xoa tat ca cac bien bo nho Matlab clf; t=-1:0.001:5; ft=t.*(heaviside(t)-heaviside(t-2))+(6-2*t).*(heaviside(t-2)-heaviside(t-3)) plot(t,ft) title('Tin hieu F(t)'); % ten tieu de cua tin hieu xlabel('Thoi gian lien tuc t') %truc x ylabel('Bien do') %truc y grid on; b Vẽ lại đồ thị phương trình tổng hợp lên hình 4.8 -10- Phần II: Lý thuyết tín hiệu Hình 4.8 -11- ... t=-1:0.001:5; ft=t.*(heaviside(t)-heaviside(t -2) )+(6 -2* t).*(heaviside(t -2) -heaviside(t-3)) plot(t,ft) title( ''Tin hieu F(t)''); % ten tieu de cua tin hieu xlabel(''Thoi gian lien tuc t'') %truc x ylabel(''Bien... thời gian hình 4. 4 -5- Phần II: Lý thuyết tín hiệu Hình 4. 4 c Hãy sửa lại chương trình Matlab để mơ tín hiệu xung vng xung tam giác dịch chuyển phía trái đơn vị clear; clf; t= -2: 0.001 :2; x=1*rectpuls(t+1,1)... t= -2: 0.001 :2; % tạo vector thoi gian x=1*rectpuls(t,1) y=1*tripuls(t ,2) subplot (2, 1,1) plot(t,x) title(''Xung vuong co ban''); xlabel(''Thoi gian lien tuc t''); ylabel(''Bien do''); subplot (2, 1 ,2)