BÁO CÁO THÍ NGHIỆM MÔN HỌC
THÔNG TIN SỐ
Mã học phần : ET3250
Sinh viên thực hiện:
Mã số sinh viên:
Lớp: Điện tử truyền thông 07-K56
Nhóm thí nghiệm:
1
Bài 1: QUÁ TRÌNH NGẪU NHIÊN CỦA TÍN HIỆU
Bài 1.1:
%
% BAI 1: QUA TRINH NGAU NHIEN CUA TIN HIEU
% Bai 11: Viet chuong trinh ve ham phan phoi chuan
%
x=-5:0.1:5; % Khoi tao mang x gom 100 phan tu
Px=(1/sqrt(2*pi)*exp(-x.^2/2)); %Gia tri cua ham P(x)
plot(x,Px,'-o');
xlabel('x'); % Tieu de cua truc Ox
ylabel('P(x)'); % Tieu de cua truc Oy
title('Ham phan bo xac suat cua phan phoi chuan');
Figure:
2
Bài 1.2: Tạo 1 quá trình ngẫu nhiên với phân bố chuẩn với độ dài 100000 phần tử
%
%Bai 1.2: Qua trinh ngau nhien voi phan phoi chuan voi do dai 100000 ptu
%

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 5, cach deu 0.1

[inx xqt] = lquan(xt,-1,1,randint(1,1,3)+2);
plot(t,xt,'b',t,xqt,'r');
grid on;
title('Do thi cua tin hieu xt va xqt ');
xlabel('t');
legend('x(t)','xqt(t)');
hold off;
Figure:
4
Bài 3: TẠP ÂM LƯỢNG TỬ TRONG KỸ THUẬT LƯỢNG TỬ HÓA TUYẾN TÍNH
Bài 3.1:
%
% Bai 3.1: Luong tu hoa tuyen tin mot t/h ngay nhien x gom 1000 mau phan bo
% deu trong dai -1:1. So sanh voi cong thuc ly thuyet, ty so SNqR trong
% truong hop t/h vao co phan bo deu dc tin nhu sau: SNqR=6.02n
%
N=1000;
x= 2*rand(1,N)-1; % x phân bo deu tu -1 den 1
nbit= 1:10; % So bit luon tu tu 1 den 10
SNqR=zeros(size(nbit)); % Khoi tao mang SNqR chua ket qua
SNqR_lt= 6.02*nbit; % Mang chua SNqR tinh theo ly thuyet
Ps=sum(x.^2)/N; % Cong suat tin hieu x
% Tinh gia tri SNqR ung voi moi gia tri cua n
for i=1:length(nbit) % Size(n,2) tra ve so cot cua n
[xi xq]=lquan(x,-1,1,nbit(i)); %LTH x voi so bit nbit(i) luu kq vao xq
eq=x-xq; % Tinh sai so
Pq=sum(eq.^2)/N; % Cong suat tap am luong tu Pq
SNqR(i)=10*log10(Ps/Pq); % Tinh SNqR(i) ra deciben
end
plot(nbit,SNqR,'b',nbit,SNqR_lt,'r-');

hold off;
6
Bài 4: MẬT ĐỘ PHỔ NĂNG LƯỢNG VÀ HÀM TỰ TƯƠNG QUAN CỦA TÍN HIỆU
Bài 4.1:
%
% Bai 4.1: Tinh va ve do thi ham mat do pho nang luong cua tin hieu ngau
% nhien x voi do dai L=500 mau
%
L=500;
x=randn(1,L);
acorr_x = xcorr(x); % Ham tu tuong quan cua x
n=-(L-1):1:(L-1);
plot(n,acorr_x);
title('Ham tu tuong quan');
xlabel('n');
ylabel('acorr_x');
hold on;
% Tin hieu co bien do tang dan
x=linspace(-1,1,L);
acorr_x=xcorr(x);
n=-(L-1):1:(L-1);
plot(n,acorr_x,'k ');
hold on;
% Tin hieu since
x=sin(linspace(-10,10,L));
acorr_x=xcorr(x);
n=-(L-1):1:(L-1);
plot(n,acorr_x,'r.');
hold on;
7

