imcinit.m

类神经网路─MATLAB的应用(范例程式)

% ------------------------------>  IMCINIT.M  <------------------------------

% ----------      Switches       -----------
regty      ='imc';           % Controller type (imc, pid, none)
refty      ='siggener';      % Reference signal (siggener/<var. name>)
simul      ='simulink';      % Control object spec. (simulink/matlab/nnet)
if exist('simulink')~=4,
  simul      ='matlab';      % Simulink not present
end


% ------    General Initializations  -------
Ts = 0.20;                   % Sampling period (in seconds)
samples = 200 ;              % Number of samples to be simulated


% --  System to be Controlled (SIMULINK) --
integrator= 'ode45';         % Name of dif. eq. solver (f. ex. ode45 or ode15s)
sim_model = 'spm1';          % Name of SIMULINK model


% ---  System to be Controlled (MATLAB)  --
mat_model = 'springm';       % Name of MATLAB model
model_out = 'smout';         % Output equation (function of the states)
x0        = [0;0];           % Initial states
 

% ----- Inverse Network Specification ------
% The file must contain the variables:
% NN, NetDefi, W1i, W2i
% (i.e. regressor structure, architecture definition, and weight matrices)
nninv = 'inverse3';          % File name


% ----- Forward Network Specification ------
% Only used if simul='nnet' (no Simulink or Matlab model available)
% The file must contain: NN, NetDeff, W1f, W2f
nnforw = 'forward';          % File name


% ---------------- Filter ------------------
Am = [1 -1.4 0.49];          % Filter denominator
Bm = [0.06 0.03];            % Filter numerator
Am=1;Bm=1;


% ----------- Reference signal -------------
dc      = 0;                 % DC-level
sq_amp  = 1;                 % Amplitude of square signals (row vector)
sq_freq = 0.1;               % Frequency of square signals (column vector)
sin_amp = [0];               % Amplitude of sine signals  (row vector)
sin_freq= [0]';              % Frequency of sine signals   (column vector)
Nvar  = 0';                  % Variance of white noise signal


% -----  Linear Controller Parameters  ----- 
K=8;                         % PID parameters
Td=0.8;                      % PID
alf=0.1;                     % PID
Wi=0.2;                      % PID (1/Ti)


% ------ Specify data vectors to plot ------
% Notice that all strings in plot_a and plot_b, resp. MUST have the same length
plot_a(1,:) = 'ref_data    ';
plot_a(2,:) = 'y_data      ';
%plot_a(3,:) = 'ym_data     ';
plot_b(1,:) = 'u_data';