Throughput performance of OFDMA

2 次查看(过去 30 天)
for the below mentioned code i want to plot the number of subcarriers in x axis and throughput on yaxis.but for me its showing the number of users in x axis.can anyone help me to solve this issue.
Bmax=2000000;
noise=1e-9;
p_fix=0.01;
N_UE=[ 10 20 30 40 50];
N_SC=[ 60 70 80 90 100];
for r = 1:length(N_SC)
for t= 1 :length(N_UE)
Xmax=1; % x-coordinate for base station
Ymax=1; % y-coordinate for base station
radius=5; %coverage radius
xsd=ones(1,N_UE(t))*Xmax;
ysd=ones(1,N_UE(t))*Ymax;
rrx=radius*sqrt(rand(1,N_UE(t)));
thetarx=2*pi*rand(1,N_UE(t));
ang=0:0.01:2*pi;
xrx=xsd+rrx.*cos(thetarx); %random position of receivers-->x-coordinate
yrx=ysd+rrx.*sin(thetarx); %random position of receivers-->x-coordinate
th = 0:pi/100:2*pi;
xunit = radius * cos(th) + Xmax;
yunit = radius * sin(th) + Ymax;
ratio=0.1;
c=2*pi*ratio;
for j=1:N_UE(t)
temp=[];
for n=1:N_SC(r)
nc=(n-1)*c;
entry=(1+i*nc)/(1+nc^2);
temp=[temp entry];
end
Rchan=toeplitz(temp);
[E,L]=eig(Rchan);
SQRTL=sqrt(L);
A=E*SQRTL/sqrt(2);
GG=randn(N_SC(r), 1)+i*randn(N_SC(r),1);
H=A*GG;
HH=abs(H);
as=4;
bs=10;
Lsigma=as+(bs-as)*rand(1,N_UE(t));
Lshadow=10.^(randn(1,N_UE(t)).*(Lsigma/10));
ae=3;
be=5;
Lexp=ae+(be-ae)*rand(1,N_UE(t));
for k=1:N_SC(r)
for j=1:N_UE(t)
dist2(j)=sqrt((xrx(j)-xsd(j))^2+(yrx(j)-ysd(j))^2);
G(j,k)=(F(j,k)*Lshadow(j))/dist2(j).^Lexp(j);
end
end
siz = size(G);
idx = sub2ind(siz, randi([1,N_UE(t)], 1, siz(2)), 1:siz(2));
C = zeros(siz);
C(idx) = G(idx)
throughput =((Bmax.*log(1+((p_fix).*C))./noise));
overall_throughput = sum(sum((Bmax.*log(1+((p_fix).*C))./noise)));
output(r,t)=overall_throughput;
figure
plot(throughput,'b','linewidth',2);
legend(['No.of users=' num2str(N_UE(t))])
xlabel('No of subcarriers')
ylabel('throughput')
axis square
grid on
end
end
figure
plot(output(:,1), 'color', 'r' ,'linewidth',2);hold on
plot(output(:,2), 'color', 'b' ,'linewidth',2);hold on
plot(output(:,3), 'color', 'm' ,'linewidth',2);hold on
plot(output(:,4), 'color', 'y' ,'linewidth',2);hold on
plot(output(:,5), 'color', 'g' ,'linewidth',2);hold on
xlabel('No of subcarriers [x10]')
ylabel('Overall Throughput')
legend( 'UE1','UE2','UE3','UE4','UE5','Location','best')
grid on
title('Overall throughput performance')

回答(1 个)

Krishan Srimal
Krishan Srimal 2018-1-29
This is gives me an error. what would be the reason pls?
Attempt to execute SCRIPT throughput as a function: C:\Users\user\Documents\throughput.m
Error in throughput (line 4) plot(throughput,'b','linewidth',2);
  3 个评论
Ubaid Ullah
Ubaid Ullah 2019-2-3
Undefined function or variable F
Error in example (line 48)
G(j,k)=(F(j,k)*Lshadow(j))/dist2(j).^Lexp(j);
Walter Roberson
Walter Roberson 2019-2-3
Ubaid Ullah you are correct, that code does not define F or initialize F, or document what the purpose of F is. The author of the code has not posted anything for about 6 months, so I do not know if they will return to clarify it.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Stability Analysis 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by