式中　 η ——风机全压效率
　　　 N e ——风机有效功率， kW
　　　 N ——风机轴功率， kW
　　　 Q ——风机体积流量， m³/h
　　　 p ——风机全压， Pa 。

3.2 数据换算

　　标准状态是指通风机进口处空气的压力为 101.325kPa(1atm 或 760mmHg) ，温度为 20℃，相对湿度为 50 ％的气体状况 [3] 。风机性能测试试验的大气条件是随地点和时间的改变而变化的，而为了使性能曲线具有可比性，风机的性能曲线都必须在标准状态下给定。所以须先将试验数据换算成为标准状态下的数据。

　　改变大气压力 p_b 以及温度 t 时，根据相似理论 [4] ，换算公式为

3.3 曲线拟合

　　本程序中，曲线拟合的算法采用最小二乘法。对于给定数据 ( x_i , y_i ) ， i =0 ， 1… n ，有函数

　　在 MATLAB 中，可通过调用库函数 ployfit() 来实现最小二乘拟合，得到拟合系数 a_i ，再代入式 (7) ，即得拟合函数 p ( x ) ，据此绘制曲线，工作量大为减少。

3.4 程序设计

　　通过调用 MATLAB 中提供的 guide 命令 [5] ，可设计出整个程序的界面，如图 1 所示。具体的编程过程可参考文献 [5，6]。

4、实例

　　以一台离心风机试验为例，应用所编制的程序进行试验数据处理，绘制其性能曲线。试验时的风机进口处大气压力为99860Pa ，温度为33.2℃ 。试验数据如表 1 所示。

表 1 试验数据

　　从通风理论 [3] 知，风压、轴功率以及效率与风机流量的平方成正比。由此，采用最小二乘法拟合得到 3 个二次多项式来描述这 3 条特性曲线。拟合得到的性能曲线如图2～图4所示。

　　从图 2～图4可以看出，试验数据与拟合曲线吻合得很好，其精度也非常高。由程序得，功率的误差平方和 S_N =0.049％，效率的误差平方和 S_η =0.017%。

５、结论

　　利用 MATLAB 软件，用最小二乘法处理风机性能实测数据，图形处理功能强大，程序界面友好，数据处理精度高。本软件已应用于工程实际中，受到用户好评。本方法也可以用于处理泵的试验数据。

参 考 文 献

[1] 吴正人，程友良. 用 VB 语言实现泵或风机性能曲线的绘制和换算.电力情报，2002(2) .
[2] 邓建中，刘之行. 计算方法 ( 第 2 版 ) . 西安交通大学出版社， 2001 .
[3] 吕玉民，陈人宴.中国供暖通风空调设备手册：第二分册.暖通设备与空调风机.机械工业出版社，1998 .
[4] 吴达人.泵与风机西安交通大学出版社， 1989 .
[5] 薛定宇，陈阳泉. 基于 MATLAB/Simulink 的系统仿真技术与应用 . 清华大学出版社，2002 .
[6] 张志涌，徐彦琴等. MATLAB 教程－基于 6.x 版本 . 北京航空航天大学出版社 2001 .