df.m

为了下载东西

function f2 = df(F,GE,GCE,GU,A,W)
% function f2 = df(F,GE,GCE,GU,A,W)
%
% Describing function
%
% F       Fuzzy lookup table
% GE      Gain
% GCE     Gain
% GU      Gain
% A       Input amplitudes (vector)
% W       Angular frequencies (vector)
% f2      1. harmonic (complex vector)

% Jantzen 940909, not finished

if nargin == 4,
  A = logspace(-1,0,20) ;
  W = logspace(-1,1,20) ;
end ;

Step = 2*pi/50  ;
Wt = 0:Step:(2*pi) ;
f2 = zeros(length(A),length(W)) ;

for M=1:length(W),
  Sincos = [sin(Wt'), W(M)*cos(Wt')] ;
  for N=1:length(A),
    Gmat = [GE 0; 0 GCE] ;
    U = GU * ftab(A(N)*Sincos*Gmat,F) ;
    % Rectangular integration ;
    Y = U' * Sincos * Step/pi ;
    Y = Y(:,1) + Y(:,2) * i ;
    f2(N,M) = Y/A(N) ;
  end;
end ;