ip_09_05.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:1k
源码类别:

通讯编程文档

开发平台:

Matlab

  1. % MATLAB script for Illustrative Problem 9.5.
  2. echo on
  3. % First determine the maximal length shift-register sequences.
  4. % Assume the initial shift-register content as "00001".
  5. connections1=[1 0 1 0 0];
  6. connections2=[1 1 1 0 1];
  7. sequence1=ss_mlsrs(connections1);
  8. sequence2=ss_mlsrs(connections2);
  9. % Cyclically shift the second sequence and add it to the first one.
  10. L=2^length(connections1)-1;;
  11. for shift_amount=0:L-1,
  12.   temp=[sequence2(shift_amount+1:L) sequence2(1:shift_amount)];
  13.   gold_seq(shift_amount+1,:)=(sequence1+temp) - floor((sequence1+temp)./2).*2;
  14.   echo off ;
  15. end;
  16. echo on ;
  17. % Find the max value of the cross-correlation for these sequences.
  18. max_cross_corr=0;
  19. for i=1:L-1,
  20.   for j=i+1:L,
  21.     % equivalent sequences
  22.     c1=2*gold_seq(i,:)-1;
  23.     c2=2*gold_seq(j,:)-1;
  24.     for m=0:L-1,
  25.       shifted_c2=[c2(m+1:L) c2(1:m)];
  26.       corr=abs(sum(c1.*shifted_c2));
  27.       if (corr>max_cross_corr),
  28.         max_cross_corr=corr;
  29.       end;
  30.       echo off ; 
  31.     end;
  32.   end;
  33. end;
  34. % Note that max_cross_correlation turns out to be 9 in this example.