Bài 5.1:
%
% BAI 5: MA DUONG DAY NRZ
% Bai 5.1: Mo phong he thong truyen so dan bang tan co so su dung phuong
% phap ma hoa NRZ luong cuc
%
len=100000; % Do dai dong bit mo phong
SNR_db=0:2:8; % Tao vector SNR_db = 0 2 4 6 8
SNR= 10.^(SNR_db/10); % Doi SNR tu decibel sang lan
bsignal=randint(1,len); % Dong bit ngau nhien co do dai len
NRZ_signal=bsignal*2-1; % Bien doi dong bit 0 1 sang -1 1
N0=1./SNR; % Cong suat tap am
for i=1:length(SNR_db)
noise=sqrt(N0(i))*randn(1,len);
r_signal=NRZ_signal + noise;
NRZ_decoder=sign(r_signal);
[Num,Rat]=symerr(NRZ_decoder,NRZ_signal);
Pe(i)=Rat;
end
plot(SNR_db,Pe,'ko ');
title('Ty le bit loi');
xlabel('SNR_db');
ylabel('Pe');

Figure:
9
Bài 5.2:
%
% BAI 5: MA DUONG DAY NRZ
% Bai 5.2: Ve do thi cua ham xac suat loi ly thuyet va mo phong tren cung

% Bai 6.1: Bieu do chom sao cua tin hieu d/c QPSK va ma hoa theo ma Gray.
% Y/c: Tao tin hieu nhi phan ngau nhien chieu dai 50000.Thuc hien dieu che
% QPSK tin hieu nay
%
len=50000;
bit=randint(1,len); %tin hieu vao
for i=1:2:len; % tao vong lap, step = 2
switch bit(i)
case 0
if bit(i+1)==0
qpsk_signal((i+1)/2)=exp(j*3*pi/4);
else
qpsk_signal((i+1)/2)=exp(j*5*pi/4);
end
case 1
if bit(i+1)==0
qpsk_signal((i+1)/2)=exp(j*pi/4);
else
qpsk_signal((i+1)/2)=exp(j*7*pi/4);
end
end
end
plot(qpsk_signal,'o');
hold on;
t=0:0.01:2*pi;
plot(exp(j*t),'r '); % Ve duong tron don vi
title('Bieu do chom sao cua tin hieu dieu che QPSK');
xlabel('I');
ylabel('Q');
hold off;

ylabel('Q');
hold on;
Es=std(qpsk_signal).^2; %Tinh nang luong
% // Ti le SNR = 6 dB //
N0=Es/(10^(6/10)); % Cong suat tap am cua nhieu voi SNR = 6 dB
noise=sqrt(N0/2).*(randn(1,length(qpsk_signal))+1i.*randn(1,length(qpsk_signal)));
qpsk=qpsk_signal + noise; % Tin hieu khi co nhieu
subplot(3,1,1); plot(qpsk,'o');
title('SNR=6dB');
grid on;axis auto;
xlabel('I');
ylabel('Q');
% // Ti le SNR = 3 dB //
N1=Es/(10^(3/10)); % Cong suat tap am cua nhieu voi SNR = 3 dB
noise1=sqrt(N1/2).*(randn(1,length(qpsk_signal))+1i.*randn(1,length(qpsk_signal)));
qpsk1=qpsk_signal + noise1;
subplot(3,1,2);
plot(qpsk1,'o');
title('SNR=3dB');
grid on;
axis auto;
xlabel('I');
ylabel('Q');
% // Ti le SNR = 0 dB //
N2=Es/(10^(0/10)); % Cong suat tap am cua nhieu voi SNR = 0 dB
noise2=sqrt(N2/2).*(randn(1,length(qpsk_signal))+1i.*randn(1,length(qpsk_signal)));
qpsk2=qpsk_signal + noise2;
subplot(3,1,3);
plot(qpsk2,'o');
title('SNR=0dB');


end
BER(i)=n/len;
end
plot(SNR_db,BER,'r-');
title('Mo phong dieu che QPSK qua kenh nhieu Gauss');
xlabel('SNR_db');
ylabel('BER');
Figure:
17
Bài 8:XÁC SUẤT LỖI BIT TRONG ĐIỀU CHẾ QPSK
%
% BAI 8: XAC SUAT LOI BIT TRONG DIEU CHE QPSK
% Yeu cau: Tinh ti le loi bit ly thuyet pb cua he thong va ve do thi
% SNR=0,2,4,6,8 dB. So sanh ket qua ly thuyet nay voi ket qua mo phong o
% bai 5
%

len=100000;
x = randint(1,len);
M=4;
qpsk=qammod(x,M);
SNR_db=0:2:8; %khoi tao vecto SNR_db 0 2 4 6 8
SNR=10.^(SNR_db/10);
N0=1./SNR; %tao vecto N0
for i=1:length(SNR_db)
BERBPSK(i)=(1./2)*erfc(sqrt(SNR(i)/2));
BERQPSK(i)=(1./2)*erfc(SNR(i)/sqrt(2));
end
plot(SNR_db,BERBPSK,'b* ',SNR_db,BERQPSK,'mo ');