turbounitray.asv

根据TD-SCDMA系统的特点

s19 = x19.*(scramble'*ones(1,44));

source10 =uspreadcode18*WQ(10);
x10 = ovsf(10)' * source10;                 % spreading
s10 = x10.*(scramble'*ones(1,44));

% the second antenna spread 
spreadcode2=reshape(unitraymatrix1,1,440);% symbols frame 
uspreadcode21=spreadcode2(1:44);
uspreadcode22=spreadcode2(45:88);
uspreadcode23=spreadcode2(89:132);
uspreadcode24=spreadcode2(133:176);
uspreadcode25=spreadcode2(177:220);
uspreadcode26=spreadcode2(221:264);
uspreadcode27=spreadcode2(265:308);
uspreadcode28=spreadcode2(309:352);
uspreadcode29=spreadcode2(353:396);
uspreadcode20=spreadcode2(397:440);

source21 =uspreadcode21*WQ(1);
x21 = ovsf(1)' * source21;                 % spreading
s21 = x21.*(scramble'*ones(1,44));               


source22 = uspreadcode22*WQ(2);
x22 = ovsf(2)' * source22;                 % spreading
s22 = x22.*(scramble'*ones(1,44));        

source23 =uspreadcode23*WQ(3);
x23 = ovsf(3)' * source23;                 % spreading
s23 = x23.*(scramble'*ones(1,44));

source24 =uspreadcode24*WQ(4);
x24 = ovsf(4)' * source24;                 % spreading
s24 = x24.*(scramble'*ones(1,44));

source25 =uspreadcode25*WQ(5);
x25 = ovsf(5)' * source25;                 % spreading
s25 = x25.*(scramble'*ones(1,44));

source26 =uspreadcode26*WQ(6);
x26 = ovsf(6)' * source26;                 % spreading
s26 = x26.*(scramble'*ones(1,44));

source27 =uspreadcode27*WQ(7);
x27 = ovsf(7)' * source27;                 % spreading
s27 = x27.*(scramble'*ones(1,44));

source28 =uspreadcode28*WQ(8);
x28 = ovsf(8)' * source28;                 % spreading
s28 = x28.*(scramble'*ones(1,44));

source29 =uspreadcode29*WQ(9);
x29 = ovsf(9)' * source29;                 % spreading
s29 = x29.*(scramble'*ones(1,44));

source20 =uspreadcode20*WQ(10);
x20 = ovsf(10)' * source20;                 % spreading
s20 = x20.*(scramble'*ones(1,44));





% the first antenna despread

s1=h1*(s11+s12+s13+s14+s15+s16+s17+s18+s19+s10)+noise;



r1 = reshape(s1,16,44);


z11 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(1))*r1;
rs11 = reshape(diag(z11)/16,1,44)./WQ(1);

z12 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(2))*r1;
rs12 = reshape(diag(z12)/16,1,44)./WQ(2);
 
z13 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(3))*r1;
rs13 = reshape(diag(z13)/16,1,44)./WQ(3);

z14 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(4))*r1;
rs14 = reshape(diag(z14)/16,1,44)./WQ(4);

z15 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(5))*r1;
rs15 = reshape(diag(z15)/16,1,44)./WQ(5);

z16 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(6))*r1;
rs16 = reshape(diag(z16)/16,1,44)./WQ(6);

z17 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(7))*r1;
rs17 = reshape(diag(z17)/16,1,44)./WQ(7);

z18 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(8))*r1;
rs18 = reshape(diag(z18)/16,1,44)./WQ(8);

z19 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(9))*r1;
rs19 = reshape(diag(z19)/16,1,44)./WQ(9);

z10 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(10))*r1;
rs10 = reshape(diag(z10)/16,1,44)./WQ(10);

rdspreadcode=[rs11 rs12 rs13 rs14 rs15 rs16 rs17 rs18 rs19 rs10];
urdspreadcode=reshape(rdspreadcode,8,55);
% the second antenna despread

s2=h2*(s21+s22+s23+s24+s25+s26+s27+s28+s29+s20)+noise;


r2 = reshape(s2,16,44);


z21 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(1))*r2;
rs21 = reshape(diag(z21)/16,1,44)./WQ(1);

z22 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(2))*r2;
rs22 = reshape(diag(z22)/16,1,44)./WQ(2);
 
z23 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(3))*r2;
rs23 = reshape(diag(z23)/16,1,44)./WQ(3);

z24 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(4))*r2;
rs24 = reshape(diag(z24)/16,1,44)./WQ(4);

z25 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(5))*r2;
rs25 = reshape(diag(z25)/16,1,44)./WQ(5);

