|
林 岚* 周云波 杨玉星 林家瑞
9 X+ _6 n2 V# ~( p$ E
/ v* M8 g: B. t摘 要 本文首先介绍了当前在国内外十分流行的工程设计和系统仿真软件包Matlab的 2 Z% [0 ] O( B# U H! {
优缺点,然后展开介绍了Matcom(一种将Matlab与C++结合的工具)的功能及如何使用,最 6 u3 f6 F8 t' G' f
后对其效果进行了总结。
5 b+ i! t6 t& D; Q关键词 Matlab Visual C++ Matcom ECG
5 b* Y0 I3 n T% [, |3 c1 工程设计软件Matlab特点分析 2 A7 p! }- k% |% M( l3 ]
Matlab是MathWorks公司于1982年推出的一套高性能的数值计算和可视化软件。它集数值 $ i- _: D: {+ i- Q# [. l0 C
分析、矩阵运算、信号处理和图形显示于一体,构成了一个方便、界面友好的用户环境
( M' n( m/ r* s' i。它还包括了ToolBox(工具箱)的各类问题的求解工具,可用来求解特定学科的问题。其
- ?0 k) p v( N" M- i- }! e6 g特点是: 8 Q {+ h( r! w3 w9 d W
(1)可扩展性:Matlab最重要的特点是易于扩展,它允许用户自行建立指定功能的M文件
. ]5 p# C* F# U- P, `。对于一个从事特定领域的工程师来说,不仅可利用Matlab所提供的函数及基本工具箱
- c9 ]0 x$ g! M# U9 s8 J函数,还可方便地构造出专用的函数,从而大大扩展了其应用范围。当前支持Matlab的 % |( g; s+ m; p' D |# F6 g
商用Toolbox(工具箱)有数百种之多。而由个人开发的Toolbox则不可计数。 8 ~5 O0 }& n- @2 q% A) W8 O
(2)易学易用性:Matlab不需要用户有高深的数学知识和程序设计能力,不需要用户深刻 - a% s0 L! ^) x9 p! C0 _9 Y
了解算法及编程技巧。 ! r7 _* A& O. V9 D4 K
(3)高效性:Matlab语句功能十分强大,一条语句可完成十分复杂的任务。如fft语句可
6 i- X0 ]5 Y9 z: }3 }) L7 w$ R1 P% ]完成对指定数据的快速傅立叶变换,这相当于上百条C语言语句的功能。它大大加快了工 ' J$ Z" r5 n9 M, |0 X& k" V
程技术人员从事软件开发的效率。据MathWorks公司声称,Matlab软件中所包含的Matla T: O9 Q2 s, r+ I7 j
b源代码相当于70万行C代码。
, n" \9 y! o/ a- W; f1 `- \0 ^由于Matlab具有如此之多的特点,在欧美高等院校,Matlab已成为应用于线性代数、自 3 d" {5 e& j+ z
动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基 5 o6 F3 p. z( E4 @2 A* o' D/ h
本教学工具;在研究单位、工业部门,Matlab也被广泛用于研究和解决各种工程问题。 # b6 p+ s _) }8 C& v/ _9 s
当前在全世界有超过40万工程师和科学家使用它来分析和解决问题。
9 A8 } J3 D4 r" c4 o* V然而Matlab自身所存在的某些缺点限制了它的应用范围。
% I. e+ F4 G& h% y( h' D, u (1)Matlab是一种解释性语言,因此它的实时效率是相当差的。例如在笔者所从事的 6 V1 N* K! L5 s6 ?' d' c* X4 [
一个心电信号处理项目中,采用Matlab语言来处理5分钟的心电数据,在486/66系统机上 }+ }# I9 T% u* z/ M; s- ^, z1 J
需要25分钟,而同样的算法,用C语言实现,仅需30多秒。由此可见,Matlab的语言执行
* [( {6 }- v% t! I效率是比较低的。 . M& C8 y X1 X+ W5 f
(2)Matlab程序不能脱离其环境运行,因此它不能被用于开发商用软件。 (3)
& _9 I+ R( C1 g- T) m% KMatlab程序可以被直接看到 ?: f. d( o1 N2 Y
2 Matcom:Matlab与C++的完美结合 ; k* m @* V) d$ B' m( c8 }; P
MathTools公司推出的Matcom(可于http://www.mathtools.com处获取)可以说是以上
+ C8 f+ d- t+ p, q问题的完美解决。它可将Matlab的源代码译成同等功能的C++源码。既保持了Matlab的优
( I# t, j+ ^* P! ^+ Q( A) n4 H良算法,又保持了C++的高执行效率。
& P& \4 L1 x* _ 下面将详细介绍如何使Matlab与C++接口。 7 R2 `; b, L5 U/ d
运行环境:Matcom 3.0 Matlab 5.1 Visual C++ 5.0 中文Windows 95 PC:4 6 L0 ~7 Q" t, r# x* c2 g7 `% I: d8 ^+ T
86/66 (1)首先应设置Matcom所对应的编译器类型、数据类型。 (2)将你所需要 , e" s3 f" o, e* S
的.M文件翻译成C源码。 在这里将用一个简单的滤波器程
/ w( h' U4 E& W7 d; _) ?+ ~Matlab源程序:
/ g. K9 o1 r6 r Lin.m 6 c, X! p1 b u9 e& j
b=[1 2 1];
8 _ I R4 A# S& U: M* l; {a=[1 -1.58 0.64];n=50;
, V6 [3 L* M/ H( Xx=[1, ones(1,n-1)]; 7 S" \' r% Y# [' a+ d( A; m9 i. c" U
y=filter(b,a,x);
0 ~! r/ q7 p# Y! s/ ~1 G* Wplot(y./70,'*'); ; p' s8 w1 p( l
hold on ( x+ |) V5 S& C" M) ^! P
plot(x,'.');
+ K; u s- Y* a' [' m6 Wgrid; + i4 @+ e1 D7 i! f
下面给出了经Matcom翻译后的C++源程序。
/ i) i5 X4 h' d- y: HLin.cpp % M( F- C. }6 p) b# P) h9 @7 X0 M0 [
M b("b",0,0); 4 |7 q$ y* y7 J; C; f& J" a
M a("a",0,0); 8 w" F2 W `7 n7 s! e4 A
M n("n",0,0);
1 H0 r* R# }. d: U3 qM x("x",0,0);
* q. v8 [. y5 NM y("y",0,0); ) N0 m8 |. H6 F+ h( J8 k3 c- f
b=(BR(1.0),2.0,1.0);
" |% u6 U, d7 X) M* u* \% Ra=(BR(1.0),-1.58,0.64); 3 X, o: y& v& T7 q* j# L5 r
n=50.0;
4 P! h& ^5 x: kx=(BR(1.0),ones(1.0,n-1.0)); ; `& K6 v0 `( z1 `. a. u" B' g
y=filter(b,a,x);
# g5 ~0 j" h3 j: H# ?plot(rdivide(y,70.0),TM("*"));
" w! ~& z7 i5 i) ihold(TM("on")); # W5 z6 |4 g# c _; @& M5 m1 g
plot(x,TM(".")); ; S* V: R3 Z& C4 Z/ S% e, k
grid();
0 [6 R0 S9 N$ l, f4 d( { g(3)Visual C++环境设置。 " r4 V' h1 @8 r" x0 X) M
在Visual C++中的"Project / Settings / C, C++ / Preprocessor /Additional incl ) g' ?4 c3 ~% {) N0 A. B) K
ude subdirectories"中填入当前Matcom所处的路径.。(如"c:\matcom\")在Visual C
. [- R9 S+ _ }7 |( O: i2 u++中的工程窗口处" "Project / Add to Project / Files".加入mlibvc0.lib库文件。 $ [- T0 S: ^: J$ E( e
3 X, u/ K4 ~3 R; j8 m/ y/ E
(4)修改C++源代码。 & ^( }( I6 P0 U
注意:
& c5 e$ ?6 ?# K( n①在头文件处应包含matlib.h. #define MLIB 0 #include "matlib.h." ②在
. t$ y& O* d9 {; p Q7 S9 \* n类成员函数CMfcApp::InitInstance中调用. InitM启动对Matlab函数的调用。 BOO
8 J9 ?/ D3 e! Y7 C) f6 C" e5 b1 RL CEcgstarApp::InitInstance() { TRACE0("initM\n"); initM(s
5 X% R, v( u1 o+ j, R. Wtin,stout,sterr); ...... } ③在类构造函数CMfcApp::CMfcApp.中不要 & s( T1 _3 P2 {; R. z% ?
调用任何Matlab函数 ; A- O1 |# Y* T$ k1 ~
④在类构造函数CMfcApp::~CmfcApp中调用exitM结束对Matlab函数的调用。 CEcg
3 [, O4 N' H) l, SstarApp::~CEcgstarApp() { exitM(); } 1 A9 x; o- J% U- e5 n2 `9 I
⑤在程序中合适的位置插入经Matcom翻译过后的C++源码。 + L2 D. l# A& E! o; ~
⑥编译工程文件,就可得到你所需的结果。
. ]$ C/ H ~$ g% E/ L2 i3 结束语 0 M4 u* y0 y! V3 L6 q5 ]# r. m
笔者在从事生物医学工程ECG信号的检测与处理项目中,对采用不同方法完成同一任务在
; I0 T- m2 b( y0 z" j8 A% v& T8 s u, f. l$ Z
代码量和运行时间上进行了比较,如表1、表2所列。 表1 ECG检测核心算法代码量 全 % Z ^& L8 U% D. z: K- }4 V- u+ o
部用C++编写 全部用Matlab编写 用Matcom转换后的C编写 约3400行 约400行 约600行 * ]: I0 K; n$ q# Q
表2 ECG检测核心算法运行时间 全部用C++编写 全部用Matlab编写 用Matcom转换后的 * j+ e4 x0 B5 U' [* k+ S* D
C编写 约30秒 约25分钟 约1分钟 由此可见,采用Matcom进行程序设计,这是一种 # r! K/ _* u/ C. Z! }
时间、效率上的极好的折衷,它使得你可以将更多的时间放在项目关键问题的思考上, 0 d: } n' @7 h4 v
而不是放在算法的实现和程序的编码中,对于广大工程技术人员的开发工作是十分有益 2 U% m8 s, u4 J- g% K1 f
的。
) {& p2 e6 P0 N( |' k作者单位:华中理工大学生物工程系(武汉 430074)
# K$ T6 U5 F4 |0 E, u: U*男,24岁,硕士研究生 |