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

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

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中rank函数简介

1、rank函数：求矩阵的秩

2、用法说明

  A=[2,3,4;6,7,9;0,0,1];  rank(A);  求矩阵的秩有很多种方法，Matlab 的 rank 函数是借助矩阵奇异值分解的方法计算的

rank的c++源码实现

求矩阵的秩
输入源矩阵，返回矩阵的秩int型

int Matlab2c::rank(Matrix& a){    int i,j,k,nn,is,js,l,ll,u,v;    double q,d;    nn=a.row;    if (a.row>=a.column) nn=a.column;    k=0;    for (l=0; l<=nn-1; l++)    {         q=0.0;        for (i=l; i<=a.row-1; i++)            for (j=l; j<=a.column-1; j++)            {                 ll=i*a.column+j; d=std::fabs(a.data[ll]);                if (d>q)                 {                     q=d; is=i; js=j;                }            }            if (q+1.0==1.0) return(k);            k=k+1;            if (is!=l)            {                 for (j=l; j<=a.column-1; j++)                {                     u=l*a.column+j; v=is*a.column+j;                    d=a.data[u]; a.data[u]=a.data[v]; a.data[v]=d;                }            }            if (js!=l)            {                 for (i=l; i<=a.row-1; i++)                {                     u=i*a.column+js; v=i*a.column+l;                    d=a.data[u]; a.data[u]=a.data[v]; a.data[v]=d;                }            }            ll=l*a.column+l;            for (i=l+1; i<=a.column-1; i++)            {                 d=a(i,l)/a.data[ll];                for (j=l+1; j<=a.column-1; j++)                {                     u=i*a.column+j;                    a.data[u]=a.data[u]-d*a(l,j);                }            }    }    return(k);}

rank函数的使用测试

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