How to plot this euation?

1 次查看（过去 30 天）

Hasan Tanbir 2014-5-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/129634-how-to-plot-this-euation

Here is the code what i wrote but it does not give proper picture

clear all;
  clc;
  No=10^-6;
  Po=10;
  Rb=16000;
  Tb=1/Rb;
  Eb=Po.*Tb;
  Da=5;
  Da_db=10^(Da/10);
  Nc=32;
  Mc=4;
  L=8;
  K=1;
  M=3;
  j=0;
    Meu1=sqrt(1./(1+(No./(2.*Eb))+(2/(3.*Da_db.*Nc)).*((1+(Mc./5))*L*K-1)));
    % fact1=factorial(M-1+j);
    % fact2=factorial(j);
    % fact3=factorial(M-1);
    n=M-1+j;
    fact=factorial(n)/(factorial(j).*factorial(n-j));
    sum=0; 
    for i=1:M;
        for j=0:2;
         s(i)=fact.*(((1+Meu1)/2)^j);
         sum=sum+s(i);
         ber(i)=(((1-Meu1)/2)^M).*s(i)
        end
    end
    semilogy(Eb/No,ber)
    axis([0 10 10e-12 10e-1])
    xlabel('Eb/No,dB');
    ylabel('BER');
    grid on