conv.m

来自「在高斯信道下」· M 代码 · 共 18 行

function conv_output=conv(G,k,conv_input) 
[h,f]=size(G);
n=size(G,1);
l=length(conv_input);
j=1;
in=[zeros(1,f),conv_input];
%zeros(1,h)
if(rem(l,k)~=0) 
    error('Error,g is not of the right size.');
    
end

for i=k+1:k:length(in)-f+1
    inn=fliplr(in(i:i+f-1));
    conv_output(j:j+n-1)=rem(inn*G',2);
    j=j+n;
end