使用RTW代码生成工具是将MATLAB用于C开发流程的最好的方法（转自恒润科技）

 

使用RTW代码生成工具是将MATLAB用于C开发流程的最好的方法

-->

      嵌入式软件开发人员对MATLAB算法实现和原型设计，及在嵌入式处理器和DSP的C代码实现方面有长期的信赖。作为高级语言，MATLAB便于设计研发。相比较而言，用C编程非常适合DSP在性能、内存和处理能力方面的优化。挑战是如何将从MATLAB灵活的开发环境中的设计转换到相对约束的C编程风格。解决方法是MATLAB可以自动转换为嵌入式C代码。

      手动从MATLAB转换到C需要考虑低层代码的细节，比如数据类型分配、内存分配、对计算负载和内存的优化等。最重要的是要保证MATLAB代码和C代码的一致性。

      当MATLAB算法使用嵌入式MATLAB语言时能够将其准确地转换为C代码，使得用户能够将精力集中在设计上而不是编写C代码上。

      这篇文章概述了手动从MATLAB到C转换的困难，展示了如何使用嵌入式MATLAB进行自动转换，并提供了将MATLAB算法转换的最好方法，从而提高所生成C代码的性能。

Challenges of Translating MATLAB Concept Code into Implementation Code
MATLAB has several advantages for design exploration, such as polymorphism, matrix-based functions, and an interactive programming environment. During translation of an algorithm from MATLAB to C, however, software designers face some important constraints. For example:

MATLAB is a dynamically typed and C is a statically typed language. When writing a MATLAB program, you do not need to define data types and sizes for your variables. While this flexibility makes it easy to develop algorithms as proofs of concept, when it comes time for translation to C, the programmer must assign appropriate data types and sizes to all variables.

MATLAB is polymorphic. Functions in MATLAB can process different types of input parameters and can apply a different algorithm to each type of parameter. For example, the abs function computes the absolute value of real numbers and norm of complex numbers and can process scalars, vectors, or matrices.

>> abs(4-3i)
ans =
    5
   
>> abs([4 -3])
ans =
    4 3

This kind of flexibility is not supported in C, which assigns a single algorithm to each parameter type. To translate a polymorphic MATLAB function to C, the programmer must maintain separate function prototype for each possible parameter signature.

MATLAB is based on compact matrix notation. Most MATLAB expressions containing vectors and matrices are compact, single-line expressions similar to the corresponding mathematical formula. The equivalent C code requires iterators, such as for loops, to express the matrix operations as a sequence of scalar computations.

Automating the MATLAB to C Workflow
Automatic MATLAB to C conversion with Real-Time Workshop® addresses many of the challenges outlined in the previous section. For example, consider an algorithm depicted in the function euclidean.m. This algorithm minimizes the Euclidean distance between a column vector x and a collection of column vectors contained in the matrix cb. The function has two output variables: y, the vector in cb with the minimum distance to x, and dist, the value of the minimum distance (Figure 1). 

 
Figure 1. Source code for the euclidean.m function.

In the body of the Euclidean function, we use the MATLAB function norm to compute the distance between x and each column vector of cb. To visualize the computed distances between any pair of points, we call the plot_distances function inside the loop. We use the %#eml directive to turn on the MATLAB M-Lint code analyzer and check the function code for errors and recommend corrections.

Figure 2 shows how the plot_distances function helps us visualize the process of computing all distances in a 2D space and finding the minimum value.

 
Figure 2. Visualizing the computation of Euclidean distances by the plot_distances function. Click on image to see enlarged view.

To generate C code from this algorithm, we must use only operators and functions that are part of the Embedded MATLAB subset. Visualization functions, such as plot, line, and grid, are not supported by the Embedded MATLAB subset. When you open the euclidian.m function in the MATLAB editor, the M-Lint code analyzer identifies and reports on these unsupported lines of code.

While it makes sense to use visualization functions to debug and verify the algorithm, when we implement the algorithm as C code, we must separate the offline analysis portions of the design from the online portions involved in embedded C code generation. In our example, we can make the algorithm compliant with the Embedded MATLAB subset simply by identifying and commenting out the plot_distances function.

Now we use the emlc command in Real-Time Workshop to generate C code for the Embedded MATLAB compliant function euclidean.m.

Typical syntax for translation is

>> emlc -eg {x,cb} –report euclidean.m

The example option (following the –eg delimiter) sets the data types and dimensions of function variables by specifying an example at the function interface. The –report option opens the Embedded MATLAB compilation report with hyperlinks to C source files and header files generated from the MATLAB function. The generated C code is created in a file named euclidean.c (Figure 3).

 
Figure 3. C code generated automatically from the euclidean.m function displayed inside the Embedded MATLAB compilation report. Click on image to see enlarged view.

By using the example option, we declare the data type, size, and complexity of the variables x and cb in the function interface, enabling the Embedded MATLAB engine to assign data type and sizes automatically to all the local variable in the generated C program. The generated C code correctly maps to zero-based indexing for accessing array elements, and the vector operations are automatically mapped to scalar computations with for loops. As a result, many difficulties encountered in manual MATLAB to C conversion are eliminated through automatic translation.

