MATLAB: Anybody can help me? It can not plot ? Error is:” for” may not be aligned with ít matching END. Thanks so much!

ergodic capacityergodic capacity at d1noma

%CACULATE ERGODIC CAPACITY AT D1 IN NOMA SCHEME 
clc;
a1=0.9;
a2=0.1;
k=3;
SNRdB=0:5:30;% gama 0 (dB)
N=length(SNRdB);
SNR=10.^(SNRdB/10);
omegasd1=1;
omegaSD1=omegasd1.*SNR;
omegasr1=1;
omegasr2=1;
omegasr3=1;
omegar1d1=1;
omegar2d1=1;
omegar3d1=1;
for i=1:length(SNRdB)
    for ii=1:N
      % gama S--->Rk
      hsr1=sqrt(omegasr1/2)*(randn(1,N)+1i*randn(1,N));
      hsr1m=abs(hsr1).^2;
      gamaSR1=hsr1m.*SNR(i);    
      hsr2=sqrt(omegasr2/2)*(randn(1,N)+1i*randn(1,N));
      hsr2m=abs(hsr2).^2;
      gamaSR2=hsr2m.*SNR(i);    
      hsr3=sqrt(omegasr3/2)*(randn(1,N)+1i*randn(1,N));
      hsr3m=abs(hsr3).^2;
      gamaSR3=hsr3m.*SNR(i);
      % gama Rk-->D1
      hr1d1=sqrt(omegar1d1/2)*(randn(1,N)+1i*randn(1,N));
      hr1d1m=abs(hr1d1).^2;
      gamaR1D1=SNR(i).* hr1d1m;       
      hr2d1=sqrt(omegar2d1/2)*(randn(1,N)+1i*randn(1,N));
      hr2d1m=abs(hr2d1).^2;
      gamaR2D1=SNR(i).* hr2d1m; 
      hr3d1=sqrt(omegar3d1/2)*(randn(1,N)+1i*randn(1,N));
      hr3d1m=abs(hr3d1).^2;
      gamaR3D1=SNR(i).*hr3d1m;    
      %he so kenh truyen S->D1 
      hsd1=sqrt(omegasd1/2)*(randn(1,N)+1i*randn(1,N));
      hsd1m=abs(hsd1).^2;
      gamasd1=hsd1m .*SNR(i);  
      E1=expint(1./(a2*omegaSD1));
      E2=expint(1./omegaSD1);
      end     
      %dung luong C1 ergodic
      X1=gamasd1;
      Yk=gamaSR3;
      Z1k=gamaR3D1;
      X=(a1.*X1)./(a2.*X1+1);
      Y=(a1.*Yk.*Z1k)./(a2.*Yk.*Z1k+Yk+Z1k+1);
      MAX=max(X,Y);
      end
      C1=mean(0.5.*log2(1+MAX));
plot(SNRdB,C1,'r-');
hold on;
title('Dung luong Ergodic tai D1')

Best Answer

clc;
a1=0.9;
a2=0.1;
k=3;
SNRdB=0:5:30;% gama 0 (dB)
N=length(SNRdB);
SNR=10.^(SNRdB/10);
omegasd1=1;
omegaSD1=omegasd1.*SNR;
omegasr1=1;
omegasr2=1;
omegasr3=1;
omegar1d1=1;
omegar2d1=1;
omegar3d1=1;
for i=1:length(SNRdB)
    for ii=1:N
        % gama S--->Rk
        hsr1=sqrt(omegasr1/2)*(randn(1,N)+1i*randn(1,N));
        hsr1m=abs(hsr1).^2;
        gamaSR1=hsr1m.*SNR(i);
        hsr2=sqrt(omegasr2/2)*(randn(1,N)+1i*randn(1,N));
        hsr2m=abs(hsr2).^2;
        gamaSR2=hsr2m.*SNR(i);
        hsr3=sqrt(omegasr3/2)*(randn(1,N)+1i*randn(1,N));
        hsr3m=abs(hsr3).^2;
        gamaSR3=hsr3m.*SNR(i);
        % gama Rk-->D1
        hr1d1=sqrt(omegar1d1/2)*(randn(1,N)+1i*randn(1,N));
        hr1d1m=abs(hr1d1).^2;
        gamaR1D1=SNR(i).* hr1d1m;
        hr2d1=sqrt(omegar2d1/2)*(randn(1,N)+1i*randn(1,N));
        hr2d1m=abs(hr2d1).^2;
        gamaR2D1=SNR(i).* hr2d1m;
        hr3d1=sqrt(omegar3d1/2)*(randn(1,N)+1i*randn(1,N));
        hr3d1m=abs(hr3d1).^2;
        gamaR3D1=SNR(i).*hr3d1m;
        %he so kenh truyen S->D1
        hsd1=sqrt(omegasd1/2)*(randn(1,N)+1i*randn(1,N));
        hsd1m=abs(hsd1).^2;
        gamasd1=hsd1m .*SNR(i);
        E1=expint(1./(a2*omegaSD1));
        E2=expint(1./omegaSD1);
    end
    %dung luong C1 ergodic
    X1=gamasd1;
    Yk=gamaSR3;
    Z1k=gamaR3D1;
    X=(a1.*X1)./(a2.*X1+1);
    Y=(a1.*Yk.*Z1k)./(a2.*Yk.*Z1k+Yk+Z1k+1);
    MAX=max(X,Y);
end
C1=mean(0.5.*log2(1+MAX));
plot(SNRdB,MAX,'r-');
hold on;
title('Dung luong Ergodic tai D1')
Related Question