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

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

1、prod函数

prod函数用于求矩阵元素的积，其调用格式如下。

（1）B=prod(A)：若A为向量，则返回所有元素的积；若A为矩阵，则返回各列所有元素的积。

（2）B=prod(A,dim)：返回矩阵A中的第dim维方向的所有元素的积。

prod函数应用示例：

>> prod(1:10)             % 计算10的阶乘

ans =

     3628800

>> M = magic(3)

M =

     8     1     6

     3     5     7

     4     9     2

>> prod(M)              % 列方向求积

ans =

    96    45    84

>> prod(M,2)            % 行方向求积

ans =

    48

   105

    72

2．cumprod函数

cumprod函数用来求矩阵或向量的累积乘积，其调用格式如下。

（1）B=cumprod(A)：若输入参数A为一个向量，则返回该向量所有元素的累积乘积；若A为矩阵，则返回该矩阵列方向各元素的累积乘积。

（2）B=cumprod(A,dim)：A为矩阵，dim为指定维数，若dim=1，则表示在列方向上求累积乘积；若dim=2就代表在行方向上求累积乘积。

【例4-19】  cumprod函数应用示例。

>> cumprod(1:10)

ans =

  Columns 1 through 7

           1           2           6          24         120         720        5040

  Columns 8 through 10

       40320      362880     3628800>> A = [1 2 3; 4 5 6]

A =

     1     2     3

     4     5     6

>> cumprod(A,1)          %   列方向上求累积乘积

ans =

     1     2     3

     4    10    18

>> cumprod(A,2)          %  行方向上求累积乘积

ans =

     1     2     6

     4    20   120

cumprod函数和prod函数的区别：cumprod函数是求累积乘积，结果中含有每一步计算的结果，而prod给出的则是最终乘积的结果。

函数功能

matlab中cumsum函数通常用于计算一个数组各行的累加值。在matlab的命令窗口中输入doc
 cumsum或者help cumsum即可获得该函数的帮助信息。

调用格式及说明

格式一：B
 = cumsum(A)

这种用法返回数组不同维数的累加和。

为了便于接下来的叙述，解释一下matlab中矩阵、数组、向量的概念：



首先，matlab的是矩阵实验室的意思。也就是说matlab中的数据都被视为矩阵。

数组就是一个一行n列的矩阵，向量就是一个n行一列的矩阵。

如果A是一个向量，
 cumsum(A) 返回一个向量，该向量中第m行的元素是A中第1行到第m行的所有元素累加和；

如果A是一个矩阵，
 cumsum(A) 返回一个和A同行同列的矩阵，矩阵中第m行第n列元素是A中第1行到第m行的所有第n列元素的累加和；

如果A是一个多维数组，
 cumsum(A)只对A中第一个非奇异维进行计算。

格式二：B
 = cumsum(A,dim)

这种调用格式返回A中由标量dim所指定的维数的累加和。例如：cumsum(A,1)返回的是沿着第一维（各列）的累加和，cumsum(A,2)返回的是沿着第二维（各行）的累加和。

cumprod 与cumsum用法相同，只不过把相加换成相乘

prod 是cumprod的一种特殊用法，仅适用于向量，与向量作用于cumprod时作用相同

