Оцените ответ разомкнутого цикла объекта в модели watertank. Откройте модель.

model = 'watertank';
open_system(model);

Чтобы оценить ответ разомкнутого цикла объекта, задайте набор ввода-вывода линеаризации, который задает этот фрагмент модели с аналитическими точками. Задайте входную аналитическую точку в контроллере вывод и задайте выходную точку разомкнутого цикла на объекте вывод.

io(1)=linio('watertank/PID Controller',1,'input');
io(2)=linio('watertank/Water-Tank System',1,'openoutput');

Найдите установившуюся рабочую точку для оценки. В данном примере используйте установившуюся рабочую точку, выведенную от образцовых начальных условий.

watertank_spec = operspec(model);
opOpts = findopOptions('DisplayReport','off');
op = findop(model,watertank_spec,opOpts);

Создайте входной сигнал для оценки. В данном примере используйте сигнал sinestream, который отправляет серию отдельных синусоидальных возмущений на частотах, которые вы задаете.

input = frest.Sinestream('Frequency',logspace(-3,2,30));

Оцените частотную характеристику заданного фрагмента модели. Результатом является модель частотной характеристики, содержащая ответы на каждой из частот, заданных в сигнале sinestream.

sysest = frestimate(model,op,io,input);
size(sysest)

FRD model with 1 outputs, 1 inputs, and 30 frequency points.

Исследуйте измеренную частотную характеристику.

bode(sysest,'*')

Откройте модель watertank.

model = 'watertank';
open_system(model);

Получите линеаризацию ответа разомкнутого цикла объекта. Для этого задайте точки ввода-вывода линеаризации и найдите установившуюся рабочую точку около образцовых начальных условий. Затем линеаризуйте модель.

io(1)=linio('watertank/PID Controller',1,'input');
io(2)=linio('watertank/Water-Tank System',1,'openoutput');

watertank_spec = operspec(model);
opOpts = findopOptions('DisplayReport','off');
op = findop(model,watertank_spec,opOpts);

syslin = linearize(model,op,io);

Чтобы проверять линеаризацию, используйте те же аналитические точки и рабочую точку, чтобы оценить частотную характеристику. В данном примере используйте sinestream входной сигнал для оценки.

input = frest.Sinestream('Frequency',logspace(-3,2,20));
sysest = frestimate(model,op,io,input);

Сравните точную линеаризацию и предполагаемый ответ в частотном диапазоне с помощью Диаграммы Боде.

bode(syslin,'b-',sysest,'r*')
legend('Exact linearization','Estimation')

Оцените ответ разомкнутого цикла объекта в модели watertank. Во-первых, откройте модель.

model = 'watertank';
open_system(model);

Задайте набор ввода-вывода линеаризации, который задает объект, и найдите установившуюся рабочую точку для оценки.

io(1)=linio('watertank/PID Controller',1,'input');
io(2)=linio('watertank/Water-Tank System',1,'openoutput');

watertank_spec = operspec(model);
opOpts = findopOptions('DisplayReport','off');
op = findop(model,watertank_spec,opOpts);

Затем создайте входной сигнал для оценки и оцените частотную характеристику заданного фрагмента модели. Используйте выходной аргумент simout, чтобы хранить данные об оценке.

input = frest.Sinestream('Frequency',logspace(-3,2,10));
[sysest,simout] = frestimate(model,op,io,input);

Откройте средство просмотра результатов симуляции.

frest.simView(simout,input,sysest)

Средство просмотра показывает вам установившийся ответ времени и БПФ того ответа для всех частот в области значений, которую вы выбираете на разделе Bode Diagram средства просмотра. Эти графики могут помочь вам идентифицировать, когда ответ отклоняется от ожидаемого ответа. Для получения дополнительной информации об использовании Средства просмотра Результатов симуляции смотрите, Анализируют Предполагаемую Частотную характеристику.

Если у вас есть линейная модель системы, вы оцениваете, можно использовать модель в качестве базового ответа для сравнения в средстве просмотра. Например, можно сравнить модель, полученную точной линеаризацией к предполагаемой частотной характеристике. Используйте набор ввода-вывода линеаризации и рабочую точку, чтобы вычислить точную линеаризацию объекта watertank.

syslin = linearize(model,io,op);

Откройте Средство просмотра Результатов симуляции снова, на этот раз обеспечив syslin как входной параметр.

frest.simView(simout,input,sysest,syslin)

Раздел Bode Diagram средства просмотра включает строку, показывающую точный ответ syslin. Это представление может быть полезным, чтобы идентифицировать особые частоты, где предполагаемый ответ отклоняется от линеаризации.