📚 课程目录第3章定点数转换在代码移植中的应用

3. 定点数转换在代码移植中的应用

任务定义

在dSPACE快速控制原型（RCP）开发中，Simulink模型默认使用双精度浮点数，而量产ECU的嵌入式目标通常仅支持定点运算或需要定点实现以降低开销。本案例完成一个可独立运行的定点数转换工具：将浮点信号按指定定标（字长与小数位长）转换为Q格式定点数，并支持反向转换回浮点，同时给出转换误差。该工具用于代码移植阶段验证定标选择是否满足精度要求。

定点数转换实现

#include <stdio.h>
#include <stdint.h>
#include <math.h>

// Q格式转换：浮点转定点（四舍五入）
int32_t float_to_q15(float val, int frac_bits) {
    float scaled = val * (float)(1 << frac_bits);
    return (int32_t)lrintf(scaled);  // 四舍五入到最近整数
}

// Q格式转换：定点转浮点
float q15_to_float(int32_t q_val, int frac_bits) {
    return (float)q_val / (float)(1 << frac_bits);
}

// 检查定点数是否溢出（超出可表示范围）
int check_overflow(float val, int word_bits, int frac_bits) {
    int int_bits = word_bits - frac_bits - 1;  // 1位符号位
    float max_val = (float)((1 << int_bits) - 1) + 
                    (float)((1 << frac_bits) - 1) / (float)(1 << frac_bits);
    float min_val = -(float)(1 << int_bits);
    return (val > max_val || val < min_val);
}

int main() {
    // 测试用例：模拟dSPACE模型中的典型信号
    float test_signals[] = {0.5f, -0.25f, 1.0f, 0.1f, -0.75f};
    int num_signals = sizeof(test_signals) / sizeof(test_signals[0]);
    
    // Q15格式：16位字长，15位小数位（范围[-1, 0.99997]）
    int word_bits = 16;
    int frac_bits = 15;
    
    printf("=== Q%d.%d 定点转换测试 ===\n", word_bits - frac_bits - 1, frac_bits);
    printf("%-12s %-12s %-12s %-12s\n", "浮点输入", "定点值", "转换回浮点", "误差");
    
    for (int i = 0; i < num_signals; i++) {
        float input = test_signals[i];
        
        if (check_overflow(input, word_bits, frac_bits)) {
            printf("%-12.4f 溢出! 超出Q%d.%d范围\n", input, 
                   word_bits - frac_bits - 1, frac_bits);
            continue;
        }
        
        int32_t q_val = float_to_q15(input, frac_bits);
        float round_trip = q15_to_float(q_val, frac_bits);
        float error = round_trip - input;
        
        printf("%-12.4f %-12d %-12.4f %-12.6f\n", 
               input, q_val, round_trip, error);
    }
    
    // 验证边界值
    printf("\n=== 边界值验证 ===\n");
    float boundary[] = {0.99997f, -1.0f, 1.0f};
    for (int i = 0; i < 3; i++) {
        if (check_overflow(boundary[i], word_bits, frac_bits)) {
            printf("%-12.4f 溢出!\n", boundary[i]);
        } else {
            int32_t q = float_to_q15(boundary[i], frac_bits);
            printf("%-12.4f -> 定点 %d, 回转 %.5f\n", 
                   boundary[i], q, q15_to_float(q, frac_bits));
        }
    }
    
    return 0;
}

编译与运行

gcc -o fixed_point fixed_point.c -lm
./fixed_point

预期输出

=== Q1.15 定点转换测试 ===
浮点输入     定点值       转换回浮点   误差
0.5000       16384        0.5000       0.000000
-0.2500      -8192        -0.2500      0.000000
1.0000       溢出! 超出Q1.15范围
0.1000       3277         0.1000       0.000031
-0.7500      -24576       -0.7500      0.000000

=== 边界值验证 ===
0.99997      32767        0.99997      0.000000
-1.0000      -32768       -1.0000      0.000000
1.0000       溢出!

验证方法

运行程序后观察：

0.5、-0.25、-0.75等可精确表示的数值误差为0；
0.1f在Q15格式下产生约3.05e-5的量化误差（2^-15），符合理论预期；
1.0超出Q1.15最大表示范围0.99997，程序正确报告溢出；
边界值0.99997和-1.0可正常转换且无误差。

判定标准：若上述输出与代码逻辑一致，且误差不超过2^-15，则定点转换工具验证通过，可用于代码移植前的定标精度评估。

代码移植应用说明

在dSPACE原型到产品代码移植中，本工具用于：

定标验证：对模型中每个浮点信号，确定合适的字长与小数位长，运行转换检查误差是否在控制精度允许范围内；
溢出检测：识别超出定点表示范围的信号，提示需要调整定标或信号缩放；
代码生成对照：将转换结果与dSPACE TargetLink或Embedded Coder生成的定点代码输出对比，验证移植一致性。

实际移植时，将float_to_q15与q15_to_float替换为编译器内置定点转换或手写汇编，并将frac_bits改为对应目标芯片的Q格式即可。