z26 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(6))*r2;
rs26 = reshape(diag(z26)/16,1,44)./WQ(6);

z27 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(7))*r2;
rs27 = reshape(diag(z27)/16,1,44)./WQ(7);

z28 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(8))*r2;
rs28 = reshape(diag(z28)/16,1,44)./WQ(8);

z29 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(9))*r2;
rs29 = reshape(diag(z29)/16,1,44)./WQ(9);

z20 = ones(44, 1)*scramble.*(ones(44, 1)*ovsf(10))*r2;
rs20 = reshape(diag(z20)/16,1,44)./WQ(10);

rdspreadcode2=[rs21 rs22 rs23 rs24 rs25 rs26 rs27 rs28 rs29 rs20];
urdspreadcode2=reshape(rdspreadcode2,8,55);


for nn=1:55
    dmatrix=zeros(8,2);
    dmatrix(:,1)=urdspreadcode(:,nn);
    dmatrix(:,2)=urdspreadcode2(:,nn);
    n=0;
    for nnn=0:255
     TR=trace( dmatrix'*(smatrix^nnn*matrix0)*(smatrix^nnn*matrix0)'*dmatrix);
    % max angle
        if TR>n
         n=TR;
         length=nnn;
     end
 end
 dst=de2bi(length,8,'left-msb');
 dsource(:,nn)=dsource(:,nn)+dst';
end
rrsource1=reshape(dsource,1,440);
containmatrix(:,kk)=rrsource1.';
end  % one frame end
rsource2=reshape(containmatrix,1,1760);
linkmatrix(:,TT)=rsource2.';
end
rsourcelink=linkmatrix(1:7040);
rsource=rsourcelink(1:7038);
      
      %  unitray space-time decoded
      
      
      
      yk = demultiplex(rsource,alpha,puncture); % demultiplex to get input for decoder 1 and 2
      
% Scale the received bits      
      rec_s = 0.5*L_c*yk;

% Initialize extrinsic information      
      L_e(1:L_total) = zeros(1,L_total);
      
      for iter = 1:niter
% Decoder one
         L_a(alpha) = L_e;  % a priori info. 
         if dec_alg == 0
            L_all = logmapo(rec_s(1,:), g, L_a, 1);  % complete info.
         else   
            L_all = sova0(rec_s(1,:), g, L_a, 1);  % complete info.
         end   
         L_e = L_all - 2*rec_s(1,1:2:2*L_total) - L_a;  % extrinsic info.

% Decoder two         
         L_a = L_e(alpha);  % a priori info.
         if dec_alg == 0
            L_all = logmapo(rec_s(2,:), g, L_a, 2);  % complete info.  
         else
            L_all = sova0(rec_s(2,:), g, L_a, 2);  % complete info. 
         end
         L_e = L_all - 2*rec_s(2,1:2:2*L_total) - L_a;  % extrinsic info.
         
% Estimate the info. bits        
         xhat(alpha) = (sign(L_all)+1)/2;

% Number of bit errors in current iteration
         err(iter) = length(find(xhat(1:L_total-m)~=x));
% Count frame errors for the current iteration
         if err(iter)>0
            nferr(nEN,iter) = nferr(nEN,iter)+1;
         end   
      end	%iter
      
% Total number of bit errors for all iterations
      errs(nEN,1:niter) = errs(nEN,1:niter) + err(1:niter);

      if rem(nframe,3)==0 | nferr(nEN, niter)==ferrlim
% Bit error rate
         ber(nEN,1:niter) = errs(nEN,1:niter)/nframe/(L_total-m);
% Frame error rate
         fer(nEN,1:niter) = nferr(nEN,1:niter)/nframe;

% Display intermediate results in process  
         fprintf('************** Eb/N0 = %5.2f db **************\n', EbN0db(nEN));
         fprintf('Frame size = %d, rate 1/%d. \n', L_total, 2+puncture);
         fprintf('%d frames transmitted, %d frames in error.\n', nframe, nferr(nEN, niter));
         fprintf('Bit Error Rate (from iteration 1 to iteration %d):\n', niter);
         for i=1:niter
            fprintf('%8.4e    ', ber(nEN,i));
         end
         fprintf('\n');
         fprintf('Frame Error Rate (from iteration 1 to iteration %d):\n', niter);
         for i=1:niter
            fprintf('%8.4e    ', fer(nEN,i));
         end
         fprintf('\n');
         fprintf('***********************************************\n\n');

% Save intermediate results 
         save turbo_sys_demo EbN0db ber fer
      end
      
   end		%while
end 		%nEN

diary off