Установите математические атрибуты с фиксированной точкой

В этом примере показано, как задать математические атрибуты с фиксированной точкой в коде MATLAB ®.

Математическими атрибутами с фиксированной точкой для назначения, сложения, вычитания и умножения можно управлять с помощью fimath объект. Можно прикрепить fimath объект в fi объект, использующий setfimath. Можно удалить fimath объект из fi объект, использующий removefimath.

Можно сгенерировать код С из примеров, если у вас есть программное обеспечение MATLAB Coder™.

Установите и удалите математические атрибуты фиксированной точки

Можно изолировать операции с фиксированной точкой от глобальных и локальных fimath настройки при помощи setfimath и removefimath функций. Вы также можете вернуться из функций без fimath присоединен к выходным переменным. Это дает вам локальный контроль над математическими настройками с фиксированной точкой, не вмешиваясь в настройки в других функциях.

Код MATLAB

function y = user_written_sum(u)
    % Setup
    F = fimath('RoundingMethod','Floor',...
        'OverflowAction','Wrap',...
        'SumMode','KeepLSB',...
        'SumWordLength',32);
    u = setfimath(u,F);
    y = fi(0,true,32,get(u,'FractionLength'),F);
    % Algorithm
    for i=1:length(u)
        y(:) = y + u(i);
    end
    % Cleanup
    y = removefimath(y);
end

Выход не имеет присоединенного FIMATH

Когда вы запускаете код, fimath управляет арифметикой внутри функции, но значение возврата не присоединено fimath. Это связано с использованием setfimath и removefimath внутри функции user_written_sum.

>> u = fi(1:10,true,16,11);
>> y = user_written_sum(u)

y =
    55
          DataTypeMode: Fixed-point: binary point scaling
            Signedness: Signed
            WordLength: 32
        FractionLength: 11

Сгенерированный код C

Если у вас есть программное обеспечение MATLAB Coder, можно сгенерировать код С с помощью следующих команд.

>> u = fi(1:10,true,16,11);
>> codegen user_written_sum -args {u} -config:lib -launchreport

Функции fimath, setfimath и removefimath управляйте математикой с фиксированной точкой, но базовые данные, содержащиеся в переменных, не изменяются, и поэтому сгенерированный код C не производит никаких копий данных.

int32_T user_written_sum(const int16_T u[10])
{
  int32_T y;
  int32_T i;
  /* Setup */
  y = 0;
  /* Algorithm */
  for (i = 0; i < 10; i++) {
    y += u[i];
  }
  /* Cleanup */
  return y;
}

Несоответствующий FIMATH

Когда вы оперируете fi объекты, их fimath свойства должны быть равны, или вы получите ошибку.

>> A = fi(pi,'ProductMode','KeepLSB');
>> B = fi(2,'ProductMode','SpecifyPrecision');
>> C = A * B

Error using embedded.fi/mtimes
The embedded.fimath of both operands must be equal.

Чтобы избежать этой ошибки, можно удалить fimath из одной из переменных в выражении. В этом примере fimath удаляется из B в контексте выражения без изменения B сам, и продукт вычисляется с помощью fimath присоединен к A.

>> C = A * removefimath(B)

C =

                6.283203125

           DataTypeMode: Fixed-point: binary point scaling
             Signedness: Signed
             WordLength: 32
         FractionLength: 26

         RoundingMethod: Nearest
         OverflowAction: Saturate
            ProductMode: KeepLSB
      ProductWordLength: 32
                SumMode: FullPrecision

Изменение FIMATH на временных переменных

Если у вас есть переменные без присоединенных fimath, но вы хотите управлять конкретной операцией, тогда можно прикрепить fimath в контексте выражения без изменения переменных.

Например, продукт вычисляется с помощью fimath определяется F.

>> F = fimath('ProductMode','KeepLSB','OverflowAction','Wrap','RoundingMethod','Floor');
>> A = fi(pi);
>> B = fi(2);
>> C = A * setfimath(B,F)

C =

     6.2832

           DataTypeMode: Fixed-point: binary point scaling
             Signedness: Signed
             WordLength: 32
         FractionLength: 26

         RoundingMethod: Floor
         OverflowAction: Wrap
            ProductMode: KeepLSB
      ProductWordLength: 32
                SumMode: FullPrecision
       MaxSumWordLength: 128

Обратите внимание, что переменные B не изменяется.

>> B

