自己动手写C语言浮点数转换字符串函数

前几天，应一个小友要求，写了几个字符串转换函数帮助其进行语言学习，自觉其中的几个函数还比较满意，故发布在此，可供初学者参考。

浮点数转换字符串函数说简单也简单，说麻烦，也够烦心的，关键看你如何写了。简单点的几十行代码就行，复杂点的恐怕就的几百行代码了。如果还要考虑移植性、可读性甚至可维护性等就更麻烦的了。我一贯认为，一些事务性的项目应着重考虑移植性、可读性和可维护性等，而常用的带点系统性质的函数代码就应该以执行效率为主。

本文的浮点数转换字符串函数还是比较复杂的，基本可算得上较低层次的转换。由于我已经习惯了用BCB写C/C++代码，因此我写的浮点数转换字符串函数是80位扩展精度浮点数的，但那个小友拿回去试了一下，说他用的VC不支持80位扩展精度浮点数，虽然能定义long double变量，但实际上还是64位的，我只好又重写了一个64位双精度浮点数的，2个版本使用条件编译，这也算得上是移植性吧，呵呵。

下面是浮点数转换字符串函数的全部代码：

#define USE_EXTENDED /*************************************************************************** * 定义浮点数转换字符串结构 * ***************************************************************************/ typedef struct { SHORT exponent; /* 指数（整数位数） */ BYTE negative; /* 负数（0正，1负）*/ CHAR digits[21]; /* 十进制整数字串 */ }FloatRec; #define F_DEFDECIMALS 6 #define F_MAXDECIMALS 100 #ifdef USE_EXTENDED #define F_MAXPRECISION 19 #define F_CONEXPONENT 0x3fff typedef long double EXTENDED, *PExtended, *PEXTENDED; #include <pshpack2.h> typedef struct { UINT64 mantissa; USHORT exponent; }_Extended; #include <poppack.h> static CONST _Extended _tab0[] = { {0x8000000000000000, 0x3FFF}, /* 10**0 */ {0xA000000000000000, 0x4002}, /* 10**1 */ {0xC800000000000000, 0x4005}, /* 10**2 */ {0xFA00000000000000, 0x4008}, /* 10**3 */ {0x9C40000000000000, 0x400C}, /* 10**4 */ {0xC350000000000000, 0x400F}, /* 10**5 */ {0xF424000000000000, 0x4012}, /* 10**6 */ {0x9896800000000000, 0x4016}, /* 10**7 */ {0xBEBC200000000000, 0x4019}, /* 10**8 */ {0xEE6B280000000000, 0x401C}, /* 10**9 */ {0x9502F90000000000, 0x4020}, /* 10**10 */ {0xBA43B74000000000, 0x4023}, /* 10**11 */ {0xE8D4A51000000000, 0x4026}, /* 10**12 */ {0x9184E72A00000000, 0x402A}, /* 10**13 */ {0xB5E620F480000000, 0x402D}, /* 10**14 */ {0xE35FA931A0000000, 0x4030}, /* 10**15 */ {0x8E1BC9BF04000000, 0x4034}, /* 10**16 */ {0xB1A2BC2EC5000000, 0x4037}, /* 10**17 */ {0xDE0B6B3A76400000, 0x403A}, /* 10**18 */ {0x8AC7230489E80000, 0x403E}, /* 10**19 */ {0xAD78EBC5AC620000, 0x4041}, /* 10**20 */ {0xD8D726B7177A8000, 0x4044}, /* 10**21 */ {0x878678326EAC9000, 0x4048}, /* 10**22 */ {0xA968163F0A57B400, 0x404B}, /* 10**23 */ {0xD3C21BCECCEDA100, 0x404E}, /* 10**24 */ {0x84595161401484A0, 0x4052}, /* 10**25 */ {0xA56FA5B99019A5C8, 0x4055}, /* 10**26 */ {0xCECB8F27F4200F3A, 0x4058}, /* 10**27 */ {0x813F3978F8940984, 0x405C}, /* 10**28 */ {0xA18F07D736B90BE5, 0x405F}, /* 10**29 */ {0xC9F2C9CD04674EDF, 0x4062}, /* 10**30 */ {0xFC6F7C4045812296, 0x4065} /* 10**31 */ }; static CONST _Extended _tab1[] = { {0x9DC5ADA82B70B59E, 0x4069}, /* 10**32 */ {0xC2781F49FFCFA6D5, 0x40D3}, /* 10**64 */ {0xEFB3AB16C59B14A3, 0x413D}, /* 10**96 */ {0x93BA47C980E98CE0, 0x41A8}, /* 10**128 */ {0xB616A12B7FE617AA, 0x4212}, /* 10**160 */ {0xE070F78D3927556B, 0x427C}, /* 10**192 */ {0x8A5296FFE33CC930, 0x42E7}, /* 10**224 */ {0xAA7EEBFB9DF9DE8E, 0x4351}, /* 10**256 */ {0xD226FC195C6A2F8C, 0x43BB}, /* 10**288 */ {0x81842F29F2CCE376, 0x4426}, /* 10**320 */ {0x9FA42700DB900AD2, 0x4490}, /* 10**352 */ {0xC4C5E310AEF8AA17, 0x44FA}, /* 10**384 */ {0xF28A9C07E9B09C59, 0x4564}, /* 10**416 */ {0x957A4AE1EBF7F3D4, 0x45CF}, /* 10**448 */ {0xB83ED8DC0795A262, 0x4639} /* 10**480 */ }; static CONST _Extended _tab2[] = { {0xE319A0AEA60E91C7, 0x46A3}, /* 10**512 */ {0xC976758681750C17, 0x4D48}, /* 10**1024 */ {0xB2B8353B3993A7E4, 0x53ED}, /* 10**1536 */ {0x9E8B3B5DC53D5DE5, 0x5A92}, /* 10**2048 */ {0x8CA554C020A1F0A6, 0x6137}, /* 10**2560 */ {0xF9895D25D88B5A8B, 0x67DB}, /* 10**3072 */ {0xDD5DC8A2BF27F3F8, 0x6E80}, /* 10**3584 */ {0xC46052028A20979B, 0x7525}, /* 10**4096 */ {0xAE3511626ED559F0, 0x7BCA} /* 10**4608 */ }; static CONST EXTENDED _conPrec = 1E19; #else // USE_EXTENDED #define F_MAXPRECISION 17 #define F_CONEXPONENT 0x03ff typedef double EXTENDED, *PExtended, *PEXTENDED; static CONST UINT64 _tab0[] = { {0x3FF0000000000000}, /* 10**0 */ {0x4024000000000000}, /* 10**1 */ {0x4059000000000000}, /* 10**2 */ {0x408F400000000000}, /* 10**3 */ {0x40C3880000000000}, /* 10**4 */ {0x40F86A0000000000}, /* 10**5 */ {0x412E848000000000}, /* 10**6 */ {0x416312D000000000}, /* 10**7 */ {0x4197D78400000000}, /* 10**8 */ {0x41CDCD6500000000}, /* 10**9 */ {0x4202A05F20000000}, /* 10**10 */ {0x42374876E8000000}, /* 10**11 */ {0x426D1A94A2000000}, /* 10**12 */ {0x42A2309CE5400000}, /* 10**13 */ {0x42D6BCC41E900000}, /* 10**14 */ {0x430C6BF526340000}, /* 10**15 */ {0x4341C37937E08000}, /* 10**16 */ {0x4376345785D8A000}, /* 10**17 */ {0x43ABC16D674EC800}, /* 10**18 */ {0x43E158E460913D00}, /* 10**19 */ {0x4415AF1D78B58C40}, /* 10**20 */ {0x444B1AE4D6E2EF50}, /* 10**21 */ {0x4480F0CF064DD592}, /* 10**22 */ {0x44B52D02C7E14AF7}, /* 10**23 */ {0x44EA784379D99DB4}, /* 10**24 */ {0x45208B2A2C280291}, /* 10**25 */ {0x4554ADF4B7320335}, /* 10**26 */ {0x4589D971E4FE8402}, /* 10**27 */ {0x45C027E72F1F1281}, /* 10**28 */ {0x45F431E0FAE6D722}, /* 10**29 */ {0x46293E5939A08CEA}, /* 10**30 */ {0x465F8DEF8808B024} /* 10**31 */ }; static CONST UINT64 _tab1[] = { {0x4693B8B5B5056E17}, /* 10**32 */ {0x4D384F03E93FF9F5}, /* 10**64 */ {0x53DDF67562D8B363}, /* 10**96 */ {0x5A827748F9301D32}, /* 10**128 */ {0x6126C2D4256FFCC3}, /* 10**160 */ {0x67CC0E1EF1A724EB}, /* 10**192 */ {0x6E714A52DFFC679A}, /* 10**224 */ {0x75154FDD7F73BF3C}, /* 10**256 */ {0x7BBA44DF832B8D46}, /* 10**288 */ }; static CONST EXTENDED _conPrec = 1E17; #endif // !USE_EXTENDED static CONST UINT64 _cvttab[F_MAXPRECISION] = { #ifdef USE_EXTENDED 0xDE0B6B3A7640000, 0x16345785D8A0000, #endif 0x02386F26FC10000, 0x0038D7EA4C68000, 0x0005AF3107A4000, 0x00009184E72A000, 0x00000E8D4A51000, 0x00000174876E800, 0x0000002540BE400, 0x00000003B9ACA00, 0x000000005F5E100, 0x000000000989680, 0x0000000000F4240, 0x0000000000186A0, 0x000000000002710, 0x0000000000003E8, 0x000000000000064, 0x00000000000000A, 0x000000000000001 }; #define DECIMAL_EXP(exponent) ((((exponent - F_CONEXPONENT) * 0x4d10) >> 16) + 1) static VOID AdjFloatDigits(UINT64 value, INT precision, INT decimals, FloatRec *rec) { INT i; // value是F_MAXPRECISION位十进制整数，故从最高位开始转换为数字串 for (i = 0; value; i ++) { rec->digits[i] = (CHAR)((value / _cvttab[i]) | 0x30); value %= _cvttab[i]; } memset(rec->digits + i, 0, F_MAXPRECISION - i); // 下面对数字串作精度处理 // 如果总的精度小于0，数字串为空串，该数字转换为'0' if ((i = (rec->exponent + decimals)) < 0) { rec->exponent = rec->negative = rec->digits[0] = 0; return; } if (i > precision) i = precision; // 如果精度位数小于18，同时该位大于等于'5'，四舍五入 if (i < F_MAXPRECISION && rec->digits[i] >= '5') { do { rec->digits[i --] = 0; rec->digits[i] ++; }while (i >= 0 && rec->digits[i] > '9'); if (i < 0) { rec->digits[0] = '1'; rec->exponent ++; } } // 否则，去掉数字串尾部多余的'0' else { if (i > F_MAXPRECISION) i = F_MAXPRECISION; do rec->digits[i --] = 0; while (i >= 0 && rec->digits[i] == '0'); if (i < 0) rec->negative = 0; } } #ifdef USE_EXTENDED // 解析扩展精度浮点数为十进制字符串，并存入浮点数记录中 // 参数：浮点数指针，精度，小数位，浮点数记录 VOID FloatResolve(PEXTENDED pvalue, INT precision, INT decimals, FloatRec *rec) { INT power; EXTENDED val; // 79位：扩展精度浮点数符号位 rec->negative = ((_Extended*)pvalue)->exponent >> 15; // 64-78位：扩展精度浮点数阶码(阶码 - 0x3fff = 二进制指数) rec->exponent = ((_Extended*)pvalue)->exponent & 0x7fff; if (!rec->exponent) // *pvalue = 0 { rec->negative = rec->digits[0] = 0; return; } if (rec->exponent == 0x7fff)// *pvalue = nan or inf { if (!((*(LPBYTE)pvalue + 7) & 0x80) || !(*(PUINT64)pvalue & 0x7fffffffffffffff)) { lstrcpyA(rec->digits, "INF"); } else { rec->exponent ++; rec->negative = 0; lstrcpyA(rec->digits, "NAN"); } return; } // 阶码转换为十进制指数 rec->exponent = DECIMAL_EXP(rec->exponent); // 0-63位：扩展精度浮点数尾数转换成F_MAXPRECISION位十进制浮点整数格式 val = *pvalue; *((LPBYTE)&val + 9) &= 0x7f;// val = fabs(*pvalue) power = F_MAXPRECISION - rec->exponent; if (power > 0) // if (power > 0) val *= (10**power) { val *= *(PEXTENDED)&_tab0[power & 31]; power >>= 5; // power /= 32; if (power) { if (power & 15) val *= *(PEXTENDED)&_tab1[(power & 15) - 1]; power >>= 4; // power /= 16; if (power) val *= *(PEXTENDED)&_tab2[power - 1]; } } else if (power < 0) // if (power < 0) val /= (10**power) { power = -power; val /= *(PEXTENDED)&_tab0[power & 31]; power >>= 5; // power /= 32; if (power) { if (power & 15) val /= *(PEXTENDED)&_tab1[(power & 15) - 1]; power >>= 4; // power /= 16; if (power) val /= *(PEXTENDED)&_tab2[power - 1]; } } val += 0.5; // 四舍五入 if (val >= _conPrec) { val /= 10; rec->exponent ++; } // 调整并转换扩展精度浮点数尾数的整数部分rec->digits AdjFloatDigits(*(PUINT64)&val >> ((((_Extended*)&val)->exponent - 0x3fff) ^ 0x3f), precision, decimals, rec); } #else // USE_EXTENDED // 解析双精度浮点数为十进制字符串，并存入浮点数记录中 // 参数：浮点数指针，精度，小数位，浮点数记录 VOID FloatResolve(PEXTENDED pvalue, INT precision, INT decimals, FloatRec *rec) { INT power; EXTENDED val; // 63位：双精度浮点数符号位 rec->negative = *((LPBYTE)pvalue + 7) >> 7; // 52-62位：双精度浮点数阶码(阶码 - 0x3ff = 二进制指数) rec->exponent = (*(PUINT64)pvalue >> 52) & 0x7ff; if (!rec->exponent) // *pvalue = 0 { rec->negative = rec->digits[0] = 0; return; } if (rec->exponent == 0x7ff)// *pvalue = nan or inf { if ((*(PUINT64)pvalue & 0xfffffffffffff) == 0) { lstrcpyA(rec->digits, "INF"); } else { rec->exponent ++; rec->negative = 0; lstrcpyA(rec->digits, "NAN"); } return; } // 阶码转换为十进制指数 rec->exponent = DECIMAL_EXP(rec->exponent); // 0-51位：双精度浮点数尾数转换成F_MAXPRECISION位十进制浮点整数格式 val = *pvalue; *((LPBYTE)&val + 7) &= 0x7f;// val = fabs(*pvalue) power = F_MAXPRECISION - rec->exponent; if (power > 0) // if (power > 0) val *= (10**power) { val *= *(PEXTENDED)&_tab0[power & 31]; power >>= 5; // power /= 32; if (power) val *= *(PEXTENDED)&_tab1[power - 1]; } else if (power < 0) // if (power < 0) val /= (10**power) { power = -power; val /= *(PEXTENDED)&_tab0[power & 31]; power >>= 5; // power /= 32; if (power) val /= *(PEXTENDED)&_tab1[power - 1]; } // 16位十进制浮点整数四舍五入 val += 0.5; if (val >= _conPrec) { val /= 10; rec->exponent ++; } // 调整并转换扩展精度浮点数尾数的整数部分rec->digits // 清除52-63位，加隐藏的高位，F_MAXPRECISION=17，高位超过52位，所以左移 AdjFloatDigits(((*(PUINT64)&val & 0x000fffffffffffff) | 0x0010000000000000) << -(52 - ((*(PUINT64)&val >> 52) - 0x3ff)), precision, decimals, rec); } #endif // !USE_EXTENDED // 输出指数字符串到buffer，返回指数字符串长度 INT PutExponent(LPSTR buffer, CONST FloatRec *rec) { LPSTR p = buffer; INT e, exp = rec->digits[0]? rec->exponent - 1 : 0; *p ++ = rec->negative & 0x80? 'E' : 'e'; if (exp < 0) { exp = -exp; *p ++ = '-'; } else *p ++ = '+'; if ((e = (exp / 1000)) != 0) { *p ++ = e + 0x30; exp %= 1000; } *p ++ = exp / 100 + 0x30; exp %= 100; *(PUSHORT)p = (((exp % 10) << 8) | (exp / 10)) + 0x3030; return (INT)(p - buffer + 2); } // 浮点数转换为字符串。参数：字符串，浮点数 LPSTR FloatToStr(LPSTR str, EXTENDED value) { INT exp; FloatRec rec; LPSTR pd = rec.digits; LPSTR ps = str; // 解析浮点数，并将信息保存在rec FloatResolve(&value, 15, 9999, &rec); // 打印负数符号 if (rec.negative) *ps ++ = '-'; // NAN or INF if (*pd > '9') { memcpy(ps, pd, 4); return str; } exp = rec.exponent; // 如果十进制指数大于15或者小于-3，转换为指数形式 if (exp > 15 || exp < -3) { *ps ++ = *pd ++; if (*pd) for (*ps ++ = '.'; *pd; *ps ++ = *pd ++); ps += PutExponent(ps, &rec); *ps = 0; return str; } // 否则，转换为小数形式 if (exp <= 0) { *ps ++ = '0'; if (*pd) { for (*ps ++ = '.'; exp < 0; *ps ++ = '0', exp ++); while (*ps ++ = *pd ++); } else *ps = 0; } else { for (; exp > 0 && *pd; *ps ++ = *pd ++, exp --); if (*pd) { *ps ++ = '.'; while (*ps ++ = *pd ++); } else { memset(ps, '0', exp); ps[exp] = 0; } } return str; } 2个版本的代码加起来很长，但还有个自写的springf函数（下篇文章开始介绍）也要用到本文除FloatToSt函数外的全部代码。

