Global Index (short | long) | Local contents | Local Index (short | long)
[w,amp,phas,sdamp,sdphas] = getff(g,nu,ny)
GETFF Selects the frequency function (for own plotting). [W,AMP,PHAS] = GETFF(G,NU,NY) W: The frequency scale in radians/second AMP: The amplitude function PHAS: The phase function (in degrees) G: The frequency function in the FREQFUNC format. See also FREQFUNC. NU: The input number (noise input is counted as input # 0) (Default 1. If G contains all spectra, then default is 0) NY: The output number (default 1) If several entries in G correspond to the same input-output relation W, AMP and PHAS will have the corresponding number of columns. The standard deviation of the amplitude and the standard deviation of the phase are obtained by [W,AMP,PHAS,SD_AMP,SD_PHAS] = GETFF(G,NU,NY)
This function calls | This function is called by |
---|---|
function [w,amp,phas,sdamp,sdphas] = getff(g,nu,ny) % L. Ljung 10-2-90 % Copyright (c) 1986-98 by The MathWorks, Inc. % $Revision: 2.3 $ $Date: 1997/12/02 03:40:53 $ if nargin < 1 disp('Usage: [FREQs,AMPLITUDE,PHASE] = GETFF(G)') disp(' [FREQs,AMP,PHASE,SDamp,SDphase] = GETFF(G,INPUT_no,OUTPUT_no)') return end if nargin<3, ny=1;end if isempty(ny), ny=1;end if nargin<2, nu=1;end, if isempty(nu), nu=1;end [rg,cg]=size(g); info=g(1,:); if all(rem(info,10)==0), %then all G contain spectral information nu=0;end w=g(2:rg,find(info==(ny-1)*1000+100+nu)); amp=g(2:rg,find(info==(ny-1)*1000+nu)); sdamp=g(2:rg,find(info==(ny-1)*1000+50+nu)); phas=g(2:rg,find(info==(ny-1)*1000+20+nu)); sdphas=g(2:rg,find(info==(ny-1)*1000+70+nu));