clear
echo off
n=7; k=3;                        % Codeword and message word lengths
     m=3;                             % Number of bits per symbol
       
N=1e4;             %仿真的總包數
M=2^3;            %8進制的調制
No=0.030;           %噪聲的功率譜密度
for j=1:length(No)
    for i=1:N                      %多進制信源
        temp=rand;
        msource(i)=floor(M*temp);
    end
    [primsg,t]=decommsg(msource,k);
    msg  = gf(primsg,m);   % Two k-symbol message words
    encodout = rsenc(msg,n,k);  %%% genarate source
    channelin=reshape(encodout',1,n*t);
    channelout=pskmoto(channelin,No(j));
    codein=gf((reshape(channelout,n,t))',m);
    [codeout,cnumerr] = rsdec(codein,n,k); % Decoding failure : cnumerr(3) is -1
    receivetemp=codeout.x;
    receive=reshape(receivetemp',1,k*t);
    
    ber(j)=0;
    for i=1:N
        if(msource(i)~=receive(i))
            ber(j)=ber(j)+1;
        end
    end
    ber(j);
    p(j)=ber(j)/N;
end
ber;
p=ber/N