最近需要用到C++和Matlab的混編，記錄一下學(xué)習(xí)過程。
因為努力和真誠，有更多的客戶和我們聚集在一起，為了共同目標(biāo)，成都創(chuàng)新互聯(lián)在工作上密切配合，從創(chuàng)業(yè)型企業(yè)到如今不斷成長，要感謝客戶對我們的高要求，讓我們敢于面對挑戰(zhàn)，才有今天的進步與發(fā)展。從網(wǎng)站到重慶小程序開發(fā),軟件開發(fā),成都App定制開發(fā),十余年企業(yè)網(wǎng)站建設(shè)服務(wù)經(jīng)驗,為企業(yè)提供網(wǎng)站設(shè)計,網(wǎng)站托管運營一條龍服務(wù).為企業(yè)提供營銷型網(wǎng)站,專業(yè)公司,原創(chuàng)設(shè)計,十余年品質(zhì),值得您的信賴.
要實現(xiàn)的是調(diào)用Matlab函數(shù)，求矩陣前k個最小的特征值及其特征向量。
//C++ #include "engine.h" //使用Matlab引擎需要包含的頭文件 #includeusing 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); //將命令行窗口設(shè)為可見 char buffer[255]; //記錄調(diào)試信息，方便調(diào)試 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代碼的路徑 //調(diào)用Matlab中定義的函數(shù)“computeEigens”，文件名需與函數(shù)名一致，即“computeEigens.m” int k = 2; engEvalString(m_engine, "[eigVector,eigValue] = computeEigens(A, k);"); //存儲計算結(jié)果 engEvalString(m_engine,"save('E:\\eigVec_eigV.mat','eigVector','eigValue');"); printf("%s", buffer); mxDestroyArray(AObj); //銷毀Matlab數(shù)組 if (m_engine) //關(guān)閉Matlab引擎 { engClose(m_engine); m_engine = NULL; } return 0; }
%調(diào)用的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'最小虛部。
以上就是本文的全部內(nèi)容，希望對大家的學(xué)習(xí)有所幫助，也希望大家多多支持創(chuàng)新互聯(lián)。
重慶分公司
重慶分公司
