Несколько советов рабочего процесса симуляций
Цель этой демонстрации состоит в том, чтобы обеспечить полезные советы для выполнения нескольких симуляций с помощью параллельных инструментов симуляции. Этот пример будет работать, даже если Parallel Computing Toolbox™ не будет доступен, но симуляции запустятся в сериале. Мы будем использовать модель sldemo_suspn_3dof в этом примере.
mdl = 'sldemo_suspn_3dof';
isModelOpen = bdIsLoaded(mdl);
open_system(mdl);
Инициализация массива Simulink. Объекты SimulationInput
Обычно вы создадите массив Simulink. SimulationInput возражает для того, чтобы запустить несколько симуляций. Существует несколько способов инициализировать массив прежде, чем заполнить его с данными.
numSims = 5; Cf_sweep = Cf*linspace(.05,.95, numSims);
Метод 1: Инициализируйте массив перед циклом
in(numSims) = Simulink.SimulationInput; for idx = 1:numSims % Need to populate the model name since we get any empty array by default in(idx).ModelName = 'sldemo_suspn_3dof'; in(idx) = in(idx).setVariable('Cf', Cf_sweep(idx)); end
Метод 2: Инициализируйте массив в цикле
Обратите внимание на то, что переменная цикла idx начинает с самого большого значения так, чтобы целый массив был предварительно выделен.
for idx = numSims:-1:1 % Since we are indexing from 5 to 1, the first iteration will % initialize the array. in(idx) = Simulink.SimulationInput('sldemo_suspn_3dof'); in(idx) = in(idx).setVariable('Cf', Cf_sweep(idx)); end
Установка модели и параметров блоков
setModelParameter и setBlockParameter методы используют тот же синтаксис пары значения параметров что set_param Использование API. Это означает, что большинство значений, в которых вы передаете этим методам, должно быть символьными массивами, не их литеральным значением.
for idx = numSims:-1:1 in(idx) = Simulink.SimulationInput('sldemo_suspn_3dof'); % Incorrect in(idx) = in(idx).setModelParameter('StartTime', 5); % Correct in(idx) = in(idx).setModelParameter('StartTime', '3'); end
Установка переменных
setVariable метод ожидает, что вы передадите литеральное значение, которое вы хотите присвоить переменной. Идея состоит в том, что это тесно отражает assignin синтаксис.
for idx = numSims:-1:1 in(idx) = Simulink.SimulationInput('sldemo_suspn_3dof'); % Incorrect, Cf is expected to be a double, not a character array in(idx) = in(idx).setVariable('Cf', '2500'); % Correct, Cf is a scalar double in(idx) = in(idx).setVariable('Cf', 2500); end
Диагностирование ошибок периода выполнения
Примите, что вы случайно настроили массив Simulink. SimulationInput возражает с неправильным значением.
Mb_sweep = linspace(0, 1500, numSims); for idx = numSims:-1:1 in(idx) = Simulink.SimulationInput('sldemo_suspn_3dof'); % Accidentally set the Mass to 0 on the first iteration in(idx) = in(idx).setVariable('Mb', Mb_sweep(idx)); % Shorten the stop time in(idx) = in(idx).setModelParameter('StopTime','1'); end
Симуляция их вызовет ошибку периода выполнения
out = sim(in);
[26-Jul-2019 18:51:01] Running simulations... [26-Jul-2019 18:51:02] Completed 1 of 5 simulation runs. Run 1 has errors. [26-Jul-2019 18:51:03] Completed 2 of 5 simulation runs [26-Jul-2019 18:51:03] Completed 3 of 5 simulation runs [26-Jul-2019 18:51:04] Completed 4 of 5 simulation runs [26-Jul-2019 18:51:04] Completed 5 of 5 simulation runs Warning: Simulation(s) with indices listed below completed with errors. Please inspect the corresponding SimulationOutput to get more details about the error: [1]
К счастью, можно смотреть Simulink. Объект SimulationOutput видеть любые сообщения об ошибке, которые прибывают из симуляции.
out(1).ErrorMessage
ans =
'Derivative of state '1' in block 'sldemo_suspn_3dof/Body Dynamics/Vertical (Z) dynamics/Zdot' at time 0.0 is not finite. The simulation will be stopped. There may be a singularity in the solution. If not, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances)'
Это работает, чтобы отладить проблемы о параллельных рабочих также.
for idx = numSims:-1:1 in(idx) = Simulink.SimulationInput('sldemo_suspn_3dof'); % Accidentally set the Mass to 0 on the first iteration in(idx) = in(idx).setVariable('Mb', Mb_sweep(idx)); % Shorten the stop time in(idx) = in(idx).setModelParameter('StopTime','1'); end out = parsim(in);
[26-Jul-2019 18:51:05] Checking for availability of parallel pool... Starting parallel pool (parpool) using the 'local' profile ... Connected to the parallel pool (number of workers: 12). [26-Jul-2019 18:51:50] Starting Simulink on parallel workers... [26-Jul-2019 18:52:21] Configuring simulation cache folder on parallel workers... [26-Jul-2019 18:52:22] Loading model on parallel workers... [26-Jul-2019 18:52:43] Running simulations... [26-Jul-2019 18:52:49] Completed 1 of 5 simulation runs. Run 1 has errors. [26-Jul-2019 18:52:49] Completed 2 of 5 simulation runs [26-Jul-2019 18:52:49] Completed 3 of 5 simulation runs [26-Jul-2019 18:52:49] Completed 4 of 5 simulation runs [26-Jul-2019 18:52:49] Completed 5 of 5 simulation runs Warning: Simulation(s) with indices listed below completed with errors. Please inspect the corresponding SimulationOutput to get more details about the error: [1] [26-Jul-2019 18:52:49] Cleaning up parallel workers...
Осмотр Simulink. SimulationOutput показывает неличную производную ошибку.
out(1).ErrorMessage
ans =
'Derivative of state '1' in block 'sldemo_suspn_3dof/Body Dynamics/Vertical (Z) dynamics/Zdot' at time 0.0 is not finite. The simulation will be stopped. There may be a singularity in the solution. If not, try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances)'
applyToModel метод сконфигурирует вашу модель с настройками на SimulationInput, таким образом, можно будет отладить проблему локально.
in(1).applyToModel;
Заметьте что значение переменной Mb в изменениях базового рабочего пространства в 0, чтобы отразить значение, которое использовалось в симуляции, соответствующей первому объекту SimulationInput в in.
Закройте работников MATLAB
Наконец, закройте параллельный пул и модель, если они не были ранее открыты.
if(~isModelOpen) close_system(mdl, 0); end delete(gcp('nocreate'));
Parallel pool using the 'local' profile is shutting down.