Исследуйте файлы

Чтобы получить доступ к файлам, которые использует этот пример, нажмите Open Script.

kalmanfilter.c

kalmanfilter.c функция C, которую тестирует пример. Это оценивает положение перемещения, основанного на объектах на его прошлых положениях.

kalmanfilter.h

kalmanfilter.h заголовочный файл для kalmanfilter.c.

position.mat

position.mat содержит положения объекта.

callKalmanFilter.m

callKalmanFilter вызовы kalmanfilter при помощи coder.ceval.

function [a,b] = callKalmanFilter(position)
% Copyright 2014 - 2016 The MathWorks, Inc.

numPts = size(position,2);

a = zeros(2,numPts,'double');
b = zeros(2,numPts,'double');
y = zeros(2,1,'double');

% Main loop
for idx = 1: numPts
    z = position(:,idx);     % Get the input data
    
    % Call the initialize function
    coder.ceval('kalmanfilter_initialize');
    
    % Call the C function
    coder.ceval('kalmanfilter',z,coder.ref(y));
    
    % Call the terminate function
    coder.ceval('kalmanfilter_terminate');
    
    a(:,idx) = [z(1); z(2)];
    b(:,idx) = [y(1); y(2)];
end
end

TestKalmanFilter.m

TestKalmanFilter тесты, превышает ли ошибка между предсказанным положением и фактическим положением заданный допуск. Модульные тесты являются модульными тестами, основанными на классах. Для получения дополнительной информации смотрите Создание модульных тестов на основе классов в MATLAB.

Несмотря на то, что вы хотите протестировать MEX-функцию, модульные тесты в TestKalmanFilter вызовите исходную функцию MATLAB, от которой вы сгенерировали MEX-функцию. Когда MATLAB Coder запускает тесты, он заменяет вызовы функции MATLAB с вызовами MEX-функции. Вы не можете запустить эти тесты непосредственно в MATLAB, потому что MATLAB не распознает coder.ceval вызовы в callKalmanFilter.

classdef TestKalmanFilter < matlab.unittest.TestCase
    % Copyright 2014 - 2016 The MathWorks, Inc.
    
    methods ( Test )
        
        function SSE_LessThanTolerance( testCase )
            load position.mat;
            [z,y] = callKalmanFilter( position );
            
            tolerance = 0.001; % tolerance of 0.0001 will break
            A = z-1000*y;
            error = sum(sum(A.^2));
            
            testCase.verifyLessThanOrEqual( error, tolerance);
            
            % For debugging
            plot_kalman_filter_trajectory(z,1000*y);
        end
        
        function SampleErrorLessThanTolerance( testCase )
            load position.mat;
            [z,y] = callKalmanFilter( position );
            
            tolerance = 0.01;   % tolerance of 0.001 will break
            A = z-1000*y;

            testCase.verifyEqual(1000*y, z, 'AbsTol', tolerance);
            % For debugging
            plot_kalman_filter_trajectory(z,1000*y);
            
            [value, location] = max(A(:));
            [R,C] = ind2sub(size(A),location);
            disp(['Max value ' num2str(value) ' is located at [' num2str(R) ',' num2str(C) ']']);
        end
    end
end

run_unit_tests_kalman.m

run_unit_tests_kalman вызовы runtests запускать тесты в TestKalmanFilter.m.

% Run unit tests
% Copyright 2014 - 2016 The MathWorks, Inc.

runtests('TestKalmanFilter')

plot_kalman_filter_trajectory.m

plot_kalman_filter_trajectory строит траекторию предполагаемых и фактических положений объекта. Каждый модульный тест вызывает эту функцию.

Сгенерируйте MEX и запущенные модульные тесты в приложении MATLAB Coder

Чтобы открыть приложение MATLAB Coder, на вкладке MATLAB Toolstrip Apps, под Генерацией кода, кликают по значку приложения MATLAB Coder.

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

Модульные тесты загружают переменную position от position.mat и передайте position к callKalmanFilter. Поэтому вход к callKalmanFilter должен иметь свойства что position имеет. В рабочем пространстве MATLAB, если вы загружаете position.mat, вы видите тот position 2 310 массив типа double.

Сконфигурируйте приложение для генерации кода MEX. Задайте имена исходных и заголовочных файлов C потому что callKalmanFilter интегрирует внешний код С.

Чтобы сгенерировать MEX-функцию, нажмите Generate.

Запустите модульные тесты на сгенерированном MEX.

Когда выполнение приложения тестовый файл, это заменяет вызовы callKalmanFilter в модульном тесте с вызовами callKalmanFilter_mex. Модульные тесты работают на MEX-функции вместо исходной функции MATLAB.

Отображения приложения тест выводятся на вкладке Test Output. Передача модульных тестов.

Из графиков вы видите, что траектория предполагаемого положения сходится с траекторией фактического положения.

Запустите модульные тесты после изменения кода С

Когда вы изменяете код С, чтобы запустить модульные тесты:

Например, измените kalmanfilter.c так, чтобы значение, присвоенное y[r2] умножается на 1,1.

y[r2] += (double)d_a[r2 + (i0 << 1)] * x_est[i0] * 1.1;

Отредактируйте kalmanfilter.c за пределами приложения, потому что можно использовать приложение, чтобы отредактировать только файлы MATLAB, перечисленные в панели Исходного кода приложения.

Чтобы сгенерировать MEX-функцию для модифицированной функции, нажмите Generate.

Запускать модульные тесты:

Тесты перестали работать, потому что ошибка превышает заданный допуск.

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

Сгенерируйте MEX и запущенные модульные тесты при помощи рабочего процесса командной строки

Можно использовать рабочий процесс командной строки, чтобы запустить модульные тесты на внешнем коде С при помощи coder.runTest. Задайте тестовый файл, который запускает модульные тесты на функции MATLAB, которая вызывает ваш код С.

Сгенерируйте MEX-функцию для функции MATLAB, которая вызывает ваш код С. В данном примере сгенерируйте MEX для callKalmanFilter.

Создайте объект настройки для генерации кода MEX.

cfg = coder.config('mex');

Задайте внешний исходный код и заголовочный файл.

cfg.CustomSource = 'kalmanfilter.c';
cfg.CustomHeaderCode = '#include "kalmanfilter.h"';

Определить тип для входа к callKalmanFilter, загрузите файл положения.

load position.mat

Чтобы сгенерировать MEX-функцию, запустите codegen. Укажите что вход к callKalmanFilter имеет тот же тип как position.

codegen -config cfg callKalmanFilter -args position

Code generation successful.

Запустите модульные тесты на MEX-функции. Укажите, что тестовым файлом является run_unit_tests_kalman и что функцией является callKalmanfilter. Когда coder.runTest запускает тестовый файл, он заменяет вызовы callKalmanFilter в модульном тесте с вызовами callKalmanFilter_mex. Модульные тесты работают на MEX-функции вместо исходной функции MATLAB.

coder.runTest('run_unit_tests_kalman', 'callKalmanFilter')

Running TestKalmanFilter
Current plot held
.Current plot held
Max value 0.0010113 is located at [2,273]
.
Done TestKalmanFilter
__________


ans = 

  1x2 TestResult array with properties:

    Name
    Passed
    Failed
    Incomplete
    Duration
    Details

Totals:
   2 Passed, 0 Failed, 0 Incomplete.
   20.1959 seconds testing time.