Giải: Ta sẽ sử dụng công thức (6.26) trong sách "Mô phỏng hệ thống truyền thông" của PGS.TS. Võ Nguyễn Quốc Bảo cho kết quả lý thuyết. Chúng ta sẽ chọn số bit mô phỏng gấp 30 lần nghịch đảo của xác suất dừng lý thuyết. Do số lần thử là một số nguyên dương, nên chúng ta dùng hàm ceil() để làm tròn số.
==================== CODE MATLAB ====================
% Tham số kênh truyền và hệ thốngN = 3;
SNRdB = 0:2:20;
SNR = 10.^(SNRdB/10);
gth = 3;
lambda = 1;
% Ly thuyet
g = lambda*SNR;
OP_a = (1 - exp(-gth./g)).^N;
% Mo phong
Flen = ceil(30./OP_a);
for idx = 1:length(SNRdB)
h = sqrt(SNR(idx).*lambda/2).*(randn(N,Flen(idx)) + 1i*randn(N,Flen(idx)));
gmax = max(abs(h).^2,[],1);
OP_s(idx) = sum(gmax < gth)./Flen(idx);
end
semilogy(SNRdB,OP_a,'-',SNRdB,OP_s,'o');
xlabel('Ty le tin hieu tren nhieu trung bình theo dB');
ylabel('Xac suat dung');
legend('Ly thuyet','Mo phong');
==================== KẾT QUẢ====================