代码开头的USE_EXTENDED为编译条件，如果你的编译系统不支持80位扩展精度浮点数，可将该定义注释掉。

前面说了，由于该代码主要是学习用的，因此数据转换层次较低，涉及到的有关浮点数格式的知识，可在网上搜索到。当然，知道是一回事，而具体怎样去操作则又是另一回事了，一些关键地方，我都作了较详细的注释，相信能对初学者正确理解浮点数格式有所帮助。例如，如何在不调用C有关函数，不使用汇编而快速的求一个浮点数的绝对值？本文代码就直接对浮点数进行操作：*((LPBYTE)&val + 7) &= 0x7f;（双精度浮点数）和*((LPBYTE)&val + 9) &= 0x7f;（扩展精度浮点数），也就是直接将浮点数的最高位置零，这当然比什么if (val < 0) val = -val语句快多了。不过，如果你要说后者的可移植性好，那我就无话可说了。要说可移植性，前者也可办到的，修改一下：*((LPBYTE)&val + sizeof(val) - 1) &= 0x7f;不就行了么，除非浮点数格式规则改变。只不过本文代码主要用来学习，用显式的方式更有意义。

80位扩展精度浮点数的有效数字为19位，而64位双精度浮点数有效数字为15 - 16位，本文的解析函数FloatResolve分别用了19位和17位的最大转换精度，尽可能的多显示几位，这主要是为自写的sprintf函数（另文介绍）做准备的，本文的浮点数转换字符串函数FloatToStr只用了最大15位的精度。

下面是个很简单的调用例子：

int main(int argc, char* argv[]) { CHAR s[32]; UINT64 inf = 0x7ff0000000000000; double v = -1.230E45; INT a, b; puts(FloatToStr(s, -10000)); puts(FloatToStr(s, 123456789012345678)); puts(FloatToStr(s, 1234567890.12345678)); puts(FloatToStr(s, 0.0001234567890)); puts(FloatToStr(s, -0.00001234567890)); puts(FloatToStr(s, v)); puts(FloatToStr(s, 0)); puts(FloatToStr(s, *(double*)&inf)); system("pause"); return 0; } 下面是运行结果：

-10000
1.23456789012346e+017
1234567890.12346
0.000123456789
-1.23456789e-005
-1.23e+045
0
INF
请按任意键继续. . .

再次声明：本文代码主要供学习使用，如作其它用途，出问题慨不负责。

水平有限，错误在所难免，欢迎指正和指导。邮箱地址：maozefa@hotmail.com