欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

C++調用Matlab函數(shù)求特征值

 更新時間:2019年06月17日 11:19:37   作者:xyzZZ_  
這篇文章主要為大家詳細介紹了C++調用Matlab函數(shù)求特征值,具有一定的參考價值,感興趣的小伙伴們可以參考一下

最近需要用到C++和Matlab的混編,記錄一下學習過程。

要實現(xiàn)的是調用Matlab函數(shù),求矩陣前k個最小的特征值及其特征向量。

//C++
 
#include "engine.h" //使用Matlab引擎需要包含的頭文件
#include <iostream>
using namespace std;
int main()
{
 Engine *m_engine; //創(chuàng)建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); //創(chuàng)建Matlab的矩陣(大小3*3,實數(shù))
 
 memcpy(mxGetPr(AObj), A, 3*sizeof(double)); //將C++的數(shù)據(jù)傳入Matlab中
 
 engPutVariable(m_engine, "A", AObj); //將AObj的值賦給A
 
  engEvalString(m_engine, "cd('E:\\MatlabScripts')"); //進入Matlab代碼的路徑 
 
  //調用Matlab中定義的函數(shù)“computeEigens”,文件名需與函數(shù)名一致,即“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數(shù)組
 
 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為要求的特征值數(shù)量
 [EigenVectors,EigenValues] = eig(M);
 [sortedEigenValues, index] = sort(diag(EigenValues));
 eigValue = sortedEigenValues(1:n_Eigens);
 idx = index(1:n_Eigens);
 eigVector = EigenVectors(:,idx);
 
end

求稀疏矩陣的特征值

后來發(fā)現(xiàn)eig不能用于求解稀疏矩陣,會報“Error using eig”的錯,要改用eigs(A, k, sigma),sigma='sm'時表示求稀疏矩陣A的前k個絕對值最小的特征值及其特征向量。sigma的其他取值含義為:'lm' 絕對值最大的特征值;'sm' 絕對值最小的特征值;'la'最的大特征值;'sa'最小的特征值;'lr' 最大實部;'sr' 最小實部;'li' 最大虛部;'si'最小虛部。

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。

相關文章

  • Protocol Buffer技術深入理解(C++實例)

    Protocol Buffer技術深入理解(C++實例)

    C++實例Protocol Buffer技術詳解,感興趣的朋友可以了解下
    2013-01-01
  • C/C++語言中結構體的內存分配小例子

    C/C++語言中結構體的內存分配小例子

    當未用 #pragma 指令指定編譯器的對齊位數(shù)時,結構體按最長寬度的數(shù)據(jù)成員的寬度對齊;當使用了 #pragma 指令指定編譯器的對齊位數(shù)時,結構體按最長寬度的數(shù)據(jù)成員的寬度和 #pragma 指令指定的位數(shù)中的較小值對齊
    2013-10-10
  • C++實現(xiàn)的求解多元一次方程示例

    C++實現(xiàn)的求解多元一次方程示例

    這篇文章主要介紹了C++實現(xiàn)的求解多元一次方程,涉及C++矩陣運算相關操作技巧,需要的朋友可以參考下
    2018-01-01
  • C++中const用法小結

    C++中const用法小結

    C++ const 允許指定一個語義約束,編譯器會強制實施這個約束,允許程序員告訴編譯器某值是保持不變的。如果在編程中確實有某個值保持不變,就應該明確使用const,這樣可以獲得編譯器的幫助。
    2016-04-04
  • C語言中`||`的短路機制詳解

    C語言中`||`的短路機制詳解

    在C語言中,邏輯或運算符(||)是一種常用的邏輯運算符,用于組合多個條件表達式,C語言中的邏輯或運算符具有短路機制,這是一種非常重要的概念,本文將深入解釋C語言中的||短路機制以及其在編程中的應用,感興趣的朋友跟隨小編一起看看吧
    2024-01-01
  • C語言超全面覆蓋操作符知識點

    C語言超全面覆蓋操作符知識點

    C?語言提供了豐富的操作符,有:算術操作符,移位操作符,位操作符,賦值操作符,單目操作符,關系操作符,邏輯操作符,條件操作符等。讓我們通讀本篇來詳細了解吧
    2022-06-06
  • c/c++ 奇技淫巧(一些c語言的技巧)

    c/c++ 奇技淫巧(一些c語言的技巧)

    這篇文章主要介紹了c/c++ 奇技淫巧,需要的朋友可以參考下
    2017-03-03
  • C++實現(xiàn)浮點數(shù)精確加法

    C++實現(xiàn)浮點數(shù)精確加法

    這篇文章主要為大家詳細介紹了C++實現(xiàn)浮點數(shù)精確加法,文中示例代碼介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下
    2020-05-05
  • Qt在線安裝加速的實現(xiàn)

    Qt在線安裝加速的實現(xiàn)

    本文主要介紹了Qt在線安裝加速的實現(xiàn),文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧
    2023-02-02
  • C++中的new/delete、構造/析構函數(shù)、dynamic_cast分析

    C++中的new/delete、構造/析構函數(shù)、dynamic_cast分析

    這篇文章主要介紹了C++中的new/delete、構造/析構函數(shù)、dynamic_cast分析 本文通過實例代碼給大家介紹的非常詳細,具有一定的參考借鑒價值,需要的朋友可以參考下
    2019-05-05

最新評論