Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 12 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
12
Dung lượng
1 MB
Nội dung
TRƯỜNG ĐẠI HỌC CƠNG NGHIỆP THÀNH PHỐ HỒ CHÍ MINH KHOA CƠNG NGHỆ ĐIỆN TỬ MƠN HỌC:THÍ NGHIỆM XỬ LÝ SỐ TÍN HIỆU BÁO CÁO THỰC HÀNH BÀI 1: LÀM QUEN MATLAB VÀ ỨNG DỤNG TRONG DSP Giảng viên: Lê Văn Hùng Sinh viên: Phạm Thị Anh Thư -19517171 VD1 n=1:10 n= 10 VD2 n=1:10 y=sin(2*pi/10*n) plot(n,y) y= 0.5878 0.9511 0.5878 -0.0000 0.9511 0.5878 0.0000 -0.5878 -0.9511 -0.9511 - VD3 clc clear all close all % T=2; K=100; bn=2; xt=sinsum(T,K,bn);%khi buoc nhay bang thi day la Fourier cua ham la xung vuong t = 0:0.01:T; plot(t,xt); legend('chuoi Fourier cua ham xung vuong') %pause bn=1; xtt=sinsum(T,K,bn);%khi buoc nhay bang thi day la Fourier cua ham la xung rang cua plot(t,xt,'r',t,xtt,'b'); legend('chuoi Fourier cua ham xung vuong','chuoi Fourier cua ham xung vuong') VD4 n=[ -3 2] m=[4 7] n' -3 n=[ -3 2] m=[4 7] n' m*n' ans = 10 n=[ -3 2] m=[4 7] m'*n -12 10 -15 10 12 -18 12 14 -21 14 VD5 n=[ -3 2] m=[4 7] 2.^n ans = 2.0000 4.0000 0.1250 4.0000 1.0000 -0.6667 1.0000 n=[ -3 2] m=[4 7] 2./n ans = 2.0000 n=[ -3 2] m=[4 7] m.*n ans = 10 -18 14 BT1 n = -10:11 x = 0.6.^(abs(n)) k = -500:500; w = (pi/500)*k; X=x*(exp(-j*pi/500)) ^ (n'*k); magX = abs(X); angX = angle(X); plot(w,X); Columns 13 through 22 0.3600 0.2160 0.0060 0.0036 0.1296 0.0778 0.0467 0.0280 0.0168 0.0101 BT2 a Tìm B= [A a], giải thích kết quả: vector hàng MATLAB tạo cách đóng tập hợp phần tử dấu ngoặc vng, sử dụng dấu cách để tách phần tử A = [1 2] a = [4 0 1] B= [A a] B= 4 0 b Tìm C=[a,A], giải thích kết quả: vector hàng MATLAB tạo cách đóng tập hợp phần tử dấu ngoặc vuông, sử dụng dấu phẩy để tách phần tử A = [1 2] a = [4 0 1] C=[a,A] C= 0 1 c Tìm D=B(2:4), giải thích kết quả tạo vector B sau tham chiếu phần tử từ đến cách viết (2:4) A = [1 2] a = [4 0 1] D=B(2:4) D= d Tìm E=[a;A], giải thích kết quả: Vector cột MATLAB tạo cách đóng tập hợp phần tử dấu ngoặc vuông, sử dụng dấu chấm phẩy để tách phần tử A = [1 2] a = [4 0 1] E=[a;A] E= 0 1 e Tìm F ma trận E có kích thước 2x3 A = [1 2] a = [4 0 1] F=(2:3) F= f Tìm G vec tơ có thành phần chuyển vị A:A’ vecto chuyển vị A A = [1 2] a = [4 0 1] G=A' g Tìm H vec tơ có thành phần tích thành phần tương ứng A a: vecto a nhân với vetơ chuyển vị A A = [1 2] a = [4 0 1] H =(a*A') H= 54 h Nhập vec tơ I có thành phần sau [2 2 2 2 A = [1 2] a = [4 0 1] I=[2.^1 2.^2 2.^4 2.^5 2.^3 2.^9 2.^4 2.^2] I= 16 32 512 16 BT3 a function [n,x]=tinhieusin(N,k,m) t=0:k/N:m; n=1:1:length(t); x=sin(2*pi*t) end clc clear all close all N=100; k=2; m=2; [n,x]=tinhieusin(N,k,m); plot(n,x,'r'); legend('tin hieu sin') title(' DO THI TIN HIEU SIN') 10 B clc clear all close all k=3; N=200; m=3; [n,x]=tinhieusin(N,k,m); t=0:k/N:m; legend('tin hieu vuong') y=sign(x); subplot(2,1,1); plot(n,x,'r'); legend('tin hieu sin') subplot(2,1,2); plot(y,'b'); 11 C clc clear all close all k=3; N=200; m=3; [n,x]=tinhieusin(N,k,m); t=0:k/N:m; y=sign(x); plot(n,x,'r',n,y,'b'); legend('tin hieu sin','tin hieu vuong') 12 ...VD1 n =1: 10 n= 10 VD2 n =1: 10 y=sin(2*pi /10 *n) plot(n,y) y= 0.5878 0.9 511 0.5878 -0.0000 0.9 511 0.5878 0.0000 -0.5878 -0.9 511 -0.9 511 - VD3 clc clear all close all % T=2; K =10 0; bn=2;... -3 2] m=[4 7] n' m*n' ans = 10 n=[ -3 2] m=[4 7] m'*n -12 10 -15 10 12 -18 12 14 - 21 14 VD5 n=[ -3 2] m=[4 7] 2.^n ans = 2.0000 4.0000 0 .12 50 4.0000 1. 0000 -0.6667 1. 0000 n=[ -3 2] m=[4 7] 2./n... ans = 10 -18 14 BT1 n = -10 :11 x = 0.6.^(abs(n)) k = -500:500; w = (pi/500)*k; X=x*(exp(-j*pi/500)) ^ (n'*k); magX = abs(X); angX = angle(X); plot(w,X); Columns 13 through 22 0.3600 0. 216 0 0.0060