Báo cáo thí nghiệm Thông tin số - ĐTVT Bách Khoa Hà Nội k57

Tài liệu Báo cáo thí nghiệm thông tin số dành cho các sinh viên Điện tử viễn thông tham khảo. Hy vọng tài liệu hữu ích phục vụ nhu cầu học tập và nghiên cứu của các bạn.
Nội dung trích xuất từ tài liệu:
Báo cáo thí nghiệm Thông tin số - ĐTVT Bách Khoa Hà Nội k57 Thông tin số Hà Nội 2014 Bao cao thı́ nghiem Thong tin số ́ ́ ̣̂ ̂ Mã học phần : ET3250 Họ tên sinh viên : Lê Hồng Thái Mã số sinh viên : 20112207 Lớp : ĐTTT-03 K56 Ngày làm thí nghiệm : 26/02 /2014 Ngày nộp báo cáo : 12/3/2014 Thông tin số Hà Nội 2014 Bài 1: QUÁ TRÌNH NgẪU nhiên của tín hiệu 1.1 code x = -5:0.01:5; Px = exp(-x.^2/2)/sqrt(2*pi); plot(x,Px,'r'); title('Tin hieu ngau nhien Phan phoi chuan'); xlabel('x'); ylabel('Pe'); 1.2 code len = 100000; % Do dai cua qua trinh ngau nhien x = randn(1,len); % Tao qua trinh ngau nhien (phan phoi chuan) step = 0.05; % khoang cach lay mau 0.05 k = -5:step:5; % khoang xet -5 -> 5, cach deu 0.1 Px = hist(x,k)/len/step; % tinh xac suat ngau nhien Thông tin số Hà Nội 2014 stem (k, Px); % ve do thi roi rac cua QT ngau nhien Px_lithuyet = exp(-k.^2/2)/ sqrt(2*pi); % tinh xac suat phan bo theo ly thuyet hold on; % ve 2 do thi tren cung truc toa do plot (k, Px_lithuyet, 'r' ); % do thi ly thuyet title (' phan bo xac suat Gauss '); % tieu de xlabel(' x ' ); % truc hoanh ylabel(' Px '); % truc tung legend('Ly thuyet', 'Mo phong' ); % chu thich hold off; % tat che do 2 do thi. Bài 2: Lượng Tử Hóa Tuyến tính Hàm: lquan function[indx qy ] = lquan(x,xmin,xmax,nbit) nlevel = 2^nbit ; % so muc luong tu hoa q = (xmax - xmin )/nlevel; % Buoc luong tu [ indx qy ] = quantiz ( x,xmin + q:q:xmax- q, xmin +q/ 2:q:xmax -q/2); 2.1 code Thông tin số Hà Nội 2014 a= [ -0.2 1.9 3.2 -2.5]; % vidu ham lquan [ indx qy ] = lquan (a,-4,4,2); % vidu ham lquan xs= rand(1,5) *2-1; % lay 5 mau bang lenh random [xi xq] = lquan (xs,-1,1,3); % su dung ham lquan 2.2 code t = 0:0.01:20; xt = sin (rand() + t ).* cos (rand() *t); % Tin hieu vao [inx xqt] = lquan (xt, -1, 1,randint(1,1,3) +2 ); % xqt tin hieu da duoc luong tu hoa hold on; % mo ve 2 do thi plot(t,xt); % do thi tin hieu vao plot(t, xqt,'r'); % do thi tin hieu da luong tu hoa title ('Luong Tu Hoa Tuyen Tinh'); % tieu de xlabel('t'); % nhan truc x ylabel('xt xqt'); % nhan truc y legend('Tin hieu dau','Tin hieu qua Luong Tu'); % giai thich tung duong tren do thi hold off; Bài 3: tập âm lượng tử trong kĩ thuật Lượng Tử Hóa tuyến tính Thông tin số Hà Nội 2014 3.1 code N =1000; x = 2*rand(1,N) -1; % x phan bo deu tu -1 den 1 nbit = 1 :10; % so bit luong tu hoa tu 1 den 10 SNqR = zeros (size (nbit)); % khoi tao mang SNqR chua ket qua SNqR_lt = 6.02 *nbit; % mang SNqR_lythuyet tinh theo Ly thuyet Ps = sum (x.^2)/N; % Cong suat tin hieu x % tinh gia tri SNqR ung voi moi gia tri nbit n = [1 2 3 4 5 6 7 8 9 10]; for i=1:size(n,2) % tra ve so cot n [xi xq] = lquan (x,-1,1,nbit(i) ); % Luong tu hoa voi ket qua dua vao xq eq = x - xq; % Tinh sai so Pq= sum (eq.^2)/N; % Cong suat tap am luong tu SNqR(i) = 10 .* log (Ps ./Pq) ; % ti so nhieu %end end plot (nbit,SNqR ,'r',nbit,SNqR_lt); % do thi ly thuyet title('Luong tu hoa Tuyen Tinh'); xlabel('nbit'); ylabel('SNqR'); legend('SNqR mo Phong', 'SNqR Ly Thuyet') 3.2 code Thông tin số Hà Nội 2014 N =1000; x = sin(linspace(1,5,N)); % x phan bo deu tu -1 den 1 nbit = 1 :10; % so bit luong tu hoa tu 1 den 10 SNqR = zeros (size (nbit)); % khoi tao mang SNqR chua ket qua SNqR_lt = 6.02 *nbit; % mang SNqR_lythuyet tinh theo Ly thuyet Ps = sum (x.^2)/N; % Cong suat tin hieu x % tinh gia tri SNqR ung voi moi gia tri nbit n = [1 2 3 4 5 6 7 8 9 10]; for i=1:size(n,2) % tra ve so cot n [xi xq] = lquan (x,-1,1,nbit(i) ); % Luong tu hoa voi ket qua dua vao xq eq = x - xq; % Tinh sai so Pq= sum (eq.^2)/N; % Cong suat tap am luong tu SNqR(i) = 10 .* log (Ps ./Pq) ; % ti so nhieu %end end plot (nbit,SNqR ,'r',nbit,SNqR_lt); % do thi ly thuyet title('Luong tu hoa Tuyen Tinh'); xlabel('nbit'); ylabel('SNqR'); legend('SNqR mo Phong', 'SNqR Ly Thuyet') Bài 4: Mật độ phổ năng lượng và hàm tu tương quan của tín hiệu 4.1 code Thông tin số Hà Nội 2014 L =500; x1 = randn(1,L); % Tin hieu ngau gien roi rac phan phoi Chuan x2 = linspace(-1,1,L); % Tin hieu co bien do tang dan x3 = sin(linspace(-10,10,L)); % Tin hieu sin cua t/h co bien do tang dan [acorr_x1 n1] = xcorr(x1); % Ham tu tuong quan cua t/h x1 [acorr_x2 n2] = xcorr(x2); % Ham tu tuong quan cua t/h x2 [acorr_x3 n3] = xcorr(x3); % Ham tu tuong quan cua t/h x3 plot( n1, acorr_x1,'r',n2, acorr_x2,'g',n3, acorr_x3); title('Ham tuong Quan'); % Tieu de xlabel('n'); % Truc hoanh ylabel('Rxx'); % Truc tung legend('x1 Ngau Nhien', 'x2 Bien Do Tang Dan Deu', ' x3 Tin Hieu sin(cua TH bien do tang dan)'); grid on % bat luoi do thi 4.2 code L = 50; ...