Exponenta Banner

Обнаружение переполнений

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

Необходимые условия

Чтобы завершить этот пример, необходимо установить следующие продукты:

  • MATLAB®

  • Файлы MATLAB Coder™

  • Fixed-Point Designer™

В локальной папке с возможностью записи создайте функцию overflow.

function y = overflow(b,x,reset)
    if nargin<3, reset = true; end
    persistent z p
    if isempty(z) || reset
        p = 0;
        z = zeros(size(b));
    end
    [y,z,p] = fir_filter(b,x,z,p);
end
function [y,z,p] = fir_filter(b,x,z,p)
    y = zeros(size(x));
    nx = length(x);
    nb = length(b);
    for n = 1:nx
        p=p+1; if p>nb, p=1; end
        z(p) = x(n);        
        acc = 0;
        k = p;
        for j=1:nb
            acc = acc + b(j)*z(k);
            k=k-1; if k<1, k=nb; end
        end        
        y(n) = acc;
    end
end

Создайте тестовый файл, overflow_test.m для выполнения overflow алгоритм. 

function overflow_test
    % The filter coefficients were computed using the FIR1 function from
    % Signal Processing Toolbox.
    %   b = fir1(11,0.25);
    b = [-0.004465461051254
         -0.004324228005260
         +0.012676739550326
         +0.074351188907780
         +0.172173206073645
         +0.249588554524763
         +0.249588554524763
         +0.172173206073645
         +0.074351188907780
         +0.012676739550326
         -0.004324228005260
         -0.004465461051254]';
    
    % Input signal
    nx = 256;
    t = linspace(0,10*pi,nx)';

    % Impulse
    x_impulse = zeros(nx,1); x_impulse(1) = 1;

    % Max Gain
    % The maximum gain of a filter will occur when the inputs line up with the
    % signs of the filter's impulse response.
    x_max_gain = sign(b)';
    x_max_gain = repmat(x_max_gain,ceil(nx/length(b)),1);
    x_max_gain = x_max_gain(1:nx);


    % Sums of sines
    f0=0.1; f1=2;
    x_sines = sin(2*pi*t*f0) + 0.1*sin(2*pi*t*f1);

    % Chirp
    f_chirp = 1/16;                  % Target frequency
    x_chirp = sin(pi*f_chirp*t.^2);  % Linear chirp

    x = [x_impulse, x_max_gain, x_sines, x_chirp];
    titles = {'Impulse', 'Max gain', 'Sum of sines', 'Chirp'};
    y = zeros(size(x));

    for i=1:size(x,2)
        reset = true;
        y(:,i) = overflow(b,x(:,i),reset);
    end

    test_plot(1,titles,t,x,y)

end
function test_plot(fig,titles,t,x,y1)
    figure(fig)
    clf
    sub_plot = 1;
    font_size = 10;
    for i=1:size(x,2)
        subplot(4,1,sub_plot)
        sub_plot = sub_plot+1;
        plot(t,x(:,i),'c',t,y1(:,i),'k')
        axis('tight')
        xlabel('t','FontSize',font_size);
        title(titles{i},'FontSize',font_size);
        ax = gca;
        ax.FontSize = 10;
    end
    figure(gcf)
end

Создайте coder.FixptConfig объект, fixptcfg, с настройками по умолчанию.

fixptcfg = coder.config('fixpt');

Установите имя испытательного стенда. В этом примере имя функции испытательного стенда overflow_test.

fixptcfg.TestBenchName = 'overflow_test';

Установите размер слова по умолчанию равную 16.

fixptcfg.DefaultWordLength = 16;

Включите обнаружение переполнения.

fixptcfg.TestNumerics = true;
fixptcfg.DetectFixptOverflows = true;

Установите fimath Product mode и Sum mode на KeepLSB. Эти настройки моделируют поведение целочисленных операций в языке C.

fixptcfg.fimath = 'fimath( ''RoundingMethod'', ''Floor'', ''OverflowAction'', ''Wrap'', ''ProductMode'', ''KeepLSB'', ''SumMode'', ''KeepLSB'')';

Создайте объект строения генерации кода, чтобы сгенерировать автономную статическую библиотеку C.

cfg = coder.config('lib');

Преобразуйте функцию MATLAB с плавающей точкой, overflow, в код С с фиксированной точкой. Вам не нужно задавать типы входа для codegen команда, поскольку выводит типы из тестового файла.

codegen -float2fixed fixptcfg -config cfg overflow

Фаза тестирования цифр сообщает о переполнении.

Overflow error in expression 'acc + b( j )*z( k )'. Percentage of Current Range = 104%.

Определите, переполнено ли сложение или умножение в этом выражении. Установите fimath ProductMode для FullPrecision чтобы умножение не переполнено, а затем запустить codegen снова команда.

fixptcfg.fimath = 'fimath(''RoundingMethod'', ''Floor'', ''OverflowAction'', ''Wrap'', ''ProductMode'', ''FullPrecision'', ''SumMode'', ''KeepLSB'')';
codegen -float2fixed fixptcfg -config cfg overflow

Фаза тестирования цифр все еще сообщает о переполнении, что указывает на то, что это сложение в выражении, которое переполнено.

Документация MATLAB Coder

Поддержка