1. Trang chủ
  2. » Luận Văn - Báo Cáo

Bao cao thi nghiem thong tin so

13 220 0

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

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ÁO CÁO THÍ NGHIỆM MƠN HỌC THƠNG TIN SỐ 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');

Ngày đăng: 31/03/2018, 18:43

TỪ KHÓA LIÊN QUAN

w