mastrovito.m

对做射频方面的朋友很有用的

clc
clear all
%% Q-Matrix
m = 4; % Number of bits per symbol(Integer between 3 and 16)
pp = 19;
p = de2bi(pp); % ,'left-msb'
q(1,:) = p(1:end-1);
% q(1,:) = p(2:end);
for i=2:m-1
    q(i,1) = q(i-1,m);
    for j=2:m
        q(i,j) = q(i-1,j-1) + q(i-1,m) & q(1,j);             
    end
end
a = 7; b = 3; primpoly = 285;
%% convert input to binary form
a_bi = de2bi(a,m,'left-msb');
b_bi = de2bi(b,m,'left-msb');
for i=1:m
    f(i,1)=a_bi(i);
    for j=2:m
        if (i-j>0)
            aij = a_bi(i-j);
        else
            aij = 0;
        end
        for t=1:j-1
            (q(j-t,i) & a_bi(n-t))
        end
        f(i,j)=aij+qa;
    end
end