eye 单位矩阵
zeros 全零矩阵
ones 全1矩阵
rand 均匀分布随机阵
genmarkov 生成随机Markov矩阵
linspace 线性等分向量
logspace 对数等分向量
logm 矩阵对数运算
cumprod 矩阵元素累计乘
cumsum 矩阵元素累计和
toeplitz Toeplitz矩阵
disp 显示矩阵和文字内容
length 确定向量的长度
size 确定矩阵的维数
diag 创建对角矩阵或抽取对角向量
find 找出非零元素1的下标
matrix 矩阵变维
rot90 矩阵逆时针旋转90度
sub2ind 全下标转换为单下标
tril 抽取下三角阵
triu 抽取上三角阵
conj 共轭矩阵
companion 伴随矩阵
det 行列式的值
norm 矩阵或向量范数
nnz 矩阵中非零元素的个数
null 清空向量或矩阵中的某个元素
orth 正交基
rank 矩阵秩
trace 矩阵迹
cond 矩阵条件数
inv 矩阵的逆
rcond 逆矩阵条件数
lu LU分解或高斯消元法
pinv 伪逆
qr QR分解
givens Givens变换
linsolve 求解线性方程
lyap Lyapunov方程
hess Hessenberg矩阵
poly 特征多项式
schur Schur分解
expm 矩阵指数
expm1 矩阵指数的Pade逼近
expm2 用泰勒级数求矩阵指数
expm3 通过特征值和特征向量求矩阵指数
funm 计算一般矩阵函数
logm 矩阵对数
sqrtm 矩阵平方根

spec 矩阵特征值
gspec 矩阵束特征值
bdiag 块矩阵，广义特征向量
eigenmar- 正则化Markov特征
kov 向量
pbig 特征空间投影
svd 奇异值分解
sva 奇异值分解近似

cumprod 元素累计积
cumsum 元素累计和
hist 统计频数直方图
max 最大值
min 最小值
mean 平均值
median 中值
prod 元素积
sort 由大到小排序
std 标准差
sum 元素和
trapz 梯形数值积分
corr 求相关系数或方差

sparse 稀疏矩阵
adj2sp 邻接矩阵转换为稀疏矩阵
full 稀疏矩阵转换为全矩阵
mtlb_sparse 将scilab稀疏矩阵转换为matlab稀疏矩阵格式
sp2adj 将稀疏矩阵转换为邻接矩阵
speye 稀疏矩阵方式单位矩阵
sprand 稀疏矩阵方式随机矩阵
spzeros 稀疏矩阵方式全零阵
lufact 稀疏矩阵LU分解
lusolve 稀疏矩阵方程求解
spchol 稀疏矩阵Cholesky分解

关于稀疏矩阵的Matlab命令集，供查阅参考。
基本稀疏矩阵
spdiags :生成稀疏带状矩阵
speye :单位稀疏矩阵
sprand :随机稀疏矩阵
sprandn :正态分布的随机稀疏矩阵
sprandsym:生成稀疏对称随机矩阵

满阵和稀疏矩阵的转换
find :寻找非零元素下标和值
full :稀疏矩阵转化为满阵
sparse :生成稀疏矩阵
spconvert:载入稀疏矩阵

稀疏矩阵的非零元素操作
nnz :非零元素个数
nonzeros :矩阵中的非零元素
nzmax :为非零元素分配的存储空间数
spalloc :稀疏矩阵存储空间
spfun :稀疏矩阵中非零元素的函数计算
spones :非零元素全部用1替换

稀疏矩阵的可视化
spy :稀疏矩阵的图形表示

排序算法
colmmd :进行列的最小度排序
colperm :基于非零算法排序
dmperm :Dulmage-Mendelsohn分解
randperm:随机置换
symmmd :对称最小度排序
symrcm :反向Cuthill-McKee排序

范数、条件数和秩
condest :I范数矩阵条件数
normest :2范数估计

线性方程的稀疏系统
bicg :双共扼梯度法
bicgstab:双共扼梯度稳定法
cgs :二次共扼梯度法
cholinc :不完全Cholesky分解
cholupdate:Cholesky分解的秩1修正
gmres :广义最小残差法
luinc :不完全的LU分解
pcg :预处理共扼梯度法
qmr :Quasi_Minimal残差法

稀疏矩阵的特征值和奇异值
eigs :少数特征值和特征向量
svds :少数奇异值

杂项函数
spparms :设置稀疏矩阵程序的参数排序算法
colmmd :进行列的最小度排序
colperm :基于非零算法排序
dmperm :Dulmage-Mendelsohn分解
randperm:随机置换
symmmd :对称最小度排序
symrcm :反向Cuthill-McKee排序
转载于:https://www.cnblogs.com/haoyul/p/6905723.html