Создайте нелинейную серую модель. Модель представляет собой линейный двигатель постоянного тока с одним входом (напряжением) и двумя выходами (угловое положение и угловая скорость). Структура модели определяется файлом dcmotor_m.m.

FileName = 'dcmotor_m'; 
Order = [2 1 2];
Parameters = [0.24365;0.24964];  
nlgr = idnlgrey(FileName,Order,Parameters);  
nlgr = setinit(nlgr, 'Fixed', false(2,1)); % set initial states free

Данные загрузки для оценки исходного состояния.

load(fullfile(matlabroot,'toolbox','ident',...
     'iddemos','data','dcmotordata'));
z = iddata(y,u,0.1);

Оцените начальные состояния так, чтобы отклик модели с использованием оцененных состояний X0 и измеренного входного сигнала u был как можно ближе к измеренному выходному сигналу y.

X0 = findstates(nlgr,z,Inf);

Оценка idss моделируют его таким образом, чтобы отклик оцененной модели как можно ближе соответствовал выходному сигналу оценочных данных.

Загрузить данные образца.

load iddata1 z1;

Оцените линейную модель на основе данных.

model = ssest(z1,2);

Оценка значения начальных состояний для наилучшего соответствия измеренному выходному сигналу z1.y.

x0est = findstates(model,z1,Inf);

Моделирование модели.

opt = simOptions('InitialCondition',x0est);
sim(model,z1(:,[],:),opt);

Оценивать исходные состояния модели выборочно, фиксируя первое состояние и позволяя оценивать второе состояние модели.

Создайте нелинейную серую модель.

FileName = 'dcmotor_m'; 
Order = [2 1 2];
Parameters = [0.24365;0.24964];  
nlgr = idnlgrey(FileName,Order,Parameters);

Модель представляет собой линейный двигатель постоянного тока с одним входом (напряжением) и двумя выходами (угловое положение и угловая скорость). Структура модели определяется dcmotor_m.m файл.

Загрузите оценочные данные.

load(fullfile(matlabroot,'toolbox','ident',...
     'iddemos','data','dcmotordata'));
z = iddata(y,u,0.1);

Удерживайте первое состояние фиксированным на нуле и оцените значение второго.

x0spec = idpar('x0',[0;0]);
x0spec.Free(1) = false;
opt = findstatesOptions;
opt.InitialState = x0spec;
[X0,Report] = findstates(nlgr,z,Inf,opt)

X0 = 2×1

         0
    0.0061

Report = 
         Status: 'Estimated by simulation error minimization'
         Method: 'lsqnonlin'
     Covariance: [2x2 double]
       DataUsed: [1x1 struct]
    Termination: [1x1 struct]

Создайте нелинейную серую модель.

FileName = 'dcmotor_m'; 
Order = [2 1 2];
Parameters = [0.24365;0.24964];  
nlgr = idnlgrey(FileName,Order,Parameters);

Модель представляет собой линейный двигатель постоянного тока с одним входом (напряжением) и двумя выходами (угловое положение и угловая скорость). Структура модели определяется dcmotor_m.m файл.

Загрузите оценочные данные.

load(fullfile(matlabroot,'toolbox','ident',...
     'iddemos','data','dcmotordata'));
z = iddata(y,u,0.1);

Укажите начальное предположение для начальных состояний.

x0spec = idpar('x0',[10;10]);

x0spec.Free по умолчанию имеет значение true

Оценка начальных состояний

opt = findstatesOptions;
opt.InitialState = x0spec;
x0 = findstates(nlgr,z,Inf,opt)

x0 = 2×1

    0.0362
   -0.1322

Создайте нелинейную серую модель.

FileName = 'dcmotor_m'; 
Order = [2 1 2];
Parameters = [0.24365;0.24964];  
nlgr = idnlgrey(FileName,Order,Parameters);  
set(nlgr, 'InputName','Voltage','OutputName', ...
		{'Angular position','Angular velocity'});

Модель представляет собой линейный двигатель постоянного тока с одним входом (напряжением) и двумя выходами (угловое положение и угловая скорость). Структура модели определяется dcmotor_m.m файл.

Загрузите оценочные данные.

load(fullfile(matlabroot,'toolbox','ident',...
     'iddemos','data','dcmotordata'));
z = iddata(y,u,0.1,'Name','DC-motor',...
     'InputName','Voltage','OutputName',...
     {'Angular position','Angular velocity'});

Создайте набор данных из трех экспериментов.

z3 = merge(z,z,z);

Выберите эксперимент для оценки начальных состояний:

Фиксированные начальные состояния имеют нулевые значения.

x0spec = idpar('x0',zeros(2,3));
x0spec.Free(1,2) = false;
x0spec.Free(2,[2 3]) = false;
opt = findstatesOptions;
opt.InitialState = x0spec;

Оценка начальных состояний

[X0,EstInfo] = findstates(nlgr,z3,Inf,opt);