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;
}