Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 13 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
13
Dung lượng
4,35 MB
Nội dung
BÁOCÁOTHÍNGHIỆM MƠN HỌC THƠNGTINSỐ Sinh viên thực hiện: Lớp: Mã sinh viên: BÀI SỐ 1: Q TRÌNH NGẪU NHIÊN CỦA TÍN HIỆU Bài 1.1 Vẽ hàm phân bố xác suất phân phối chuẩn Mã nguồn: x=-5:.1:5; % Khoi tao mang x gom Px=(1/sqrt(2*pi))*exp(-x.^2/2); plot(x,Px); % Ve thi title('Ham phan bo xac suat cua phan thi xlabel('x'); % Ten truc hoanh ylabel('P(x)'); % Ten truc tung Figure: 100 phan tu % Gia tri cua ham P(x) ham phan phoi chuan phoi chuan'); % tieu de cua Bài 1.2 Vẽ hàm phân bố xác suất phân phối chuẩn lý thuyết mô Mã nguồn: len=100000; % Do dai cua qua trinh ngau nhien x=randn(1,len); % Tao qua trinh ngau nhien theo phan phoi chuan step=0.1; % step = 0.1 k=-5:step:5; % khoang xet tu -5 den cach deu 0.1 px=hist(x,k)/len/step; % Xac dinh so vector stem(k,px); % Ve thi roi rac cua Px mo phong Px_lythuyet=exp(-k.^2/2)/sqrt(2*pi); % Tinh P(x) ly thuyet hold on; plot(k,Px_lythuyet); % Ve thi Px ly thuyet title('Phan bo xac suat Gauss'); % tieu de cua thi xlabel('x'); % Ten truc hoanh ylabel('P(x)'); % Ten truc tung legend('Mo phong','Ly thuyet'); % Ky hieu hold off; Figure: BÀI SỐ 2: LƯỢNG TỬ HĨA TUYẾN TÍNH Bài 2.1 Lượng tử hóa tín hiệu Mã nguồn: xs = rand(1,5)*2-1; % tin hieu ngau nhien gom mau [xi xq] = lquan(xs,-1,1,3); % luong tu hoa tin hieu Kết mô phỏng: xs = -0.6484 0.3862 -0.1578 0.7446 0.5684 xi = 6 xq = -0.6250 0.3750 -0.1250 0.6250 0.6250 Bài 2.2 Cho tín hiệu lượng tử x(t) ngẫu nhiên có biên độ -1 đến xét khoảng 0=0 if imag(r_signal((j+1)/2))>=0 r_bsignal(j) = 1; r_bsignal(j+1) = 0; else r_bsignal(j) = 1; r_bsignal(j+1) = 1; end else if imag(r_signal((j+1)/2))>=0 r_bsignal(j) = 0; r_bsignal(j+1) = 0; else r_bsignal(j) = 0; r_bsignal(j+1) = 1; end end end [n,BER(i)] = biterr(r_bsignal,bsignal);% ty le loi bit BER mo phng end % Tinh ty le loi bit BER ly thuyet Pb = 1/2*erfc(1/sqrt(2).*sqrt(SNR)); % ty le loi bit BER ly thuyet plot(SNR_db,Pb,'r'); % Ve thi Pb ly thuyet hold on; plot(SNR_db,BER,'o') % Ve thi Pb mo phong title('Ty le loi bit BER voi ty le tin hieu tren tap am SNR') xlabel('SNR_db (dB)'); ylabel('BER'); legend('ly thuyet','mo phong'); hold off; Figure: BÀI SỐ 8: MÔ PHỎNG ĐIỀU CHẾ M-QAM QUA KÊNH NHIỄU GAUSS Bài 8.1 Tính tỷ lệ lỗi bit BER kỹ thuật điều chế M-QAM với M=16,64,128 với mức độ nhiễu kênh có SNR[dB] = 0,1,…,25 Vẽ đồ thị đường BER trường hợp điều chế Mã nguồn: n_sym = 50000; % So ky tu dieu che M = [16 64 128]; % So symbol ky hieu SNR_db = 0:25; % Tao vector SNR = - 25 Decibel BER = zeros(length(M),length(SNR_db)); % BER de luu ti le loi bit for k = 1:size(M,2) % size(M,2) la so cot cua M s_stream = randi([0 M(k)-1],1,n_sym); % Tao dong bieu tuong dai n_sym s_mod=qammod(s_stream,M(k),0,'BIN'); % Dieu che M-QAM %Tinh ty le loi BER for r = 1:size(SNR_db,2) s_mod_awgn = awgn(s_mod,SNR_db(r),'measured'); % Tin hieu qua nhieu AWGN s_demod = qamdemod(s_mod_awgn,M(k),0,'BIN');% Giai dieu che M-QAM [n BER(k,r)]= biterr(s_stream,s_demod); % Tinh ti le loi bit end end semilogy(SNR_db,BER(1,:),'bo:'); % Ve thi BER ung voi M = 16 hold on; semilogy(SNR_db,BER(2,:),'ro:'); % Ve thi BER ung voi M = 64 semilogy(SNR_db,BER(3,:),'go:'); % Ve thi BER ung voi M = 128 title('Do thi BER cua M-QAM'); xlabel('SNR_dB'); ylabel('BER'); legend('16-QAM','64-QAM','128-QAM'); hold off; Figure: ... Luong tu hoa tin hieu x_uni [indx2 x_qsin]=lquan(x_sin,-1,1,nbit(i)); % Luong tu hoa tin hieu x_sin eq_uni=x_uni-x_quni; % tinh sai so luong tu hoa x_uni eq_sin=x_sin-x_qsin; % tinh sai so luong... % Tin hieu ngau nhien co bien -1 den [inx xqt]=lquan(xt,-1,1,randi(3)+1); % Tinh xqt plot(t,xt,'b',t,xqt,'r'); % Ve thi xt va xq grid on; title('Luong tu hoa tin hieu x(t)'); % tieu de cua thi. .. SNqR_sin(i)=10*log10(Ps_sin/Pq_sin); % tinh SNR_db cua x_sin end % ve thi SNqR_uni, SNqR_sin, SNqR_lt plot(nbit,SNqR_uni,'m:',nbit,SNqR_sin,'b-',nbit,SNqR_lt,'g '); title('Ty so tin hieu tren tap am SNqR'); % tieu de cua thi xlabel('nbit');