C++调用Matlab函数求特征值

 更新时间:2019年06月17日 11:19:37   作者:xyzZZ_  
这篇文章主要为大家详细介绍了C++调用Matlab函数求特征值,具有一定的参考价值,感兴趣的小伙伴们可以参考一下

最近需要用到C++和Matlab的混编,记录一下学习过程。

要实现的是调用Matlab函数,求矩阵前k个最小的特征值及其特征向量。

//C++
 
#include "engine.h" //使用Matlab引擎需要包含的头文件
#include <iostream>
using namespace std;
int main()
{
 Engine *m_engine; //创建Matlab引擎
 m_engine = NULL; //初始化引擎
 
 if((!m_engine && !(m_engine = engOpen(NULL)))) //打开引擎,此时会打开一个Matlab命令行窗口
 {
  return -1;
 }
 engSetVisible(m_engine,1); //将命令行窗口设为可见
 
 char buffer[255]; //记录调试信息,方便调试
 engOutputBuffer(m_engine, buffer, 255); 
 
 
 double A[3][3] = {-1,1,0,
   -4,3,0,
    1,0,2}; 
 mxArray* AObj = mxCreateDoubleMatrix(3, 3, mxREAL); //创建Matlab的矩阵(大小3*3,实数)
 
 memcpy(mxGetPr(AObj), A, 3*sizeof(double)); //将C++的数据传入Matlab中
 
 engPutVariable(m_engine, "A", AObj); //将AObj的值赋给A
 
  engEvalString(m_engine, "cd('E:\\MatlabScripts')"); //进入Matlab代码的路径 
 
  //调用Matlab中定义的函数“computeEigens”,文件名需与函数名一致,即“computeEigens.m”
  int k = 2;
 engEvalString(m_engine, "[eigVector,eigValue] = computeEigens(A, k);"); 
 
  //存储计算结果
 engEvalString(m_engine,"save('E:\\eigVec_eigV.mat','eigVector','eigValue');");
 
 printf("%s", buffer);
 
 mxDestroyArray(AObj); //销毁Matlab数组
 
 if (m_engine) //关闭Matlab引擎
 {
 engClose(m_engine);
 m_engine = NULL;
 }
 
 return 0;
 
}
%调用的Matlab代码
%E:\MatlabScripts\computeEigens.m
 
function [ eigVector, eigValue ] = computeEigens( M, n_Eigens ) %n_Eigens为要求的特征值数量
 [EigenVectors,EigenValues] = eig(M);
 [sortedEigenValues, index] = sort(diag(EigenValues));
 eigValue = sortedEigenValues(1:n_Eigens);
 idx = index(1:n_Eigens);
 eigVector = EigenVectors(:,idx);
 
end

求稀疏矩阵的特征值

后来发现eig不能用于求解稀疏矩阵,会报“Error using eig”的错,要改用eigs(A, k, sigma),sigma='sm'时表示求稀疏矩阵A的前k个绝对值最小的特征值及其特征向量。sigma的其他取值含义为:'lm' 绝对值最大的特征值;'sm' 绝对值最小的特征值;'la'最的大特征值;'sa'最小的特征值;'lr' 最大实部;'sr' 最小实部;'li' 最大虚部;'si'最小虚部。

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持脚本之家。

相关文章

  • Qt读取和写入配置(ini)文件

    Qt读取和写入配置(ini)文件

    ini文件在windows系统中可以存储需要持久保存的配置信息,本文主要介绍了Qt读取和写入配置(ini)文件,具有一定的参考价值,感兴趣的可以了解一下
    2023-09-09
  • Linux下实现C++操作Mysql数据库

    Linux下实现C++操作Mysql数据库

    由于工作需要抽出一周的时间来研究C/C++访问各种数据库的方法,并打算封装一套数据库操作类,现在奉上最简单的一部分:在Linux下访问MySQL数据库。
    2017-05-05
  • 在C++中加载TorchScript模型的方法

    在C++中加载TorchScript模型的方法

    这篇文章主要介绍了在C++中加载TorchScript模型的方法,本文通过实例代码给大家介绍的非常详细,具有一定的参考借鉴价值,需要的朋友可以参考下
    2020-02-02
  • C++ string替换单个指定字符为其它字符问题

    C++ string替换单个指定字符为其它字符问题

    这篇文章主要介绍了C++ string替换单个指定字符为其它字符问题,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2023-06-06
  • Qt实战案例之如何利用QProcess类实现启动进程

    Qt实战案例之如何利用QProcess类实现启动进程

    这篇文章主要介绍了Qt实战案例之如何利用QProcess类实现启动进程,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下
    2022-02-02
  • cocos2d-x学习笔记之CCLayer、CCLayerColor、CCLayerGradient、CCLayerMultiplex场景层介绍

    cocos2d-x学习笔记之CCLayer、CCLayerColor、CCLayerGradient、CCLayerMu

    这篇文章主要介绍了cocos2d-x学习笔记之CCLayer、CCLayerColor、CCLayerGradient、CCLayerMultiplex场景层介绍,需要的朋友可以参考下
    2014-09-09
  • C++中rapidjson将嵌套map转为嵌套json的讲解

    C++中rapidjson将嵌套map转为嵌套json的讲解

    今天小编就为大家分享一篇关于C++中rapidjson将嵌套map转为嵌套json的讲解,小编觉得内容挺不错的,现在分享给大家,具有很好的参考价值,需要的朋友一起跟随小编来看看吧
    2019-04-04
  • 简单讲解c++ vector

    简单讲解c++ vector

    这篇文章主要介绍了c++ vector的相关资料,帮助大家更好的理解和学习c++,感兴趣的朋友可以了解下
    2020-09-09
  • 八皇后问题的相关C++代码解答示例

    八皇后问题的相关C++代码解答示例

    这篇文章主要介绍了八皇后问题的相关C++代码解答示例,文中包括ACM竞赛的八皇后相关知识的练习实例,需要的朋友可以参考下
    2015-08-08
  • 关于C++的.cpp文件运行全过程

    关于C++的.cpp文件运行全过程

    这篇文章主要介绍了C++的.cpp文件运行全过程,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2023-02-02

最新评论