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

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

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

1、cumsum函数：

矩阵按列累加求和

2、用法说明

B =cumsum(A)返回矩阵按列累加求和的结果。
如果A是一个向量，将返回一个包含A各元素累加求和的结果，一个同等大小的向量。
如果A是一个矩阵，将返回一个同等大小的矩阵，包含A每一列向量累加求和的结果。

cumsum的c++源码实现

到返回各维度顺序累加和，输入为行向量，返回累加求和
输入：矩阵a
输出：矩阵a各列各维度的累加值。
提示：返回值大小与输入矩阵大小相等。

Matrix Matlab2c::cumsum(Matrix& a) 
{
	int i,j;
	Matrix p(a.row,a.column);
	double sum=0;
	if(a.row==1)  //若数据矩阵为行向量
	{
		for (j=0;j<a.column;j++)
		{
			sum+=a.data[j];
			p.data[j]=sum;
		}
		return p;
	}
	for (j=0;j<p.column;j++)
	{
		sum=0;
		for (i=0;i<p.row;i++)
		{
			sum+=(a(i,j));
			p(i,j)=sum;
		}
	}
	return p;
}

CMatrix Matlab2c::cumsum(CMatrix& a) 
{
	CMatrix p(a.row,a.column);
	Complex sum=0;
	if(a.row==1)  //若数据矩阵为行向量
	{
		for (int j=0;j<a.column;j++)
		{
			sum+=a.data[j];
			p.data[j]=sum;
		}
		return p;
	}
	for (int j=0;j<p.column;j++)
	{
		sum=0;
		for (int i=0;i<p.row;i++)
		{
			sum+=(a(i,j));
			p(i,j)=sum;
		}
	}
	return p;
}

cumsum函数的使用测试

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


int main()
{
	double a[]={1,2,3,4,5,6};
	Matrix aa=Matrix(1,6,a);
	Matrix bb=Matrix(2,3,a);
	Matrix cc=Matlab2c::cumsum(aa);
	cout<<cc.toString()<<endl;
	Matrix dd=Matlab2c::cumsum(bb);
	cout<<dd.toString()<<endl;


	system("pause");
	return 0;
}