pptmfq1.m

阵列信号处理的工具箱

function pptmfq1(signals, noTgt, snapshot, smplPoints)

%PPTMFQ1 Plots the error function Q for the PTMF method.
%
%--------
%Synopsis:
%  function pptmfq(signals, noTgt, snapshot, smplPoints)
%
%Output and Input:
%  signals  (RxRadarSigT) : The 'true' signals.
%  noTgt  (RealScalarT) : No. of targets. Have to be supplied.
%  snapshot (RealScalarT) : Plot for this snapshot.
%  smplPoints (Vector of DoaT)  : Plot for these doas. Not implemented yet.
%
%Description
%  Plots the error function Q for the standard antenna signal model:
%   x(t) = A(Theta)*b(t) + n(t)
%  The prediction is xhat(t) = A(Thetahat)*bhat, bhat(t)=pinv(A(Thetahat))*x(t)
%  The error function Q = (x(t)-xhat(t))'*(x(t)-xhat(t))
%  x(t) is the signal supplied in the 'signals' input parameter. 'snapshot'
%  chooses the time t. Theta = [theta1, theta2, ...] angles to different
%  targets.
%
%Known Bugs:
%  This function (DBT release 1.11) seems to give the negative values
%  of the angles.
%
%See Also:
%  doapar1, perrfun, dbtex4



% *****************************************************************************
%   *  DBT, A Matlab Toolbox for Radar Signal Processing  *
% (c) FOA 1994-2000. See the file dbtright.m for copyright notice.
%
%
%  Start        : 970317 Svante Bj鰎klund (svabj).
%  Latest change: $Date: 2000/10/16 15:21:23 $ $Author: svabj $.
%  $Revision: 1.13 $
% *****************************************************************************

% ----------------------------------------------------------------------- %
% Pick out fields from input parameters.
% ----------------------------------------------------------------------- %

chkdtype(signals, 'RxCompSigT','RxRadarSigT')
X   = getm(signals.signals, ':',1,':',1,1)';
T   = size(X,2);    % Number of snapshots == number of iterations.

ant   = signals.antenna;

chkdtype(ant, 'AntDefT')
K   = ant.noElem;
D   = ant.distElem;
%lambda  = ant.lambda;
waveform = signals.waveform;
lambda  = waveform.wavelength;


if (nargin >= 4)
  error('DBT-Error: smplPoints is not implemented yet.')
end
if (nargin < 4)
  smplPoints = [];
end
if isempty(smplPoints)
  smplPoints = d2r(-90:0.5:90);
end


if (noTgt == 1)
  maxth2 = 1; % Maximum number of angle steps for the second angle.
elseif (noTgt == 2)
  maxth2 = 30;  % Maximum number of angle steps for the second angle.
else
  error('DBT-Error: Only possible to plot Q for one or two angles.')
end%if


%-------------------------------------------------------------------------- %
% Calculate the Q-function.
% ------------------------------------------------------------------------- %

% Qniva2.m Ritar graf p