ambiguity.m

来自「这是一个关于模糊函数的matlab源代码」· M 代码 · 共 34 行

function Z = ambiguity(Ut, normalized)
  % FOR GOOD RESULTS  the length of the input vector (Ut)should be more
  % than 7 
  % Example - a rectangular pulse shuld be [1 1 1 1 1 1 1]
  % 


  if length(Ut)<7
      warning ('for good visual result the length of the input vector (Ut)should be more than 7 ');
  end
      
  if nargin == 1
	normalize = 'false';
  end

  A=xcorril(Ut);
  B=fft(A);
  Z=abs(B);
  clear A B;
  halfLength=floor(length(Z)/2);
  A1=Z(1:halfLength,:);
  A2=Z(halfLength+1:length(Z),:);
  X=cat(1,A2,A1);
  Z=X;
  clear A1 A2 X;
  
  k = strcmp(normalize,'true');
  if k == 1 
      m = max (Z);
      mm = max (m);
      Z=Z./mm;
  end