Design Patterns for a MATLAB to C Workflow
In an embedded system, the size and data type of each variable must be set before implementation. In addition, if the performance requirements are not met, the algorithm’s memory and computational footprint must be optimized. The following sections examine design patterns that use supported Embedded MATLAB features to ensure that the generated C code adheres to these requirements.

Accommodating Changes in Variable Dimensions
In the MATLAB language, all data can vary in size. Embedded MATLAB supports variable-sized arrays and matrices with known upper bounds. This means you can accommodate the use of variable-sized data for embedded implementations by using buffers of a constant maximum size and by addressing subportions of the constant-size buffers. Within your Embedded MATLAB functions you can define variable-size inputs, outputs, and local variables, with known upper bounds. For inputs and outputs, you must specify the upper bounds explicitly at the function interface. For local data, Embedded MATLAB uses in-depth analysis to calculate the upper bounds at compile time. However, when the analysis fails to detect an accurate upper bound, you must specify them explicitly for local variables.

We update the euclidean.m function to accommodate changes in dimensions over which we compute the distances.We want to compute only the distance between first N elements of a given vector x with the first N elements of every column vector contained in the matrix cb. The resulting function, euclidean_varsize.m, will have a third input argument, N (Figure 4).

 
Figure 4. Input argument for euclidean_varsize.m

The compilation of this function will result in errors because we have not yet specified an upper bound for the value of N. As a result, the local variable 1:N will have no specified upper bound. To impose the upper bound we can constrain the value of the parameter N in the first line of the function by using the assert function with relational operators (Figure 5).

 
Figure 5. Using the assert function with relational operators.

The function varsize_example.m shows another common pattern, one where an array such as y is first initialized and then grows in size based on a condition related to the value of an input or local variable:

 
The compilation of this function will again result in errors since we have not specified an upper bound for the variable y. To accommodate this type of size change for the local variable y, we can specify the upper bound using the eml.varsize function for all instances of that local variable. In this example we constrain a maximum dimension of 1-by-8 for the variable y. In the upper branch of the if statement, the variable y has a dimension of 1-by-6 and in the lower branch, a dimension of 1-by-7. The resulting function will be

 
Optimizing for Memory or Computational Complexity
Another common refinement is to optimize the generated C code for memory and complexity. In our euclidean.m algorithm, to compute the distance between two points, norm takes the square root of the squared values of each element of a given vector. Because computing the square root is computationally expensive, we can update our function by computing only the sum of the squared elements, without any loss in the intended behavior. The resulting function, which has a much lower computational load, uses the sum function supported by the Embedded MATLAB language subset (Figure 6).

 
Figure 6. Optimized euclidean.m function.

It may be desirable to reduce memory footprint of the generated C code. In some cases, the initialization of new variables in your Embedded MATLAB function may produce redundant copies in the generated C code. Although Embedded MATLAB technology eliminates many copies automatically, you can eliminate data copies that are not automatically handled by declaring uninitialized variables using the eml.nullcopy function. In Figure 7, the variable y is initialized with such a construction.

 
Figure 7. Initializing the variable y using the eml.nullcopy function.

Using Fixed-Point and Native Integer Data Types
By default, MATLAB uses 64-bit double-precision numerical representation for variables created in the workspace. As convenient as this choice is for design exploration, it is not memory-efficient for real-time processing of many common signals, such as image or audio signals represented natively with word lengths of 8 or 16 bits. To handle these types of signals and to implement your MATLAB algorithm on target processors with limited word lengths, you must convert the design to a fixed-point or integer-based representation. You can use Fixed-Point Toolbox™ to create fixed-point variables and perform fixed-point computations. Since the Embedded MATLAB language subset supports the fixed-point data object (fi), by using Real-Time Workshop you can generate pure integer C code from your Embedded MATLAB code. This usually involves modifying your original MATLAB function to declare variables based on integer or fixed-point representations.

Our euclidean_optimized.m function can process integer data types or fixed-point data types as its input variables. To generate C code we only need to compile the same function with integer or fixed-point variables in the example option of the emlc command. The command syntax for input variables of 16-bit signed integer type, for example, will be

>> emlc -eg {int16(x),int16(cb)} –report euclidean_optimized.m
The resulting generated C code contains only integer C data types, and can be readily compiled into fixed-point processors (Figure 8). 

Figure 8. C code for euclidean_optimized.m compiled with integer input variables. Computations are purely integer-based. Click on image to see enlarged view.

A Common Language and Development Environment
Automatic translation of MATLAB to C with the Embedded MATLAB subset eliminates the need to produce, maintain, and verify hand written C code. Design iterations become easier, as you stay within the MATLAB environment and take advantage of its interactive debugging and visualization capabilities. Many desirable features of MATLAB programs, such as matrix-based operations, polymorphism, variable-size data, and fixed-point numerical representations, are automatically translated to C code, enabling you to focus on improving your design rather than maintaining multiple copies of the source code written in different languages.

Products Used

• MATLAB®
• Fixed-Point Toolbox™
• Real-Time Workshop®
• Signal Processing Toolbox™

For More Information

• Convert MATLAB Code to Embedded C Using Embedded MATLAB
• From MATLAB to Embedded C. DSP Design Line