B =

           DataTypeMode: Fixed-point: binary point scaling
             Signedness: Signed
             WordLength: 16
         FractionLength: 13

Удаление конфликта FIMATH в цикле

Можно вычислить продукты и суммы, чтобы соответствовать аккумулятору DSP с округлением пола и обернуть переполнение, а также использовать ближайшее округление и насыщение переполнения на выходе. Чтобы избежать несоответствия fimath ошибки, можно удалить fimath на переменном выходе, когда она используется в расчет с другими переменными.

Код MATLAB

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

function [y,z] = setfimath_removefimath_in_a_loop(b,a,x,z)
    % Setup
    F_floor = fimath('RoundingMethod','Floor',...
           'OverflowAction','Wrap',...
           'ProductMode','KeepLSB',...
           'ProductWordLength',32,...
           'SumMode','KeepLSB',...
           'SumWordLength',40);
    F_nearest = fimath('RoundingMethod','Nearest',...
        'OverflowAction','Wrap');
    % Set fimaths that are local to this function
    b = setfimath(b,F_floor);
    a = setfimath(a,F_floor);
    x = setfimath(x,F_floor);
    z = setfimath(z,F_floor);
    % Create y with nearest rounding
    y = coder.nullcopy(fi(zeros(size(x)),true,16,14,F_nearest));
    % Algorithm
    for j=1:length(x)
        % Nearest assignment into y
        y(j) =  b(1)*x(j) + z(1);
        % Remove y's fimath conflict with other fimaths
        z(1) = (b(2)*x(j) + z(2)) - a(2) * removefimath(y(j));
        z(2) =  b(3)*x(j)         - a(3) * removefimath(y(j));
    end
    % Cleanup: Remove fimath from outputs
    y = removefimath(y);
    z = removefimath(z);
end

Инструкции по генерации кода

Если у вас есть программное обеспечение MATLAB Coder, можно сгенерировать код С с заданными аппаратными характеристиками с помощью следующих команд.

N = 256;
t = 1:N;
xstep = [ones(N/2,1);-ones(N/2,1)];
num = [0.0299545822080925  0.0599091644161849  0.0299545822080925];
den = [1                  -1.4542435862515900  0.5740619150839550];

b = fi(num,true,16);
a = fi(den,true,16);
x = fi(xstep,true,16,15);
zi = fi(zeros(2,1),true,16,14);

B = coder.Constant(b);
A = coder.Constant(a);

config_obj = coder.config('lib');
config_obj.GenerateReport = true;
config_obj.LaunchReport = true;
config_obj.TargetLang = 'C';
config_obj.GenerateComments = true;
config_obj.GenCodeOnly = true;
config_obj.HardwareImplementation.ProdBitPerChar=8;
config_obj.HardwareImplementation.ProdBitPerShort=16;
config_obj.HardwareImplementation.ProdBitPerInt=32;
config_obj.HardwareImplementation.ProdBitPerLong=40;

codegen -config config_obj setfimath_removefimath_in_a_loop -args {B,A,x,zi} -launchreport

Сгенерированный код C

void setfimath_removefimath_in_a_loop(const int16_T x[256], int16_T z[2],
  int16_T y[256])
{
  int32_T j;
  int40_T i0;
  int16_T b_y;

  /* Setup */
  /* Set fimaths that are local to this function */
  /* Create y with nearest rounding */
  /* Algorithm */
  for (j = 0; j < 256; j++) {
    /* Nearest assignment into y */
    i0 = 15705 * x[j] + ((int40_T)z[0] << 20);
    b_y = (int16_T)((int32_T)(i0 >> 20) + ((i0 & 524288L) != 0L));

    /* Remove y's fimath conflict with other fimaths */
    z[0] = (int16_T)(((31410 * x[j] + ((int40_T)z[1] << 20)) - ((int40_T)(-23826
      * b_y) << 6)) >> 20);
    z[1] = (int16_T)((15705 * x[j] - ((int40_T)(9405 * b_y) << 6)) >> 20);
    y[j] = b_y;
  }

  /* Cleanup: Remove fimath from outputs */
}

Полиморфный код

Можно записать код MATLAB, который можно использовать как для типов с плавающей точкой, так и для фиксированных точек, используя setfimath и removefimath.

function y = user_written_function(u)
    % Setup
    F = fimath('RoundingMethod','Floor',...
        'OverflowAction','Wrap',...
        'SumMode','KeepLSB');
    u = setfimath(u,F);
    % Algorithm
    y = u + u;
    % Cleanup
    y = removefimath(y);
