Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

zhdgzhdg 发表于 2006-4-1 07:36

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程序 的编译问题》中举到了另一个例子。有兴趣可以参考，但“鱼”并不重要，“渔”才重要

浪尖发表于 2006-4-1 08:39

楼主高手

FIEforever 发表于 2009-4-18 01:14

入门者慢慢研究

dbx12358 发表于 2009-9-26 16:17

楼主厉害~~

页: [1]

声振论坛's Archiver

m文件转c/c++ ，“不可能”变“可能”

入门者慢慢研究