- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Matlab与Visual c 溷合编程
Matlab 与Visual c++接口编程
以广义特征值为例
1,Matlab 程序
function [v,d]=qrfact(a,b)
[v,d]=eig(a,b);
end
2,Matlab 程序编译
deploytool
创建一个工程qrfact.prj,
将qrfat.m 加入qrfact.prj
编译
编译完成后在该项目目录的qrfact\src 下会出现qrfact.dll,qrfact.h 和qrfact.lib 等文件。
3,Vc 设置
首先在Visual c++下建立一个MFC 项目cmatlab (笔者利用Visual Studio2008 )
将qrfact.dll 拷贝至c:\windows\system32 下,将qrfact.h 和qrfact.lib 拷贝至cmatlab 项目目
录下。
在打开 Visual Studio2008 ,在“项目-选项-项目和解决方案”中,选择“包含文件”中加
入matlab 的include 路径,笔者为“C:\ Program Files\ MATLAB\ R2009a\extern\include”;在“项
目-选项-项目和解决方案”中,选择“库文件”中加入matlab 静态库的路径“C:\ Program
Files\ MATLAB\ R2009a\extern\lib\win32\lcc”。可根据自己系统的实际情况选择,笔者为32 位
系统,故选择win32,如果为64 位系统则选win64.
在“项目-属性”按下图设置,在附加依赖项中加入mclmcrrt.lib,mclmcr.lib 和qrfact.lib 三个
静态库。
4 ,源代码实现
4.1 打开项目文件夹中的 targetver.h 头文件,把其中的 0X0600 修改为 0X0501,即:
#define WINVER 0x0600 === #define WINVER 0x0501
#define _WIN32_WINNT 0x0600 === #define _WIN32_WINNT 0x0501
4.2 在*.cpp 文件中加入头文件#include “qrfact.h”
void CcmatlabDoc ::OnQrfactorization ()
{
// TODO: 在此添加命令处理程序代码
if (pmesh.VertexNumber=0)
return ;
if(!qrfactInitialize ())
return ;
//计算图的临接矩阵,可根据实际情况,自己编写,也可以指定广义特征值的两个矩阵a,b;
int vtxnumber=pmesh.VertexNumber ;
double *pneighbor ;
pneighbor=new double[vtxnumber*vtxnumber];
memset (pneighbor,0,sizeof(double)*vtxnumber*vtxnumber);
int i,j,vid0,vid1;
for (i=0;ipmesh.EdgeNumber ;i++)
{
vid0=pmesh.pedgeset [i].vertexid0;
vid1=pmesh.pedgeset [i].vertexid1;
pneighbor[vid0*vtxnumber+vid1]=-1.0;
pneighbor[vid1*vtxnumber+vid0]=-1.0;
}
double sum;
for (i=0;ivtxnumber ;i++)
{
sum=0.0;
for (j=0;jvtxnumber ;j++)
{
sum+=pneighbor [i*vtxnumber+j];
}
pneighbor[i*vtxnumber+i]=-sum;
}
double *parea ;
parea=new double[vtxnumber*vtxnumber];
memset (pa
文档评论(0)