picmud.m

PIC呀

function [bits_out,S_out,r]=picmud(bits_in,bits_add,k)
m=length(bits_in);
n=length(bits_add(1,:));
p=floor(m/n);
%bits_finalout=zeros(k,p);  %picmud最终输出的信号
bits_out=zeros(k,p);       %picmud没有经过判定，输出的信号
S_out=zeros(k,p);          %对输入信号进行扩频得到的有多址干扰的信号
r=zeros(k,k);              %相关系数矩阵
for i=1:k
    for j=1:k
        r(i,j)=bits_add(i,:)*transpose(bits_add(j,:))/n;  %归一化
    end
end
for i=1:k
    S_out(i,:)=decdma(bits_in,bits_add(i,:));
end
for i=1:k
    bits_mid=zeros(1,p);   %picmud的中间信号
    for x=1:k
        bits_mid(1,:)=bits_mid(1,:)+r(i,x)*S_out(x,:);
    end
    bits_out(i,:)=2*S_out(i,:)-bits_mid(1,:);
end
%以上结果有误差，需进一步的判定