encode.m

来自「采用互补序列降低ＯＦＤＭ的ＰＡＲ方法」· M 代码 · 共 31 行

%function par=encode(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;
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);
    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));