clear all
clc
arph_set=[0.1,0.5,1];
beta_set=[0.1,0.5,1];
P=20;%transmit power
N0=1;%spatial power spectral density N0
z0=4;% user location
M=(1:13);% define the variation range of the number of small LIS-Units M
A=16;% the surface area of LIS
N_=zeros(length(arph_set),length(M));
for i=1:length(arph_set)
arph=arph_set(i);
beta=beta_set(i);
% upper-bounded of noise density
N_(i,:)=N0+1./M.^(2.*beta).*4.^(beta-1).*P.*arph.*A.^(2.*beta)./((beta+1).*z0^2.*pi.^(beta+1));
end
figure(1)
hold on
sizeType={‘^’,’<','*','o','s'};
Color = linspecer(size(N_,1));
for i=1:size(N_,1)
plot(M,N_(i,:),'-','LineWidth',1.5,'color',Color(i,:),'marker',sizeType{i});
hold on
end
xlabel('M')
ylabel('N_s [dB]')
legend('w. HWI,\alpha=\beta=0.1','w. HWI,\alpha=\beta=0.5','w. HWI,\alpha=\beta=1','Location','northeast');
grid on