encode.m

来自「ofdm中互补序列降低par的ber程序」· M 代码 · 共 33 行

function modul=encode(inf,len);
G=[1 0 0 3 0 3 3 2;0 1 0 1 0 1 0 1;0 0 1 1 0 0 1 1;0 0 0 0 1 1 1 1];
b=[0 0 0 2 0 0 2 0];
%inf=[1 3 2 0];
%len=800;
modul=zeros(1,2*len);
par=zeros(1,len);
%inf=source(len);
for k=1:len/4
    tt=inf((4*k-3):4*k);
    tt=tt*G+b;
    tt=mod(tt,4);
    for j=1:8
        if tt(j)==0
            tt(j)=1+i;
        elseif tt(j)==1
            tt(j)=-1+i;
        elseif tt(j)==2
            tt(j)=-1-i;
        else tt(j)=1-i;
        end
    end
    out=ifft(tt,8);
    modul(8*k-7:8*k)=out;
    av=mean(out.*(conj(out)));
    m=max(out.*(conj(out)));
    p=m/av;
    par(k)=10*log10(p);
end
%x=[1:len/4];
%plot(x,par(x));