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

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

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

1、cumprod函数：

矩阵按列累积相乘

2、用法说明

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

cumprod的c++源码实现

矩阵各维度累积连乘
返回各维度累乘积，输入为矩阵，返回按列累积连乘，输入为向量，返回按序列累积连乘结果

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

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

cumprod函数的使用测试

#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::cumprod(aa);
	cout<<cc.toString()<<endl;
	Matrix dd=Matlab2c::cumprod(bb);
	cout<<dd.toString()<<endl;


	system("pause");
	return 0;
}