Bài 1
1.1V hàm phân bó c a bi n ng u nhiên µ v i giá tr trung bình và đẽ ủ ế ẫ ớ ị ộ
l ch chu n cho tr c mệ ẩ ướ µ=0;σµ=1
Ch ng trình:ươ
function tn1
x=-4:0.05:4;
p=exp(-x.^2/2)/(sqrt(2*pi));
check=trapz(x,p);
plot(x,p);
xlabel('x','Fontsize',12);
ylabel('p(x)','fontsize',12);
K t qu th c hi n:ế ả ự ệ
1.2V hàm phân b xác su t c a bi n ng u nhiên Gauss v i tr trungẽ ố ấ ủ ế ẫ ớ ị
bình và đ l ch chu n c a 1 quá trình ng u nhiên t t oộ ệ ẩ ủ ẫ ự ạ
Ch ng trình:ươ
function tn1_2
clear;
n=1000000;
x=-4:0.05:4;
p=exp(-x.^2/2)/(sqrt(2*pi));
check=trapz(x,p);
plot(x,p,'r');
hold on;
y=randn(1,n);
m=mean(y);
variance=std(y)^2;
x2=-4:0.05:4;
c=hist(y,x2);
stem(x2,c/n/(x2(2)-x2(1)));
title('Gaussian Distributed pdf');
xlabel('x');
ylabel('p(x)');
legend('theoretical','experimetal');
hold off;
K t qu th c hi n :ế ả ự ệ
Nh n xét :Đ th c a hàm phân b Gauss nh n đ c gi ng v i đ th bài 1-1ậ ồ ị ủ ố ậ ượ ố ớ ồ ị ở
1.3Xây d ng hàm t t ng quan cho quá trình ng u nhiên bài 1.2ự ự ươ ẫ ở
Nh n xét: Hàm t t ng quan c a quá trình ng u nhiên ví d 1.2 gi ng v iậ ự ươ ủ ẫ ở ụ ố ớ
hàm m t đ ph năng l ng c a chính nóậ ộ ổ ượ ủ
Bài 3
3.1 V bi u đ chòm sao c a tín hi u đi u ch QPSKẽ ể ồ ủ ệ ề ế
Ch ng trình:ươ
clear;
x=round(rand(1,10000));
plot(x);
for i=1:2:length(x)
if x(i)==0&x(i+1)==0
s((i+1)/2)=exp(j*pi/4);
elseif x(i)==0&x(i+1)==1
s((i+1)/2)=exp(j*3*pi/4);
elseif x(i)==1&x(i+1)==1
s((i+1)/2)=exp(j*5*pi/4);
elseif x(i)==1&x(i+1)==0
s((i+1)/2)=exp(j*7*pi/4);
end
end
save ex5p1_res s x;
plot(s,'*');
hold on;
t=0:0.01:2*pi;
plot(exp(j*t),('r--'));
xlabel('\phi(t)');
ylabel('s_m');
title('the complex signel-space diagram for 4-QPSK');
K t qu th c hi n:ế ả ự ệ
3.2 Cho tín hi u đi u ch đ c qua kênh nhi u tr ng (AWGN) v bi u đệ ề ế ượ ễ ắ ẽ ể ồ
chòm sao.
Ch ng trình:ươ
clear;
load ex5p1_res;
es=var(s);
eb=es/2;
snr_db=6;
n_0=eb/10^(snr_db/10);
n=sqrt(n_0/2)*(randn(size(s))+j*randn(size(s)));
r=s+n;
plot(r,'.');
hold on;
plot(s,'.');
hold on;
t=0:0.01:2*pi;
plot('exp(j*t),]r--');
legend('s_m','s');
xlabel('I');
ylabel('Q');
title('the complex signal diagram of 4-QPSK');
K t qu th c hi n :ế ả ự ệ
TH1: SNR=6db
TH2: SNR=3dB