end

Входные параметры фиксированной точки

Когда функция вызывается с входами с фиксированной точкой, то fimath F используется для арифметики, и выход не присоединен fimath.

>> u = fi(pi/8,true,16,15,'RoundingMethod','Convergent');
>> y = user_written_function(u)

y =

             0.785400390625

           DataTypeMode: Fixed-point: binary point scaling
             Signedness: Signed
             WordLength: 32
         FractionLength: 15

Сгенерированный код C для фиксированной точки

Если у вас есть программное обеспечение MATLAB Coder, можно сгенерировать код С с помощью следующих команд.

>> u = fi(pi/8,true,16,15,'RoundingMethod','Convergent');
>> codegen user_written_function -args {u} -config:lib -launchreport

int32_T user_written_function(int16_T u)
{
  /* Setup */
  /* Algorithm */
  /* Cleanup */
  return u + u;
}

Двойные входы

Начиная с setfimath и removefimath являются транзитными для типов с плавающей точкой, user_written_function пример также работает с типами с плавающей точкой.

function y = user_written_function(u)
    % Setup
    F = fimath('RoundingMethod','Floor',...
        'OverflowAction','Wrap',...
        'SumMode','KeepLSB');
    u = setfimath(u,F);
    % Algorithm
    y = u + u;
    % Cleanup
    y = removefimath(y);
end

Сгенерированный код C для Double

При компиляции с вход с плавающей точкой вы получаете следующий сгенерированный код C.

>> codegen user_written_function -args {0} -config:lib -launchreport

real_T user_written_function(real_T u)
{
  return u + u;
}

Где находится real_T тип определяется как double:

typedef double real_T;

Больше полиморфного кода

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

function y = user_written_sum_polymorphic(u)
    % Setup
    F = fimath('RoundingMethod','Floor',...
        'OverflowAction','Wrap',...
        'SumMode','KeepLSB',...
        'SumWordLength',32);

     u = setfimath(u,F);

     if isfi(u)
         y = fi(0,true,32,get(u,'FractionLength'),F);
     else
         y = zeros(1,1,class(u));
     end

     % Algorithm
     for i=1:length(u)
         y(:) = y + u(i);
     end

     % Cleanup
     y = removefimath(y);

end

Система фиксированных Сгенерированных кодов C

Если у вас есть программное обеспечение MATLAB Coder, можно сгенерировать код С с фиксированной точкой с помощью следующих команд.

>> u = fi(1:10,true,16,11);
>> codegen user_written_sum_polymorphic -args {u} -config:lib -launchreport

int32_T user_written_sum_polymorphic(const int16_T u[10])
{
  int32_T y;
  int32_T i;

  /* Setup */
  y = 0;

  /* Algorithm */
  for (i = 0; i < 10; i++) {
    y += u[i];
  }

  /* Cleanup */
  return y;
}

Система с плавающей точкой Сгенерированного кода C

Если у вас есть программное обеспечение MATLAB Coder, можно сгенерировать код С с плавающей точкой с помощью следующих команд.

>> u = 1:10;
>> codegen user_written_sum_polymorphic -args {u} -config:lib -launchreport

real_T user_written_sum_polymorphic(const real_T u[10])
{
  real_T y;
  int32_T i;

  /* Setup */
  y = 0.0;

  /* Algorithm */
  for (i = 0; i < 10; i++) {
    y += u[i];
  }

  /* Cleanup */
  return y;
}

Где находится real_T тип определяется как double:

typedef double real_T;

SETFIMATH для целых типов

Следуя установленному шаблону обработки встроенных целых чисел, таких как fi объекты, setfimath преобразует целочисленный вход в эквивалентный fi с приложенными fimath.

>> u = int8(5);
>> codegen user_written_u_plus_u -args {u} -config:lib -launchreport

function y = user_written_u_plus_u(u)
    % Setup
    F = fimath('RoundingMethod','Floor',...
        'OverflowAction','Wrap',...
        'SumMode','KeepLSB',...
        'SumWordLength',32);
    u = setfimath(u,F);
    % Algorithm
    y = u + u;
    % Cleanup
    y = removefimath(y);
end

Тип выхода задан как fimath чтобы быть 32-битным.

int32_T user_written_u_plus_u(int8_T u)
{
  /* Setup */
  /* Algorithm */
  /* Cleanup */
  return u + u;
}

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