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;

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 :

s_m=exp(j*theta_m);

for i=1:length(s)
d=abs(s_m-ns(i));
md=min(abs(s_m-ns(i)));
if md==d(1);
r(2*i-1)=0;
r(2*i)=0;
elseif md==d(2);
r(2*i-1)=0;
r(2*i)=1;
elseif md==d(3);
r(2*i-1)=1;
r(2*i)=1;
elseif md==d(4);
r(2*i-1)=1;
r(2*i)=0;
end
end
c=0;
for i=1:length(x)
if r(i)~=x(i);
c=c+1;
end
end
y=c;
Kết quả thực hiện:
4.2 So sánh tỷ lệ bit lỗi thu được bằng mô phỏng so với lý thuyết
 Chương trình:
clear;