clear all%function [p1, p2] = fund_mac(K,ct)
M=6; K=3; snrdb=[0 : 5 : 30];
Rate1=1; Rate2=4; ct=5000;
for k = 1: length(snrdb)
snr = 10^(snrdb(k)/10); ep1 = 2^Rate1-1; ep2 = 2^Rate2-1; Ptarget = 1- exp(-2*ep1/snr);
a2 = (1+ep1/snr/log(1-Ptarget))/(1+ep1); a1=1-a2;
sum1=0; sum2=0; sum3=0; sum4=0; sum5=0; sum6=0;sum7=0; sum8=0; sum9=0;
Ptarget1 = 1- exp(-2*ep1/snr);% Ptarget2 = 1- exp(-2*ep2/snr); Ptarget3 = 1- exp(-2*ep3/snr);
a2 = (1+ep1/snr/log(1-Ptarget1))/(1+ep1); a1=1-a2;
%b2 = (1+ep2/snr/log(1-Ptarget2))/(1+ep2); b1=1-b2;
%c2 = (1+ep3/snr/log(1-Ptarget3))/(1+ep3); c1=1-c2;
for i = 1 :ct
H1=complex(sqrt(0.5)*randn(K,M),sqrt(0.5)*randn(K,M));
H2=complex(sqrt(0.5)*randn(K,M),sqrt(0.5)*randn(K,M));
[Qx,Rx] = qr(H2');
R2=Rx(1:K,:); V2 = Qx(:,1:K);
x = abs(diag(R2)).^2;
%H11 = H1*V2';
%[Q1, R1] = qr(H11);
if log2(1+a1*x(1)/(a2*x(1)+1/snr))<Rate1 %first layer: user 1 outage
sum1 = sum1 +1;
sum2 = sum2+1; sum3 = sum3+1;
else
if log2(1+a2*x(1)*snr)<Rate2 %first layer: user 2 outage
sum1 = sum1 +1;
sum2 = sum2+1; sum3 = sum3+1;
else
if log2(1+a1*x(2)/(a2*x(2)+1/snr))<Rate1 %second layer: user 1 outage
sum2 = sum2 +1;
sum3 = sum3+1;
else
if log2(1+a2*x(2)*snr)<Rate2 %second layer: user 2 outage
sum2 = sum2+1;
sum3 = sum3+1;
else
if log2(1+a1*x(3)/(a2*x(3)+1/snr))<Rate1 %third layer: user 1 outage
sum3 = sum3 +1;
else
if log2(1+a2*x(3)*snr)<Rate2 %third layer: user 2 outage
sum3 = sum3+1;
end
end
end
end
end
end
%%%%%%%%%% OMA
if log2(1+x(1)*snr)/2<Rate2 %first layer: user 2 outage
sum4 = sum4 +1;
sum5 = sum5+1; sum6 = sum6+1;
else
if log2(1+x(2)*snr)/2<Rate2 %second layer: user 2 outage
sum5 = sum5+1; sum6 = sum6+1;
else
if log2(1+x(3)*snr)/2<Rate2 %third layer: user 2 outage
sum6 = sum6+1;
end
end
end
end
p1(k) = sum1/ct; p2(k) = sum2/ct; p3(k) = sum3/ct; p4(k)=sum4/ct; p5(k) = sum5/ct; p6(k)=sum6/ct; p7(k)=sum7/ct; p8(k) = sum8/ct; p9(k)=sum9/ct;
end
semilogy(snrdb,p1,'-o',snrdb,p2,'-s',snrdb,p3,'-x',snrdb,p4,'-o',snrdb,p5,'-*',snrdb,p6,'-d')%,snrdb,p7,'-o',snrdb,p8,'-*',snrdb,p9,'-d')
jghjhj.rar （22个子文件） 
QR_NOMA_u2_policy_II.m 2KB figure_comp1.fig 16KB figure_comp_oma2.fig 16KB QR_NOMA_u2_policy_I.m 3KB test1.m 1KB NOMA_compare_policy_I_OMA.m 3KB figure_comp1_v2.fig 16KB figure_comp_oma1.fig 16KB data_com1.mat 747B data_policy_I_u1_1.mat 856B QR_NOMA_u1_policyI.m 2KB data_policy_I_u1_2.mat 871B data_policy_II_u2_2_4.mat 503B data_policy_II_u1_2.mat 824B data_policy_I_u2_1.mat 1KB QR_NOMA_u1_policyII.m 2KB policy.eps 17KB data_policy_II_u1_1.mat 817B data_com_oma2.mat 723B data_policy_II_u2_2_3.mat 495B data_policy_II_u2_2_1.mat 505B data_policy_II_u2_2_2.mat 481B
局外狗
- 粉丝: 78
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
VIP享7折，此内容立减3.57元 开通VIP
评论1