qam64ys.m

来自「对数字信号进行分析」· M 代码 · 共 85 行

function [qam64ys1,qam64ys5]=qam64ys(ldpc1,tb)
runum1=length(ldpc1);
runum2=fix(runum1/6);
qam64ys2=reshape(ldpc1,runum2,6);
qam64ys3=zeros(runum2,3);
qam64ys4=zeros(runum2,3);
qam64ys3=qam64ys2(:,1:3);
qam64ys4=qam64ys2(:,4:6);
qam64ys1=zeros(1,runum2);
qam64ys5=zeros(1,runum2);
w1=0;
ii=zeros(1,3);
for i=1:1:runum2
        ii=qam64ys3(i,:);       
        if ii(1,1)==0
            if ii(1,2)==0
                if ii(1,3)==0
                    w1=-7;
                else
                    w1=-5;
                end
            else
                if ii(1,3)==0
                    w1=-1;
                else
                    w1=-3;
                end
            end
        else
            if ii(1,2)==0
                if ii(1,3)==0
                    w1=7;
                else
                    w1=5;
                end
            else
                if ii(1,3)==0
                    w1=1;
                else
                    w1=3;
                end
            end
        end
       qam64ys1(1,i)=w1;
end
t1=tb:tb:tb*runum2;
% figure
% stairs(t1,qam64ys1)

for i=1:1:runum2
        ii=qam64ys4(i,:);       
        if ii(1,1)==0
            if ii(1,2)==0
                if ii(1,3)==0
                    w1=-7;
                else
                    w1=-5;
                end
            else
                if ii(1,3)==0
                    w1=-1;
                else
                    w1=-3;
                end
            end
        else
            if ii(1,2)==0
                if ii(1,3)==0
                    w1=7;
                else
                    w1=5;
                end
            else
                if ii(1,3)==0
                    w1=1;
                else
                    w1=3;
                end
            end
        end
       qam64ys5(1,i)=w1;
end
% figure
% stairs(t1,qam64ys5)