发信人: fork (撒哈拉沙漠的沙), 信区: Matlab
标 题: m文件转c/c++ ，“不可能”变“可能”
发信站: BBS 哈工大紫丁香站 (Sun Aug 1 15:42:01 2004)
m文件转c/c++ ，“不可能”变“可能”
当你用“mcc -B sglcpp 文件名”成功将自己的一个程序编译为c语言的程序后是不是
觉得很兴奋呢，以为大功告成了？让我来给你泼一盆冷水吧。试试下面这个程序：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function fork4()
I = imread('lena.tif');
imshow(I)
BW = roipoly;%该函数让你在图像上选取一个多边形
imshow(BW)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
用 mcc -B sglcpp fork4 编译，成功。（虽然编译器给出了很多警告，但可执行程序
最终还是生成了）下面来运行那个自动生成的可执行程序试试......阿哦，不能选取，出
现如下错误：
Exception! File: handler.cpp, Line: 73
Output argument 'a' was not assigned during call to 'xlim'.
此路不通？！
roipoly 是一个经常要用到的函数，如果凡是有它的程序都不能转成VC的话，那就很
让我们这些搞图像处理的为难了。能不能想法使它能够工作呢？试一试：从错误提示中了
解到问题有可能出在xlim函数中，好的我们搜索matlab的安装路径，在“C:\MATLAB6p5toolbox\matlab\graph3d”里找到了它，我把它的关键代码贴出来：
........................................
if nargin == 0
a = get(gca,'xlim');
else
if length(arg1)==1 & ishandle(arg1) & strcmp(get(arg1, 'type'), 'axes')
ax = arg1;
if nargin==2
val = arg2;
else
a = get(ax,'xlim');
return
end
else
if nargin==2
error('Wrong number of arguments')
else
ax = gca;
val = arg1;
end
end
if isstr(val)
if(strcmp(val,'mode'))
a = get(ax,'xlimmode');
else
set(ax,'xlimmode',val);
end
else
set(ax,'xlim',val);
end
end
........................................
错误提示是“a（返回值）没有被赋值”就返回了。那么我们这样做：将xlim.m拷出来
放到我们的工作目录下（这一步很重要，因为下面我们要修改xlim.m程序，如果直接在源
文件中修改将会影响到以后所有基于该函数的程序，放到我们的工作目录中修改则只会影响
到当前的程序），将xlim.m中的一段改为如下：
.........................
if isstr(val)
if(strcmp(val,'mode'))
a = get(ax,'xlimmode');
else
set(ax,'xlimmode',val);
disp('here a')
end
else
set(ax,'xlim',val);
disp('here b')
end
.........................
这么做的目的是要探察a是在哪里不被赋值就返回了。在matlab中运行fork4发现命令
行中只显示'here a'，这说明在fork4的运行过程中流程经过了该位置，但a未被赋值。为了
使编译后的程序能够运行，我在该处手工的给a赋值，于是该段代码变为：
.........................
if isstr(val)
if(strcmp(val,'mode'))
a = get(ax,'xlimmode');
else
set(ax,'xlimmode',val);
a=1;
end
else
set(ax,'xlim',val);
end
.........................
再编译，运行你会发现xlim.m中没有问题了，但ylim.m中又有问题，于是仿照前法，
将ylim.m拷出来（一定记得不要在原位置对文件进行修改）相应代码修改：
.........................
if isstr(val)
if(strcmp(val,'mode'))
a = get(ax,'ylimmode');
else
set(ax,'ylimmode',val);
a = 1;
end
else
set(ax,'ylim',val);
end
.........................
再编译，运行，OK！！！
好了以后凡是涉及 roipoly 函数都可以这样来做。
最后我想说明一下这样做为什么可行。
“a没有被赋值就被返回”正说明a没有必要在该情况下进行赋值，所以它干脆不赋值
就返回了，这在matlab中运行没有问题，但转成C语言后可能因为一些特殊的原因，要求函
数必须在函数体内对返回值进行赋值，因此就发生了运行时错误。这个错误并非是逻辑或数
值上的错误，而只是违反了一些“形式上的规范”，所以我们人为的给a赋一个值遵守了这个规范就能顺利通过了。那该赋什么值呢？正如前面所说，既然a根本就没有必要进行赋值，那么给它赋什么值都无所谓了。所以随便给它赋什么值都行。
很多原本看似“Mission Imposible”的“m文件转c/c++文件，VC编译”的任务其实都
可以通过修改matlab自身的m文件来实现，在另一篇文章《分析：解决m程序转化成cpp程序
的编译问题》中举到了另一个例子。有兴趣可以参考，但“鱼”并不重要，“渔”才重要
。
