苏勇,[email protected],2022年12月
本文在于记录使用Simulink建模仿真,到最终部署到微控制器上的过程。从Simulink算法模型到目标硬件平台可以支持的可执行文件,大体经过Matlab的各种模型描述语言,编译成通用C语言,再在目标硬件平台的编译环境中编译生成最终的可执行文件。市面上的相关的学术研究和实践经验,大多是使用了MATLAB已经官方认证集成的BSP,例如:Arduino、TI DSP、STM32、树莓派等,在这些工具包的加持下,MATLAB配合具体的目标平台的集成开发环境,可以一通到底。但实际上,MATLAB只要在将Simulink工具描述的算法模型转化成C语言,完全可以由开发者自行进一步适配的具体平台,而不需要依赖于具体的官方认证的设备支持包。更进一步,适用于Arm兼容的微控制器,Matlab已经集成了一个通用的模板,在将Matlab算法语言转化成C语言的过程中特别考虑在嵌入式平台上执行的代码优化。
本文详细演绎了从Simulink创建模型,仿真验证,之后再生成C源码部署到plus-f5270开发板的全过程。 其中描述的方法和操作步骤,为后续部署更多仿真模型奠定了基础。
从Simulink中以符号图形和算法框图表示的仿真模型,转化成可编译的C源码,再到可以运行在嵌入式平台的可执行程序,经历了多个阶段的处理,在每个阶段需要对应的工具,将上一阶段的输出文件作为本阶段的输入,经过处理后再输出给后续的环节。
图x完整展示了MATLAB从算法模型到可执行程序的各个环节:
使用MATLAB、Simulink和Stateflow工具可以在MATLAB中创建模型,此时的模型可能是以逻辑图符号的方式表示,保存成model.mdl文件
RealTime Workshop(RTW)是MATLAB中的一个工具,是用来做仿真的软件开发环境,可以使用描述模型符号及其逻辑关系的model.mdl文件作为输入,经过本地编译,形成描述模型的model.rtw文件。RTW文件是以ASCII码的形式存放,开发者可以直接读其中的内容。
通过TLC目标语言编译器生成的C源码文件model.c文件,可以直接加入到嵌入式系统的集成开发环境当中(例如Keil)参与编译。但实际上,这里生成的C源码文件主要是描述算法的实现,是应用程序的框架,仍需要开发者在具体嵌入式硬件平台上,结合实际应用,手动向其中嵌入具体的输入输出函数。
MATLAB的工具生成算法执行的源码框架,配合目标嵌入式硬件平台的适配移植,最终在嵌入式硬件平台的编译环境中完成编译链接,最终生成可以在目标嵌入式硬件平台上运行的可执行文件。后续可以使用目标嵌入式硬件平台的一般操作过程完成下载和调试。
MathWorks ®
Mathworks公司的MATLAB软件环境中有3个代码生成工具:Matlab Coder、Simulink Coder和Embedded Coder。MATLAB可以将M语言转化成目标代码,Simulink Coder可以将Simulink中的仿真模型转化成目标代码,而Embedded Coder 依赖于MATLAB Coder以及Simulink Coder的,通过进一步优化,主要生成用于部署在嵌入式平台上的源码。
考虑最终将生成的C代码部署到嵌入式平台上,从MATLAB程序和从Simulink模型生成源码的工具流有一些区别,如图x所示。
此处以实现一个流水灯的模型为例,说明基于模型设计的全过程。
在MATLAB中启动Simulink后,可通过工具栏的“建模”->“模型设置”,激活“配置参数”对话框。如图x所示。
在“配置参数”对话框中,首先配置“求解器”。求解器是Simulink用数值方法结算模型切片状态的引擎,对应于Simulink的仿真过程的行为。这里考虑将来是在嵌入式系统中作为周期任务更新求解过程,设定在Simulink的仿真行为同嵌入式系统类似,为“离散”的“固定步长”,计算步长的时间为10ms,对应于配置对话框中的0.01s。如图x所示。在实际的嵌入式系统平台运行求解器时,将使用一个硬件定时器,以10ms的周期触发中断,执行求解器的程序。
在“配置参数”对话框的“硬件实现”标签中,选择即将生成代码的目标平台。考虑到未来将会在ARM微控制器的平台上部署(但不在Simulink的支持设备清单中),这里将“设备供应商”选择为“ARM Compatible”,选择“设备”类型为“ARM Cortex-M”。特别注意将“Code Generation system target file”选择成为“ert.tlc”,这是最终生成C语言的配置文件。
实际上,在配置TLC时,配置参数对话框会自动切换到“代码生成”页面,在这里选择使用etc.tlc文件。如图x所示。
从图x中可以看到,Simulink提供了多种TLC配置文件可选。这里选择的ert.tlc,是专门适用于Embedded Coder生成嵌入式系统C代码的。
配置好项目之后,保存模型文件到文件系统中。从图x中可以看到,Simulink保存的模型文件可以有两种后缀名,一种是“.slx”,另一种是“.mdl”,这里,我们使用默认的“.slx”文件后缀名,保存成“leds_model.slx”文件。
在plus-f5270开发板上设计了4个可编程的LED灯,有原理图如图x所示。
这里要设计一个逻辑,用一个计数器驱动小灯的亮暗状态,当计数器输出0时,控制仅LED1亮,当计数器输出1时,控制仅LED2亮,以此类推。计数器从0到3周期计数,小灯也依次闪烁。
在“库浏览器”中选择需要的符号元件,拖放到绘图工作区中,例如,常量“Constant”、等于号“Equal”、有限计数器“Counter Limited”、输出端点“Out”等,并可用示波器“Scope”和显示器“Display”查看仿真输出结果。使用连线,将符号的输入输出信号连接在一起,形成一个有输入输出的完整系统的功能框图。如图x所示。
从图x中可以看出,计数器的输出值驱动4个LED灯的亮暗状态。但默认情况下,计数器的计数节奏是同求解器的计算步长保持一致的,如之前设定求解器的计算步长是10ms,这就意味着这些小灯会闪烁得非常快,难以被人眼观察到。这里希望调整计数器的计数节奏为1m,需要调整计数器的计数节奏。双击计数器的图标,可以激活计数器模块的参数配置对话框,配置其中的“采样时间”值为1。如图x所示
此时,使用计数器控制小灯闪烁的流水灯模型就已经搭建完成。
在仿真的工具栏中,点击“步进”即可逐步运行,点击“运行”可以全速仿真运行。如图x所示。
启动仿真后,可以通过系统框图中的示波器模块查看仿真结果,如图x所示。
从图x中显示的示波器输出波形可以看到,4个输出通道以1s为控制节奏,依次输出高电平。验证仿真结果达到预期设计。
经过仿真验证的系统设计框图,可以直接转化成嵌入式系统可以可以执行的C程序源码。在调用RTW编译模型生成C源码之前,必须再次确认几个配置要点:
然后,试着生成C源码。在Simulink的界面中,在工具栏“APP”页面中选择“Embedded Coder”激活“C代码”标签页,然后在“C代码”标签页中其中选择“生成待代码”。如图x所示。
此时,Simulink的目标语言编译器TLC将启动编译,并生成C源码程序文件。如图x所示。
至此,在Simulink中调用Embedded Coder,已经将创建的leds_model模型转化成了C源码,保存在一系列源码文件中。此时,可以在预先设定的当前工作目录中,找到新建保存生成源码文件的文件夹“leds_model_ert_rtw”,顾名思义:leds_model是模型的名字;ert是Embedded Real-Time Target,对应源码生成器的配置选项;rtw是Real-Time Workshop,这是源码生成器的工具包,在MATLAB的早期版本中独立作为一个工具包,现在已经被集成在源码生成器内部。
这里准备了一个plus-f5270开发板上可以运行的工程,plus-f5270_hello_world_mdk,作为部署模型的基础。
首先,将工程的目录名改为“plus-f5270_leds_model_mdk” ,便于标识该工程专用于部署leds_model模型。在该工程根目录下创建“model”目录,将由Embedded Coder生成的模型C源码文件全部复制到其中。此时,工程目录的源码结构如图x所示。
然后,在Keil工程中添加模型源码,先试着编译一下,验证模型文件本身源码的完整性。如图x所示。
经编译验证可知,模型源码本身是闭包无误的,不会为基础工程引入额外的错误,也不需要专门针对嵌入式编译环境调整任何源码。在Keil中编译模型文件需要注意:
/*
* File: ert_main.c
*
* Code generated for Simulink model 'leds'.
*
* Model version : 1.4
* Simulink Coder version : 9.8 (R2022b) 13-May-2022
* C/C++ source code generated on : Mon Dec 19 18:09:56 2022
*
* Target selection: ert.tlc
* Embedded hardware selection: ARM Compatible->ARM Cortex-M
* Code generation objectives: Unspecified
* Validation result: Not run
*/
#include
#include /* This example main program uses printf/fflush */
#include "leds_model.h" /* Model header file */
/*
* Associating rt_OneStep with a real-time clock or interrupt service routine
* is what makes the generated code "real-time". The function rt_OneStep is
* always associated with the base rate of the model. Subrates are managed
* by the base rate from inside the generated code. Enabling/disabling
* interrupts and floating point context switches are target specific. This
* example code indicates where these should take place relative to executing
* the generated code step function. Overrun behavior should be tailored to
* your application needs. This example simply sets an error status in the
* real-time model and returns from rt_OneStep.
*/
void rt_OneStep(void);
void rt_OneStep(void)
{
static boolean_T OverrunFlag = false;
/* Disable interrupts here */
/* Check for overrun */
if (OverrunFlag) {
rtmSetErrorStatus(leds_model_M, "Overrun");
return;
}
OverrunFlag = true;
/* Save FPU context here (if necessary) */
/* Re-enable timer or interrupt here */
/* Set model inputs here */
/* Step the model */
leds_model_step();
/* Get model outputs here */
/* Indicate task complete */
OverrunFlag = false;
/* Disable interrupts here */
/* Restore FPU context here (if necessary) */
/* Enable interrupts here */
}
/*
* The example main function illustrates what is required by your
* application code to initialize, execute, and terminate the generated code.
* Attaching rt_OneStep to a real-time clock is target specific. This example
* illustrates how you do this relative to initializing the model.
*/
int_T main(int_T argc, const char *argv[])
{
/* Unused arguments */
(void)(argc);
(void)(argv);
/* Initialize model */
leds_model_initialize();
/* Attach rt_OneStep to a timer or interrupt service routine with
* period 0.001 seconds (base rate of the model) here.
* The call syntax for rt_OneStep is
*
* rt_OneStep();
*/
printf("Warning: The simulation will run forever. "
"Generated ERT main won't simulate model step behavior. "
"To change this behavior select the 'MAT-file logging' option.\n");
fflush((NULL));
while (rtmGetErrorStatus(leds_model_M) == (NULL)) {
/* Perform application tasks here */
}
/* Terminate model */
leds_model_terminate();
return 0;
}
/*
* File trailer for generated code.
*
* [EOF]
*/
在plus-f5270工程的boards目录下创建model_port.c文件,集中存放基于plus-f5270开发板适配leds_model模型的配置源码。在model_port.c文件中创建model_init()函数,在顶级调用main()函数中将调用model_init()函数初始化并启动模型求解器。同时,在model_port.c文件中创建model_input()和model_output()函数,用于同步模型同硬件的绑定关系。
Note:在模型中使用的输入和输出,都是保存在内存中的数据。在微控制器芯片内部的电路系统中,通过访问地址映射空间中的寄存器,控制的状态和行为。在面向嵌入式平台设计的仿真系统中,生成的C源码,是不会直接控制具体嵌入式平台的电路系统的。这个比较好理解,因为仿真系统是同一份,但运行在不同的嵌入式平台上,对应的寄存器是不同的。这就需要嵌入式平台的开发者自行适配代码,将电路系统中状态量和控制量,与仿真模型的输入和输出进行同步。
在plus-f5270开发板上,选用硬件定时器TIM4产生周期为10ms的中断服务程序,在中断服务中周期调用模型生成的求解器处理函数,执行模型。
在model_init()函数中配置硬件定时器TIM4的代码如下:
/* model_port.c */
#include "board_init.h"
#include
#include /* This example main program uses printf/fflush */
#include "leds_model.h" /* Model header file */
extern void leds_model_initialize(void);
extern void leds_model_step(void);
void model_init(void)
{
/* prepare the models. */
leds_model_initialize();
/* setup timer. */
TIM_Init_Type tim_init;
tim_init.ClockFreqHz = BOARD_MODEL_TIMER_CLKSRC_HZ;
tim_init.StepFreqHz = 1000000; /* 1mhz per step. */
tim_init.Period = 10000; /* 100hz for a round. */
tim_init.EnablePreloadPeriod = true;
tim_init.PeriodMode = TIM_PeriodMode_Continuous;
tim_init.CountMode = TIM_CountMode_Increasing;
TIM_Init(BOARD_MODEL_TIMER_PORT, &tim_init);
/* enable interrupt. */
NVIC_EnableIRQ(BOARD_MODEL_TIMER_NVIC_IRQn);
TIM_EnableInterrupts(BOARD_MODEL_TIMER_PORT, TIM_INT_UPDATE_PERIOD, true);
/* start timer. */
TIM_Start(BOARD_MODEL_TIMER_PORT);
}
/* bind the model input to hardware. */
void model_input(void)
{
}
/* bind the model output to hardware. */
void model_output(void)
{
}
/* hardware timer interrupt for model solver. */
void BOARD_MODEL_TIMER_NVIC_IRQHandler(void)
{
uint32_t flags = TIM_GetInterruptStatus(BOARD_MODEL_TIMER_PORT);
if (flags & TIM_STATUS_UPDATE_PERIOD)
{
model_input();
leds_model_step(); /* call the model solver step by step. */
model_output();
}
TIM_ClearInterruptStatus(BOARD_MODEL_TIMER_PORT, flags);
}
/* EOF. */
当在model_init()函数中调用TIM_Start()函数之后,硬件定时器开始工作,自动周期触发定时器中断服务,由硬件中断服务程序自动周期调用模型生成的求解器程序leds_model_step()。在定时器中断服务程序中调用求解器之前,需要手动执行在适配过程中创建的model_input()函数,将嵌入式平台电子系统中的状态量同步到模型内部的输入变量中;在调用求解器之后,需要手动执行在适配过程中创建的model_output()函数,将模型求解器更新的输出变量的状态同步到嵌入式平台电子系统中,控制电路完成相关的操作。
具体在本例中使用定时器控制LED灯,仅需要在model_output()函数中实现根据模型输出状态控制LED灯的亮暗状态。但嵌入式平台进行编程时需要注意,使用GPIO控制LED灯之前,需要对将要使用到的GPIO硬件资源预先初始化,配置好对应的外设时钟、引脚复用功能,以及初始状态等。
在model_port.c文件中的model_init()函数中,添加配置GPIO引脚控制LED灯全部熄灭,在model_output()函数中,根据模型中表示LED灯状态的变量leds_model_Y.ledn(n=1, 2, 3, 4)配置到控制LED灯的GPIO引脚输出信号上。
void model_init(void)
{
...
/* setup leds. */
GPIO_WriteBit(BOARD_LED1_GPIO_PORT, BOARD_LED1_GPIO_PIN, 1u);
GPIO_WriteBit(BOARD_LED2_GPIO_PORT, BOARD_LED2_GPIO_PIN, 1u);
GPIO_WriteBit(BOARD_LED3_GPIO_PORT, BOARD_LED3_GPIO_PIN, 1u);
GPIO_WriteBit(BOARD_LED4_GPIO_PORT, BOARD_LED4_GPIO_PIN, 1u);
...
}
/* bind the model output to hardware. */
void model_output(void)
{
GPIO_WriteBit(BOARD_LED1_GPIO_PORT, BOARD_LED1_GPIO_PIN, (leds_model_Y.led1?0u:1u) );
GPIO_WriteBit(BOARD_LED2_GPIO_PORT, BOARD_LED2_GPIO_PIN, (leds_model_Y.led2?0u:1u) );
GPIO_WriteBit(BOARD_LED3_GPIO_PORT, BOARD_LED3_GPIO_PIN, (leds_model_Y.led3?0u:1u) );
GPIO_WriteBit(BOARD_LED4_GPIO_PORT, BOARD_LED4_GPIO_PIN, (leds_model_Y.led4?0u:1u) );
}
在board_init.h文件中定义LED和硬件定时器的具体资源,并添加对model_init()函数的声明。board_init.h文件也是main.c中包含硬件相关资源的统一入口,在适配模型的过程中所需要的其他软硬件相关的组件,也都可以在board_init.h文件的头部声明引用。
#ifndef __BOARD_INIT_H__
#define __BOARD_INIT_H__
#include
#include
#include "hal_common.h"
#include "hal_rcc.h"
#include "hal_uart.h"
#include "hal_tim.h"
#include "hal_gpio.h"
#include "clock_init.h"
#include "pin_init.h"
/* DEBUG UART. */
#define BOARD_DEBUG_UART_PORT UART1
#define BOARD_DEBUG_UART_BAUDRATE 115200u
#define BOARD_DEBUG_UART_FREQ CLOCK_APB2_FREQ
/* leds. */
#define BOARD_LED1_GPIO_PORT GPIOI
#define BOARD_LED1_GPIO_PIN GPIO_PIN_0
#define BOARD_LED2_GPIO_PORT GPIOD
#define BOARD_LED2_GPIO_PIN GPIO_PIN_2
#define BOARD_LED3_GPIO_PORT GPIOB
#define BOARD_LED3_GPIO_PIN GPIO_PIN_14
#define BOARD_LED4_GPIO_PORT GPIOC
#define BOARD_LED4_GPIO_PIN GPIO_PIN_9
/* model timer. */
#define BOARD_MODEL_TIMER_PORT (TIM_Type *)TIM4
#define BOARD_MODEL_TIMER_CLKSRC_HZ CLOCK_APB1_FREQ
#define BOARD_MODEL_TIMER_NVIC_IRQn TIM4_IRQn
#define BOARD_MODEL_TIMER_NVIC_IRQHandler TIM4_IRQHandler
void BOARD_Init(void);
void model_init(void);
#endif /* __BOARD_INIT_H__ */
在clock_init.c文件中的BOARD_InitBoardClocks()函数中添加代码,启用GPIOB、GPIOC、GPIOD、GPIOI,以及TIM4外设模块的总线时钟。
void BOARD_InitBootClocks(void)
{
CLOCK_ResetToDefault();
CLOCK_BootToHSE120MHz();
/* UART1. */
RCC_EnableAPB2Periphs(RCC_APB2_PERIPH_UART1, true);
RCC_ResetAPB2Periphs(RCC_APB2_PERIPH_UART1);
/* GPIOB. */
RCC_EnableAHB1Periphs(RCC_AHB1_PERIPH_GPIOB, true);
RCC_ResetAHB1Periphs(RCC_AHB1_PERIPH_GPIOB);
/* GPIOC. */
RCC_EnableAHB1Periphs(RCC_AHB1_PERIPH_GPIOC, true);
RCC_ResetAHB1Periphs(RCC_AHB1_PERIPH_GPIOC);
/* GPIOD. */
RCC_EnableAHB1Periphs(RCC_AHB1_PERIPH_GPIOD, true);
RCC_ResetAHB1Periphs(RCC_AHB1_PERIPH_GPIOD);
/* GPIOI. */
RCC_EnableAHB1Periphs(RCC_AHB1_PERIPH_GPIOI, true);
RCC_ResetAHB1Periphs(RCC_AHB1_PERIPH_GPIOI);
/* TIM4. */
RCC_EnableAPB1Periphs(RCC_APB1_PERIPH_TIM4, true);
RCC_ResetAPB1Periphs(RCC_APB1_PERIPH_TIM4);
}
在pin_init.c文件中的BOARD_InitPins()函数中添加代码,配置LED灯对应的PI0、PD2、PB14、PC9引脚的复用功能为GPIO输出。
void BOARD_InitPins(void)
{
/* PB7 - UART1_TX. */
GPIO_Init_Type gpio_init;
gpio_init.Pins = GPIO_PIN_6;
gpio_init.PinMode = GPIO_PinMode_AF_PushPull; //GPIO_PinMode_AF_PushPull
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOB, &gpio_init);
GPIO_PinAFConf(GPIOB, gpio_init.Pins, GPIO_AF_7);
/* PB6 - UART1_RX. */
gpio_init.Pins = GPIO_PIN_7;
gpio_init.PinMode = GPIO_PinMode_In_Floating; //GPIO_PinMode_In_Floating
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOB, &gpio_init);
GPIO_PinAFConf(GPIOB, gpio_init.Pins, GPIO_AF_7);
/* PI0 - led1. */
gpio_init.Pins = GPIO_PIN_0;
gpio_init.PinMode = GPIO_PinMode_Out_PushPull;
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOI, &gpio_init);
GPIO_PinAFConf(GPIOI, gpio_init.Pins, GPIO_AF_15);
/* PD2 - led2. */
gpio_init.Pins = GPIO_PIN_2;
gpio_init.PinMode = GPIO_PinMode_Out_PushPull;
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOD, &gpio_init);
GPIO_PinAFConf(GPIOD, gpio_init.Pins, GPIO_AF_15);
/* PB14 - led3. */
gpio_init.Pins = GPIO_PIN_14;
gpio_init.PinMode = GPIO_PinMode_Out_PushPull;
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOB, &gpio_init);
GPIO_PinAFConf(GPIOB, gpio_init.Pins, GPIO_AF_15);
/* PC9 - led4. */
gpio_init.Pins = GPIO_PIN_9;
gpio_init.PinMode = GPIO_PinMode_Out_PushPull;
gpio_init.Speed = GPIO_Speed_50MHz;
GPIO_Init(GPIOC, &gpio_init);
GPIO_PinAFConf(GPIOC, gpio_init.Pins, GPIO_AF_15);
}
模型的求解器是由硬件定时器中断自动触发执行的,除了模型的初始化函数mode_init()之外,不需要嵌入式平台的开发者在main() 函数中人为调用任何模型生成的函数,并且也不影响原有程序的正常执行。实际上,本例工程中的main.c,在原有hello_world工程的main.c中,仅添加了一行model_init()函数的调用。
int main(void)
{
uint8_t ch;
BOARD_Init();
printf("hello, world\r\n");
/* init & start the model solver. */
model_init();
while (1)
{
ch = getchar();
putchar(ch);
}
}
在嵌入式平台上编译工程和下载的操作同一般嵌入式开发无异。在Keil中编译,下载可执行程序到plus-f5270开发板中。
下载可执行程序到到开发板上,运行,可以观察到开发板上的4个LED小灯以1s为节奏演示流水灯的效果。如图x所示。