clear all clc SNR_dB = 40; Ps = 3; alpha = 0.1:0.1:0.8; eta = 0.8; R = 1; M = 2; N = 4; xR = 0.5; PL = 1; %LV_MP22(SNR_dB,Ps,alpha,eta,R,M,N,xR,PL)%tinh OP LV_LT22(SNR_dB,Ps,alpha,eta,R,M,N,xR,PL)%tinh OP Mô phỏng function LV_MP22(SNR_dB,Ps,alpha,eta,R,M,N,xR,PL)%tinh OP %PL = 1; Nt = 10; Pth = 20; LSR = (xR.PL); LRD = ((1xR).PL); %gamath = 2(2R(1alpha))1; SNRW = 10.(SNR_dB.10); sobit = 105; N0 = Ps.SNRW; %fhi = 0.2exp(2etaalpha.(1alpha)) %fhi = 2etaalpha.(1alpha); SNRRDD = 0; OP_MP = zeros(1,length(alpha)); for aa = 1 : length(alpha) alpha1=alpha(aa); gamath = 2.(2.R.(1alpha1))1; fhi = 2etaalpha1.(1alpha1); aa for num = 1 : sobit for n = 1:N h1 = sqrt(LSR)6.2(randn + jrandn); end hSRR = max(abs(h1).2,,1); for m = 1:M h2 = sqrt(LRD)6.2(randn + jrandn); hRDD = max(abs(h2).2,,1); if (PshSRR