matlab2c使用c++实现matlab函数系列教程-inv函数

全栈工程师开发手册 （作者：栾鹏）

matlab2c动态链接库下载
matlab库函数大全
matlab2c基础教程
matlab2c开发全解教程

matlab2c调用方法：

1、下载动态链接库
2、将Matlab2c.dll拷贝到exe同目录下
3、将Matlab2c.h、Matlab2c.lib放到项目头文件目录下
4、在cpp文件中引入下面的代码

#include "Matlab2c.h"#pragma comment(lib,"Matlab2c.lib")  using namespace Matlab2c;

---

matlab中inv函数简介

1、inv函数：求矩阵的逆
2、用法说明

 A=inv(B)，其中B是输入的可逆矩阵，输出A就是B的逆矩阵，逆矩阵满足性质  AB=BA=E (E是单位阵)。如果输入的是不可逆矩阵会弹出警告，并返回inf

inv的c++源码实现

求矩阵的逆
输入：为可逆方矩阵 输出：矩阵的逆。不可逆时返回null

void swap(double *a,double *b);  //声明子程序Matrix Matlab2c::inv(Matrix& a){    if (a.row!=a.column)    {        throw_logic_error("求逆矩阵不为方阵");        return NULL;    }    int *is,*js,i,j,k;    int n=a.row;    double temp,fmax;    Matrix p(a.row,a.column);    for(i=0;i<n;i++)    {         for(j=0;j<n;j++)            p(i,j)=a(i,j);    }    is=new int[n];    js=new int[n];    for(k=0;k<n;k++)    {        fmax=0.0;        for(i=k;i<n;i++)            for(j=k;j<n;j++)            {                 temp=std::fabs(p(i,j));//找最大值                if(temp>fmax)                {                     fmax=temp;                    is[k]=i;js[k]=j;                }            }            if((fmax+1.0)==1.0)            {                  delete[] is;                delete[] js;                throw_logic_error("不存在逆矩阵");                return NULL;            }            if((i=is[k])!=k)                for(j=0;j<n;j++)                    swap(p.data+k*n+j,p.data+i*n+j);//交换指针            if((j=js[k])!=k)                for(i=0;i<n;i++)                    swap(p.data+i*n+k,p.data+i*n+j);  //交换指针            p(k,k)=1.0/p(k,k);            for(j=0;j<n;j++)                if(j!=k)                    p(k,j)*=p(k,k);            for(i=0;i<n;i++)                if(i!=k)                    for(j=0;j<n;j++)                        if(j!=k)                            p(i,j)=p(i,j)-p(i,k)*p(k,j);            for(i=0;i<n;i++)                if(i!=k)                    p(i,k)*=-p(k,k);    }    for(k=n-1;k>=0;k--)    {        if((j=js[k])!=k)            for(i=0;i<n;i++)                swap((p.data+j*n+i),(p.data+k*n+i));        if((i=is[k])!=k)            for(j=0;j<n;j++)                swap((p.data+j*n+i),(p.data+j*n+k));    }    delete[] is;    delete[] js;    return p;}void swap(double *a,double *b){      double c;    c=*a;    *a=*b;    *b=c;}

inv函数的使用测试

#include "Matlab2c.h"#pragma comment(lib,"Matlab2c.lib")  using namespace Matlab2c;int main(){    double a[]={1,2,3,2,2,1,3,4,3};    Matrix aa(3,3,a);    Matrix bb = Matlab2c::inv(aa);    cout<<bb.toString()<<endl;    system("pause");    return 0;}