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.')endif (nargin < 4)  smplPoints = [];endif isempty(smplPoints)  smplPoints = d2r(-90:0.5:90);endif (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